diff options
Diffstat (limited to 'man/saem.Rd')
| -rw-r--r-- | man/saem.Rd | 16 |
1 files changed, 13 insertions, 3 deletions
diff --git a/man/saem.Rd b/man/saem.Rd index 2b9199dd..d7b04691 100644 --- a/man/saem.Rd +++ b/man/saem.Rd @@ -82,6 +82,13 @@ automatic choice is not desired} default, uncorrelated random effects are specified for all degradation parameters.} +\item{covariates}{A data frame with covariate data for use in +'covariate_models', with dataset names as row names.} + +\item{covariate_models}{A list containing linear model formulas with one explanatory +variable, i.e. of the type 'parameter ~ covariate'. Covariates must be available +in the 'covariates' data frame.} + \item{no_random_effect}{Character vector of degradation parameters for which there should be no variability over the groups. Only used if the covariance model is not explicitly specified.} @@ -142,10 +149,13 @@ f_mmkin_parent <- mmkin(c("SFO", "FOMC", "DFOP"), ds, quiet = TRUE) f_saem_sfo <- saem(f_mmkin_parent["SFO", ]) f_saem_fomc <- saem(f_mmkin_parent["FOMC", ]) f_saem_dfop <- saem(f_mmkin_parent["DFOP", ]) +anova(f_saem_sfo, f_saem_fomc, f_saem_dfop) +anova(f_saem_sfo, f_saem_dfop, test = TRUE) illparms(f_saem_dfop) -update(f_saem_dfop, covariance.model = diag(c(1, 1, 1, 0))) -AIC(f_saem_dfop) +f_saem_dfop_red <- update(f_saem_dfop, no_random_effect = "g_qlogis") +anova(f_saem_dfop, f_saem_dfop_red, test = TRUE) +anova(f_saem_sfo, f_saem_fomc, f_saem_dfop) # The returned saem.mmkin object contains an SaemixObject, therefore we can use # functions from saemix library(saemix) @@ -157,7 +167,7 @@ plot(f_saem_fomc$so, plot.type = "vpc") f_mmkin_parent_tc <- update(f_mmkin_parent, error_model = "tc") f_saem_fomc_tc <- saem(f_mmkin_parent_tc["FOMC", ]) -compare.saemix(f_saem_fomc$so, f_saem_fomc_tc$so) +anova(f_saem_fomc, f_saem_fomc_tc, test = TRUE) sfo_sfo <- mkinmod(parent = mkinsub("SFO", "A1"), A1 = mkinsub("SFO")) |