From c29e27b9bf5f5361db44e28b06da7b8a1e636e85 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Sat, 9 Apr 2022 06:21:23 +0200 Subject: Improvements to mean_degparms() and plot.mixed.mmkin() - New argument 'default_log_parms' for mean_degparms() - 'plot.mixed.mmkin': Pass the frame argument also to residual plots, take the 'default_log_parms' argument for 'mean_degparms' used for constructing approximate population curves, plot population curve last to avoid that it is covered by data --- man/dimethenamid_2018.Rd | 10 +++++++--- 1 file changed, 7 insertions(+), 3 deletions(-) (limited to 'man/dimethenamid_2018.Rd') diff --git a/man/dimethenamid_2018.Rd b/man/dimethenamid_2018.Rd index 6c28ab7b..20d2dc1b 100644 --- a/man/dimethenamid_2018.Rd +++ b/man/dimethenamid_2018.Rd @@ -62,12 +62,16 @@ plot(mixed(dmta_sfo_sfo3p_tc), test_log_parms = FALSE) # look more plausible, but the truth is likely to be in # between these variants plot(mixed(dmta_sfo_sfo3p_tc), test_log_parms = TRUE) -# Therefore we use nonlinear mixed-effects models +# We can also specify a default value for the failing +# log parameters, to mimic FOCUS guidance +plot(mixed(dmta_sfo_sfo3p_tc), test_log_parms = TRUE, + default_log_parms = log(2)/1000) +# As these attempts are not satisfying, we use nonlinear mixed-effects models # f_dmta_nlme_tc <- nlme(dmta_sfo_sfo3p_tc) # nlme reaches maxIter = 50 without convergence f_dmta_saem_tc <- saem(dmta_sfo_sfo3p_tc) # I am commenting out the convergence plot as rendering them -# with pkgdown fails (at least without further tweaks to the +# with pkgdown fails (at least without further tweaks to the # graphics device used) #saemix::plot(f_dmta_saem_tc$so, plot.type = "convergence") summary(f_dmta_saem_tc) @@ -78,7 +82,7 @@ summary(f_dmta_saem_tc) # covariance.model = diag(c(0, rep(1, 7)))) # saemix::plot(f_dmta_saem_tc_2$so, plot.type = "convergence") # This does not perform better judged by AIC and BIC -saemix::compare.saemix(f_dmta_saem_tc$so, f_dmta_saem_tc_2$so) +# saemix::compare.saemix(f_dmta_saem_tc$so, f_dmta_saem_tc_2$so) } } \keyword{datasets} -- cgit v1.2.1