aboutsummaryrefslogtreecommitdiff
path: root/docs/dev/reference/summary.saem.mmkin.html
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2023-02-13 05:19:08 +0100
committerJohannes Ranke <jranke@uni-bremen.de>2023-02-13 05:19:08 +0100
commit8d1a84ac2190538ed3bac53a303064e281595868 (patch)
treeacb894d85ab7ec87c4911c355a5264a77e08e34b /docs/dev/reference/summary.saem.mmkin.html
parent51d63256a7b3020ee11931d61b4db97b9ded02c0 (diff)
parent4200e566ad2600f56bc3987669aeab88582139eb (diff)
Merge branch 'main' into custom_lsoda_call
Diffstat (limited to 'docs/dev/reference/summary.saem.mmkin.html')
-rw-r--r--docs/dev/reference/summary.saem.mmkin.html551
1 files changed, 305 insertions, 246 deletions
diff --git a/docs/dev/reference/summary.saem.mmkin.html b/docs/dev/reference/summary.saem.mmkin.html
index d8a41356..3b5869f1 100644
--- a/docs/dev/reference/summary.saem.mmkin.html
+++ b/docs/dev/reference/summary.saem.mmkin.html
@@ -21,7 +21,7 @@ endpoints such as formation fractions and DT50 values. Optionally
</button>
<span class="navbar-brand">
<a class="navbar-link" href="../index.html">mkin</a>
- <span class="version label label-info" data-toggle="tooltip" data-placement="bottom" title="In-development version">1.2.0</span>
+ <span class="version label label-info" data-toggle="tooltip" data-placement="bottom" title="In-development version">1.2.2</span>
</span>
</div>
@@ -63,7 +63,10 @@ endpoints such as formation fractions and DT50 values. Optionally
<a href="../articles/web_only/NAFTA_examples.html">Example evaluation of NAFTA SOP Attachment examples</a>
</li>
<li>
- <a href="../articles/web_only/benchmarks.html">Some benchmark timings</a>
+ <a href="../articles/web_only/benchmarks.html">Benchmark timings for mkin</a>
+ </li>
+ <li>
+ <a href="../articles/web_only/saem_benchmarks.html">Benchmark timings for saem.mmkin</a>
</li>
</ul></li>
<li>
@@ -99,7 +102,7 @@ endpoints such as formation fractions and DT50 values. Optionally
<div id="ref-usage">
<div class="sourceCode"><pre class="sourceCode r"><code><span><span class="co"># S3 method for saem.mmkin</span></span>
-<span><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/summary-methods.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">object</span>, data <span class="op">=</span> <span class="cn">FALSE</span>, verbose <span class="op">=</span> <span class="cn">FALSE</span>, distimes <span class="op">=</span> <span class="cn">TRUE</span>, <span class="va">...</span><span class="op">)</span></span>
+<span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">object</span>, data <span class="op">=</span> <span class="cn">FALSE</span>, verbose <span class="op">=</span> <span class="cn">FALSE</span>, distimes <span class="op">=</span> <span class="cn">TRUE</span>, <span class="va">...</span><span class="op">)</span></span>
<span></span>
<span><span class="co"># S3 method for summary.saem.mmkin</span></span>
<span><span class="fu"><a href="https://rdrr.io/r/base/print.html" class="external-link">print</a></span><span class="op">(</span><span class="va">x</span>, digits <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/Extremes.html" class="external-link">max</a></span><span class="op">(</span><span class="fl">3</span>, <span class="fu"><a href="https://rdrr.io/r/base/options.html" class="external-link">getOption</a></span><span class="op">(</span><span class="st">"digits"</span><span class="op">)</span> <span class="op">-</span> <span class="fl">3</span><span class="op">)</span>, verbose <span class="op">=</span> <span class="va">x</span><span class="op">$</span><span class="va">verbose</span>, <span class="va">...</span><span class="op">)</span></span></code></pre></div>
@@ -243,56 +246,58 @@ saemix authors for the parts inherited from saemix.</p>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Likelihood computed by importance sampling</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> AIC BIC logLik</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> 828.1 822.7 -400.1</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 810.8 805.4 -391.4</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Fitted parameters:</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> estimate lower upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 100.74378 97.81291 103.67465</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k_m1 -4.06168 -4.17104 -3.95231</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_qlogis -0.92584 -1.31273 -0.53894</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 -2.81914 -3.60206 -2.03623</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 -3.63916 -4.32672 -2.95161</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis -0.02927 -1.15247 1.09394</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.86164 0.67928 1.04400</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.07973 0.06437 0.09509</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.parent_0 0.73313 -7.46512 8.93137</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k_m1 0.06488 -0.06041 0.19017</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.f_parent_qlogis 0.41955 0.15206 0.68705</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 0.81750 0.29140 1.34361</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 0.75265 0.27939 1.22590</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.g_qlogis 0.34411 -1.70964 2.39786</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> estimate lower upper</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 100.86947 97.81542 103.92353</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k_m1 -4.06947 -4.16944 -3.96950</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_qlogis -0.93256 -1.34200 -0.52312</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 -2.37017 -2.72660 -2.01375</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 -4.06264 -4.21344 -3.91184</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis -0.02174 -0.45898 0.41549</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.87598 0.67275 1.07922</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.07949 0.06389 0.09509</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.parent_0 0.19170 -30.36286 30.74626</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k_m1 0.01883 -0.28736 0.32502</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.f_parent_qlogis 0.44300 0.16391 0.72209</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 0.35320 0.09661 0.60978</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 0.13707 0.02359 0.25056</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.g_qlogis 0.37478 0.04490 0.70467</span>
<span class="r-in"><span><span class="fu"><a href="illparms.html">illparms</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span><span class="op">)</span></span></span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [1] "sd(parent_0)" "sd(log_k_m1)" "sd(g_qlogis)"</span>
-<span class="r-in"><span><span class="va">f_saem_dfop_sfo_2</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/stats/update.html" class="external-link">update</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>, covariance.model <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/diag.html" class="external-link">diag</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">0</span>, <span class="fl">0</span>, <span class="fl">1</span>, <span class="fl">1</span>, <span class="fl">1</span>, <span class="fl">0</span><span class="op">)</span><span class="op">)</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [1] "sd(parent_0)" "sd(log_k_m1)"</span>
+<span class="r-in"><span><span class="va">f_saem_dfop_sfo_2</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/stats/update.html" class="external-link">update</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>,</span></span>
+<span class="r-in"><span> no_random_effect <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"parent_0"</span>, <span class="st">"log_k_m1"</span><span class="op">)</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="illparms.html">illparms</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo_2</span><span class="op">)</span></span></span>
<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/pkg/nlme/man/intervals.html" class="external-link">intervals</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo_2</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Approximate 95% confidence intervals</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Fixed effects:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> lower est. upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 97.54844979 100.46239264 103.37633550</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> k_m1 0.01575805 0.01729111 0.01897331</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_to_m1 0.21014925 0.28626877 0.37680664</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> k1 0.02651112 0.05601399 0.11834908</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> k2 0.01326524 0.02649799 0.05293107</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> g 0.31467778 0.51297098 0.70726363</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 98.36731429 101.42508066 104.48284703</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k_m1 0.01513234 0.01670094 0.01843214</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_to_m1 0.20221431 0.27608850 0.36461630</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k1 0.06915073 0.09759718 0.13774560</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k2 0.01487068 0.01740389 0.02036863</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> g 0.37365671 0.48384821 0.59563299</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Random effects:</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> lower est. upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> sd(f_parent_qlogis) 0.1658367 0.4471180 0.7283993</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> sd(log_k1) 0.2768757 0.7929203 1.3089649</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> sd(log_k2) 0.2693629 0.7566116 1.2438602</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> lower est. upper</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> sd(f_parent_qlogis) 0.16439770 0.4427585 0.7211193</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> sd(log_k1) 0.08304243 0.3345213 0.5860002</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> sd(log_k2) 0.03146410 0.1490210 0.2665779</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> sd(g_qlogis) 0.06216385 0.4023430 0.7425221</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> lower est. upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.70273100 0.88750764 1.07228428</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.06781347 0.08328016 0.09874685</span>
-<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/summary-methods.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo_2</span>, data <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.67696663 0.87777355 1.07858048</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.06363957 0.07878001 0.09392044</span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/summary.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo_2</span>, data <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span> saemix version used for fitting: 3.2 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> mkin version used for pre-fitting: 1.2.0 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> mkin version used for pre-fitting: 1.2.2 </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> R version used for fitting: 4.2.2 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Date of fit: Tue Nov 1 14:12:50 2022 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Date of summary: Tue Nov 1 14:12:50 2022 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Date of fit: Thu Dec 15 14:47:14 2022 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Date of summary: Thu Dec 15 14:47:14 2022 </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Equations:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *</span>
@@ -307,12 +312,12 @@ saemix authors for the parts inherited from saemix.</p>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Model predictions using solution type analytical </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Fitted in 25.006 s</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Fitted in 9.623 s</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Using 300, 100 iterations and 10 chains</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Variance model: Two-component variance function </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Mean of starting values for individual parameters:</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Starting values for degradation parameters:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 log_k_m1 f_parent_qlogis log_k1 log_k2 </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> 101.65645 -4.05368 -0.94311 -2.35943 -4.07006 </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis </span>
@@ -321,237 +326,291 @@ saemix authors for the parts inherited from saemix.</p>
<span class="r-out co"><span class="r-pr">#&gt;</span> Fixed degradation parameter values:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> None</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Starting values for random effects (square root of initial entries in omega):</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 log_k_m1 f_parent_qlogis log_k1 log_k2 g_qlogis</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 6.742 0.0000 0.0000 0.0000 0.0000 0.000</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k_m1 0.000 0.2236 0.0000 0.0000 0.0000 0.000</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_qlogis 0.000 0.0000 0.5572 0.0000 0.0000 0.000</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 0.000 0.0000 0.0000 0.8031 0.0000 0.000</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 0.000 0.0000 0.0000 0.0000 0.2931 0.000</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis 0.000 0.0000 0.0000 0.0000 0.0000 0.807</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Starting values for error model parameters:</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 b.1 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 1 1 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Results:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Likelihood computed by importance sampling</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> AIC BIC logLik</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> 825.6 821.3 -401.8</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> AIC BIC logLik</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 807 802.3 -391.5</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Optimised parameters:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 100.46239 97.54845 103.37634</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k_m1 -4.05756 -4.15040 -3.96472</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_qlogis -0.91358 -1.32403 -0.50312</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 -2.88215 -3.63019 -2.13412</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 -3.63069 -4.32261 -2.93876</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis 0.05190 -0.77834 0.88213</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.88751 0.70273 1.07228</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.08328 0.06781 0.09875</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.f_parent_qlogis 0.44712 0.16584 0.72840</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 0.79292 0.27688 1.30896</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 0.75661 0.26936 1.24386</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 101.42508 98.36731 104.48285</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k_m1 -4.09229 -4.19092 -3.99366</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_qlogis -0.96395 -1.37251 -0.55538</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 -2.32691 -2.67147 -1.98235</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 -4.05106 -4.20836 -3.89376</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis -0.06463 -0.51656 0.38730</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.87777 0.67697 1.07858</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.07878 0.06364 0.09392</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.f_parent_qlogis 0.44276 0.16440 0.72112</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 0.33452 0.08304 0.58600</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 0.14902 0.03146 0.26658</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.g_qlogis 0.40234 0.06216 0.74252</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Correlation: </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> parnt_0 lg_k_m1 f_prnt_ log_k1 log_k2 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k_m1 -0.4102 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_qlogis -0.2113 0.2439 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 0.1308 -0.1305 -0.0504 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 -0.0383 0.0592 0.0151 0.0001 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis -0.0029 -0.0118 0.0131 -0.2547 -0.1942</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k_m1 -0.4693 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_qlogis -0.2378 0.2595 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 0.1720 -0.1593 -0.0669 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 0.0179 0.0594 0.0035 0.1995 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis 0.1073 -0.1060 -0.0322 -0.2299 -0.3168</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Random effects:</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.f_parent_qlogis 0.4471 0.1658 0.7284</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 0.7929 0.2769 1.3090</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 0.7566 0.2694 1.2439</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.f_parent_qlogis 0.4428 0.16440 0.7211</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 0.3345 0.08304 0.5860</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 0.1490 0.03146 0.2666</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.g_qlogis 0.4023 0.06216 0.7425</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Variance model:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.88751 0.70273 1.07228</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.08328 0.06781 0.09875</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.87777 0.67697 1.07858</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.07878 0.06364 0.09392</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Backtransformed parameters:</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 100.46239 97.54845 103.37634</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> k_m1 0.01729 0.01576 0.01897</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_to_m1 0.28627 0.21015 0.37681</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> k1 0.05601 0.02651 0.11835</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> k2 0.02650 0.01327 0.05293</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> g 0.51297 0.31468 0.70726</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 101.4251 98.36731 104.48285</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k_m1 0.0167 0.01513 0.01843</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_to_m1 0.2761 0.20221 0.36462</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k1 0.0976 0.06915 0.13775</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k2 0.0174 0.01487 0.02037</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> g 0.4838 0.37366 0.59563</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Resulting formation fractions:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> ff</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> parent_m1 0.2863</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> parent_sink 0.7137</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_m1 0.2761</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_sink 0.7239</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Estimated disappearance times:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> DT50 DT90 DT50back DT50_k1 DT50_k2</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> parent 17.44 65.15 19.61 12.37 26.16</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> m1 40.09 133.17 NA NA NA</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent 15.54 94.33 28.4 7.102 39.83</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> m1 41.50 137.87 NA NA NA</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Data:</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds name time observed predicted residual std standardized</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 0 89.8 1.005e+02 -10.662393 8.4135 -1.267301</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 0 104.1 1.005e+02 3.637607 8.4135 0.432355</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 1 88.7 9.576e+01 -7.063498 8.0244 -0.880249</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 1 95.5 9.576e+01 -0.263498 8.0244 -0.032837</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 3 81.8 8.717e+01 -5.369491 7.3135 -0.734185</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 3 94.5 8.717e+01 7.330509 7.3135 1.002320</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 7 71.5 7.274e+01 -1.238672 6.1224 -0.202319</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 7 70.3 7.274e+01 -2.438672 6.1224 -0.398322</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 14 54.2 5.418e+01 0.022691 4.5984 0.004935</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 14 49.6 5.418e+01 -4.577309 4.5984 -0.995423</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 28 31.5 3.241e+01 -0.914545 2.8416 -0.321837</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 28 28.8 3.241e+01 -3.614545 2.8416 -1.271993</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 60 12.1 1.283e+01 -0.730904 1.3891 -0.526186</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 60 13.6 1.283e+01 0.769096 1.3891 0.553681</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 90 6.2 6.128e+00 0.071981 1.0238 0.070309</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 90 8.3 6.128e+00 2.171981 1.0238 2.121538</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 120 2.2 3.022e+00 -0.822164 0.9225 -0.891230</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 120 2.4 3.022e+00 -0.622164 0.9225 -0.674429</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 1 0.3 1.163e+00 -0.863423 0.8928 -0.967116</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 1 0.2 1.163e+00 -0.963423 0.8928 -1.079126</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 3 2.2 3.233e+00 -1.032930 0.9274 -1.113734</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 3 3.0 3.233e+00 -0.232930 0.9274 -0.251152</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 7 6.5 6.495e+00 0.005314 1.0393 0.005113</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 7 5.0 6.495e+00 -1.494686 1.0393 -1.438116</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 14 10.2 1.010e+01 0.096372 1.2230 0.078801</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 14 9.5 1.010e+01 -0.603628 1.2230 -0.493572</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 28 12.2 1.269e+01 -0.492073 1.3802 -0.356526</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 28 13.4 1.269e+01 0.707927 1.3802 0.512922</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 60 11.8 1.086e+01 0.944360 1.2669 0.745420</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 60 13.2 1.086e+01 2.344360 1.2669 1.850494</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 90 6.6 7.723e+00 -1.123088 1.0961 -1.024658</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 90 9.3 7.723e+00 1.576912 1.0961 1.438708</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 120 3.5 5.184e+00 -1.683936 0.9869 -1.706219</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 120 5.4 5.184e+00 0.216064 0.9869 0.218923</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 0 118.0 1.005e+02 17.537607 8.4135 2.084469</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 0 99.8 1.005e+02 -0.662393 8.4135 -0.078730</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 1 90.2 9.566e+01 -5.456414 8.0156 -0.680727</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 1 94.6 9.566e+01 -1.056414 8.0156 -0.131795</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 3 96.1 8.702e+01 9.082833 7.3009 1.244062</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 3 78.4 8.702e+01 -8.617167 7.3009 -1.180281</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 7 77.9 7.298e+01 4.919834 6.1423 0.800981</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 7 77.7 7.298e+01 4.719834 6.1423 0.768420</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 14 56.0 5.588e+01 0.124003 4.7372 0.026176</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 14 54.7 5.588e+01 -1.175997 4.7372 -0.248245</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 28 36.6 3.719e+01 -0.587869 3.2217 -0.182474</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 28 36.8 3.719e+01 -0.387869 3.2217 -0.120394</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 60 22.1 2.013e+01 1.973728 1.8966 1.040673</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 60 24.7 2.013e+01 4.573728 1.8966 2.411556</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 90 12.4 1.259e+01 -0.185933 1.3734 -0.135379</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 90 10.8 1.259e+01 -1.785933 1.3734 -1.300347</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 120 6.8 7.981e+00 -1.180542 1.1088 -1.064723</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 120 7.9 7.981e+00 -0.080542 1.1088 -0.072640</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 1 1.3 1.306e+00 -0.006246 0.8941 -0.006986</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 3 3.7 3.589e+00 0.110879 0.9365 0.118399</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 3 4.7 3.589e+00 1.110879 0.9365 1.186217</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 7 8.1 7.062e+00 1.038045 1.0647 0.974978</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 7 7.9 7.062e+00 0.838045 1.0647 0.787129</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 14 10.1 1.065e+01 -0.553713 1.2549 -0.441227</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 14 10.3 1.065e+01 -0.353713 1.2549 -0.281857</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 28 10.7 1.284e+01 -2.144854 1.3900 -1.543111</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 28 12.2 1.284e+01 -0.644854 1.3900 -0.463939</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 60 10.7 1.082e+01 -0.115278 1.2645 -0.091165</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 60 12.5 1.082e+01 1.684722 1.2645 1.332337</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 90 9.1 8.014e+00 1.085607 1.1105 0.977610</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 90 7.4 8.014e+00 -0.614393 1.1105 -0.553272</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 120 6.1 5.736e+00 0.363593 1.0079 0.360737</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 120 4.5 5.736e+00 -1.236407 1.0079 -1.226697</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 0 106.2 1.005e+02 5.737607 8.4135 0.681955</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 0 106.9 1.005e+02 6.437607 8.4135 0.765155</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 1 107.4 9.343e+01 13.972212 7.8311 1.784188</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 1 96.1 9.343e+01 2.672212 7.8311 0.341229</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 3 79.4 8.160e+01 -2.196297 6.8531 -0.320484</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 3 82.6 8.160e+01 1.003703 6.8531 0.146460</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 7 63.9 6.464e+01 -0.737220 5.4557 -0.135129</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 7 62.4 6.464e+01 -2.237220 5.4557 -0.410072</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 14 51.0 4.772e+01 3.278433 4.0722 0.805086</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 14 47.1 4.772e+01 -0.621567 4.0722 -0.152638</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 28 36.1 3.303e+01 3.070676 2.8903 1.062400</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 28 36.6 3.303e+01 3.570676 2.8903 1.235391</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 60 20.1 1.929e+01 0.808039 1.8355 0.440235</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 60 19.8 1.929e+01 0.508039 1.8355 0.276789</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 90 11.3 1.209e+01 -0.794443 1.3425 -0.591785</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 90 10.7 1.209e+01 -1.394443 1.3425 -1.038728</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 120 8.2 7.590e+00 0.610002 1.0896 0.559843</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 120 7.3 7.590e+00 -0.289998 1.0896 -0.266152</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 0 0.8 -4.263e-14 0.800000 0.8875 0.901401</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 1 1.8 1.692e+00 0.107665 0.8986 0.119811</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 1 2.3 1.692e+00 0.607665 0.8986 0.676214</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 3 4.2 4.455e+00 -0.255347 0.9619 -0.265449</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 3 4.1 4.455e+00 -0.355347 0.9619 -0.369404</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 7 6.8 8.124e+00 -1.324338 1.1160 -1.186685</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 7 10.1 8.124e+00 1.975662 1.1160 1.770309</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 14 11.4 1.104e+01 0.361860 1.2778 0.283196</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 14 12.8 1.104e+01 1.761860 1.2778 1.378852</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 28 11.5 1.177e+01 -0.272554 1.3225 -0.206097</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 28 10.6 1.177e+01 -1.172554 1.3225 -0.886648</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 60 7.5 9.242e+00 -1.741667 1.1747 -1.482591</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 60 8.6 9.242e+00 -0.641667 1.1747 -0.546218</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 90 7.3 6.837e+00 0.463318 1.0544 0.439398</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 90 8.1 6.837e+00 1.263318 1.0544 1.198095</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 120 5.3 4.906e+00 0.394322 0.9770 0.403595</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 120 3.8 4.906e+00 -1.105678 0.9770 -1.131677</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 0 104.7 1.005e+02 4.237607 8.4135 0.503670</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 0 88.3 1.005e+02 -12.162393 8.4135 -1.445587</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 1 94.2 9.723e+01 -3.029220 8.1458 -0.371877</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 1 94.6 9.723e+01 -2.629220 8.1458 -0.322772</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 3 78.1 9.114e+01 -13.041804 7.6420 -1.706592</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 3 96.5 9.114e+01 5.358196 7.6420 0.701150</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 7 76.2 8.033e+01 -4.133084 6.7488 -0.612421</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 7 77.8 8.033e+01 -2.533084 6.7488 -0.375340</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 14 70.8 6.504e+01 5.757987 5.4889 1.049017</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 14 67.3 6.504e+01 2.257987 5.4889 0.411371</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 28 43.1 4.418e+01 -1.080806 3.7849 -0.285557</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 28 45.1 4.418e+01 0.919194 3.7849 0.242858</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 60 21.3 2.110e+01 0.200596 1.9686 0.101899</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 60 23.5 2.110e+01 2.400596 1.9686 1.219459</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 90 11.8 1.183e+01 -0.034206 1.3263 -0.025791</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 90 12.1 1.183e+01 0.265794 1.3263 0.200408</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 120 7.0 6.985e+00 0.014647 1.0612 0.013803</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 120 6.2 6.985e+00 -0.785353 1.0612 -0.740078</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 0 1.6 -1.705e-13 1.600000 0.8875 1.802801</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 1 0.9 6.803e-01 0.219655 0.8893 0.246994</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 3 3.7 1.927e+00 1.773027 0.9019 1.965880</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 3 2.0 1.927e+00 0.073027 0.9019 0.080970</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 7 3.6 4.013e+00 -0.412926 0.9483 -0.435417</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 7 3.8 4.013e+00 -0.212926 0.9483 -0.224523</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 14 7.1 6.604e+00 0.495843 1.0441 0.474896</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 14 6.6 6.604e+00 -0.004157 1.0441 -0.003981</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 28 9.5 9.077e+00 0.422700 1.1658 0.362576</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 28 9.3 9.077e+00 0.222700 1.1658 0.191024</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 60 8.3 8.818e+00 -0.518498 1.1520 -0.450099</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 60 9.0 8.818e+00 0.181502 1.1520 0.157559</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 90 6.6 6.738e+00 -0.137785 1.0500 -0.131222</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 90 7.7 6.738e+00 0.962215 1.0500 0.916383</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 120 3.7 4.794e+00 -1.093754 0.9732 -1.123914</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 120 3.5 4.794e+00 -1.293754 0.9732 -1.329429</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 0 110.4 1.005e+02 9.937607 8.4135 1.181155</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 0 112.1 1.005e+02 11.637607 8.4135 1.383212</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 1 93.5 9.372e+01 -0.215694 7.8550 -0.027460</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 1 91.0 9.372e+01 -2.715694 7.8550 -0.345730</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 3 71.0 8.226e+01 -11.257156 6.9076 -1.629667</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 3 89.7 8.226e+01 7.442844 6.9076 1.077480</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 7 60.4 6.553e+01 -5.128464 5.5289 -0.927571</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 7 59.1 6.553e+01 -6.428464 5.5289 -1.162699</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 14 56.5 4.835e+01 8.146351 4.1235 1.975572</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 14 47.0 4.835e+01 -1.353649 4.1235 -0.328273</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 28 30.2 3.300e+01 -2.803303 2.8883 -0.970586</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 28 23.9 3.300e+01 -9.103303 2.8883 -3.151832</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 60 17.0 1.891e+01 -1.905909 1.8074 -1.054506</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 60 18.7 1.891e+01 -0.205909 1.8074 -0.113926</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 90 11.3 1.172e+01 -0.423434 1.3194 -0.320923</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 90 11.9 1.172e+01 0.176566 1.3194 0.133820</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 120 9.0 7.281e+00 1.719138 1.0749 1.599402</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 120 8.1 7.281e+00 0.819138 1.0749 0.762086</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 0 0.7 -2.842e-13 0.700000 0.8875 0.788726</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 1 3.0 3.252e+00 -0.252227 0.9279 -0.271821</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 1 2.6 3.252e+00 -0.652227 0.9279 -0.702895</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 3 5.1 8.615e+00 -3.515326 1.1413 -3.080237</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 3 7.5 8.615e+00 -1.115326 1.1413 -0.977283</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 7 16.5 1.588e+01 0.619041 1.5928 0.388661</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 7 19.0 1.588e+01 3.119041 1.5928 1.958272</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 14 22.9 2.189e+01 1.014705 2.0272 0.500543</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 14 23.2 2.189e+01 1.314705 2.0272 0.648529</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 28 22.2 2.369e+01 -1.487604 2.1632 -0.687701</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 28 24.4 2.369e+01 0.712396 2.1632 0.329332</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 60 15.5 1.869e+01 -3.193942 1.7920 -1.782295</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 60 19.8 1.869e+01 1.106058 1.7920 0.617206</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 90 14.9 1.380e+01 1.103454 1.4518 0.760041</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 90 14.2 1.380e+01 0.403454 1.4518 0.277892</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 120 10.9 9.864e+00 1.035963 1.2093 0.856637</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 120 10.4 9.864e+00 0.535963 1.2093 0.443187</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds name time observed predicted residual std standardized</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 0 89.8 1.014e+02 -11.62508 8.0383 -1.44620</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 0 104.1 1.014e+02 2.67492 8.0383 0.33277</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 1 88.7 9.650e+01 -7.80311 7.6530 -1.01961</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 1 95.5 9.650e+01 -1.00311 7.6530 -0.13107</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 3 81.8 8.753e+01 -5.72638 6.9510 -0.82382</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 3 94.5 8.753e+01 6.97362 6.9510 1.00326</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 7 71.5 7.254e+01 -1.04133 5.7818 -0.18010</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 7 70.3 7.254e+01 -2.24133 5.7818 -0.38765</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 14 54.2 5.349e+01 0.71029 4.3044 0.16502</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 14 49.6 5.349e+01 -3.88971 4.3044 -0.90366</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 28 31.5 3.167e+01 -0.16616 2.6446 -0.06283</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 28 28.8 3.167e+01 -2.86616 2.6446 -1.08379</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 60 12.1 1.279e+01 -0.69287 1.3365 -0.51843</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 60 13.6 1.279e+01 0.80713 1.3365 0.60392</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 90 6.2 6.397e+00 -0.19718 1.0122 -0.19481</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 90 8.3 6.397e+00 1.90282 1.0122 1.87996</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 120 2.2 3.323e+00 -1.12320 0.9160 -1.22623</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 parent 120 2.4 3.323e+00 -0.92320 0.9160 -1.00788</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 1 0.3 1.179e+00 -0.87919 0.8827 -0.99605</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 1 0.2 1.179e+00 -0.97919 0.8827 -1.10935</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 3 2.2 3.273e+00 -1.07272 0.9149 -1.17256</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 3 3.0 3.273e+00 -0.27272 0.9149 -0.29811</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 7 6.5 6.559e+00 -0.05872 1.0186 -0.05765</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 7 5.0 6.559e+00 -1.55872 1.0186 -1.53032</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 14 10.2 1.016e+01 0.03787 1.1880 0.03188</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 14 9.5 1.016e+01 -0.66213 1.1880 -0.55734</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 28 12.2 1.268e+01 -0.47913 1.3297 -0.36032</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 28 13.4 1.268e+01 0.72087 1.3297 0.54211</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 60 11.8 1.078e+01 1.02493 1.2211 0.83936</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 60 13.2 1.078e+01 2.42493 1.2211 1.98588</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 90 6.6 7.705e+00 -1.10464 1.0672 -1.03509</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 90 9.3 7.705e+00 1.59536 1.0672 1.49491</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 120 3.5 5.236e+00 -1.73617 0.9699 -1.79010</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 1 m1 120 5.4 5.236e+00 0.16383 0.9699 0.16892</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 0 118.0 1.014e+02 16.57492 8.0383 2.06198</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 0 99.8 1.014e+02 -1.62508 8.0383 -0.20217</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 1 90.2 9.599e+01 -5.79045 7.6129 -0.76061</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 1 94.6 9.599e+01 -1.39045 7.6129 -0.18264</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 3 96.1 8.652e+01 9.57931 6.8724 1.39388</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 3 78.4 8.652e+01 -8.12069 6.8724 -1.18164</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 7 77.9 7.197e+01 5.93429 5.7370 1.03439</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 7 77.7 7.197e+01 5.73429 5.7370 0.99953</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 14 56.0 5.555e+01 0.44657 4.4637 0.10005</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 14 54.7 5.555e+01 -0.85343 4.4637 -0.19120</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 28 36.6 3.853e+01 -1.93170 3.1599 -0.61132</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 28 36.8 3.853e+01 -1.73170 3.1599 -0.54803</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 60 22.1 2.110e+01 1.00360 1.8795 0.53396</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 60 24.7 2.110e+01 3.60360 1.8795 1.91728</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 90 12.4 1.250e+01 -0.09712 1.3190 -0.07363</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 90 10.8 1.250e+01 -1.69712 1.3190 -1.28667</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 120 6.8 7.419e+00 -0.61913 1.0546 -0.58709</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 parent 120 7.9 7.419e+00 0.48087 1.0546 0.45599</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 1 1.3 1.422e+00 -0.12194 0.8849 -0.13781</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 3 3.7 3.831e+00 -0.13149 0.9282 -0.14166</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 3 4.7 3.831e+00 0.86851 0.9282 0.93567</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 7 8.1 7.292e+00 0.80812 1.0490 0.77034</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 7 7.9 7.292e+00 0.60812 1.0490 0.57969</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 14 10.1 1.055e+01 -0.45332 1.2090 -0.37495</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 14 10.3 1.055e+01 -0.25332 1.2090 -0.20953</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 28 10.7 1.230e+01 -1.59960 1.3074 -1.22347</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 28 12.2 1.230e+01 -0.09960 1.3074 -0.07618</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 60 10.7 1.065e+01 0.05342 1.2141 0.04400</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 60 12.5 1.065e+01 1.85342 1.2141 1.52661</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 90 9.1 8.196e+00 0.90368 1.0897 0.82930</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 90 7.4 8.196e+00 -0.79632 1.0897 -0.73078</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 120 6.1 5.997e+00 0.10252 0.9969 0.10284</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 2 m1 120 4.5 5.997e+00 -1.49748 0.9969 -1.50220</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 0 106.2 1.014e+02 4.77492 8.0383 0.59402</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 0 106.9 1.014e+02 5.47492 8.0383 0.68110</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 1 107.4 9.390e+01 13.49935 7.4494 1.81214</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 1 96.1 9.390e+01 2.19935 7.4494 0.29524</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 3 79.4 8.152e+01 -2.12307 6.4821 -0.32753</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 3 82.6 8.152e+01 1.07693 6.4821 0.16614</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 7 63.9 6.446e+01 -0.55834 5.1533 -0.10834</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 7 62.4 6.446e+01 -2.05834 5.1533 -0.39942</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 14 51.0 4.826e+01 2.74073 3.9019 0.70241</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 14 47.1 4.826e+01 -1.15927 3.9019 -0.29711</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 28 36.1 3.424e+01 1.86399 2.8364 0.65718</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 28 36.6 3.424e+01 2.36399 2.8364 0.83346</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 60 20.1 1.968e+01 0.42172 1.7815 0.23672</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 60 19.8 1.968e+01 0.12172 1.7815 0.06833</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 90 11.3 1.195e+01 -0.64633 1.2869 -0.50222</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 90 10.7 1.195e+01 -1.24633 1.2869 -0.96844</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 120 8.2 7.255e+00 0.94532 1.0474 0.90251</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 parent 120 7.3 7.255e+00 0.04532 1.0474 0.04327</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 0 0.8 2.956e-11 0.80000 0.8778 0.91140</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 1 1.8 1.758e+00 0.04187 0.8886 0.04712</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 1 2.3 1.758e+00 0.54187 0.8886 0.60978</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 3 4.2 4.567e+00 -0.36697 0.9486 -0.38683</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 3 4.1 4.567e+00 -0.46697 0.9486 -0.49224</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 7 6.8 8.151e+00 -1.35124 1.0876 -1.24242</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 7 10.1 8.151e+00 1.94876 1.0876 1.79182</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 14 11.4 1.083e+01 0.57098 1.2240 0.46647</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 14 12.8 1.083e+01 1.97098 1.2240 1.61022</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 28 11.5 1.147e+01 0.03175 1.2597 0.02520</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 28 10.6 1.147e+01 -0.86825 1.2597 -0.68928</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 60 7.5 9.298e+00 -1.79834 1.1433 -1.57298</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 60 8.6 9.298e+00 -0.69834 1.1433 -0.61083</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 90 7.3 7.038e+00 0.26249 1.0382 0.25283</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 90 8.1 7.038e+00 1.06249 1.0382 1.02340</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 120 5.3 5.116e+00 0.18417 0.9659 0.19068</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 3 m1 120 3.8 5.116e+00 -1.31583 0.9659 -1.36232</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 0 104.7 1.014e+02 3.27492 8.0383 0.40741</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 0 88.3 1.014e+02 -13.12508 8.0383 -1.63281</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 1 94.2 9.781e+01 -3.61183 7.7555 -0.46572</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 1 94.6 9.781e+01 -3.21183 7.7555 -0.41414</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 3 78.1 9.110e+01 -13.00467 7.2307 -1.79853</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 3 96.5 9.110e+01 5.39533 7.2307 0.74617</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 7 76.2 7.951e+01 -3.30511 6.3246 -0.52258</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 7 77.8 7.951e+01 -1.70511 6.3246 -0.26960</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 14 70.8 6.376e+01 7.03783 5.0993 1.38016</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 14 67.3 6.376e+01 3.53783 5.0993 0.69379</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 28 43.1 4.340e+01 -0.30456 3.5303 -0.08627</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 28 45.1 4.340e+01 1.69544 3.5303 0.48026</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 60 21.3 2.142e+01 -0.12077 1.9022 -0.06349</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 60 23.5 2.142e+01 2.07923 1.9022 1.09308</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 90 11.8 1.207e+01 -0.26813 1.2940 -0.20721</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 90 12.1 1.207e+01 0.03187 1.2940 0.02463</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 120 7.0 6.954e+00 0.04554 1.0347 0.04402</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 parent 120 6.2 6.954e+00 -0.75446 1.0347 -0.72914</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 0 1.6 1.990e-13 1.60000 0.8778 1.82279</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 1 0.9 7.305e-01 0.16949 0.8797 0.19267</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 3 3.7 2.051e+00 1.64896 0.8925 1.84753</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 3 2.0 2.051e+00 -0.05104 0.8925 -0.05719</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 7 3.6 4.204e+00 -0.60375 0.9382 -0.64354</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 7 3.8 4.204e+00 -0.40375 0.9382 -0.43036</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 14 7.1 6.760e+00 0.34021 1.0267 0.33137</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 14 6.6 6.760e+00 -0.15979 1.0267 -0.15563</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 28 9.5 9.011e+00 0.48856 1.1289 0.43277</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 28 9.3 9.011e+00 0.28856 1.1289 0.25561</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 60 8.3 8.611e+00 -0.31077 1.1093 -0.28014</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 60 9.0 8.611e+00 0.38923 1.1093 0.35086</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 90 6.6 6.678e+00 -0.07753 1.0233 -0.07576</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 90 7.7 6.678e+00 1.02247 1.0233 0.99915</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 120 3.7 4.847e+00 -1.14679 0.9572 -1.19804</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 4 m1 120 3.5 4.847e+00 -1.34679 0.9572 -1.40698</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 0 110.4 1.014e+02 8.97492 8.0383 1.11651</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 0 112.1 1.014e+02 10.67492 8.0383 1.32800</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 1 93.5 9.466e+01 -1.16118 7.5089 -0.15464</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 1 91.0 9.466e+01 -3.66118 7.5089 -0.48758</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 3 71.0 8.302e+01 -12.01844 6.5988 -1.82130</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 3 89.7 8.302e+01 6.68156 6.5988 1.01254</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 7 60.4 6.563e+01 -5.22574 5.2440 -0.99652</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 7 59.1 6.563e+01 -6.52574 5.2440 -1.24442</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 14 56.5 4.727e+01 9.22621 3.8263 2.41128</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 14 47.0 4.727e+01 -0.27379 3.8263 -0.07156</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 28 30.2 3.103e+01 -0.83405 2.5977 -0.32108</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 28 23.9 3.103e+01 -7.13405 2.5977 -2.74634</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 60 17.0 1.800e+01 -0.99696 1.6675 -0.59787</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 60 18.7 1.800e+01 0.70304 1.6675 0.42161</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 90 11.3 1.167e+01 -0.36809 1.2710 -0.28961</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 90 11.9 1.167e+01 0.23191 1.2710 0.18246</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 120 9.0 7.595e+00 1.40496 1.0623 1.32256</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 parent 120 8.1 7.595e+00 0.50496 1.0623 0.47535</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 0 0.7 0.000e+00 0.70000 0.8778 0.79747</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 1 3.0 3.158e+00 -0.15799 0.9123 -0.17317</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 1 2.6 3.158e+00 -0.55799 0.9123 -0.61160</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 3 5.1 8.443e+00 -3.34286 1.1013 -3.03535</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 3 7.5 8.443e+00 -0.94286 1.1013 -0.85613</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 7 16.5 1.580e+01 0.69781 1.5232 0.45811</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 7 19.0 1.580e+01 3.19781 1.5232 2.09935</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 14 22.9 2.216e+01 0.73604 1.9543 0.37663</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 14 23.2 2.216e+01 1.03604 1.9543 0.53014</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 28 22.2 2.423e+01 -2.03128 2.1011 -0.96678</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 28 24.4 2.423e+01 0.16872 2.1011 0.08030</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 60 15.5 1.876e+01 -3.25610 1.7187 -1.89455</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 60 19.8 1.876e+01 1.04390 1.7187 0.60739</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 90 14.9 1.366e+01 1.23585 1.3890 0.88976</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 90 14.2 1.366e+01 0.53585 1.3890 0.38579</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 120 10.9 9.761e+00 1.13911 1.1670 0.97613</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> ds 5 m1 120 10.4 9.761e+00 0.63911 1.1670 0.54767</span>
+<span class="r-in"><span><span class="co"># Add a correlation between random effects of g and k2</span></span></span>
+<span class="r-in"><span><span class="va">cov_model_3</span> <span class="op">&lt;-</span> <span class="va">f_saem_dfop_sfo_2</span><span class="op">$</span><span class="va">so</span><span class="op">@</span><span class="va">model</span><span class="op">@</span><span class="va">covariance.model</span></span></span>
+<span class="r-in"><span><span class="va">cov_model_3</span><span class="op">[</span><span class="st">"log_k2"</span>, <span class="st">"g_qlogis"</span><span class="op">]</span> <span class="op">&lt;-</span> <span class="fl">1</span></span></span>
+<span class="r-in"><span><span class="va">cov_model_3</span><span class="op">[</span><span class="st">"g_qlogis"</span>, <span class="st">"log_k2"</span><span class="op">]</span> <span class="op">&lt;-</span> <span class="fl">1</span></span></span>
+<span class="r-in"><span><span class="va">f_saem_dfop_sfo_3</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/stats/update.html" class="external-link">update</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>,</span></span>
+<span class="r-in"><span> covariance.model <span class="op">=</span> <span class="va">cov_model_3</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/pkg/nlme/man/intervals.html" class="external-link">intervals</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo_3</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Approximate 95% confidence intervals</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Fixed effects:</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> lower est. upper</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 98.39888363 101.48951337 104.58014311</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k_m1 0.01508704 0.01665986 0.01839665</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_parent_to_m1 0.20141557 0.27540583 0.36418131</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k1 0.07708759 0.10430866 0.14114200</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> k2 0.01476621 0.01786384 0.02161129</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> g 0.33679867 0.45083525 0.57028162</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Random effects:</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> lower est. upper</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> sd(f_parent_qlogis) 0.38085375 0.4441841 0.5075145</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> sd(log_k1) 0.04774819 0.2660384 0.4843286</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> sd(log_k2) -0.63842736 0.1977024 1.0338321</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> sd(g_qlogis) 0.22711289 0.4502227 0.6733326</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> corr(log_k2,g_qlogis) -0.83271473 -0.6176939 -0.4026730</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> lower est. upper</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 0.67347568 0.87437392 1.07527216</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> b.1 0.06393032 0.07912417 0.09431802</span>
+<span class="r-in"><span><span class="co"># The correlation does not improve the fit judged by AIC and BIC, although</span></span></span>
+<span class="r-in"><span><span class="co"># the likelihood is higher with the additional parameter</span></span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/anova.html" class="external-link">anova</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>, <span class="va">f_saem_dfop_sfo_2</span>, <span class="va">f_saem_dfop_sfo_3</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Data: 171 observations of 2 variable(s) grouped in 5 datasets</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> npar AIC BIC Lik</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_saem_dfop_sfo_2 12 806.96 802.27 -391.48</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_saem_dfop_sfo_3 13 807.99 802.91 -391.00</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> f_saem_dfop_sfo 14 810.83 805.36 -391.42</span>
<span class="r-in"><span><span class="co"># }</span></span></span>
<span class="r-in"><span></span></span>
</code></pre></div>

Contact - Imprint