diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2022-04-09 06:21:23 +0200 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2022-04-09 06:21:23 +0200 |
| commit | c29e27b9bf5f5361db44e28b06da7b8a1e636e85 (patch) | |
| tree | 4f1b44f2cc6981c8fef8a2085783ad9dc2bdb306 /R/dimethenamid_2018.R | |
| parent | 2728910b96f0ec7dd7ccd97fc6c1f6677e5e352d (diff) | |
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
Diffstat (limited to 'R/dimethenamid_2018.R')
| -rw-r--r-- | R/dimethenamid_2018.R | 10 |
1 files changed, 7 insertions, 3 deletions
diff --git a/R/dimethenamid_2018.R b/R/dimethenamid_2018.R index 2fdd1981..00ed9073 100644 --- a/R/dimethenamid_2018.R +++ b/R/dimethenamid_2018.R @@ -49,12 +49,16 @@ #' # 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) @@ -65,6 +69,6 @@ #' # 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) #' } "dimethenamid_2018" |