From 6e559ef8680053873a8da51a000f6c2f2d28d719 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Mon, 11 Apr 2022 08:53:38 +0200 Subject: Update static docs --- docs/reference/dimethenamid_2018.html | 18 +++++++++++------- 1 file changed, 11 insertions(+), 7 deletions(-) (limited to 'docs/reference/dimethenamid_2018.html') diff --git a/docs/reference/dimethenamid_2018.html b/docs/reference/dimethenamid_2018.html index 85e954ac..d2ed8647 100644 --- a/docs/reference/dimethenamid_2018.html +++ b/docs/reference/dimethenamid_2018.html @@ -186,20 +186,25 @@ specific pieces of information in the comments.

# 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) #> saemix version used for fitting: 3.0 #> mkin version used for pre-fitting: 1.1.0 #> R version used for fitting: 4.1.3 -#> Date of fit: Wed Mar 16 14:32:04 2022 -#> Date of summary: Wed Mar 16 14:32:04 2022 +#> Date of fit: Sat Apr 9 18:03:34 2022 +#> Date of summary: Sat Apr 9 18:03:34 2022 #> #> Equations: #> d_DMTA/dt = - k_DMTA * DMTA @@ -212,7 +217,7 @@ specific pieces of information in the comments.

#> #> Model predictions using solution type deSolve #> -#> Fitted in 927.963 s +#> Fitted in 787.836 s #> Using 300, 100 iterations and 9 chains #> #> Variance model: Two-component variance function @@ -300,8 +305,7 @@ specific pieces of information in the comments.

# 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) -#> Error in saemix::compare.saemix(f_dmta_saem_tc$so, f_dmta_saem_tc_2$so): object 'f_dmta_saem_tc_2' not found +# saemix::compare.saemix(f_dmta_saem_tc$so, f_dmta_saem_tc_2$so) # } -- cgit v1.2.1