aboutsummaryrefslogtreecommitdiff
path: root/man/nlmixr.mmkin.Rd
diff options
context:
space:
mode:
Diffstat (limited to 'man/nlmixr.mmkin.Rd')
-rw-r--r--man/nlmixr.mmkin.Rd39
1 files changed, 25 insertions, 14 deletions
diff --git a/man/nlmixr.mmkin.Rd b/man/nlmixr.mmkin.Rd
index 173b0d39..698c04f0 100644
--- a/man/nlmixr.mmkin.Rd
+++ b/man/nlmixr.mmkin.Rd
@@ -94,7 +94,8 @@ An dataframe suitable for use with \link[nlmixr:nlmixr]{nlmixr::nlmixr}
\description{
This function uses \code{\link[nlmixr:nlmixr]{nlmixr::nlmixr()}} as a backend for fitting nonlinear mixed
effects models created from \link{mmkin} row objects using the Stochastic Approximation
-Expectation Maximisation algorithm (SAEM).
+Expectation Maximisation algorithm (SAEM) or First Order Conditional
+Estimation with Interaction (FOCEI).
}
\details{
An mmkin row object is essentially a list of mkinfit objects that have been
@@ -110,20 +111,31 @@ f_mmkin_parent <- mmkin(c("SFO", "FOMC", "DFOP", "HS"), ds, quiet = TRUE, cores
f_mmkin_parent_tc <- mmkin(c("SFO", "FOMC", "DFOP"), ds, error_model = "tc",
cores = 1, quiet = TRUE)
-f_nlmixr_sfo_saem <- nlmixr(f_mmkin_parent["SFO", ], est = "saem")
-f_nlmixr_sfo_focei <- nlmixr(f_mmkin_parent["SFO", ], est = "focei")
+library(nlmixr)
+f_nlmixr_sfo_saem <- nlmixr(f_mmkin_parent["SFO", ], est = "saem",
+ control = saemControl(print = 0))
+f_nlmixr_sfo_focei <- nlmixr(f_mmkin_parent["SFO", ], est = "focei",
+ control = foceiControl(print = 0))
-f_nlmixr_fomc_saem <- nlmixr(f_mmkin_parent["FOMC", ], est = "saem")
-f_nlmixr_fomc_focei <- nlmixr(f_mmkin_parent["FOMC", ], est = "focei")
+f_nlmixr_fomc_saem <- nlmixr(f_mmkin_parent["FOMC", ], est = "saem",
+ control = saemControl(print = 0))
+f_nlmixr_fomc_focei <- nlmixr(f_mmkin_parent["FOMC", ], est = "focei",
+ control = foceiControl(print = 0))
-f_nlmixr_dfop_saem <- nlmixr(f_mmkin_parent["DFOP", ], est = "saem")
-f_nlmixr_dfop_focei <- nlmixr(f_mmkin_parent["DFOP", ], est = "focei")
+f_nlmixr_dfop_saem <- nlmixr(f_mmkin_parent["DFOP", ], est = "saem",
+ control = saemControl(print = 0))
+f_nlmixr_dfop_focei <- nlmixr(f_mmkin_parent["DFOP", ], est = "focei",
+ control = foceiControl(print = 0))
-f_nlmixr_hs_saem <- nlmixr(f_mmkin_parent["HS", ], est = "saem")
-f_nlmixr_hs_focei <- nlmixr(f_mmkin_parent["HS", ], est = "focei")
+f_nlmixr_hs_saem <- nlmixr(f_mmkin_parent["HS", ], est = "saem",
+ control = saemControl(print = 0))
+f_nlmixr_hs_focei <- nlmixr(f_mmkin_parent["HS", ], est = "focei",
+ control = foceiControl(print = 0))
-f_nlmixr_fomc_saem_tc <- nlmixr(f_mmkin_parent_tc["FOMC", ], est = "saem")
-f_nlmixr_fomc_focei_tc <- nlmixr(f_mmkin_parent_tc["FOMC", ], est = "focei")
+f_nlmixr_fomc_saem_tc <- nlmixr(f_mmkin_parent_tc["FOMC", ], est = "saem",
+ control = saemControl(print = 0))
+f_nlmixr_fomc_focei_tc <- nlmixr(f_mmkin_parent_tc["FOMC", ], est = "focei",
+ control = foceiControl(print = 0))
AIC(
f_nlmixr_sfo_saem$nm, f_nlmixr_sfo_focei$nm,
@@ -135,9 +147,8 @@ AIC(
AIC(nlme(f_mmkin_parent["FOMC", ]))
AIC(nlme(f_mmkin_parent["HS", ]))
-# nlme is comparable to nlmixr with focei, saem finds a better
-# solution, the two-component error model does not improve it
-plot(f_nlmixr_fomc_saem)
+# The FOCEI fit of FOMC with constant variance or the tc error model is best
+plot(f_nlmixr_fomc_saem_tc)
sfo_sfo <- mkinmod(parent = mkinsub("SFO", "A1"),
A1 = mkinsub("SFO"))

Contact - Imprint