From 9275bcb39b5ee25753ef489d334b4906401970b3 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Mon, 14 Nov 2022 21:47:45 +0100 Subject: Update online docs --- docs/dev/reference/saem.html | 516 ++++++++++++++++++++++--------------------- 1 file changed, 264 insertions(+), 252 deletions(-) (limited to 'docs/dev/reference/saem.html') diff --git a/docs/dev/reference/saem.html b/docs/dev/reference/saem.html index 8ea0ef6e..968dc68c 100644 --- a/docs/dev/reference/saem.html +++ b/docs/dev/reference/saem.html @@ -106,13 +106,14 @@ Expectation Maximisation algorithm (SAEM).

conf.level = 0.6, solution_type = "auto", covariance.model = "auto", + omega.init = "auto", covariates = NULL, covariate_models = NULL, no_random_effect = NULL, + error.init = c(3, 0.1), nbiter.saemix = c(300, 100), control = list(displayProgress = FALSE, print = FALSE, nbiter.saemix = nbiter.saemix, save = FALSE, save.graphs = FALSE), - fail_with_errors = TRUE, verbose = FALSE, quiet = FALSE, ... @@ -129,8 +130,10 @@ Expectation Maximisation algorithm (SAEM).

degparms_start = numeric(), covariance.model = "auto", no_random_effect = NULL, + omega.init = "auto", covariates = NULL, covariate_models = NULL, + error.init = numeric(), test_log_parms = FALSE, conf.level = 0.6, verbose = FALSE, @@ -189,11 +192,21 @@ automatic choice is not desired

covariance.model
-

Will be passed to saemix::SaemixModel(). Per +

Will be passed to saemix::saemixModel(). Per default, uncorrelated random effects are specified for all degradation parameters.

+
omega.init
+

Will be passed to saemix::saemixModel(). If using +mkin transformations and the default covariance model with optionally +excluded random effects, the variances of the degradation parameters +are estimated using mean_degparms, with testing of untransformed +log parameters for significant difference from zero. If not using +mkin transformations or a custom covariance model, the default +initialisation of saemix::saemixModel is used for omega.init.

+ +
covariates

A data frame with covariate data for use in 'covariate_models', with dataset names as row names.

@@ -211,6 +224,10 @@ which there should be no variability over the groups. Only used if the covariance model is not explicitly specified.

+
error.init
+

Will be passed to saemix::saemixModel().

+ +
nbiter.saemix

Convenience option to increase the number of iterations

@@ -220,11 +237,6 @@ iterations

Passed to saemix::saemix.

-
fail_with_errors
-

Should a failure to compute standard errors -from the inverse of the Fisher Information Matrix be a failure?

- -
verbose

Should we print information about created objects of type saemix::SaemixModel and saemix::SaemixData?

@@ -293,15 +305,15 @@ using mmkin.

#> Data: 90 observations of 1 variable(s) grouped in 5 datasets #> #> npar AIC BIC Lik -#> f_saem_sfo 5 624.26 622.31 -307.13 -#> f_saem_fomc 7 467.87 465.13 -226.93 -#> f_saem_dfop 9 493.98 490.47 -237.99 +#> f_saem_sfo 5 624.33 622.38 -307.17 +#> f_saem_fomc 7 467.85 465.11 -226.92 +#> f_saem_dfop 9 493.76 490.24 -237.88 anova(f_saem_sfo, f_saem_dfop, test = TRUE) #> Data: 90 observations of 1 variable(s) grouped in 5 datasets #> #> npar AIC BIC Lik Chisq Df Pr(>Chisq) -#> f_saem_sfo 5 624.26 622.31 -307.13 -#> f_saem_dfop 9 493.98 490.47 -237.99 138.28 4 < 2.2e-16 *** +#> f_saem_sfo 5 624.33 622.38 -307.17 +#> f_saem_dfop 9 493.76 490.24 -237.88 138.57 4 < 2.2e-16 *** #> --- #> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 illparms(f_saem_dfop) @@ -311,16 +323,16 @@ using mmkin.

#> Data: 90 observations of 1 variable(s) grouped in 5 datasets #> #> npar AIC BIC Lik Chisq Df Pr(>Chisq) -#> f_saem_dfop_red 8 490.64 487.52 -237.32 -#> f_saem_dfop 9 493.98 490.47 -237.99 0 1 1 +#> f_saem_dfop_red 8 488.68 485.55 -236.34 +#> f_saem_dfop 9 493.76 490.24 -237.88 0 1 1 anova(f_saem_sfo, f_saem_fomc, f_saem_dfop) #> Data: 90 observations of 1 variable(s) grouped in 5 datasets #> #> npar AIC BIC Lik -#> f_saem_sfo 5 624.26 622.31 -307.13 -#> f_saem_fomc 7 467.87 465.13 -226.93 -#> f_saem_dfop 9 493.98 490.47 -237.99 +#> f_saem_sfo 5 624.33 622.38 -307.17 +#> f_saem_fomc 7 467.85 465.11 -226.92 +#> f_saem_dfop 9 493.76 490.24 -237.88 # The returned saem.mmkin object contains an SaemixObject, therefore we can use # functions from saemix library(saemix) @@ -335,9 +347,9 @@ using mmkin.

compare.saemix(f_saem_sfo$so, f_saem_fomc$so, f_saem_dfop$so) #> Likelihoods calculated by importance sampling #> AIC BIC -#> 1 624.2598 622.3070 -#> 2 467.8664 465.1324 -#> 3 493.9811 490.4660 +#> 1 624.3316 622.3788 +#> 2 467.8472 465.1132 +#> 3 493.7592 490.2441 plot(f_saem_fomc$so, plot.type = "convergence") plot(f_saem_fomc$so, plot.type = "individual.fit") @@ -354,9 +366,9 @@ using mmkin.

anova(f_saem_fomc, f_saem_fomc_tc, test = TRUE) #> Data: 90 observations of 1 variable(s) grouped in 5 datasets #> -#> npar AIC BIC Lik Chisq Df Pr(>Chisq) -#> f_saem_fomc 7 467.87 465.13 -226.93 -#> f_saem_fomc_tc 8 469.80 466.68 -226.90 0.0645 1 0.7995 +#> npar AIC BIC Lik Chisq Df Pr(>Chisq) +#> f_saem_fomc 7 467.85 465.11 -226.92 +#> f_saem_fomc_tc 8 469.83 466.71 -226.92 0.015 1 0.9027 sfo_sfo <- mkinmod(parent = mkinsub("SFO", "A1"), A1 = mkinsub("SFO")) @@ -393,32 +405,32 @@ using mmkin.

#> 170 observations of 2 variable(s) grouped in 5 datasets #> #> Likelihood computed by importance sampling -#> AIC BIC logLik -#> 842 836.9 -408 +#> AIC BIC logLik +#> 839.2 834.1 -406.6 #> #> Fitted parameters: -#> estimate lower upper -#> parent_0 93.7701 91.1458 96.3945 -#> log_k_A1 -5.8116 -7.5998 -4.0234 -#> f_parent_qlogis -0.9608 -1.3654 -0.5562 -#> log_k1 -2.5841 -3.6876 -1.4805 -#> log_k2 -3.5228 -5.3254 -1.7203 -#> g_qlogis -0.1027 -0.8719 0.6665 -#> a.1 1.8856 1.6676 2.1037 -#> SD.parent_0 2.7682 0.7668 4.7695 -#> SD.log_k_A1 1.7447 0.4047 3.0848 -#> SD.f_parent_qlogis 0.4525 0.1620 0.7431 -#> SD.log_k1 1.2423 0.4560 2.0285 -#> SD.log_k2 2.0390 0.7601 3.3180 -#> SD.g_qlogis 0.4439 -0.3069 1.1947 +#> estimate lower upper +#> parent_0 93.70402 91.04104 96.3670 +#> log_k_A1 -5.83760 -7.66452 -4.0107 +#> f_parent_qlogis -0.95718 -1.35955 -0.5548 +#> log_k1 -2.35514 -3.39402 -1.3163 +#> log_k2 -3.79634 -5.64009 -1.9526 +#> g_qlogis -0.02108 -0.66463 0.6225 +#> a.1 1.88191 1.66491 2.0989 +#> SD.parent_0 2.81628 0.78922 4.8433 +#> SD.log_k_A1 1.78751 0.42105 3.1540 +#> SD.f_parent_qlogis 0.45016 0.16116 0.7391 +#> SD.log_k1 1.06923 0.31676 1.8217 +#> SD.log_k2 2.03768 0.70938 3.3660 +#> SD.g_qlogis 0.44024 -0.09262 0.9731 plot(f_saem_dfop_sfo) summary(f_saem_dfop_sfo, data = TRUE) #> saemix version used for fitting: 3.2 #> mkin version used for pre-fitting: 1.2.0 #> R version used for fitting: 4.2.2 -#> Date of fit: Tue Nov 1 14:12:07 2022 -#> Date of summary: Tue Nov 1 14:12:07 2022 +#> Date of fit: Mon Nov 14 21:06:30 2022 +#> Date of summary: Mon Nov 14 21:06:30 2022 #> #> Equations: #> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -433,7 +445,7 @@ using mmkin.

#> #> Model predictions using solution type analytical #> -#> Fitted in 8.45 s +#> Fitted in 9.625 s #> Using 300, 100 iterations and 10 chains #> #> Variance model: Constant variance @@ -450,237 +462,237 @@ using mmkin.

#> Results: #> #> Likelihood computed by importance sampling -#> AIC BIC logLik -#> 842 836.9 -408 +#> AIC BIC logLik +#> 839.2 834.1 -406.6 #> #> Optimised parameters: -#> est. lower upper -#> parent_0 93.7701 91.1458 96.3945 -#> log_k_A1 -5.8116 -7.5998 -4.0234 -#> f_parent_qlogis -0.9608 -1.3654 -0.5562 -#> log_k1 -2.5841 -3.6876 -1.4805 -#> log_k2 -3.5228 -5.3254 -1.7203 -#> g_qlogis -0.1027 -0.8719 0.6665 -#> a.1 1.8856 1.6676 2.1037 -#> SD.parent_0 2.7682 0.7668 4.7695 -#> SD.log_k_A1 1.7447 0.4047 3.0848 -#> SD.f_parent_qlogis 0.4525 0.1620 0.7431 -#> SD.log_k1 1.2423 0.4560 2.0285 -#> SD.log_k2 2.0390 0.7601 3.3180 -#> SD.g_qlogis 0.4439 -0.3069 1.1947 +#> est. lower upper +#> parent_0 93.70402 91.04104 96.3670 +#> log_k_A1 -5.83760 -7.66452 -4.0107 +#> f_parent_qlogis -0.95718 -1.35955 -0.5548 +#> log_k1 -2.35514 -3.39402 -1.3163 +#> log_k2 -3.79634 -5.64009 -1.9526 +#> g_qlogis -0.02108 -0.66463 0.6225 +#> a.1 1.88191 1.66491 2.0989 +#> SD.parent_0 2.81628 0.78922 4.8433 +#> SD.log_k_A1 1.78751 0.42105 3.1540 +#> SD.f_parent_qlogis 0.45016 0.16116 0.7391 +#> SD.log_k1 1.06923 0.31676 1.8217 +#> SD.log_k2 2.03768 0.70938 3.3660 +#> SD.g_qlogis 0.44024 -0.09262 0.9731 #> #> Correlation: #> parnt_0 lg_k_A1 f_prnt_ log_k1 log_k2 -#> log_k_A1 -0.0160 -#> f_parent_qlogis -0.0263 0.0612 -#> log_k1 0.0100 -0.0014 -0.0033 -#> log_k2 0.0131 0.0050 -0.0011 0.0071 -#> g_qlogis -0.0419 -0.0199 0.0026 -0.0765 -0.0707 +#> log_k_A1 -0.0147 +#> f_parent_qlogis -0.0269 0.0573 +#> log_k1 0.0263 -0.0011 -0.0040 +#> log_k2 0.0020 0.0065 -0.0002 -0.0776 +#> g_qlogis -0.0248 -0.0180 -0.0004 -0.0903 -0.0603 #> #> Random effects: -#> est. lower upper -#> SD.parent_0 2.7682 0.7668 4.7695 -#> SD.log_k_A1 1.7447 0.4047 3.0848 -#> SD.f_parent_qlogis 0.4525 0.1620 0.7431 -#> SD.log_k1 1.2423 0.4560 2.0285 -#> SD.log_k2 2.0390 0.7601 3.3180 -#> SD.g_qlogis 0.4439 -0.3069 1.1947 +#> est. lower upper +#> SD.parent_0 2.8163 0.78922 4.8433 +#> SD.log_k_A1 1.7875 0.42105 3.1540 +#> SD.f_parent_qlogis 0.4502 0.16116 0.7391 +#> SD.log_k1 1.0692 0.31676 1.8217 +#> SD.log_k2 2.0377 0.70938 3.3660 +#> SD.g_qlogis 0.4402 -0.09262 0.9731 #> #> Variance model: #> est. lower upper -#> a.1 1.886 1.668 2.104 +#> a.1 1.882 1.665 2.099 #> #> Backtransformed parameters: #> est. lower upper -#> parent_0 93.770115 9.115e+01 96.39447 -#> k_A1 0.002993 5.005e-04 0.01789 -#> f_parent_to_A1 0.276720 2.034e-01 0.36443 -#> k1 0.075467 2.503e-02 0.22753 -#> k2 0.029516 4.867e-03 0.17902 -#> g 0.474353 2.949e-01 0.66073 +#> parent_0 93.704015 9.104e+01 96.36699 +#> k_A1 0.002916 4.692e-04 0.01812 +#> f_parent_to_A1 0.277443 2.043e-01 0.36475 +#> k1 0.094880 3.357e-02 0.26813 +#> k2 0.022453 3.553e-03 0.14191 +#> g 0.494731 3.397e-01 0.65078 #> #> Resulting formation fractions: #> ff -#> parent_A1 0.2767 -#> parent_sink 0.7233 +#> parent_A1 0.2774 +#> parent_sink 0.7226 #> #> Estimated disappearance times: -#> DT50 DT90 DT50back DT50_k1 DT50_k2 -#> parent 14.56 58.26 17.54 9.185 23.48 -#> A1 231.62 769.41 NA NA NA +#> DT50 DT90 DT50back DT50_k1 DT50_k2 +#> parent 14.0 72.38 21.79 7.306 30.87 +#> A1 237.7 789.68 NA NA NA #> #> Data: #> ds name time observed predicted residual std standardized -#> Dataset 6 parent 0 97.2 95.78623 1.41377 1.886 0.749758 -#> Dataset 6 parent 0 96.4 95.78623 0.61377 1.886 0.325498 -#> Dataset 6 parent 3 71.1 71.34666 -0.24666 1.886 -0.130812 -#> Dataset 6 parent 3 69.2 71.34666 -2.14666 1.886 -1.138429 -#> Dataset 6 parent 6 58.1 56.49768 1.60232 1.886 0.849749 -#> Dataset 6 parent 6 56.6 56.49768 0.10232 1.886 0.054262 -#> Dataset 6 parent 10 44.4 44.53511 -0.13511 1.886 -0.071650 -#> Dataset 6 parent 10 43.4 44.53511 -1.13511 1.886 -0.601974 -#> Dataset 6 parent 20 33.3 29.77451 3.52549 1.886 1.869656 -#> Dataset 6 parent 20 29.2 29.77451 -0.57451 1.886 -0.304675 -#> Dataset 6 parent 34 17.6 19.32540 -1.72540 1.886 -0.915023 -#> Dataset 6 parent 34 18.0 19.32540 -1.32540 1.886 -0.702894 -#> Dataset 6 parent 55 10.5 10.42781 0.07219 1.886 0.038282 -#> Dataset 6 parent 55 9.3 10.42781 -1.12781 1.886 -0.598107 -#> Dataset 6 parent 90 4.5 3.74190 0.75810 1.886 0.402037 -#> Dataset 6 parent 90 4.7 3.74190 0.95810 1.886 0.508102 -#> Dataset 6 parent 112 3.0 1.96485 1.03515 1.886 0.548966 -#> Dataset 6 parent 112 3.4 1.96485 1.43515 1.886 0.761096 -#> Dataset 6 parent 132 2.3 1.09395 1.20605 1.886 0.639596 -#> Dataset 6 parent 132 2.7 1.09395 1.60605 1.886 0.851726 -#> Dataset 6 A1 3 4.3 4.72702 -0.42702 1.886 -0.226458 -#> Dataset 6 A1 3 4.6 4.72702 -0.12702 1.886 -0.067361 -#> Dataset 6 A1 6 7.0 7.51314 -0.51314 1.886 -0.272128 -#> Dataset 6 A1 6 7.2 7.51314 -0.31314 1.886 -0.166063 -#> Dataset 6 A1 10 8.2 9.63719 -1.43719 1.886 -0.762179 -#> Dataset 6 A1 10 8.0 9.63719 -1.63719 1.886 -0.868244 -#> Dataset 6 A1 20 11.0 11.84931 -0.84931 1.886 -0.450409 -#> Dataset 6 A1 20 13.7 11.84931 1.85069 1.886 0.981468 -#> Dataset 6 A1 34 11.5 12.82336 -1.32336 1.886 -0.701808 -#> Dataset 6 A1 34 12.7 12.82336 -0.12336 1.886 -0.065418 -#> Dataset 6 A1 55 14.9 12.89456 2.00544 1.886 1.063533 -#> Dataset 6 A1 55 14.5 12.89456 1.60544 1.886 0.851403 -#> Dataset 6 A1 90 12.1 11.55919 0.54081 1.886 0.286806 -#> Dataset 6 A1 90 12.3 11.55919 0.74081 1.886 0.392871 -#> Dataset 6 A1 112 9.9 10.42334 -0.52334 1.886 -0.277539 -#> Dataset 6 A1 112 10.2 10.42334 -0.22334 1.886 -0.118442 -#> Dataset 6 A1 132 8.8 9.37987 -0.57987 1.886 -0.307519 -#> Dataset 6 A1 132 7.8 9.37987 -1.57987 1.886 -0.837844 -#> Dataset 7 parent 0 93.6 90.95702 2.64298 1.886 1.401639 -#> Dataset 7 parent 0 92.3 90.95702 1.34298 1.886 0.712217 -#> Dataset 7 parent 3 87.0 84.77506 2.22494 1.886 1.179942 -#> Dataset 7 parent 3 82.2 84.77506 -2.57506 1.886 -1.365616 -#> Dataset 7 parent 7 74.0 77.60962 -3.60962 1.886 -1.914268 -#> Dataset 7 parent 7 73.9 77.60962 -3.70962 1.886 -1.967301 -#> Dataset 7 parent 14 64.2 67.50646 -3.30646 1.886 -1.753499 -#> Dataset 7 parent 14 69.5 67.50646 1.99354 1.886 1.057221 -#> Dataset 7 parent 30 54.0 52.48909 1.51091 1.886 0.801271 -#> Dataset 7 parent 30 54.6 52.48909 2.11091 1.886 1.119465 -#> Dataset 7 parent 60 41.1 39.54372 1.55628 1.886 0.825335 -#> Dataset 7 parent 60 38.4 39.54372 -1.14372 1.886 -0.606542 -#> Dataset 7 parent 90 32.5 33.87968 -1.37968 1.886 -0.731676 -#> Dataset 7 parent 90 35.5 33.87968 1.62032 1.886 0.859298 -#> Dataset 7 parent 120 28.1 30.41071 -2.31071 1.886 -1.225427 -#> Dataset 7 parent 120 29.0 30.41071 -1.41071 1.886 -0.748135 -#> Dataset 7 parent 180 26.5 25.36386 1.13614 1.886 0.602524 -#> Dataset 7 parent 180 27.6 25.36386 2.23614 1.886 1.185881 -#> Dataset 7 A1 3 3.9 2.74863 1.15137 1.886 0.610600 -#> Dataset 7 A1 3 3.1 2.74863 0.35137 1.886 0.186341 -#> Dataset 7 A1 7 6.9 5.92686 0.97314 1.886 0.516081 -#> Dataset 7 A1 7 6.6 5.92686 0.67314 1.886 0.356983 -#> Dataset 7 A1 14 10.4 10.38800 0.01200 1.886 0.006362 -#> Dataset 7 A1 14 8.3 10.38800 -2.08800 1.886 -1.107320 -#> Dataset 7 A1 30 14.4 16.93529 -2.53529 1.886 -1.344524 -#> Dataset 7 A1 30 13.7 16.93529 -3.23529 1.886 -1.715751 -#> Dataset 7 A1 60 22.1 22.33044 -0.23044 1.886 -0.122209 -#> Dataset 7 A1 60 22.3 22.33044 -0.03044 1.886 -0.016144 -#> Dataset 7 A1 90 27.5 24.42300 3.07700 1.886 1.631809 -#> Dataset 7 A1 90 25.4 24.42300 0.97700 1.886 0.518127 -#> Dataset 7 A1 120 28.0 25.51140 2.48860 1.886 1.319768 -#> Dataset 7 A1 120 26.6 25.51140 1.08860 1.886 0.577313 -#> Dataset 7 A1 180 25.8 26.80282 -1.00282 1.886 -0.531818 -#> Dataset 7 A1 180 25.3 26.80282 -1.50282 1.886 -0.796981 -#> Dataset 8 parent 0 91.9 91.08733 0.81267 1.886 0.430980 -#> Dataset 8 parent 0 90.8 91.08733 -0.28733 1.886 -0.152377 -#> Dataset 8 parent 1 64.9 67.55332 -2.65332 1.886 -1.407123 -#> Dataset 8 parent 1 66.2 67.55332 -1.35332 1.886 -0.717701 -#> Dataset 8 parent 3 43.5 41.65811 1.84189 1.886 0.976800 -#> Dataset 8 parent 3 44.1 41.65811 2.44189 1.886 1.294994 -#> Dataset 8 parent 8 18.3 19.65773 -1.35773 1.886 -0.720038 -#> Dataset 8 parent 8 18.1 19.65773 -1.55773 1.886 -0.826103 -#> Dataset 8 parent 14 10.2 10.65118 -0.45118 1.886 -0.239269 -#> Dataset 8 parent 14 10.8 10.65118 0.14882 1.886 0.078925 -#> Dataset 8 parent 27 4.9 3.11694 1.78306 1.886 0.945601 -#> Dataset 8 parent 27 3.3 3.11694 0.18306 1.886 0.097082 -#> Dataset 8 parent 48 1.6 0.43165 1.16835 1.886 0.619603 -#> Dataset 8 parent 48 1.5 0.43165 1.06835 1.886 0.566570 -#> Dataset 8 parent 70 1.1 0.05441 1.04559 1.886 0.554503 -#> Dataset 8 parent 70 0.9 0.05441 0.84559 1.886 0.448438 -#> Dataset 8 A1 1 9.6 7.66431 1.93569 1.886 1.026546 -#> Dataset 8 A1 1 7.7 7.66431 0.03569 1.886 0.018930 -#> Dataset 8 A1 3 15.0 15.57948 -0.57948 1.886 -0.307311 -#> Dataset 8 A1 3 15.1 15.57948 -0.47948 1.886 -0.254279 -#> Dataset 8 A1 8 21.2 20.38988 0.81012 1.886 0.429625 -#> Dataset 8 A1 8 21.1 20.38988 0.71012 1.886 0.376593 -#> Dataset 8 A1 14 19.7 20.16439 -0.46439 1.886 -0.246276 -#> Dataset 8 A1 14 18.9 20.16439 -1.26439 1.886 -0.670535 -#> Dataset 8 A1 27 17.5 16.40918 1.09082 1.886 0.578489 -#> Dataset 8 A1 27 15.9 16.40918 -0.50918 1.886 -0.270030 -#> Dataset 8 A1 48 9.5 10.12011 -0.62011 1.886 -0.328861 -#> Dataset 8 A1 48 9.8 10.12011 -0.32011 1.886 -0.169764 -#> Dataset 8 A1 70 6.2 5.79080 0.40920 1.886 0.217011 -#> Dataset 8 A1 70 6.1 5.79080 0.30920 1.886 0.163979 -#> Dataset 9 parent 0 99.8 97.38786 2.41214 1.886 1.279218 -#> Dataset 9 parent 0 98.3 97.38786 0.91214 1.886 0.483731 -#> Dataset 9 parent 1 77.1 79.25431 -2.15431 1.886 -1.142481 -#> Dataset 9 parent 1 77.2 79.25431 -2.05431 1.886 -1.089449 -#> Dataset 9 parent 3 59.0 55.69866 3.30134 1.886 1.750781 -#> Dataset 9 parent 3 58.1 55.69866 2.40134 1.886 1.273489 -#> Dataset 9 parent 8 27.4 31.64893 -4.24893 1.886 -2.253314 -#> Dataset 9 parent 8 29.2 31.64893 -2.44893 1.886 -1.298729 -#> Dataset 9 parent 14 19.1 22.57316 -3.47316 1.886 -1.841901 -#> Dataset 9 parent 14 29.6 22.57316 7.02684 1.886 3.726507 -#> Dataset 9 parent 27 10.1 14.11345 -4.01345 1.886 -2.128430 -#> Dataset 9 parent 27 18.2 14.11345 4.08655 1.886 2.167199 -#> Dataset 9 parent 48 4.5 6.95586 -2.45586 1.886 -1.302400 -#> Dataset 9 parent 48 9.1 6.95586 2.14414 1.886 1.137093 -#> Dataset 9 parent 70 2.3 3.31753 -1.01753 1.886 -0.539619 -#> Dataset 9 parent 70 2.9 3.31753 -0.41753 1.886 -0.221424 -#> Dataset 9 parent 91 2.0 1.63642 0.36358 1.886 0.192816 -#> Dataset 9 parent 91 1.8 1.63642 0.16358 1.886 0.086751 -#> Dataset 9 parent 120 2.0 0.61667 1.38333 1.886 0.733614 -#> Dataset 9 parent 120 2.2 0.61667 1.58333 1.886 0.839679 -#> Dataset 9 A1 1 4.2 3.67247 0.52753 1.886 0.279763 -#> Dataset 9 A1 1 3.9 3.67247 0.22753 1.886 0.120666 -#> Dataset 9 A1 3 7.4 8.36240 -0.96240 1.886 -0.510385 -#> Dataset 9 A1 3 7.9 8.36240 -0.46240 1.886 -0.245223 -#> Dataset 9 A1 8 14.5 12.80590 1.69410 1.886 0.898422 -#> Dataset 9 A1 8 13.7 12.80590 0.89410 1.886 0.474162 -#> Dataset 9 A1 14 14.2 13.99625 0.20375 1.886 0.108053 -#> Dataset 9 A1 14 12.2 13.99625 -1.79625 1.886 -0.952596 -#> Dataset 9 A1 27 13.7 14.22730 -0.52730 1.886 -0.279641 -#> Dataset 9 A1 27 13.2 14.22730 -1.02730 1.886 -0.544803 -#> Dataset 9 A1 48 13.6 13.33713 0.26287 1.886 0.139406 -#> Dataset 9 A1 48 15.4 13.33713 2.06287 1.886 1.093991 -#> Dataset 9 A1 70 10.4 11.84008 -1.44008 1.886 -0.763708 -#> Dataset 9 A1 70 11.6 11.84008 -0.24008 1.886 -0.127318 -#> Dataset 9 A1 91 10.0 10.30732 -0.30732 1.886 -0.162980 -#> Dataset 9 A1 91 9.5 10.30732 -0.80732 1.886 -0.428142 -#> Dataset 9 A1 120 9.1 8.33981 0.76019 1.886 0.403149 -#> Dataset 9 A1 120 9.0 8.33981 0.66019 1.886 0.350117 -#> Dataset 10 parent 0 96.1 93.70349 2.39651 1.886 1.270926 -#> Dataset 10 parent 0 94.3 93.70349 0.59651 1.886 0.316342 -#> Dataset 10 parent 8 73.9 77.86253 -3.96253 1.886 -2.101429 -#> Dataset 10 parent 8 73.9 77.86253 -3.96253 1.886 -2.101429 -#> Dataset 10 parent 14 69.4 70.18665 -0.78665 1.886 -0.417182 -#> Dataset 10 parent 14 73.1 70.18665 2.91335 1.886 1.545019 -#> Dataset 10 parent 21 65.6 64.03245 1.56755 1.886 0.831308 -#> Dataset 10 parent 21 65.3 64.03245 1.26755 1.886 0.672210 -#> Dataset 10 parent 41 55.9 54.71491 1.18509 1.886 0.628480 -#> Dataset 10 parent 41 54.4 54.71491 -0.31491 1.886 -0.167007 -#> Dataset 10 parent 63 47.0 49.63436 -2.63436 1.886 -1.397065 -#> Dataset 10 parent 63 49.3 49.63436 -0.33436 1.886 -0.177319 -#> Dataset 10 parent 91 44.7 45.08853 -0.38853 1.886 -0.206049 -#> Dataset 10 parent 91 46.7 45.08853 1.61147 1.886 0.854600 -#> Dataset 10 parent 120 42.1 41.07653 1.02347 1.886 0.542772 -#> Dataset 10 parent 120 41.3 41.07653 0.22347 1.886 0.118513 -#> Dataset 10 A1 8 3.3 4.08295 -0.78295 1.886 -0.415218 -#> Dataset 10 A1 8 3.4 4.08295 -0.68295 1.886 -0.362186 -#> Dataset 10 A1 14 3.9 6.04367 -2.14367 1.886 -1.136841 -#> Dataset 10 A1 14 2.9 6.04367 -3.14367 1.886 -1.667165 -#> Dataset 10 A1 21 6.4 7.59693 -1.19693 1.886 -0.634761 -#> Dataset 10 A1 21 7.2 7.59693 -0.39693 1.886 -0.210502 -#> Dataset 10 A1 41 9.1 9.86436 -0.76436 1.886 -0.405361 -#> Dataset 10 A1 41 8.5 9.86436 -1.36436 1.886 -0.723555 -#> Dataset 10 A1 63 11.7 10.99397 0.70603 1.886 0.374425 -#> Dataset 10 A1 63 12.0 10.99397 1.00603 1.886 0.533522 -#> Dataset 10 A1 91 13.3 11.91274 1.38726 1.886 0.735696 -#> Dataset 10 A1 91 13.2 11.91274 1.28726 1.886 0.682663 -#> Dataset 10 A1 120 14.3 12.66519 1.63481 1.886 0.866981 -#> Dataset 10 A1 120 12.1 12.66519 -0.56519 1.886 -0.299733 +#> Dataset 6 parent 0 97.2 95.70025 1.49975 1.882 0.79693 +#> Dataset 6 parent 0 96.4 95.70025 0.69975 1.882 0.37183 +#> Dataset 6 parent 3 71.1 71.44670 -0.34670 1.882 -0.18423 +#> Dataset 6 parent 3 69.2 71.44670 -2.24670 1.882 -1.19384 +#> Dataset 6 parent 6 58.1 56.59283 1.50717 1.882 0.80087 +#> Dataset 6 parent 6 56.6 56.59283 0.00717 1.882 0.00381 +#> Dataset 6 parent 10 44.4 44.56648 -0.16648 1.882 -0.08847 +#> Dataset 6 parent 10 43.4 44.56648 -1.16648 1.882 -0.61984 +#> Dataset 6 parent 20 33.3 29.76020 3.53980 1.882 1.88096 +#> Dataset 6 parent 20 29.2 29.76020 -0.56020 1.882 -0.29767 +#> Dataset 6 parent 34 17.6 19.39208 -1.79208 1.882 -0.95226 +#> Dataset 6 parent 34 18.0 19.39208 -1.39208 1.882 -0.73971 +#> Dataset 6 parent 55 10.5 10.55761 -0.05761 1.882 -0.03061 +#> Dataset 6 parent 55 9.3 10.55761 -1.25761 1.882 -0.66826 +#> Dataset 6 parent 90 4.5 3.84742 0.65258 1.882 0.34676 +#> Dataset 6 parent 90 4.7 3.84742 0.85258 1.882 0.45304 +#> Dataset 6 parent 112 3.0 2.03997 0.96003 1.882 0.51013 +#> Dataset 6 parent 112 3.4 2.03997 1.36003 1.882 0.72268 +#> Dataset 6 parent 132 2.3 1.14585 1.15415 1.882 0.61328 +#> Dataset 6 parent 132 2.7 1.14585 1.55415 1.882 0.82583 +#> Dataset 6 A1 3 4.3 4.86054 -0.56054 1.882 -0.29786 +#> Dataset 6 A1 3 4.6 4.86054 -0.26054 1.882 -0.13844 +#> Dataset 6 A1 6 7.0 7.74179 -0.74179 1.882 -0.39417 +#> Dataset 6 A1 6 7.2 7.74179 -0.54179 1.882 -0.28789 +#> Dataset 6 A1 10 8.2 9.94048 -1.74048 1.882 -0.92485 +#> Dataset 6 A1 10 8.0 9.94048 -1.94048 1.882 -1.03112 +#> Dataset 6 A1 20 11.0 12.19109 -1.19109 1.882 -0.63291 +#> Dataset 6 A1 20 13.7 12.19109 1.50891 1.882 0.80180 +#> Dataset 6 A1 34 11.5 13.10706 -1.60706 1.882 -0.85395 +#> Dataset 6 A1 34 12.7 13.10706 -0.40706 1.882 -0.21630 +#> Dataset 6 A1 55 14.9 13.06131 1.83869 1.882 0.97703 +#> Dataset 6 A1 55 14.5 13.06131 1.43869 1.882 0.76448 +#> Dataset 6 A1 90 12.1 11.54495 0.55505 1.882 0.29494 +#> Dataset 6 A1 90 12.3 11.54495 0.75505 1.882 0.40122 +#> Dataset 6 A1 112 9.9 10.31533 -0.41533 1.882 -0.22070 +#> Dataset 6 A1 112 10.2 10.31533 -0.11533 1.882 -0.06128 +#> Dataset 6 A1 132 8.8 9.20222 -0.40222 1.882 -0.21373 +#> Dataset 6 A1 132 7.8 9.20222 -1.40222 1.882 -0.74510 +#> Dataset 7 parent 0 93.6 90.82357 2.77643 1.882 1.47532 +#> Dataset 7 parent 0 92.3 90.82357 1.47643 1.882 0.78453 +#> Dataset 7 parent 3 87.0 84.73448 2.26552 1.882 1.20384 +#> Dataset 7 parent 3 82.2 84.73448 -2.53448 1.882 -1.34675 +#> Dataset 7 parent 7 74.0 77.65013 -3.65013 1.882 -1.93958 +#> Dataset 7 parent 7 73.9 77.65013 -3.75013 1.882 -1.99272 +#> Dataset 7 parent 14 64.2 67.60639 -3.40639 1.882 -1.81007 +#> Dataset 7 parent 14 69.5 67.60639 1.89361 1.882 1.00621 +#> Dataset 7 parent 30 54.0 52.53663 1.46337 1.882 0.77760 +#> Dataset 7 parent 30 54.6 52.53663 2.06337 1.882 1.09642 +#> Dataset 7 parent 60 41.1 39.42728 1.67272 1.882 0.88884 +#> Dataset 7 parent 60 38.4 39.42728 -1.02728 1.882 -0.54587 +#> Dataset 7 parent 90 32.5 33.76360 -1.26360 1.882 -0.67144 +#> Dataset 7 parent 90 35.5 33.76360 1.73640 1.882 0.92268 +#> Dataset 7 parent 120 28.1 30.39975 -2.29975 1.882 -1.22203 +#> Dataset 7 parent 120 29.0 30.39975 -1.39975 1.882 -0.74379 +#> Dataset 7 parent 180 26.5 25.62379 0.87621 1.882 0.46559 +#> Dataset 7 parent 180 27.6 25.62379 1.97621 1.882 1.05010 +#> Dataset 7 A1 3 3.9 2.70005 1.19995 1.882 0.63762 +#> Dataset 7 A1 3 3.1 2.70005 0.39995 1.882 0.21252 +#> Dataset 7 A1 7 6.9 5.83475 1.06525 1.882 0.56605 +#> Dataset 7 A1 7 6.6 5.83475 0.76525 1.882 0.40663 +#> Dataset 7 A1 14 10.4 10.26142 0.13858 1.882 0.07364 +#> Dataset 7 A1 14 8.3 10.26142 -1.96142 1.882 -1.04225 +#> Dataset 7 A1 30 14.4 16.82999 -2.42999 1.882 -1.29123 +#> Dataset 7 A1 30 13.7 16.82999 -3.12999 1.882 -1.66319 +#> Dataset 7 A1 60 22.1 22.32486 -0.22486 1.882 -0.11949 +#> Dataset 7 A1 60 22.3 22.32486 -0.02486 1.882 -0.01321 +#> Dataset 7 A1 90 27.5 24.45927 3.04073 1.882 1.61576 +#> Dataset 7 A1 90 25.4 24.45927 0.94073 1.882 0.49988 +#> Dataset 7 A1 120 28.0 25.54862 2.45138 1.882 1.30260 +#> Dataset 7 A1 120 26.6 25.54862 1.05138 1.882 0.55868 +#> Dataset 7 A1 180 25.8 26.82277 -1.02277 1.882 -0.54347 +#> Dataset 7 A1 180 25.3 26.82277 -1.52277 1.882 -0.80916 +#> Dataset 8 parent 0 91.9 91.16791 0.73209 1.882 0.38901 +#> Dataset 8 parent 0 90.8 91.16791 -0.36791 1.882 -0.19550 +#> Dataset 8 parent 1 64.9 67.58358 -2.68358 1.882 -1.42598 +#> Dataset 8 parent 1 66.2 67.58358 -1.38358 1.882 -0.73520 +#> Dataset 8 parent 3 43.5 41.62086 1.87914 1.882 0.99853 +#> Dataset 8 parent 3 44.1 41.62086 2.47914 1.882 1.31735 +#> Dataset 8 parent 8 18.3 19.60116 -1.30116 1.882 -0.69140 +#> Dataset 8 parent 8 18.1 19.60116 -1.50116 1.882 -0.79768 +#> Dataset 8 parent 14 10.2 10.63101 -0.43101 1.882 -0.22903 +#> Dataset 8 parent 14 10.8 10.63101 0.16899 1.882 0.08980 +#> Dataset 8 parent 27 4.9 3.12435 1.77565 1.882 0.94354 +#> Dataset 8 parent 27 3.3 3.12435 0.17565 1.882 0.09334 +#> Dataset 8 parent 48 1.6 0.43578 1.16422 1.882 0.61864 +#> Dataset 8 parent 48 1.5 0.43578 1.06422 1.882 0.56550 +#> Dataset 8 parent 70 1.1 0.05534 1.04466 1.882 0.55510 +#> Dataset 8 parent 70 0.9 0.05534 0.84466 1.882 0.44883 +#> Dataset 8 A1 1 9.6 7.63450 1.96550 1.882 1.04442 +#> Dataset 8 A1 1 7.7 7.63450 0.06550 1.882 0.03481 +#> Dataset 8 A1 3 15.0 15.52593 -0.52593 1.882 -0.27947 +#> Dataset 8 A1 3 15.1 15.52593 -0.42593 1.882 -0.22633 +#> Dataset 8 A1 8 21.2 20.32192 0.87808 1.882 0.46659 +#> Dataset 8 A1 8 21.1 20.32192 0.77808 1.882 0.41345 +#> Dataset 8 A1 14 19.7 20.09721 -0.39721 1.882 -0.21107 +#> Dataset 8 A1 14 18.9 20.09721 -1.19721 1.882 -0.63617 +#> Dataset 8 A1 27 17.5 16.37477 1.12523 1.882 0.59792 +#> Dataset 8 A1 27 15.9 16.37477 -0.47477 1.882 -0.25228 +#> Dataset 8 A1 48 9.5 10.13141 -0.63141 1.882 -0.33551 +#> Dataset 8 A1 48 9.8 10.13141 -0.33141 1.882 -0.17610 +#> Dataset 8 A1 70 6.2 5.81827 0.38173 1.882 0.20284 +#> Dataset 8 A1 70 6.1 5.81827 0.28173 1.882 0.14970 +#> Dataset 9 parent 0 99.8 97.48728 2.31272 1.882 1.22892 +#> Dataset 9 parent 0 98.3 97.48728 0.81272 1.882 0.43186 +#> Dataset 9 parent 1 77.1 79.29476 -2.19476 1.882 -1.16624 +#> Dataset 9 parent 1 77.2 79.29476 -2.09476 1.882 -1.11310 +#> Dataset 9 parent 3 59.0 55.67060 3.32940 1.882 1.76915 +#> Dataset 9 parent 3 58.1 55.67060 2.42940 1.882 1.29092 +#> Dataset 9 parent 8 27.4 31.57871 -4.17871 1.882 -2.22046 +#> Dataset 9 parent 8 29.2 31.57871 -2.37871 1.882 -1.26398 +#> Dataset 9 parent 14 19.1 22.51546 -3.41546 1.882 -1.81489 +#> Dataset 9 parent 14 29.6 22.51546 7.08454 1.882 3.76454 +#> Dataset 9 parent 27 10.1 14.09074 -3.99074 1.882 -2.12057 +#> Dataset 9 parent 27 18.2 14.09074 4.10926 1.882 2.18355 +#> Dataset 9 parent 48 4.5 6.95747 -2.45747 1.882 -1.30584 +#> Dataset 9 parent 48 9.1 6.95747 2.14253 1.882 1.13848 +#> Dataset 9 parent 70 2.3 3.32472 -1.02472 1.882 -0.54451 +#> Dataset 9 parent 70 2.9 3.32472 -0.42472 1.882 -0.22569 +#> Dataset 9 parent 91 2.0 1.64300 0.35700 1.882 0.18970 +#> Dataset 9 parent 91 1.8 1.64300 0.15700 1.882 0.08343 +#> Dataset 9 parent 120 2.0 0.62073 1.37927 1.882 0.73291 +#> Dataset 9 parent 120 2.2 0.62073 1.57927 1.882 0.83918 +#> Dataset 9 A1 1 4.2 3.64568 0.55432 1.882 0.29455 +#> Dataset 9 A1 1 3.9 3.64568 0.25432 1.882 0.13514 +#> Dataset 9 A1 3 7.4 8.30173 -0.90173 1.882 -0.47916 +#> Dataset 9 A1 3 7.9 8.30173 -0.40173 1.882 -0.21347 +#> Dataset 9 A1 8 14.5 12.71589 1.78411 1.882 0.94803 +#> Dataset 9 A1 8 13.7 12.71589 0.98411 1.882 0.52293 +#> Dataset 9 A1 14 14.2 13.90452 0.29548 1.882 0.15701 +#> Dataset 9 A1 14 12.2 13.90452 -1.70452 1.882 -0.90574 +#> Dataset 9 A1 27 13.7 14.15523 -0.45523 1.882 -0.24190 +#> Dataset 9 A1 27 13.2 14.15523 -0.95523 1.882 -0.50759 +#> Dataset 9 A1 48 13.6 13.31038 0.28962 1.882 0.15389 +#> Dataset 9 A1 48 15.4 13.31038 2.08962 1.882 1.11037 +#> Dataset 9 A1 70 10.4 11.85965 -1.45965 1.882 -0.77562 +#> Dataset 9 A1 70 11.6 11.85965 -0.25965 1.882 -0.13797 +#> Dataset 9 A1 91 10.0 10.36294 -0.36294 1.882 -0.19286 +#> Dataset 9 A1 91 9.5 10.36294 -0.86294 1.882 -0.45855 +#> Dataset 9 A1 120 9.1 8.43003 0.66997 1.882 0.35601 +#> Dataset 9 A1 120 9.0 8.43003 0.56997 1.882 0.30287 +#> Dataset 10 parent 0 96.1 93.95603 2.14397 1.882 1.13925 +#> Dataset 10 parent 0 94.3 93.95603 0.34397 1.882 0.18278 +#> Dataset 10 parent 8 73.9 77.70592 -3.80592 1.882 -2.02237 +#> Dataset 10 parent 8 73.9 77.70592 -3.80592 1.882 -2.02237 +#> Dataset 10 parent 14 69.4 70.04570 -0.64570 1.882 -0.34311 +#> Dataset 10 parent 14 73.1 70.04570 3.05430 1.882 1.62298 +#> Dataset 10 parent 21 65.6 64.01710 1.58290 1.882 0.84111 +#> Dataset 10 parent 21 65.3 64.01710 1.28290 1.882 0.68170 +#> Dataset 10 parent 41 55.9 54.98434 0.91566 1.882 0.48656 +#> Dataset 10 parent 41 54.4 54.98434 -0.58434 1.882 -0.31050 +#> Dataset 10 parent 63 47.0 49.87137 -2.87137 1.882 -1.52577 +#> Dataset 10 parent 63 49.3 49.87137 -0.57137 1.882 -0.30361 +#> Dataset 10 parent 91 44.7 45.06727 -0.36727 1.882 -0.19516 +#> Dataset 10 parent 91 46.7 45.06727 1.63273 1.882 0.86759 +#> Dataset 10 parent 120 42.1 40.76402 1.33598 1.882 0.70991 +#> Dataset 10 parent 120 41.3 40.76402 0.53598 1.882 0.28481 +#> Dataset 10 A1 8 3.3 4.14599 -0.84599 1.882 -0.44954 +#> Dataset 10 A1 8 3.4 4.14599 -0.74599 1.882 -0.39640 +#> Dataset 10 A1 14 3.9 6.08478 -2.18478 1.882 -1.16093 +#> Dataset 10 A1 14 2.9 6.08478 -3.18478 1.882 -1.69231 +#> Dataset 10 A1 21 6.4 7.59411 -1.19411 1.882 -0.63452 +#> Dataset 10 A1 21 7.2 7.59411 -0.39411 1.882 -0.20942 +#> Dataset 10 A1 41 9.1 9.78292 -0.68292 1.882 -0.36289 +#> Dataset 10 A1 41 8.5 9.78292 -1.28292 1.882 -0.68171 +#> Dataset 10 A1 63 11.7 10.93274 0.76726 1.882 0.40770 +#> Dataset 10 A1 63 12.0 10.93274 1.06726 1.882 0.56711 +#> Dataset 10 A1 91 13.3 11.93986 1.36014 1.882 0.72274 +#> Dataset 10 A1 91 13.2 11.93986 1.26014 1.882 0.66961 +#> Dataset 10 A1 120 14.3 12.79238 1.50762 1.882 0.80111 +#> Dataset 10 A1 120 12.1 12.79238 -0.69238 1.882 -0.36791 # The following takes about 6 minutes #f_saem_dfop_sfo_deSolve <- saem(f_mmkin["DFOP-SFO", ], solution_type = "deSolve", -- cgit v1.2.1