Update static docsv1.2.5

author: Johannes Ranke <jranke@uni-bremen.de> 2023-08-09 17:59:42 +0200
committer: Johannes Ranke <jranke@uni-bremen.de> 2023-08-09 17:59:42 +0200
commit: 9abab1e2d4385039b01ad3dc0d9c5966bbe94fee (patch)
tree: cc9f23c3abe70e289b40c2cb03a4cda4ed157122 /docs
parent: f21373a16d368c56234246417a32bdfb95589cf9 (diff)
22 files changed, 2645 insertions, 73 deletions
diff --git a/docs/404.html b/docs/404.html
index b9ab1f25..932d8345 100644
--- a/docs/404.html
+++ b/docs/404.html
@@ -32,7 +32,7 @@
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="https://pkgdown.jrwb.de/mkin/index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
diff --git a/docs/articles/FOCUS_L.html b/docs/articles/FOCUS_L.html
index dabb16f2..7b5acf17 100644
--- a/docs/articles/FOCUS_L.html
+++ b/docs/articles/FOCUS_L.html
@@ -33,7 +33,7 @@
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="../index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
@@ -135,7 +135,7 @@ to L3</h1>
 Ranke</h4>
             
             <h4 data-toc-skip class="date">Last change 18 May 2023
-(rebuilt 2023-05-19)</h4>
+(rebuilt 2023-08-09)</h4>
       
       <small class="dont-index">Source: <a href="https://github.com/jranke/mkin/blob/HEAD/vignettes/FOCUS_L.rmd" class="external-link"><code>vignettes/FOCUS_L.rmd</code></a></small>
       <div class="hidden name"><code>FOCUS_L.rmd</code></div>
@@ -168,18 +168,18 @@ model fit. This covers the numerical analysis given in the FOCUS
 report.</p>
 <div class="sourceCode" id="cb2"><pre class="downlit sourceCode r">
 <code class="sourceCode R"><span><span class="va">m.L1.SFO</span> <span class="op">&lt;-</span> <span class="fu"><a href="../reference/mkinfit.html">mkinfit</a></span><span class="op">(</span><span class="st">"SFO"</span>, <span class="va">FOCUS_2006_L1_mkin</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</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">m.L1.SFO</span><span class="op">)</span></span></code></pre></div>
-<pre><code><span><span class="co">## mkin version used for fitting:    1.2.4 </span></span>
-<span><span class="co">## R version used for fitting:       4.3.0 </span></span>
-<span><span class="co">## Date of fit:     Fri May 19 09:20:25 2023 </span></span>
-<span><span class="co">## Date of summary: Fri May 19 09:20:25 2023 </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">m.L1.SFO</span><span class="op">)</span></span></code></pre></div>
+<pre><code><span><span class="co">## mkin version used for fitting:    1.2.5 </span></span>
+<span><span class="co">## R version used for fitting:       4.3.1 </span></span>
+<span><span class="co">## Date of fit:     Wed Aug  9 17:55:39 2023 </span></span>
+<span><span class="co">## Date of summary: Wed Aug  9 17:55:39 2023 </span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Equations:</span></span>
 <span><span class="co">## d_parent/dt = - k_parent * parent</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Model predictions using solution type analytical </span></span>
 <span><span class="co">## </span></span>
-<span><span class="co">## Fitted using 133 model solutions performed in 0.011 s</span></span>
+<span><span class="co">## Fitted using 133 model solutions performed in 0.031 s</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Error model: Constant variance </span></span>
 <span><span class="co">## </span></span>
@@ -256,7 +256,7 @@ report.</p>
 <p>A plot of the fit is obtained with the plot function for mkinfit
 objects.</p>
 <div class="sourceCode" id="cb4"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L1.SFO</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span>, main <span class="op">=</span> <span class="st">"FOCUS L1 - SFO"</span><span class="op">)</span></span></code></pre></div>
+<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L1.SFO</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span>, main <span class="op">=</span> <span class="st">"FOCUS L1 - SFO"</span><span class="op">)</span></span></code></pre></div>
 <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-4-1.png" width="576"></p>
 <p>The residual plot can be easily obtained by</p>
 <div class="sourceCode" id="cb5"><pre class="downlit sourceCode r">
@@ -268,25 +268,25 @@ objects.</p>
 <pre><code><span><span class="co">## Warning in mkinfit("FOMC", FOCUS_2006_L1_mkin, quiet = TRUE): Optimisation did not converge:</span></span>
 <span><span class="co">## false convergence (8)</span></span></code></pre>
 <div class="sourceCode" id="cb8"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L1.FOMC</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span>, main <span class="op">=</span> <span class="st">"FOCUS L1 - FOMC"</span><span class="op">)</span></span></code></pre></div>
+<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L1.FOMC</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span>, main <span class="op">=</span> <span class="st">"FOCUS L1 - FOMC"</span><span class="op">)</span></span></code></pre></div>
 <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-6-1.png" width="576"></p>
 <div class="sourceCode" id="cb9"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><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">m.L1.FOMC</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
+<code class="sourceCode R"><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">m.L1.FOMC</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
 <pre><code><span><span class="co">## Warning in sqrt(diag(covar)): NaNs produced</span></span></code></pre>
 <pre><code><span><span class="co">## Warning in sqrt(1/diag(V)): NaNs produced</span></span></code></pre>
 <pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(.) had 0 or NA entries; non-finite result</span></span>
 <span><span class="co">## is doubtful</span></span></code></pre>
-<pre><code><span><span class="co">## mkin version used for fitting:    1.2.4 </span></span>
-<span><span class="co">## R version used for fitting:       4.3.0 </span></span>
-<span><span class="co">## Date of fit:     Fri May 19 09:20:25 2023 </span></span>
-<span><span class="co">## Date of summary: Fri May 19 09:20:25 2023 </span></span>
+<pre><code><span><span class="co">## mkin version used for fitting:    1.2.5 </span></span>
+<span><span class="co">## R version used for fitting:       4.3.1 </span></span>
+<span><span class="co">## Date of fit:     Wed Aug  9 17:55:39 2023 </span></span>
+<span><span class="co">## Date of summary: Wed Aug  9 17:55:39 2023 </span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Equations:</span></span>
 <span><span class="co">## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Model predictions using solution type analytical </span></span>
 <span><span class="co">## </span></span>
-<span><span class="co">## Fitted using 342 model solutions performed in 0.021 s</span></span>
+<span><span class="co">## Fitted using 342 model solutions performed in 0.07 s</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Error model: Constant variance </span></span>
 <span><span class="co">## </span></span>
@@ -398,7 +398,7 @@ residual plot can be obtained simply by adding the argument
 <code>show_residuals</code> to the plot command.</p>
 <div class="sourceCode" id="cb15"><pre class="downlit sourceCode r">
 <code class="sourceCode R"><span><span class="va">m.L2.SFO</span> <span class="op">&lt;-</span> <span class="fu"><a href="../reference/mkinfit.html">mkinfit</a></span><span class="op">(</span><span class="st">"SFO"</span>, <span class="va">FOCUS_2006_L2_mkin</span>, quiet<span class="op">=</span><span class="cn">TRUE</span><span class="op">)</span></span>
-<span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L2.SFO</span>, show_residuals <span class="op">=</span> <span class="cn">TRUE</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span>,</span>
+<span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L2.SFO</span>, show_residuals <span class="op">=</span> <span class="cn">TRUE</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span>,</span>
 <span>     main <span class="op">=</span> <span class="st">"FOCUS L2 - SFO"</span><span class="op">)</span></span></code></pre></div>
 <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-8-1.png" width="672"></p>
 <p>The <span class="math inline">\(\chi^2\)</span> error level of 14%
@@ -422,22 +422,22 @@ kinetics.</p>
 <p>For comparison, the FOMC model is fitted as well, and the <span class="math inline">\(\chi^2\)</span> error level is checked.</p>
 <div class="sourceCode" id="cb16"><pre class="downlit sourceCode r">
 <code class="sourceCode R"><span><span class="va">m.L2.FOMC</span> <span class="op">&lt;-</span> <span class="fu"><a href="../reference/mkinfit.html">mkinfit</a></span><span class="op">(</span><span class="st">"FOMC"</span>, <span class="va">FOCUS_2006_L2_mkin</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span>
-<span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L2.FOMC</span>, show_residuals <span class="op">=</span> <span class="cn">TRUE</span>,</span>
+<span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L2.FOMC</span>, show_residuals <span class="op">=</span> <span class="cn">TRUE</span>,</span>
 <span>     main <span class="op">=</span> <span class="st">"FOCUS L2 - FOMC"</span><span class="op">)</span></span></code></pre></div>
 <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-9-1.png" width="672"></p>
 <div class="sourceCode" id="cb17"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><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">m.L2.FOMC</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
-<pre><code><span><span class="co">## mkin version used for fitting:    1.2.4 </span></span>
-<span><span class="co">## R version used for fitting:       4.3.0 </span></span>
-<span><span class="co">## Date of fit:     Fri May 19 09:20:25 2023 </span></span>
-<span><span class="co">## Date of summary: Fri May 19 09:20:25 2023 </span></span>
+<code class="sourceCode R"><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">m.L2.FOMC</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
+<pre><code><span><span class="co">## mkin version used for fitting:    1.2.5 </span></span>
+<span><span class="co">## R version used for fitting:       4.3.1 </span></span>
+<span><span class="co">## Date of fit:     Wed Aug  9 17:55:40 2023 </span></span>
+<span><span class="co">## Date of summary: Wed Aug  9 17:55:40 2023 </span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Equations:</span></span>
 <span><span class="co">## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Model predictions using solution type analytical </span></span>
 <span><span class="co">## </span></span>
-<span><span class="co">## Fitted using 239 model solutions performed in 0.013 s</span></span>
+<span><span class="co">## Fitted using 239 model solutions performed in 0.044 s</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Error model: Constant variance </span></span>
 <span><span class="co">## </span></span>
@@ -506,15 +506,15 @@ the data.</p>
 <p>Fitting the four parameter DFOP model further reduces the <span class="math inline">\(\chi^2\)</span> error level.</p>
 <div class="sourceCode" id="cb19"><pre class="downlit sourceCode r">
 <code class="sourceCode R"><span><span class="va">m.L2.DFOP</span> <span class="op">&lt;-</span> <span class="fu"><a href="../reference/mkinfit.html">mkinfit</a></span><span class="op">(</span><span class="st">"DFOP"</span>, <span class="va">FOCUS_2006_L2_mkin</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span>
-<span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L2.DFOP</span>, show_residuals <span class="op">=</span> <span class="cn">TRUE</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span>,</span>
+<span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">m.L2.DFOP</span>, show_residuals <span class="op">=</span> <span class="cn">TRUE</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span>,</span>
 <span>     main <span class="op">=</span> <span class="st">"FOCUS L2 - DFOP"</span><span class="op">)</span></span></code></pre></div>
 <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-10-1.png" width="672"></p>
 <div class="sourceCode" id="cb20"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><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">m.L2.DFOP</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
-<pre><code><span><span class="co">## mkin version used for fitting:    1.2.4 </span></span>
-<span><span class="co">## R version used for fitting:       4.3.0 </span></span>
-<span><span class="co">## Date of fit:     Fri May 19 09:20:25 2023 </span></span>
-<span><span class="co">## Date of summary: Fri May 19 09:20:25 2023 </span></span>
+<code class="sourceCode R"><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">m.L2.DFOP</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
+<pre><code><span><span class="co">## mkin version used for fitting:    1.2.5 </span></span>
+<span><span class="co">## R version used for fitting:       4.3.1 </span></span>
+<span><span class="co">## Date of fit:     Wed Aug  9 17:55:40 2023 </span></span>
+<span><span class="co">## Date of summary: Wed Aug  9 17:55:40 2023 </span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Equations:</span></span>
 <span><span class="co">## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *</span></span>
@@ -523,7 +523,7 @@ the data.</p>
 <span><span class="co">## </span></span>
 <span><span class="co">## Model predictions using solution type analytical </span></span>
 <span><span class="co">## </span></span>
-<span><span class="co">## Fitted using 581 model solutions performed in 0.038 s</span></span>
+<span><span class="co">## Fitted using 581 model solutions performed in 0.119 s</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Error model: Constant variance </span></span>
 <span><span class="co">## </span></span>
@@ -611,7 +611,7 @@ only the L3 dataset prepared above.</p>
 <code class="sourceCode R"><span><span class="co"># Only use one core here, not to offend the CRAN checks</span></span>
 <span><span class="va">mm.L3</span> <span class="op">&lt;-</span> <span class="fu"><a href="../reference/mmkin.html">mmkin</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="st">"SFO"</span>, <span class="st">"FOMC"</span>, <span class="st">"DFOP"</span><span class="op">)</span>, cores <span class="op">=</span> <span class="fl">1</span>,</span>
 <span>               <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="st">"FOCUS L3"</span> <span class="op">=</span> <span class="va">FOCUS_2006_L3_mkin</span><span class="op">)</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span>
-<span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mm.L3</span><span class="op">)</span></span></code></pre></div>
+<span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mm.L3</span><span class="op">)</span></span></code></pre></div>
 <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-12-1.png" width="700"></p>
 <p>The <span class="math inline">\(\chi^2\)</span> error level of 21% as
 well as the plot suggest that the SFO model does not fit very well. The
@@ -627,11 +627,11 @@ as a row index and datasets as a column index.</p>
 using square brackets for indexing which will result in the use of the
 summary and plot functions working on mkinfit objects.</p>
 <div class="sourceCode" id="cb24"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><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">mm.L3</span><span class="op">[[</span><span class="st">"DFOP"</span>, <span class="fl">1</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div>
-<pre><code><span><span class="co">## mkin version used for fitting:    1.2.4 </span></span>
-<span><span class="co">## R version used for fitting:       4.3.0 </span></span>
-<span><span class="co">## Date of fit:     Fri May 19 09:20:26 2023 </span></span>
-<span><span class="co">## Date of summary: Fri May 19 09:20:26 2023 </span></span>
+<code class="sourceCode R"><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">mm.L3</span><span class="op">[[</span><span class="st">"DFOP"</span>, <span class="fl">1</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div>
+<pre><code><span><span class="co">## mkin version used for fitting:    1.2.5 </span></span>
+<span><span class="co">## R version used for fitting:       4.3.1 </span></span>
+<span><span class="co">## Date of fit:     Wed Aug  9 17:55:41 2023 </span></span>
+<span><span class="co">## Date of summary: Wed Aug  9 17:55:41 2023 </span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Equations:</span></span>
 <span><span class="co">## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *</span></span>
@@ -640,7 +640,7 @@ summary and plot functions working on mkinfit objects.</p>
 <span><span class="co">## </span></span>
 <span><span class="co">## Model predictions using solution type analytical </span></span>
 <span><span class="co">## </span></span>
-<span><span class="co">## Fitted using 376 model solutions performed in 0.021 s</span></span>
+<span><span class="co">## Fitted using 376 model solutions performed in 0.075 s</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Error model: Constant variance </span></span>
 <span><span class="co">## </span></span>
@@ -715,7 +715,7 @@ summary and plot functions working on mkinfit objects.</p>
 <span><span class="co">##    91   parent     15.0     15.18 -0.18181</span></span>
 <span><span class="co">##   120   parent     12.0     10.19  1.81395</span></span></code></pre>
 <div class="sourceCode" id="cb26"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mm.L3</span><span class="op">[[</span><span class="st">"DFOP"</span>, <span class="fl">1</span><span class="op">]</span><span class="op">]</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
+<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mm.L3</span><span class="op">[[</span><span class="st">"DFOP"</span>, <span class="fl">1</span><span class="op">]</span><span class="op">]</span>, show_errmin <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
 <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-13-1.png" width="700"></p>
 <p>Here, a look to the model plot, the confidence intervals of the
 parameters and the correlation matrix suggest that the parameter
@@ -746,7 +746,7 @@ below:</p>
 <span><span class="va">mm.L4</span> <span class="op">&lt;-</span> <span class="fu"><a href="../reference/mmkin.html">mmkin</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="st">"SFO"</span>, <span class="st">"FOMC"</span><span class="op">)</span>, cores <span class="op">=</span> <span class="fl">1</span>,</span>
 <span>               <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="st">"FOCUS L4"</span> <span class="op">=</span> <span class="va">FOCUS_2006_L4_mkin</span><span class="op">)</span>,</span>
 <span>               quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span>
-<span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mm.L4</span><span class="op">)</span></span></code></pre></div>
+<span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">mm.L4</span><span class="op">)</span></span></code></pre></div>
 <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-15-1.png" width="700"></p>
 <p>The <span class="math inline">\(\chi^2\)</span> error level of 3.3%
 as well as the plot suggest that the SFO model fits very well. The error
@@ -754,18 +754,18 @@ level at which the <span class="math inline">\(\chi^2\)</span> test
 passes is slightly lower for the FOMC model. However, the difference
 appears negligible.</p>
 <div class="sourceCode" id="cb29"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><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">mm.L4</span><span class="op">[[</span><span class="st">"SFO"</span>, <span class="fl">1</span><span class="op">]</span><span class="op">]</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
-<pre><code><span><span class="co">## mkin version used for fitting:    1.2.4 </span></span>
-<span><span class="co">## R version used for fitting:       4.3.0 </span></span>
-<span><span class="co">## Date of fit:     Fri May 19 09:20:26 2023 </span></span>
-<span><span class="co">## Date of summary: Fri May 19 09:20:26 2023 </span></span>
+<code class="sourceCode R"><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">mm.L4</span><span class="op">[[</span><span class="st">"SFO"</span>, <span class="fl">1</span><span class="op">]</span><span class="op">]</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
+<pre><code><span><span class="co">## mkin version used for fitting:    1.2.5 </span></span>
+<span><span class="co">## R version used for fitting:       4.3.1 </span></span>
+<span><span class="co">## Date of fit:     Wed Aug  9 17:55:42 2023 </span></span>
+<span><span class="co">## Date of summary: Wed Aug  9 17:55:42 2023 </span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Equations:</span></span>
 <span><span class="co">## d_parent/dt = - k_parent * parent</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Model predictions using solution type analytical </span></span>
 <span><span class="co">## </span></span>
-<span><span class="co">## Fitted using 142 model solutions performed in 0.008 s</span></span>
+<span><span class="co">## Fitted using 142 model solutions performed in 0.027 s</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Error model: Constant variance </span></span>
 <span><span class="co">## </span></span>
@@ -819,18 +819,18 @@ appears negligible.</p>
 <span><span class="co">##        DT50 DT90</span></span>
 <span><span class="co">## parent  106  352</span></span></code></pre>
 <div class="sourceCode" id="cb31"><pre class="downlit sourceCode r">
-<code class="sourceCode R"><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">mm.L4</span><span class="op">[[</span><span class="st">"FOMC"</span>, <span class="fl">1</span><span class="op">]</span><span class="op">]</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
-<pre><code><span><span class="co">## mkin version used for fitting:    1.2.4 </span></span>
-<span><span class="co">## R version used for fitting:       4.3.0 </span></span>
-<span><span class="co">## Date of fit:     Fri May 19 09:20:26 2023 </span></span>
-<span><span class="co">## Date of summary: Fri May 19 09:20:26 2023 </span></span>
+<code class="sourceCode R"><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">mm.L4</span><span class="op">[[</span><span class="st">"FOMC"</span>, <span class="fl">1</span><span class="op">]</span><span class="op">]</span>, data <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></code></pre></div>
+<pre><code><span><span class="co">## mkin version used for fitting:    1.2.5 </span></span>
+<span><span class="co">## R version used for fitting:       4.3.1 </span></span>
+<span><span class="co">## Date of fit:     Wed Aug  9 17:55:42 2023 </span></span>
+<span><span class="co">## Date of summary: Wed Aug  9 17:55:42 2023 </span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Equations:</span></span>
 <span><span class="co">## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Model predictions using solution type analytical </span></span>
 <span><span class="co">## </span></span>
-<span><span class="co">## Fitted using 224 model solutions performed in 0.012 s</span></span>
+<span><span class="co">## Fitted using 224 model solutions performed in 0.04 s</span></span>
 <span><span class="co">## </span></span>
 <span><span class="co">## Error model: Constant variance </span></span>
 <span><span class="co">## </span></span>
diff --git a/docs/articles/index.html b/docs/articles/index.html
index 7060b377..ac6236f1 100644
--- a/docs/articles/index.html
+++ b/docs/articles/index.html
@@ -17,7 +17,7 @@
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="../index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
@@ -119,6 +119,8 @@
         <dd>
         </dd><dt><a href="prebuilt/2022_dmta_pathway.html">Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</a></dt>
         <dd>
+        </dd><dt><a href="prebuilt/2023_mesotrione_parent.html">Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione</a></dt>
+        <dd>
         </dd><dt><a href="twa.html">Calculation of time weighted average concentrations with mkin</a></dt>
         <dd>
         </dd><dt><a href="web_only/FOCUS_Z.html">Example evaluation of FOCUS dataset Z</a></dt>
diff --git a/docs/articles/mkin.html b/docs/articles/mkin.html
index 6213916b..3bd0d4d3 100644
--- a/docs/articles/mkin.html
+++ b/docs/articles/mkin.html
@@ -33,7 +33,7 @@
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="../index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
@@ -134,7 +134,7 @@
 Ranke</h4>
             
             <h4 data-toc-skip class="date">Last change 18 May 2023
-(rebuilt 2023-05-19)</h4>
+(rebuilt 2023-08-09)</h4>
       
       <small class="dont-index">Source: <a href="https://github.com/jranke/mkin/blob/HEAD/vignettes/mkin.rmd" class="external-link"><code>vignettes/mkin.rmd</code></a></small>
       <div class="hidden name"><code>mkin.rmd</code></div>
diff --git a/docs/articles/prebuilt/2023_mesotrione_parent.html b/docs/articles/prebuilt/2023_mesotrione_parent.html
new file mode 100644
index 00000000..b233fc3c
--- /dev/null
+++ b/docs/articles/prebuilt/2023_mesotrione_parent.html
@@ -0,0 +1,2560 @@
+<!DOCTYPE html>
+<!-- Generated by pkgdown: do not edit by hand --><html lang="en">
+<head>
+<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
+<meta charset="utf-8">
+<meta http-equiv="X-UA-Compatible" content="IE=edge">
+<meta name="viewport" content="width=device-width, initial-scale=1.0">
+<title>Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione • mkin</title>
+<!-- jquery --><script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.4.1/jquery.min.js" integrity="sha256-CSXorXvZcTkaix6Yvo6HppcZGetbYMGWSFlBw8HfCJo=" crossorigin="anonymous"></script><!-- Bootstrap --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/css/bootstrap.min.css" integrity="sha256-bZLfwXAP04zRMK2BjiO8iu9pf4FbLqX6zitd+tIvLhE=" crossorigin="anonymous">
+<script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/js/bootstrap.min.js" integrity="sha256-nuL8/2cJ5NDSSwnKD8VqreErSWHtnEP9E7AySL+1ev4=" crossorigin="anonymous"></script><!-- bootstrap-toc --><link rel="stylesheet" href="../../bootstrap-toc.css">
+<script src="../../bootstrap-toc.js"></script><!-- Font Awesome icons --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/all.min.css" integrity="sha256-mmgLkCYLUQbXn0B1SRqzHar6dCnv9oZFPEC1g1cwlkk=" crossorigin="anonymous">
+<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous">
+<!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- pkgdown --><link href="../../pkgdown.css" rel="stylesheet">
+<script src="../../pkgdown.js"></script><meta property="og:title" content="Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione">
+<meta property="og:description" content="mkin">
+<!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]>
+<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
+<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
+<![endif]-->
+</head>
+<body data-spy="scroll" data-target="#toc">
+    
+
+    <div class="container template-article">
+      <header><div class="navbar navbar-default navbar-fixed-top" role="navigation">
+  <div class="container">
+    <div class="navbar-header">
+      <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar" aria-expanded="false">
+        <span class="sr-only">Toggle navigation</span>
+        <span class="icon-bar"></span>
+        <span class="icon-bar"></span>
+        <span class="icon-bar"></span>
+      </button>
+      <span class="navbar-brand">
+        <a class="navbar-link" href="../../index.html">mkin</a>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
+      </span>
+    </div>
+
+    <div id="navbar" class="navbar-collapse collapse">
+      <ul class="nav navbar-nav">
+<li>
+  <a href="../../reference/index.html">Reference</a>
+</li>
+<li class="dropdown">
+  <a href="#" class="dropdown-toggle" data-toggle="dropdown" role="button" data-bs-toggle="dropdown" aria-expanded="false">
+    Articles
+     
+    <span class="caret"></span>
+  </a>
+  <ul class="dropdown-menu" role="menu">
+<li>
+      <a href="../../articles/mkin.html">Introduction to mkin</a>
+    </li>
+    <li class="divider">
+    </li>
+<li class="dropdown-header">Example evaluations with (generalised) nonlinear least squares</li>
+    <li>
+      <a href="../../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a>
+    </li>
+    <li>
+      <a href="../../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a>
+    </li>
+    <li>
+      <a href="../../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a>
+    </li>
+    <li class="divider">
+    </li>
+<li class="dropdown-header">Example evaluations with hierarchical models (nonlinear mixed-effects models)</li>
+    <li>
+      <a href="../../articles/prebuilt/2022_dmta_parent.html">Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P</a>
+    </li>
+    <li>
+      <a href="../../articles/prebuilt/2022_dmta_pathway.html">Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</a>
+    </li>
+    <li>
+      <a href="../../articles/prebuilt/2022_cyan_pathway.html">Testing hierarchical pathway kinetics with residue data on cyantraniliprole</a>
+    </li>
+    <li>
+      <a href="../../articles/web_only/dimethenamid_2018.html">Comparison of saemix and nlme evaluations of dimethenamid data from 2018</a>
+    </li>
+    <li>
+      <a href="../../articles/web_only/multistart.html">Short demo of the multistart method</a>
+    </li>
+    <li class="divider">
+    </li>
+<li class="dropdown-header">Performance</li>
+    <li>
+      <a href="../../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a>
+    </li>
+    <li>
+      <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>
+    <li class="divider">
+    </li>
+<li class="dropdown-header">Miscellaneous</li>
+    <li>
+      <a href="../../articles/twa.html">Calculation of time weighted average concentrations with mkin</a>
+    </li>
+    <li>
+      <a href="../../articles/web_only/NAFTA_examples.html">Example evaluation of NAFTA SOP Attachment examples</a>
+    </li>
+  </ul>
+</li>
+<li>
+  <a href="../../news/index.html">News</a>
+</li>
+      </ul>
+<ul class="nav navbar-nav navbar-right">
+<li>
+  <a href="https://github.com/jranke/mkin/" class="external-link">
+    <span class="fab fa-github fa-lg"></span>
+     
+  </a>
+</li>
+      </ul>
+</div>
+<!--/.nav-collapse -->
+  </div>
+<!--/.container -->
+</div>
+<!--/.navbar -->
+
+      
+
+      </header><div class="row">
+  <div class="col-md-9 contents">
+    <div class="page-header toc-ignore">
+      <h1 data-toc-skip>Testing covariate modelling in hierarchical
+parent degradation kinetics with residue data on mesotrione</h1>
+                        <h4 data-toc-skip class="author">Johannes
+Ranke</h4>
+            
+            <h4 data-toc-skip class="date">Last change on 4 August 2023,
+last compiled on 9 August 2023</h4>
+      
+      <small class="dont-index">Source: <a href="https://github.com/jranke/mkin/blob/HEAD/vignettes/prebuilt/2023_mesotrione_parent.rmd" class="external-link"><code>vignettes/prebuilt/2023_mesotrione_parent.rmd</code></a></small>
+      <div class="hidden name"><code>2023_mesotrione_parent.rmd</code></div>
+
+    </div>
+
+    
+    
+<div class="section level2">
+<h2 id="introduction">Introduction<a class="anchor" aria-label="anchor" href="#introduction"></a>
+</h2>
+<p>The purpose of this document is to test demonstrate how nonlinear
+hierarchical models (NLHM) based on the parent degradation models SFO,
+FOMC, DFOP and HS can be fitted with the mkin package, also considering
+the influence of covariates like soil pH on different degradation
+parameters. Because in some other case studies, the SFORB
+parameterisation of biexponential decline has shown some advantages over
+the DFOP parameterisation, SFORB was included in the list of tested
+models as well.</p>
+<p>The mkin package is used in version 1.2.5, which is contains the
+functions that were used for the evaluations. The <code>saemix</code>
+package is used as a backend for fitting the NLHM, but is also loaded to
+make the convergence plot function available.</p>
+<p>This document is processed with the <code>knitr</code> package, which
+also provides the <code>kable</code> function that is used to improve
+the display of tabular data in R markdown documents. For parallel
+processing, the <code>parallel</code> package is used.</p>
+<div class="sourceCode" id="cb1"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="kw"><a href="https://rdrr.io/r/base/library.html" class="external-link">library</a></span><span class="op">(</span><span class="va"><a href="https://pkgdown.jrwb.de/mkin/">mkin</a></span><span class="op">)</span></span>
+<span><span class="kw"><a href="https://rdrr.io/r/base/library.html" class="external-link">library</a></span><span class="op">(</span><span class="va"><a href="https://yihui.org/knitr/" class="external-link">knitr</a></span><span class="op">)</span></span>
+<span><span class="kw"><a href="https://rdrr.io/r/base/library.html" class="external-link">library</a></span><span class="op">(</span><span class="va">saemix</span><span class="op">)</span></span>
+<span><span class="kw"><a href="https://rdrr.io/r/base/library.html" class="external-link">library</a></span><span class="op">(</span><span class="va">parallel</span><span class="op">)</span></span>
+<span><span class="va">n_cores</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/parallel/detectCores.html" class="external-link">detectCores</a></span><span class="op">(</span><span class="op">)</span></span>
+<span><span class="kw">if</span> <span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/Sys.info.html" class="external-link">Sys.info</a></span><span class="op">(</span><span class="op">)</span><span class="op">[</span><span class="st">"sysname"</span><span class="op">]</span> <span class="op">==</span> <span class="st">"Windows"</span><span class="op">)</span> <span class="op">{</span></span>
+<span>  <span class="va">cl</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/parallel/makeCluster.html" class="external-link">makePSOCKcluster</a></span><span class="op">(</span><span class="va">n_cores</span><span class="op">)</span></span>
+<span><span class="op">}</span> <span class="kw">else</span> <span class="op">{</span></span>
+<span>  <span class="va">cl</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/parallel/makeCluster.html" class="external-link">makeForkCluster</a></span><span class="op">(</span><span class="va">n_cores</span><span class="op">)</span></span>
+<span><span class="op">}</span></span></code></pre></div>
+<div class="section level3">
+<h3 id="test-data">Test data<a class="anchor" aria-label="anchor" href="#test-data"></a>
+</h3>
+<div class="sourceCode" id="cb2"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">data_file</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/system.file.html" class="external-link">system.file</a></span><span class="op">(</span></span>
+<span>  <span class="st">"testdata"</span>, <span class="st">"mesotrione_soil_efsa_2016.xlsx"</span>, package <span class="op">=</span> <span class="st">"mkin"</span><span class="op">)</span></span>
+<span><span class="va">meso_ds</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/read_spreadsheet.html">read_spreadsheet</a></span><span class="op">(</span><span class="va">data_file</span>, parent_only <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
+<p>The following tables show the covariate data and the 18 datasets that
+were read in from the spreadsheet file.</p>
+<div class="sourceCode" id="cb3"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">pH</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/attr.html" class="external-link">attr</a></span><span class="op">(</span><span class="va">meso_ds</span>, <span class="st">"covariates"</span><span class="op">)</span></span>
+<span><span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="va">pH</span>, caption <span class="op">=</span> <span class="st">"Covariate data"</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<caption>Covariate data</caption>
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">pH</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">Richmond</td>
+<td align="right">6.2</td>
+</tr>
+<tr class="even">
+<td align="left">Richmond 2</td>
+<td align="right">6.2</td>
+</tr>
+<tr class="odd">
+<td align="left">ERTC</td>
+<td align="right">6.4</td>
+</tr>
+<tr class="even">
+<td align="left">Toulouse</td>
+<td align="right">7.7</td>
+</tr>
+<tr class="odd">
+<td align="left">Picket Piece</td>
+<td align="right">7.1</td>
+</tr>
+<tr class="even">
+<td align="left">721</td>
+<td align="right">5.6</td>
+</tr>
+<tr class="odd">
+<td align="left">722</td>
+<td align="right">5.7</td>
+</tr>
+<tr class="even">
+<td align="left">723</td>
+<td align="right">5.4</td>
+</tr>
+<tr class="odd">
+<td align="left">724</td>
+<td align="right">4.8</td>
+</tr>
+<tr class="even">
+<td align="left">725</td>
+<td align="right">5.8</td>
+</tr>
+<tr class="odd">
+<td align="left">727</td>
+<td align="right">5.1</td>
+</tr>
+<tr class="even">
+<td align="left">728</td>
+<td align="right">5.9</td>
+</tr>
+<tr class="odd">
+<td align="left">729</td>
+<td align="right">5.6</td>
+</tr>
+<tr class="even">
+<td align="left">730</td>
+<td align="right">5.3</td>
+</tr>
+<tr class="odd">
+<td align="left">731</td>
+<td align="right">6.1</td>
+</tr>
+<tr class="even">
+<td align="left">732</td>
+<td align="right">5.0</td>
+</tr>
+<tr class="odd">
+<td align="left">741</td>
+<td align="right">5.7</td>
+</tr>
+<tr class="even">
+<td align="left">742</td>
+<td align="right">7.2</td>
+</tr>
+</tbody>
+</table>
+<div class="sourceCode" id="cb4"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="kw">for</span> <span class="op">(</span><span class="va">ds_name</span> <span class="kw">in</span> <span class="fu"><a href="https://rdrr.io/r/base/names.html" class="external-link">names</a></span><span class="op">(</span><span class="va">meso_ds</span><span class="op">)</span><span class="op">)</span> <span class="op">{</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>
+<span>    <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="fu"><a href="../../reference/mkin_long_to_wide.html">mkin_long_to_wide</a></span><span class="op">(</span><span class="va">meso_ds</span><span class="op">[[</span><span class="va">ds_name</span><span class="op">]</span><span class="op">]</span><span class="op">)</span>,</span>
+<span>      caption <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/paste.html" class="external-link">paste</a></span><span class="op">(</span><span class="st">"Dataset"</span>, <span class="va">ds_name</span><span class="op">)</span>,</span>
+<span>      booktabs <span class="op">=</span> <span class="cn">TRUE</span>, row.names <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span><span class="op">)</span></span>
+<span><span class="op">}</span></span></code></pre></div>
+<table class="table">
+<caption>Dataset Richmond</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.000000</td>
+<td align="right">91.00</td>
+</tr>
+<tr class="even">
+<td align="right">1.179050</td>
+<td align="right">86.70</td>
+</tr>
+<tr class="odd">
+<td align="right">3.537149</td>
+<td align="right">73.60</td>
+</tr>
+<tr class="even">
+<td align="right">7.074299</td>
+<td align="right">61.50</td>
+</tr>
+<tr class="odd">
+<td align="right">10.611448</td>
+<td align="right">55.70</td>
+</tr>
+<tr class="even">
+<td align="right">15.327647</td>
+<td align="right">47.70</td>
+</tr>
+<tr class="odd">
+<td align="right">17.685747</td>
+<td align="right">39.50</td>
+</tr>
+<tr class="even">
+<td align="right">24.760046</td>
+<td align="right">29.80</td>
+</tr>
+<tr class="odd">
+<td align="right">35.371494</td>
+<td align="right">19.60</td>
+</tr>
+<tr class="even">
+<td align="right">68.384889</td>
+<td align="right">5.67</td>
+</tr>
+<tr class="odd">
+<td align="right">0.000000</td>
+<td align="right">97.90</td>
+</tr>
+<tr class="even">
+<td align="right">1.179050</td>
+<td align="right">96.40</td>
+</tr>
+<tr class="odd">
+<td align="right">3.537149</td>
+<td align="right">89.10</td>
+</tr>
+<tr class="even">
+<td align="right">7.074299</td>
+<td align="right">74.40</td>
+</tr>
+<tr class="odd">
+<td align="right">10.611448</td>
+<td align="right">57.40</td>
+</tr>
+<tr class="even">
+<td align="right">15.327647</td>
+<td align="right">46.30</td>
+</tr>
+<tr class="odd">
+<td align="right">18.864797</td>
+<td align="right">35.50</td>
+</tr>
+<tr class="even">
+<td align="right">27.118146</td>
+<td align="right">27.20</td>
+</tr>
+<tr class="odd">
+<td align="right">35.371494</td>
+<td align="right">19.10</td>
+</tr>
+<tr class="even">
+<td align="right">74.280138</td>
+<td align="right">6.50</td>
+</tr>
+<tr class="odd">
+<td align="right">108.472582</td>
+<td align="right">3.40</td>
+</tr>
+<tr class="even">
+<td align="right">142.665027</td>
+<td align="right">2.20</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset Richmond 2</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.000000</td>
+<td align="right">96.0</td>
+</tr>
+<tr class="even">
+<td align="right">2.422004</td>
+<td align="right">82.4</td>
+</tr>
+<tr class="odd">
+<td align="right">5.651343</td>
+<td align="right">71.2</td>
+</tr>
+<tr class="even">
+<td align="right">8.073348</td>
+<td align="right">53.1</td>
+</tr>
+<tr class="odd">
+<td align="right">11.302687</td>
+<td align="right">48.5</td>
+</tr>
+<tr class="even">
+<td align="right">16.954030</td>
+<td align="right">33.4</td>
+</tr>
+<tr class="odd">
+<td align="right">22.605373</td>
+<td align="right">24.2</td>
+</tr>
+<tr class="even">
+<td align="right">45.210746</td>
+<td align="right">11.9</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset ERTC</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.000000</td>
+<td align="right">99.9</td>
+</tr>
+<tr class="even">
+<td align="right">2.755193</td>
+<td align="right">80.0</td>
+</tr>
+<tr class="odd">
+<td align="right">6.428782</td>
+<td align="right">42.1</td>
+</tr>
+<tr class="even">
+<td align="right">9.183975</td>
+<td align="right">50.1</td>
+</tr>
+<tr class="odd">
+<td align="right">12.857565</td>
+<td align="right">28.4</td>
+</tr>
+<tr class="even">
+<td align="right">19.286347</td>
+<td align="right">39.8</td>
+</tr>
+<tr class="odd">
+<td align="right">25.715130</td>
+<td align="right">29.9</td>
+</tr>
+<tr class="even">
+<td align="right">51.430259</td>
+<td align="right">2.5</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset Toulouse</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.000000</td>
+<td align="right">96.8</td>
+</tr>
+<tr class="even">
+<td align="right">2.897983</td>
+<td align="right">63.3</td>
+</tr>
+<tr class="odd">
+<td align="right">6.761960</td>
+<td align="right">22.3</td>
+</tr>
+<tr class="even">
+<td align="right">9.659942</td>
+<td align="right">16.6</td>
+</tr>
+<tr class="odd">
+<td align="right">13.523919</td>
+<td align="right">16.1</td>
+</tr>
+<tr class="even">
+<td align="right">20.285879</td>
+<td align="right">17.2</td>
+</tr>
+<tr class="odd">
+<td align="right">27.047838</td>
+<td align="right">1.8</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset Picket Piece</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.000000</td>
+<td align="right">102.0</td>
+</tr>
+<tr class="even">
+<td align="right">2.841195</td>
+<td align="right">73.7</td>
+</tr>
+<tr class="odd">
+<td align="right">6.629454</td>
+<td align="right">35.5</td>
+</tr>
+<tr class="even">
+<td align="right">9.470649</td>
+<td align="right">31.8</td>
+</tr>
+<tr class="odd">
+<td align="right">13.258909</td>
+<td align="right">18.0</td>
+</tr>
+<tr class="even">
+<td align="right">19.888364</td>
+<td align="right">3.7</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 721</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">86.4</td>
+</tr>
+<tr class="even">
+<td align="right">11.24366</td>
+<td align="right">61.4</td>
+</tr>
+<tr class="odd">
+<td align="right">22.48733</td>
+<td align="right">49.8</td>
+</tr>
+<tr class="even">
+<td align="right">33.73099</td>
+<td align="right">41.0</td>
+</tr>
+<tr class="odd">
+<td align="right">44.97466</td>
+<td align="right">35.1</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 722</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">90.3</td>
+</tr>
+<tr class="even">
+<td align="right">11.24366</td>
+<td align="right">52.1</td>
+</tr>
+<tr class="odd">
+<td align="right">22.48733</td>
+<td align="right">37.4</td>
+</tr>
+<tr class="even">
+<td align="right">33.73099</td>
+<td align="right">21.2</td>
+</tr>
+<tr class="odd">
+<td align="right">44.97466</td>
+<td align="right">14.3</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 723</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">89.3</td>
+</tr>
+<tr class="even">
+<td align="right">11.24366</td>
+<td align="right">70.8</td>
+</tr>
+<tr class="odd">
+<td align="right">22.48733</td>
+<td align="right">51.1</td>
+</tr>
+<tr class="even">
+<td align="right">33.73099</td>
+<td align="right">42.7</td>
+</tr>
+<tr class="odd">
+<td align="right">44.97466</td>
+<td align="right">26.7</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 724</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.000000</td>
+<td align="right">89.4</td>
+</tr>
+<tr class="even">
+<td align="right">9.008208</td>
+<td align="right">65.2</td>
+</tr>
+<tr class="odd">
+<td align="right">18.016415</td>
+<td align="right">55.8</td>
+</tr>
+<tr class="even">
+<td align="right">27.024623</td>
+<td align="right">46.0</td>
+</tr>
+<tr class="odd">
+<td align="right">36.032831</td>
+<td align="right">41.7</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 725</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">89.0</td>
+</tr>
+<tr class="even">
+<td align="right">10.99058</td>
+<td align="right">35.4</td>
+</tr>
+<tr class="odd">
+<td align="right">21.98116</td>
+<td align="right">18.6</td>
+</tr>
+<tr class="even">
+<td align="right">32.97174</td>
+<td align="right">11.6</td>
+</tr>
+<tr class="odd">
+<td align="right">43.96232</td>
+<td align="right">7.6</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 727</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">91.3</td>
+</tr>
+<tr class="even">
+<td align="right">10.96104</td>
+<td align="right">63.2</td>
+</tr>
+<tr class="odd">
+<td align="right">21.92209</td>
+<td align="right">51.1</td>
+</tr>
+<tr class="even">
+<td align="right">32.88313</td>
+<td align="right">42.0</td>
+</tr>
+<tr class="odd">
+<td align="right">43.84417</td>
+<td align="right">40.8</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 728</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">91.8</td>
+</tr>
+<tr class="even">
+<td align="right">11.24366</td>
+<td align="right">43.6</td>
+</tr>
+<tr class="odd">
+<td align="right">22.48733</td>
+<td align="right">22.0</td>
+</tr>
+<tr class="even">
+<td align="right">33.73099</td>
+<td align="right">15.9</td>
+</tr>
+<tr class="odd">
+<td align="right">44.97466</td>
+<td align="right">8.8</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 729</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">91.6</td>
+</tr>
+<tr class="even">
+<td align="right">11.24366</td>
+<td align="right">60.5</td>
+</tr>
+<tr class="odd">
+<td align="right">22.48733</td>
+<td align="right">43.5</td>
+</tr>
+<tr class="even">
+<td align="right">33.73099</td>
+<td align="right">28.4</td>
+</tr>
+<tr class="odd">
+<td align="right">44.97466</td>
+<td align="right">20.5</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 730</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">92.7</td>
+</tr>
+<tr class="even">
+<td align="right">11.07446</td>
+<td align="right">58.9</td>
+</tr>
+<tr class="odd">
+<td align="right">22.14893</td>
+<td align="right">44.0</td>
+</tr>
+<tr class="even">
+<td align="right">33.22339</td>
+<td align="right">46.0</td>
+</tr>
+<tr class="odd">
+<td align="right">44.29785</td>
+<td align="right">29.3</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 731</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">92.1</td>
+</tr>
+<tr class="even">
+<td align="right">11.24366</td>
+<td align="right">64.4</td>
+</tr>
+<tr class="odd">
+<td align="right">22.48733</td>
+<td align="right">45.3</td>
+</tr>
+<tr class="even">
+<td align="right">33.73099</td>
+<td align="right">33.6</td>
+</tr>
+<tr class="odd">
+<td align="right">44.97466</td>
+<td align="right">23.5</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 732</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">90.3</td>
+</tr>
+<tr class="even">
+<td align="right">11.24366</td>
+<td align="right">58.2</td>
+</tr>
+<tr class="odd">
+<td align="right">22.48733</td>
+<td align="right">40.1</td>
+</tr>
+<tr class="even">
+<td align="right">33.73099</td>
+<td align="right">33.1</td>
+</tr>
+<tr class="odd">
+<td align="right">44.97466</td>
+<td align="right">25.8</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 741</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">90.3</td>
+</tr>
+<tr class="even">
+<td align="right">10.84712</td>
+<td align="right">68.7</td>
+</tr>
+<tr class="odd">
+<td align="right">21.69424</td>
+<td align="right">58.0</td>
+</tr>
+<tr class="even">
+<td align="right">32.54136</td>
+<td align="right">52.2</td>
+</tr>
+<tr class="odd">
+<td align="right">43.38848</td>
+<td align="right">48.0</td>
+</tr>
+</tbody>
+</table>
+<table class="table">
+<caption>Dataset 742</caption>
+<thead><tr class="header">
+<th align="right">time</th>
+<th align="right">meso</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="right">0.00000</td>
+<td align="right">92.0</td>
+</tr>
+<tr class="even">
+<td align="right">11.24366</td>
+<td align="right">60.9</td>
+</tr>
+<tr class="odd">
+<td align="right">22.48733</td>
+<td align="right">36.2</td>
+</tr>
+<tr class="even">
+<td align="right">33.73099</td>
+<td align="right">18.3</td>
+</tr>
+<tr class="odd">
+<td align="right">44.97466</td>
+<td align="right">8.7</td>
+</tr>
+</tbody>
+</table>
+</div>
+</div>
+<div class="section level2">
+<h2 id="separate-evaluations">Separate evaluations<a class="anchor" aria-label="anchor" href="#separate-evaluations"></a>
+</h2>
+<p>In order to obtain suitable starting parameters for the NLHM fits,
+separate fits of the five models to the data for each soil are generated
+using the <code>mmkin</code> function from the mkin package. In a first
+step, constant variance is assumed. Convergence is checked with the
+<code>status</code> function.</p>
+<div class="sourceCode" id="cb5"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">deg_mods</span> <span class="op">&lt;-</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">"SFO"</span>, <span class="st">"FOMC"</span>, <span class="st">"DFOP"</span>, <span class="st">"SFORB"</span>, <span class="st">"HS"</span><span class="op">)</span></span>
+<span><span class="va">f_sep_const</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/mmkin.html">mmkin</a></span><span class="op">(</span></span>
+<span>  <span class="va">deg_mods</span>,</span>
+<span>  <span class="va">meso_ds</span>,</span>
+<span>  error_model <span class="op">=</span> <span class="st">"const"</span>,</span>
+<span>  cluster <span class="op">=</span> <span class="va">cl</span>,</span>
+<span>  quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
+<div class="sourceCode" id="cb6"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/status.html">status</a></span><span class="op">(</span><span class="va">f_sep_const</span><span class="op">[</span>, <span class="fl">1</span><span class="op">:</span><span class="fl">5</span><span class="op">]</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="left">Richmond</th>
+<th align="left">Richmond 2</th>
+<th align="left">ERTC</th>
+<th align="left">Toulouse</th>
+<th align="left">Picket Piece</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">HS</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+</tbody>
+</table>
+<div class="sourceCode" id="cb7"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/status.html">status</a></span><span class="op">(</span><span class="va">f_sep_const</span><span class="op">[</span>, <span class="fl">6</span><span class="op">:</span><span class="fl">18</span><span class="op">]</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<colgroup>
+<col width="10%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+</colgroup>
+<thead><tr class="header">
+<th align="left"></th>
+<th align="left">721</th>
+<th align="left">722</th>
+<th align="left">723</th>
+<th align="left">724</th>
+<th align="left">725</th>
+<th align="left">727</th>
+<th align="left">728</th>
+<th align="left">729</th>
+<th align="left">730</th>
+<th align="left">731</th>
+<th align="left">732</th>
+<th align="left">741</th>
+<th align="left">742</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">HS</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+</tbody>
+</table>
+<p>In the tables above, OK indicates convergence and C indicates failure
+to converge. Most separate fits with constant variance converged, with
+the exception of two FOMC fits, one SFORB fit and one HS fit.</p>
+<div class="sourceCode" id="cb8"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">f_sep_tc</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_sep_const</span>, error_model <span class="op">=</span> <span class="st">"tc"</span><span class="op">)</span></span></code></pre></div>
+<div class="sourceCode" id="cb9"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/status.html">status</a></span><span class="op">(</span><span class="va">f_sep_tc</span><span class="op">[</span>, <span class="fl">1</span><span class="op">:</span><span class="fl">5</span><span class="op">]</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="left">Richmond</th>
+<th align="left">Richmond 2</th>
+<th align="left">ERTC</th>
+<th align="left">Toulouse</th>
+<th align="left">Picket Piece</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">HS</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+</tbody>
+</table>
+<div class="sourceCode" id="cb10"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/status.html">status</a></span><span class="op">(</span><span class="va">f_sep_tc</span><span class="op">[</span>, <span class="fl">6</span><span class="op">:</span><span class="fl">18</span><span class="op">]</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<colgroup>
+<col width="10%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+<col width="6%">
+</colgroup>
+<thead><tr class="header">
+<th align="left"></th>
+<th align="left">721</th>
+<th align="left">722</th>
+<th align="left">723</th>
+<th align="left">724</th>
+<th align="left">725</th>
+<th align="left">727</th>
+<th align="left">728</th>
+<th align="left">729</th>
+<th align="left">730</th>
+<th align="left">731</th>
+<th align="left">732</th>
+<th align="left">741</th>
+<th align="left">742</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">HS</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">C</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+</tbody>
+</table>
+<p>With the two-component error model, the set of fits that did not
+converge is larger, with convergence problems appearing for a number of
+non-SFO fits.</p>
+</div>
+<div class="section level2">
+<h2 id="hierarchical-model-fits-without-covariate-effect">Hierarchical model fits without covariate effect<a class="anchor" aria-label="anchor" href="#hierarchical-model-fits-without-covariate-effect"></a>
+</h2>
+<p>The following code fits hierarchical kinetic models for the ten
+combinations of the five different degradation models with the two
+different error models in parallel.</p>
+<div class="sourceCode" id="cb11"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">f_saem_1</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/mhmkin.html">mhmkin</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">f_sep_const</span>, <span class="va">f_sep_tc</span><span class="op">)</span>, cluster <span class="op">=</span> <span class="va">cl</span><span class="op">)</span></span>
+<span><span class="fu"><a href="../../reference/status.html">status</a></span><span class="op">(</span><span class="va">f_saem_1</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="left">const</th>
+<th align="left">tc</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">HS</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+</tbody>
+</table>
+<p>All fits terminate without errors (status OK).</p>
+<div class="sourceCode" id="cb12"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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_1</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">1</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">npar</th>
+<th align="right">AIC</th>
+<th align="right">BIC</th>
+<th align="right">Lik</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO const</td>
+<td align="right">5</td>
+<td align="right">800.0</td>
+<td align="right">804.5</td>
+<td align="right">-395.0</td>
+</tr>
+<tr class="even">
+<td align="left">SFO tc</td>
+<td align="right">6</td>
+<td align="right">801.9</td>
+<td align="right">807.2</td>
+<td align="right">-394.9</td>
+</tr>
+<tr class="odd">
+<td align="left">FOMC const</td>
+<td align="right">7</td>
+<td align="right">787.4</td>
+<td align="right">793.6</td>
+<td align="right">-386.7</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC tc</td>
+<td align="right">8</td>
+<td align="right">788.9</td>
+<td align="right">796.1</td>
+<td align="right">-386.5</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP const</td>
+<td align="right">9</td>
+<td align="right">787.6</td>
+<td align="right">795.6</td>
+<td align="right">-384.8</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB const</td>
+<td align="right">9</td>
+<td align="right">787.4</td>
+<td align="right">795.4</td>
+<td align="right">-384.7</td>
+</tr>
+<tr class="odd">
+<td align="left">HS const</td>
+<td align="right">9</td>
+<td align="right">781.9</td>
+<td align="right">789.9</td>
+<td align="right">-382.0</td>
+</tr>
+<tr class="even">
+<td align="left">DFOP tc</td>
+<td align="right">10</td>
+<td align="right">787.4</td>
+<td align="right">796.3</td>
+<td align="right">-383.7</td>
+</tr>
+<tr class="odd">
+<td align="left">SFORB tc</td>
+<td align="right">10</td>
+<td align="right">795.8</td>
+<td align="right">804.7</td>
+<td align="right">-387.9</td>
+</tr>
+<tr class="even">
+<td align="left">HS tc</td>
+<td align="right">10</td>
+<td align="right">783.7</td>
+<td align="right">792.7</td>
+<td align="right">-381.9</td>
+</tr>
+</tbody>
+</table>
+<p>The model comparisons show that the fits with constant variance are
+consistently preferable to the corresponding fits with two-component
+error for these data. This is confirmed by the fact that the parameter
+<code>b.1</code> (the relative standard deviation in the fits obtained
+with the saemix package), is ill-defined in all fits.</p>
+<div class="sourceCode" id="cb13"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">f_saem_1</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<colgroup>
+<col width="6%">
+<col width="44%">
+<col width="49%">
+</colgroup>
+<thead><tr class="header">
+<th align="left"></th>
+<th align="left">const</th>
+<th align="left">tc</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO</td>
+<td align="left">sd(meso_0)</td>
+<td align="left">sd(meso_0), b.1</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC</td>
+<td align="left">sd(meso_0), sd(log_beta)</td>
+<td align="left">sd(meso_0), sd(log_beta), b.1</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP</td>
+<td align="left">sd(meso_0), sd(log_k1)</td>
+<td align="left">sd(meso_0), sd(g_qlogis), b.1</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB</td>
+<td align="left">sd(meso_free_0), sd(log_k_meso_free_bound)</td>
+<td align="left">sd(meso_free_0), sd(log_k_meso_free_bound), b.1</td>
+</tr>
+<tr class="odd">
+<td align="left">HS</td>
+<td align="left">sd(meso_0)</td>
+<td align="left">sd(meso_0), b.1</td>
+</tr>
+</tbody>
+</table>
+<p>For obtaining fits with only well-defined random effects, we update
+the set of fits, excluding random effects that were ill-defined
+according to the <code>illparms</code> function.</p>
+<div class="sourceCode" id="cb14"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">f_saem_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_1</span>, no_random_effect <span class="op">=</span> <span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">f_saem_1</span><span class="op">)</span><span class="op">)</span></span>
+<span><span class="fu"><a href="../../reference/status.html">status</a></span><span class="op">(</span><span class="va">f_saem_2</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="left">const</th>
+<th align="left">tc</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+<tr class="odd">
+<td align="left">HS</td>
+<td align="left">OK</td>
+<td align="left">OK</td>
+</tr>
+</tbody>
+</table>
+<p>The updated fits terminate without errors.</p>
+<div class="sourceCode" id="cb15"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">f_saem_2</span><span class="op">)</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="left">const</th>
+<th align="left">tc</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">SFO</td>
+<td align="left"></td>
+<td align="left">b.1</td>
+</tr>
+<tr class="even">
+<td align="left">FOMC</td>
+<td align="left"></td>
+<td align="left">b.1</td>
+</tr>
+<tr class="odd">
+<td align="left">DFOP</td>
+<td align="left"></td>
+<td align="left">b.1</td>
+</tr>
+<tr class="even">
+<td align="left">SFORB</td>
+<td align="left"></td>
+<td align="left">b.1</td>
+</tr>
+<tr class="odd">
+<td align="left">HS</td>
+<td align="left"></td>
+<td align="left"></td>
+</tr>
+</tbody>
+</table>
+<p>No ill-defined errors remain in the fits with constant variance.</p>
+</div>
+<div class="section level2">
+<h2 id="hierarchical-model-fits-with-covariate-effect">Hierarchical model fits with covariate effect<a class="anchor" aria-label="anchor" href="#hierarchical-model-fits-with-covariate-effect"></a>
+</h2>
+<p>In the following sections, hierarchical fits including a model for
+the influence of pH on selected degradation parameters are shown for all
+parent models. Constant variance is selected as the error model based on
+the fits without covariate effects. Random effects that were ill-defined
+in the fits without pH influence are excluded. A potential influence of
+the soil pH is only included for parameters with a well-defined random
+effect, because experience has shown that only for such parameters a
+significant pH effect could be found.</p>
+<div class="section level3">
+<h3 id="sfo">SFO<a class="anchor" aria-label="anchor" href="#sfo"></a>
+</h3>
+<div class="sourceCode" id="cb16"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">sfo_pH</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/saem.html">saem</a></span><span class="op">(</span><span class="va">f_sep_const</span><span class="op">[</span><span class="st">"SFO"</span>, <span class="op">]</span>, no_random_effect <span class="op">=</span> <span class="st">"meso_0"</span>, covariates <span class="op">=</span> <span class="va">pH</span>,</span>
+<span>  covariate_models <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">log_k_meso</span> <span class="op">~</span> <span class="va">pH</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<div class="sourceCode" id="cb17"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">sfo_pH</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_0</td>
+<td align="right">91.35</td>
+<td align="right">89.27</td>
+<td align="right">93.43</td>
+</tr>
+<tr class="even">
+<td align="left">log_k_meso</td>
+<td align="right">-6.66</td>
+<td align="right">-7.97</td>
+<td align="right">-5.35</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_k_meso)</td>
+<td align="right">0.59</td>
+<td align="right">0.37</td>
+<td align="right">0.81</td>
+</tr>
+<tr class="even">
+<td align="left">a.1</td>
+<td align="right">5.48</td>
+<td align="right">4.71</td>
+<td align="right">6.24</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.log_k_meso</td>
+<td align="right">0.35</td>
+<td align="right">0.23</td>
+<td align="right">0.47</td>
+</tr>
+</tbody>
+</table>
+<p>The parameter showing the pH influence in the above table is
+<code>beta_pH(log_k_meso)</code>. Its confidence interval does not
+include zero, indicating that the influence of soil pH on the log of the
+degradation rate constant is significantly greater than zero.</p>
+<div class="sourceCode" id="cb18"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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_2</span><span class="op">[[</span><span class="st">"SFO"</span>, <span class="st">"const"</span><span class="op">]</span><span class="op">]</span>, <span class="va">sfo_pH</span>, test <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
+<pre><code>Data: 116 observations of 1 variable(s) grouped in 18 datasets
+
+                           npar    AIC    BIC     Lik  Chisq Df Pr(&gt;Chisq)    
+f_saem_2[["SFO", "const"]]    4 797.56 801.12 -394.78                         
+sfo_pH                        5 783.09 787.54 -386.54 16.473  1  4.934e-05 ***
+---
+Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1</code></pre>
+<p>The comparison with the SFO fit without covariate effect confirms
+that considering the soil pH improves the model, both by comparison of
+AIC and BIC and by the likelihood ratio test.</p>
+<div class="sourceCode" id="cb20"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">sfo_pH</span><span class="op">)</span></span></code></pre></div>
+<p><img src="2023_mesotrione_parent_files/figure-html/unnamed-chunk-8-1.png" width="700" style="display: block; margin: auto;"></p>
+<p>Endpoints for a model with covariates are by default calculated for
+the median of the covariate values. This quantile can be adapted, or a
+specific covariate value can be given as shown below.</p>
+<div class="sourceCode" id="cb21"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">sfo_pH</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+      pH
+50% 5.75
+
+$distimes
+         DT50     DT90
+meso 18.52069 61.52441</code></pre>
+<div class="sourceCode" id="cb23"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">sfo_pH</span>, covariate_quantile <span class="op">=</span> <span class="fl">0.9</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+      pH
+90% 7.13
+
+$distimes
+         DT50     DT90
+meso 8.237019 27.36278</code></pre>
+<div class="sourceCode" id="cb25"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">sfo_pH</span>, covariates <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>pH <span class="op">=</span> <span class="fl">7.0</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+     pH
+User  7
+
+$distimes
+        DT50    DT90
+meso 8.89035 29.5331</code></pre>
+</div>
+<div class="section level3">
+<h3 id="fomc">FOMC<a class="anchor" aria-label="anchor" href="#fomc"></a>
+</h3>
+<div class="sourceCode" id="cb27"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">fomc_pH</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/saem.html">saem</a></span><span class="op">(</span><span class="va">f_sep_const</span><span class="op">[</span><span class="st">"FOMC"</span>, <span class="op">]</span>, no_random_effect <span class="op">=</span> <span class="st">"meso_0"</span>, covariates <span class="op">=</span> <span class="va">pH</span>,</span>
+<span>  covariate_models <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">log_alpha</span> <span class="op">~</span> <span class="va">pH</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<div class="sourceCode" id="cb28"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">fomc_pH</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_0</td>
+<td align="right">92.84</td>
+<td align="right">90.75</td>
+<td align="right">94.93</td>
+</tr>
+<tr class="even">
+<td align="left">log_alpha</td>
+<td align="right">-2.21</td>
+<td align="right">-3.49</td>
+<td align="right">-0.92</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_alpha)</td>
+<td align="right">0.58</td>
+<td align="right">0.37</td>
+<td align="right">0.79</td>
+</tr>
+<tr class="even">
+<td align="left">log_beta</td>
+<td align="right">4.21</td>
+<td align="right">3.44</td>
+<td align="right">4.99</td>
+</tr>
+<tr class="odd">
+<td align="left">a.1</td>
+<td align="right">5.03</td>
+<td align="right">4.32</td>
+<td align="right">5.73</td>
+</tr>
+<tr class="even">
+<td align="left">SD.log_alpha</td>
+<td align="right">0.00</td>
+<td align="right">-23.77</td>
+<td align="right">23.78</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.log_beta</td>
+<td align="right">0.37</td>
+<td align="right">0.01</td>
+<td align="right">0.74</td>
+</tr>
+</tbody>
+</table>
+<p>As in the case of SFO, the confidence interval of the slope parameter
+(here <code>beta_pH(log_alpha)</code>) quantifying the influence of soil
+pH does not include zero, and the model comparison clearly indicates
+that the model with covariate influence is preferable. However, the
+random effect for <code>alpha</code> is not well-defined any more after
+inclusion of the covariate effect (the confidence interval of
+<code>SD.log_alpha</code> includes zero).</p>
+<div class="sourceCode" id="cb29"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">fomc_pH</span><span class="op">)</span></span></code></pre></div>
+<pre><code>[1] "sd(log_alpha)"</code></pre>
+<p>Therefore, the model is updated without this random effect, and no
+ill-defined parameters remain.</p>
+<div class="sourceCode" id="cb31"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">fomc_pH_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">fomc_pH</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">"meso_0"</span>, <span class="st">"log_alpha"</span><span class="op">)</span><span class="op">)</span></span>
+<span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">fomc_pH_2</span><span class="op">)</span></span></code></pre></div>
+<div class="sourceCode" id="cb32"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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_2</span><span class="op">[[</span><span class="st">"FOMC"</span>, <span class="st">"const"</span><span class="op">]</span><span class="op">]</span>, <span class="va">fomc_pH</span>, <span class="va">fomc_pH_2</span>, test <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
+<pre><code>Data: 116 observations of 1 variable(s) grouped in 18 datasets
+
+                            npar    AIC    BIC     Lik  Chisq Df Pr(&gt;Chisq)    
+f_saem_2[["FOMC", "const"]]    5 783.25 787.71 -386.63                         
+fomc_pH_2                      6 767.49 772.83 -377.75 17.762  1  2.503e-05 ***
+fomc_pH                        7 770.07 776.30 -378.04  0.000  1          1    
+---
+Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1</code></pre>
+<p>Model comparison indicates that including pH dependence significantly
+improves the fit, and that the reduced model with covariate influence
+results in the most preferable FOMC fit.</p>
+<div class="sourceCode" id="cb34"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">fomc_pH_2</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_0</td>
+<td align="right">93.05</td>
+<td align="right">90.98</td>
+<td align="right">95.13</td>
+</tr>
+<tr class="even">
+<td align="left">log_alpha</td>
+<td align="right">-2.91</td>
+<td align="right">-4.18</td>
+<td align="right">-1.63</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_alpha)</td>
+<td align="right">0.66</td>
+<td align="right">0.44</td>
+<td align="right">0.87</td>
+</tr>
+<tr class="even">
+<td align="left">log_beta</td>
+<td align="right">3.95</td>
+<td align="right">3.29</td>
+<td align="right">4.62</td>
+</tr>
+<tr class="odd">
+<td align="left">a.1</td>
+<td align="right">4.98</td>
+<td align="right">4.28</td>
+<td align="right">5.68</td>
+</tr>
+<tr class="even">
+<td align="left">SD.log_beta</td>
+<td align="right">0.40</td>
+<td align="right">0.26</td>
+<td align="right">0.54</td>
+</tr>
+</tbody>
+</table>
+<div class="sourceCode" id="cb35"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">fomc_pH_2</span><span class="op">)</span></span></code></pre></div>
+<p><img src="2023_mesotrione_parent_files/figure-html/unnamed-chunk-14-1.png" width="700" style="display: block; margin: auto;"></p>
+<div class="sourceCode" id="cb36"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">fomc_pH_2</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+      pH
+50% 5.75
+
+$distimes
+         DT50     DT90 DT50back
+meso 17.30248 82.91343 24.95943</code></pre>
+<div class="sourceCode" id="cb38"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">fomc_pH_2</span>, covariates <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>pH <span class="op">=</span> <span class="fl">7</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+     pH
+User  7
+
+$distimes
+         DT50     DT90 DT50back
+meso 6.986239 27.02927 8.136621</code></pre>
+</div>
+<div class="section level3">
+<h3 id="dfop">DFOP<a class="anchor" aria-label="anchor" href="#dfop"></a>
+</h3>
+<p>In the DFOP fits without covariate effects, random effects for two
+degradation parameters (<code>k2</code> and <code>g</code>) were
+identifiable.</p>
+<div class="sourceCode" id="cb40"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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_2</span><span class="op">[[</span><span class="st">"DFOP"</span>, <span class="st">"const"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_0</td>
+<td align="right">93.61</td>
+<td align="right">91.58</td>
+<td align="right">95.63</td>
+</tr>
+<tr class="even">
+<td align="left">log_k1</td>
+<td align="right">-1.53</td>
+<td align="right">-2.27</td>
+<td align="right">-0.79</td>
+</tr>
+<tr class="odd">
+<td align="left">log_k2</td>
+<td align="right">-3.42</td>
+<td align="right">-3.73</td>
+<td align="right">-3.11</td>
+</tr>
+<tr class="even">
+<td align="left">g_qlogis</td>
+<td align="right">-1.67</td>
+<td align="right">-2.57</td>
+<td align="right">-0.77</td>
+</tr>
+<tr class="odd">
+<td align="left">a.1</td>
+<td align="right">4.74</td>
+<td align="right">4.02</td>
+<td align="right">5.45</td>
+</tr>
+<tr class="even">
+<td align="left">SD.log_k2</td>
+<td align="right">0.60</td>
+<td align="right">0.38</td>
+<td align="right">0.81</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.g_qlogis</td>
+<td align="right">0.94</td>
+<td align="right">0.33</td>
+<td align="right">1.54</td>
+</tr>
+</tbody>
+</table>
+<p>A fit with pH dependent degradation parameters was obtained by
+excluding the same random effects as in the refined DFOP fit without
+covariate influence, and including covariate models for the two
+identifiable parameters <code>k2</code> and <code>g</code>.</p>
+<div class="sourceCode" id="cb41"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">dfop_pH</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/saem.html">saem</a></span><span class="op">(</span><span class="va">f_sep_const</span><span class="op">[</span><span class="st">"DFOP"</span>, <span class="op">]</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">"meso_0"</span>, <span class="st">"log_k1"</span><span class="op">)</span>,</span>
+<span>  covariates <span class="op">=</span> <span class="va">pH</span>,</span>
+<span>  covariate_models <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">log_k2</span> <span class="op">~</span> <span class="va">pH</span>, <span class="va">g_qlogis</span> <span class="op">~</span> <span class="va">pH</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<p>The corresponding parameters for the influence of soil pH are
+<code>beta_pH(log_k2)</code> for the influence of soil pH on
+<code>k2</code>, and <code>beta_pH(g_qlogis)</code> for its influence on
+<code>g</code>.</p>
+<div class="sourceCode" id="cb42"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">dfop_pH</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_0</td>
+<td align="right">92.84</td>
+<td align="right">90.85</td>
+<td align="right">94.84</td>
+</tr>
+<tr class="even">
+<td align="left">log_k1</td>
+<td align="right">-2.82</td>
+<td align="right">-3.09</td>
+<td align="right">-2.54</td>
+</tr>
+<tr class="odd">
+<td align="left">log_k2</td>
+<td align="right">-11.48</td>
+<td align="right">-15.32</td>
+<td align="right">-7.64</td>
+</tr>
+<tr class="even">
+<td align="left">beta_pH(log_k2)</td>
+<td align="right">1.31</td>
+<td align="right">0.69</td>
+<td align="right">1.92</td>
+</tr>
+<tr class="odd">
+<td align="left">g_qlogis</td>
+<td align="right">3.13</td>
+<td align="right">0.47</td>
+<td align="right">5.80</td>
+</tr>
+<tr class="even">
+<td align="left">beta_pH(g_qlogis)</td>
+<td align="right">-0.57</td>
+<td align="right">-1.04</td>
+<td align="right">-0.09</td>
+</tr>
+<tr class="odd">
+<td align="left">a.1</td>
+<td align="right">4.96</td>
+<td align="right">4.26</td>
+<td align="right">5.65</td>
+</tr>
+<tr class="even">
+<td align="left">SD.log_k2</td>
+<td align="right">0.76</td>
+<td align="right">0.47</td>
+<td align="right">1.05</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.g_qlogis</td>
+<td align="right">0.01</td>
+<td align="right">-9.96</td>
+<td align="right">9.97</td>
+</tr>
+</tbody>
+</table>
+<div class="sourceCode" id="cb43"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">dfop_pH</span><span class="op">)</span></span></code></pre></div>
+<pre><code>[1] "sd(g_qlogis)"</code></pre>
+<p>Confidence intervals for neither of them include zero, indicating a
+significant difference from zero. However, the random effect for
+<code>g</code> is now ill-defined. The fit is updated without this
+ill-defined random effect.</p>
+<div class="sourceCode" id="cb45"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">dfop_pH_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">dfop_pH</span>,</span>
+<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">"meso_0"</span>, <span class="st">"log_k1"</span>, <span class="st">"g_qlogis"</span><span class="op">)</span><span class="op">)</span></span>
+<span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">dfop_pH_2</span><span class="op">)</span></span></code></pre></div>
+<pre><code>[1] "beta_pH(g_qlogis)"</code></pre>
+<p>Now, the slope parameter for the pH effect on <code>g</code> is
+ill-defined. Therefore, another attempt is made without the
+corresponding covariate model.</p>
+<div class="sourceCode" id="cb47"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">dfop_pH_3</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/saem.html">saem</a></span><span class="op">(</span><span class="va">f_sep_const</span><span class="op">[</span><span class="st">"DFOP"</span>, <span class="op">]</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">"meso_0"</span>, <span class="st">"log_k1"</span><span class="op">)</span>,</span>
+<span>  covariates <span class="op">=</span> <span class="va">pH</span>,</span>
+<span>  covariate_models <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">log_k2</span> <span class="op">~</span> <span class="va">pH</span><span class="op">)</span><span class="op">)</span></span>
+<span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">dfop_pH_3</span><span class="op">)</span></span></code></pre></div>
+<pre><code>[1] "sd(g_qlogis)"</code></pre>
+<p>As the random effect for <code>g</code> is again ill-defined, the fit
+is repeated without it.</p>
+<div class="sourceCode" id="cb49"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">dfop_pH_4</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">dfop_pH_3</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">"meso_0"</span>, <span class="st">"log_k1"</span>, <span class="st">"g_qlogis"</span><span class="op">)</span><span class="op">)</span></span>
+<span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">dfop_pH_4</span><span class="op">)</span></span></code></pre></div>
+<p>While no ill-defined parameters remain, model comparison suggests
+that the previous model <code>dfop_pH_2</code> with two pH dependent
+parameters is preferable, based on information criteria as well as based
+on the likelihood ratio test.</p>
+<div class="sourceCode" id="cb50"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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_2</span><span class="op">[[</span><span class="st">"DFOP"</span>, <span class="st">"const"</span><span class="op">]</span><span class="op">]</span>, <span class="va">dfop_pH</span>, <span class="va">dfop_pH_2</span>, <span class="va">dfop_pH_3</span>, <span class="va">dfop_pH_4</span><span class="op">)</span></span></code></pre></div>
+<pre><code>Data: 116 observations of 1 variable(s) grouped in 18 datasets
+
+                            npar    AIC    BIC     Lik
+f_saem_2[["DFOP", "const"]]    7 782.94 789.18 -384.47
+dfop_pH_4                      7 767.35 773.58 -376.68
+dfop_pH_2                      8 765.14 772.26 -374.57
+dfop_pH_3                      8 769.00 776.12 -376.50
+dfop_pH                        9 769.10 777.11 -375.55</code></pre>
+<div class="sourceCode" id="cb52"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">dfop_pH_2</span>, <span class="va">dfop_pH_4</span>, test <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
+<pre><code>Data: 116 observations of 1 variable(s) grouped in 18 datasets
+
+          npar    AIC    BIC     Lik  Chisq Df Pr(&gt;Chisq)  
+dfop_pH_4    7 767.35 773.58 -376.68                       
+dfop_pH_2    8 765.14 772.26 -374.57 4.2153  1    0.04006 *
+---
+Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1</code></pre>
+<p>When focussing on parameter identifiability using the test if the
+confidence interval includes zero, <code>dfop_pH_4</code> would still be
+the preferred model. However, it should be kept in mind that parameter
+confidence intervals are constructed using a simple linearisation of the
+likelihood. As the confidence interval of the random effect for
+<code>g</code> only marginally includes zero, it is suggested that this
+is acceptable, and that <code>dfop_pH_2</code> can be considered the
+most preferable model.</p>
+<div class="sourceCode" id="cb54"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">dfop_pH_2</span><span class="op">)</span></span></code></pre></div>
+<p><img src="2023_mesotrione_parent_files/figure-html/unnamed-chunk-19-1.png" width="700" style="display: block; margin: auto;"></p>
+<div class="sourceCode" id="cb55"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">dfop_pH_2</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+      pH
+50% 5.75
+
+$distimes
+         DT50     DT90 DT50back  DT50_k1  DT50_k2
+meso 18.36876 73.51841 22.13125 4.191901 23.98672</code></pre>
+<div class="sourceCode" id="cb57"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">dfop_pH_2</span>, covariates <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>pH <span class="op">=</span> <span class="fl">7</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+     pH
+User  7
+
+$distimes
+         DT50     DT90 DT50back  DT50_k1  DT50_k2
+meso 8.346428 28.34437 8.532507 4.191901 8.753618</code></pre>
+</div>
+<div class="section level3">
+<h3 id="sforb">SFORB<a class="anchor" aria-label="anchor" href="#sforb"></a>
+</h3>
+<div class="sourceCode" id="cb59"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">sforb_pH</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/saem.html">saem</a></span><span class="op">(</span><span class="va">f_sep_const</span><span class="op">[</span><span class="st">"SFORB"</span>, <span class="op">]</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">"meso_free_0"</span>, <span class="st">"log_k_meso_free_bound"</span><span class="op">)</span>,</span>
+<span>  covariates <span class="op">=</span> <span class="va">pH</span>,</span>
+<span>  covariate_models <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">log_k_meso_free</span> <span class="op">~</span> <span class="va">pH</span>, <span class="va">log_k_meso_bound_free</span> <span class="op">~</span> <span class="va">pH</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<div class="sourceCode" id="cb60"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">sforb_pH</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_free_0</td>
+<td align="right">93.42</td>
+<td align="right">91.32</td>
+<td align="right">95.52</td>
+</tr>
+<tr class="even">
+<td align="left">log_k_meso_free</td>
+<td align="right">-5.37</td>
+<td align="right">-6.94</td>
+<td align="right">-3.81</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_k_meso_free)</td>
+<td align="right">0.42</td>
+<td align="right">0.18</td>
+<td align="right">0.67</td>
+</tr>
+<tr class="even">
+<td align="left">log_k_meso_free_bound</td>
+<td align="right">-3.49</td>
+<td align="right">-4.92</td>
+<td align="right">-2.05</td>
+</tr>
+<tr class="odd">
+<td align="left">log_k_meso_bound_free</td>
+<td align="right">-9.98</td>
+<td align="right">-19.22</td>
+<td align="right">-0.74</td>
+</tr>
+<tr class="even">
+<td align="left">beta_pH(log_k_meso_bound_free)</td>
+<td align="right">1.23</td>
+<td align="right">-0.21</td>
+<td align="right">2.67</td>
+</tr>
+<tr class="odd">
+<td align="left">a.1</td>
+<td align="right">4.90</td>
+<td align="right">4.18</td>
+<td align="right">5.63</td>
+</tr>
+<tr class="even">
+<td align="left">SD.log_k_meso_free</td>
+<td align="right">0.35</td>
+<td align="right">0.23</td>
+<td align="right">0.47</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.log_k_meso_bound_free</td>
+<td align="right">0.13</td>
+<td align="right">-1.95</td>
+<td align="right">2.20</td>
+</tr>
+</tbody>
+</table>
+<p>The confidence interval of
+<code>beta_pH(log_k_meso_bound_free)</code> includes zero, indicating
+that the influence of soil pH on <code>k_meso_bound_free</code> cannot
+reliably be quantified. Also, the confidence interval for the random
+effect on this parameter (<code>SD.log_k_meso_bound_free</code>)
+includes zero.</p>
+<p>Using the <code>illparms</code> function, these ill-defined
+parameters can be found more conveniently.</p>
+<div class="sourceCode" id="cb61"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">sforb_pH</span><span class="op">)</span></span></code></pre></div>
+<pre><code>[1] "sd(log_k_meso_bound_free)"      "beta_pH(log_k_meso_bound_free)"</code></pre>
+<p>To remove the ill-defined parameters, a second variant of the SFORB
+model with pH influence is fitted. No ill-defined parameters remain.</p>
+<div class="sourceCode" id="cb63"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">sforb_pH_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">sforb_pH</span>,</span>
+<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">"meso_free_0"</span>, <span class="st">"log_k_meso_free_bound"</span>, <span class="st">"log_k_meso_bound_free"</span><span class="op">)</span>,</span>
+<span>  covariate_models <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">log_k_meso_free</span> <span class="op">~</span> <span class="va">pH</span><span class="op">)</span><span class="op">)</span></span>
+<span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">sforb_pH_2</span><span class="op">)</span></span></code></pre></div>
+<p>The model comparison of the SFORB fits includes the refined model
+without covariate effect, and both versions of the SFORB fit with
+covariate effect.</p>
+<div class="sourceCode" id="cb64"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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_2</span><span class="op">[[</span><span class="st">"SFORB"</span>, <span class="st">"const"</span><span class="op">]</span><span class="op">]</span>, <span class="va">sforb_pH</span>, <span class="va">sforb_pH_2</span>, test <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
+<pre><code>Data: 116 observations of 1 variable(s) grouped in 18 datasets
+
+                             npar    AIC    BIC     Lik   Chisq Df Pr(&gt;Chisq)  
+f_saem_2[["SFORB", "const"]]    7 783.40 789.63 -384.70                        
+sforb_pH_2                      7 770.94 777.17 -378.47 12.4616  0             
+sforb_pH                        9 768.81 776.83 -375.41  6.1258  2    0.04675 *
+---
+Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1</code></pre>
+<p>The first model including pH influence is preferable based on
+information criteria and the likelihood ratio test. However, as it is
+not fully identifiable, the second model is selected.</p>
+<div class="sourceCode" id="cb66"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">sforb_pH_2</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_free_0</td>
+<td align="right">93.32</td>
+<td align="right">91.16</td>
+<td align="right">95.48</td>
+</tr>
+<tr class="even">
+<td align="left">log_k_meso_free</td>
+<td align="right">-6.15</td>
+<td align="right">-7.43</td>
+<td align="right">-4.86</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_k_meso_free)</td>
+<td align="right">0.54</td>
+<td align="right">0.33</td>
+<td align="right">0.75</td>
+</tr>
+<tr class="even">
+<td align="left">log_k_meso_free_bound</td>
+<td align="right">-3.80</td>
+<td align="right">-5.20</td>
+<td align="right">-2.40</td>
+</tr>
+<tr class="odd">
+<td align="left">log_k_meso_bound_free</td>
+<td align="right">-2.95</td>
+<td align="right">-4.26</td>
+<td align="right">-1.64</td>
+</tr>
+<tr class="even">
+<td align="left">a.1</td>
+<td align="right">5.08</td>
+<td align="right">4.38</td>
+<td align="right">5.79</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.log_k_meso_free</td>
+<td align="right">0.33</td>
+<td align="right">0.22</td>
+<td align="right">0.45</td>
+</tr>
+</tbody>
+</table>
+<div class="sourceCode" id="cb67"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">sforb_pH_2</span><span class="op">)</span></span></code></pre></div>
+<p><img src="2023_mesotrione_parent_files/figure-html/unnamed-chunk-25-1.png" width="700" style="display: block; margin: auto;"></p>
+<div class="sourceCode" id="cb68"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">sforb_pH_2</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+      pH
+50% 5.75
+
+$ff
+meso_free 
+        1 
+
+$SFORB
+   meso_b1    meso_b2     meso_g 
+0.09735824 0.02631699 0.31602120 
+
+$distimes
+         DT50     DT90 DT50back DT50_meso_b1 DT50_meso_b2
+meso 16.86549 73.15824 22.02282     7.119554     26.33839</code></pre>
+<div class="sourceCode" id="cb70"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">sforb_pH_2</span>, covariates <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>pH <span class="op">=</span> <span class="fl">7</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+     pH
+User  7
+
+$ff
+meso_free 
+        1 
+
+$SFORB
+   meso_b1    meso_b2     meso_g 
+0.13315233 0.03795988 0.61186191 
+
+$distimes
+         DT50     DT90 DT50back DT50_meso_b1 DT50_meso_b2
+meso 7.932495 36.93311 11.11797     5.205671        18.26</code></pre>
+</div>
+<div class="section level3">
+<h3 id="hs">HS<a class="anchor" aria-label="anchor" href="#hs"></a>
+</h3>
+<div class="sourceCode" id="cb72"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">hs_pH</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/saem.html">saem</a></span><span class="op">(</span><span class="va">f_sep_const</span><span class="op">[</span><span class="st">"HS"</span>, <span class="op">]</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">"meso_0"</span><span class="op">)</span>,</span>
+<span>  covariates <span class="op">=</span> <span class="va">pH</span>,</span>
+<span>  covariate_models <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">log_k1</span> <span class="op">~</span> <span class="va">pH</span>, <span class="va">log_k2</span> <span class="op">~</span> <span class="va">pH</span>, <span class="va">log_tb</span> <span class="op">~</span> <span class="va">pH</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<div class="sourceCode" id="cb73"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">hs_pH</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_0</td>
+<td align="right">93.33</td>
+<td align="right">91.47</td>
+<td align="right">95.19</td>
+</tr>
+<tr class="even">
+<td align="left">log_k1</td>
+<td align="right">-5.81</td>
+<td align="right">-7.27</td>
+<td align="right">-4.36</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_k1)</td>
+<td align="right">0.47</td>
+<td align="right">0.23</td>
+<td align="right">0.72</td>
+</tr>
+<tr class="even">
+<td align="left">log_k2</td>
+<td align="right">-6.80</td>
+<td align="right">-8.76</td>
+<td align="right">-4.83</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_k2)</td>
+<td align="right">0.54</td>
+<td align="right">0.21</td>
+<td align="right">0.87</td>
+</tr>
+<tr class="even">
+<td align="left">log_tb</td>
+<td align="right">3.25</td>
+<td align="right">1.25</td>
+<td align="right">5.25</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_tb)</td>
+<td align="right">-0.10</td>
+<td align="right">-0.43</td>
+<td align="right">0.23</td>
+</tr>
+<tr class="even">
+<td align="left">a.1</td>
+<td align="right">4.49</td>
+<td align="right">3.78</td>
+<td align="right">5.21</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.log_k1</td>
+<td align="right">0.37</td>
+<td align="right">0.24</td>
+<td align="right">0.51</td>
+</tr>
+<tr class="even">
+<td align="left">SD.log_k2</td>
+<td align="right">0.29</td>
+<td align="right">0.10</td>
+<td align="right">0.48</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.log_tb</td>
+<td align="right">0.25</td>
+<td align="right">-0.07</td>
+<td align="right">0.57</td>
+</tr>
+</tbody>
+</table>
+<div class="sourceCode" id="cb74"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">hs_pH</span><span class="op">)</span></span></code></pre></div>
+<pre><code>[1] "sd(log_tb)"      "beta_pH(log_tb)"</code></pre>
+<p>According to the output of the <code>illparms</code> function, the
+random effect on the break time <code>tb</code> cannot reliably be
+quantified, neither can the influence of soil pH on <code>tb</code>. The
+fit is repeated without the corresponding covariate model, and no
+ill-defined parameters remain.</p>
+<div class="sourceCode" id="cb76"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="va">hs_pH_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">hs_pH</span>, covariate_models <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span><span class="va">log_k1</span> <span class="op">~</span> <span class="va">pH</span>, <span class="va">log_k2</span> <span class="op">~</span> <span class="va">pH</span><span class="op">)</span><span class="op">)</span></span>
+<span><span class="fu"><a href="../../reference/illparms.html">illparms</a></span><span class="op">(</span><span class="va">hs_pH_2</span><span class="op">)</span></span></code></pre></div>
+<p>Model comparison confirms that this model is preferable to the fit
+without covariate influence, and also to the first version with
+covariate influence.</p>
+<div class="sourceCode" id="cb77"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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_2</span><span class="op">[[</span><span class="st">"HS"</span>, <span class="st">"const"</span><span class="op">]</span><span class="op">]</span>, <span class="va">hs_pH</span>, <span class="va">hs_pH_2</span>, test <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></code></pre></div>
+<pre><code>Data: 116 observations of 1 variable(s) grouped in 18 datasets
+
+                          npar    AIC    BIC     Lik  Chisq Df Pr(&gt;Chisq)    
+f_saem_2[["HS", "const"]]    8 780.08 787.20 -382.04                         
+hs_pH_2                     10 766.47 775.37 -373.23 17.606  2  0.0001503 ***
+hs_pH                       11 769.80 779.59 -373.90  0.000  1  1.0000000    
+---
+Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1</code></pre>
+<div class="sourceCode" id="cb79"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">hs_pH_2</span><span class="op">)</span><span class="op">$</span><span class="va">confint_trans</span> <span class="op">|&gt;</span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>digits <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div>
+<table class="table">
+<thead><tr class="header">
+<th align="left"></th>
+<th align="right">est.</th>
+<th align="right">lower</th>
+<th align="right">upper</th>
+</tr></thead>
+<tbody>
+<tr class="odd">
+<td align="left">meso_0</td>
+<td align="right">93.33</td>
+<td align="right">91.50</td>
+<td align="right">95.15</td>
+</tr>
+<tr class="even">
+<td align="left">log_k1</td>
+<td align="right">-5.68</td>
+<td align="right">-7.09</td>
+<td align="right">-4.27</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_k1)</td>
+<td align="right">0.46</td>
+<td align="right">0.22</td>
+<td align="right">0.69</td>
+</tr>
+<tr class="even">
+<td align="left">log_k2</td>
+<td align="right">-6.61</td>
+<td align="right">-8.34</td>
+<td align="right">-4.88</td>
+</tr>
+<tr class="odd">
+<td align="left">beta_pH(log_k2)</td>
+<td align="right">0.50</td>
+<td align="right">0.21</td>
+<td align="right">0.79</td>
+</tr>
+<tr class="even">
+<td align="left">log_tb</td>
+<td align="right">2.70</td>
+<td align="right">2.33</td>
+<td align="right">3.08</td>
+</tr>
+<tr class="odd">
+<td align="left">a.1</td>
+<td align="right">4.45</td>
+<td align="right">3.74</td>
+<td align="right">5.16</td>
+</tr>
+<tr class="even">
+<td align="left">SD.log_k1</td>
+<td align="right">0.36</td>
+<td align="right">0.22</td>
+<td align="right">0.49</td>
+</tr>
+<tr class="odd">
+<td align="left">SD.log_k2</td>
+<td align="right">0.23</td>
+<td align="right">0.02</td>
+<td align="right">0.43</td>
+</tr>
+<tr class="even">
+<td align="left">SD.log_tb</td>
+<td align="right">0.55</td>
+<td align="right">0.25</td>
+<td align="right">0.85</td>
+</tr>
+</tbody>
+</table>
+<div class="sourceCode" id="cb80"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">hs_pH_2</span><span class="op">)</span></span></code></pre></div>
+<p><img src="2023_mesotrione_parent_files/figure-html/unnamed-chunk-30-1.png" width="700" style="display: block; margin: auto;"></p>
+<div class="sourceCode" id="cb81"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">hs_pH_2</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+      pH
+50% 5.75
+
+$distimes
+         DT50     DT90 DT50back  DT50_k1  DT50_k2
+meso 14.68725 82.45287 24.82079 14.68725 29.29299</code></pre>
+<div class="sourceCode" id="cb83"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">hs_pH_2</span>, covariates <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>pH <span class="op">=</span> <span class="fl">7</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+     pH
+User  7
+
+$distimes
+         DT50     DT90 DT50back  DT50_k1  DT50_k2
+meso 8.298536 38.85371 11.69613 8.298536 15.71561</code></pre>
+</div>
+<div class="section level3">
+<h3 id="comparison-across-parent-models">Comparison across parent models<a class="anchor" aria-label="anchor" href="#comparison-across-parent-models"></a>
+</h3>
+<p>After model reduction for all models with pH influence, they are
+compared with each other.</p>
+<div class="sourceCode" id="cb85"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><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">sfo_pH</span>, <span class="va">fomc_pH_2</span>, <span class="va">dfop_pH_2</span>, <span class="va">dfop_pH_4</span>, <span class="va">sforb_pH_2</span>, <span class="va">hs_pH_2</span><span class="op">)</span></span></code></pre></div>
+<pre><code>Data: 116 observations of 1 variable(s) grouped in 18 datasets
+
+           npar    AIC    BIC     Lik
+sfo_pH        5 783.09 787.54 -386.54
+fomc_pH_2     6 767.49 772.83 -377.75
+dfop_pH_4     7 767.35 773.58 -376.68
+sforb_pH_2    7 770.94 777.17 -378.47
+dfop_pH_2     8 765.14 772.26 -374.57
+hs_pH_2      10 766.47 775.37 -373.23</code></pre>
+<p>The DFOP model with pH influence on <code>k2</code> and
+<code>g</code> and a random effect only on <code>k2</code> is finally
+selected as the best fit.</p>
+<p>The endpoints resulting from this model are listed below. Please
+refer to the Appendix for a detailed listing.</p>
+<div class="sourceCode" id="cb87"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">dfop_pH_2</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+      pH
+50% 5.75
+
+$distimes
+         DT50     DT90 DT50back  DT50_k1  DT50_k2
+meso 18.36876 73.51841 22.13125 4.191901 23.98672</code></pre>
+<div class="sourceCode" id="cb89"><pre class="downlit sourceCode r">
+<code class="sourceCode R"><span><span class="fu"><a href="../../reference/endpoints.html">endpoints</a></span><span class="op">(</span><span class="va">dfop_pH_2</span>, covariates <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>pH <span class="op">=</span> <span class="fl">7</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
+<pre><code>$covariates
+     pH
+User  7
+
+$distimes
+         DT50     DT90 DT50back  DT50_k1  DT50_k2
+meso 8.346428 28.34437 8.532507 4.191901 8.753618</code></pre>
+</div>
+</div>
+<div class="section level2">
+<h2 id="conclusions">Conclusions<a class="anchor" aria-label="anchor" href="#conclusions"></a>
+</h2>
+<p>These evaluations demonstrate that covariate effects can be included
+for all types of parent degradation models. These models can then be
+further refined to make them fully identifiable.</p>
+</div>
+<div class="section level2">
+<h2 id="appendix">Appendix<a class="anchor" aria-label="anchor" href="#appendix"></a>
+</h2>
+<div class="section level3">
+<h3 id="hierarchical-fit-listings">Hierarchical fit listings<a class="anchor" aria-label="anchor" href="#hierarchical-fit-listings"></a>
+</h3>
+<div class="section level4">
+<h4 id="fits-without-covariate-effects">Fits without covariate effects<a class="anchor" aria-label="anchor" href="#fits-without-covariate-effects"></a>
+</h4>
+
+</div>
+<div class="section level4">
+<h4 id="fits-with-covariate-effects">Fits with covariate effects<a class="anchor" aria-label="anchor" href="#fits-with-covariate-effects"></a>
+</h4>
+
+</div>
+</div>
+<div class="section level3">
+<h3 id="session-info">Session info<a class="anchor" aria-label="anchor" href="#session-info"></a>
+</h3>
+<pre><code>R version 4.3.1 (2023-06-16)
+Platform: x86_64-pc-linux-gnu (64-bit)
+Running under: Debian GNU/Linux 12 (bookworm)
+
+Matrix products: default
+BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.11.0 
+LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.11.0
+
+locale:
+ [1] LC_CTYPE=de_DE.UTF-8       LC_NUMERIC=C              
+ [3] LC_TIME=de_DE.UTF-8        LC_COLLATE=de_DE.UTF-8    
+ [5] LC_MONETARY=de_DE.UTF-8    LC_MESSAGES=de_DE.UTF-8   
+ [7] LC_PAPER=de_DE.UTF-8       LC_NAME=C                 
+ [9] LC_ADDRESS=C               LC_TELEPHONE=C            
+[11] LC_MEASUREMENT=de_DE.UTF-8 LC_IDENTIFICATION=C       
+
+time zone: Europe/Berlin
+tzcode source: system (glibc)
+
+attached base packages:
+[1] parallel  stats     graphics  grDevices utils     datasets  methods  
+[8] base     
+
+other attached packages:
+[1] saemix_3.2 npde_3.3   knitr_1.43 mkin_1.2.5
+
+loaded via a namespace (and not attached):
+ [1] sass_0.4.6        utf8_1.2.3        generics_0.1.3    stringi_1.7.12   
+ [5] lattice_0.20-45   digest_0.6.31     magrittr_2.0.3    evaluate_0.21    
+ [9] grid_4.3.1        fastmap_1.1.1     cellranger_1.1.0  rprojroot_2.0.3  
+[13] jsonlite_1.8.5    mclust_6.0.0      gridExtra_2.3     purrr_1.0.1      
+[17] fansi_1.0.4       scales_1.2.1      codetools_0.2-19  textshaping_0.3.6
+[21] jquerylib_0.1.4   cli_3.6.1         rlang_1.1.1       munsell_0.5.0    
+[25] cachem_1.0.8      yaml_2.3.7        tools_4.3.1       memoise_2.0.1    
+[29] dplyr_1.1.2       colorspace_2.1-0  ggplot2_3.4.2     vctrs_0.6.2      
+[33] R6_2.5.1          zoo_1.8-12        lifecycle_1.0.3   stringr_1.5.0    
+[37] fs_1.6.2          ragg_1.2.5        pkgconfig_2.0.3   desc_1.4.2       
+[41] pkgdown_2.0.7     bslib_0.4.2       pillar_1.9.0      gtable_0.3.3     
+[45] glue_1.6.2        systemfonts_1.0.4 highr_0.10        xfun_0.39        
+[49] tibble_3.2.1      lmtest_0.9-40     tidyselect_1.2.0  htmltools_0.5.5  
+[53] nlme_3.1-162      rmarkdown_2.22    compiler_4.3.1    readxl_1.4.2     </code></pre>
+</div>
+<div class="section level3">
+<h3 id="hardware-info">Hardware info<a class="anchor" aria-label="anchor" href="#hardware-info"></a>
+</h3>
+<pre><code>CPU model: Intel(R) Core(TM) i7-4710MQ CPU @ 2.50GHz</code></pre>
+<pre><code>MemTotal:       12165632 kB</code></pre>
+</div>
+</div>
+  </div>
+
+  <div class="col-md-3 hidden-xs hidden-sm" id="pkgdown-sidebar">
+
+        <nav id="toc" data-toggle="toc"><h2 data-toc-skip>Contents</h2>
+    </nav>
+</div>
+
+</div>
+
+
+
+      <footer><div class="copyright">
+  <p></p>
+<p>Developed by Johannes Ranke.</p>
+</div>
+
+<div class="pkgdown">
+  <p></p>
+<p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.0.7.</p>
+</div>
+
+      </footer>
+</div>
+
+  
+
+
+  
+
+  </body>
+</html>
diff --git a/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-14-1.png b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-14-1.png
new file mode 100644
index 00000000..863a48bd
--- /dev/null
+++ b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-14-1.png
diff --git a/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-19-1.png b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-19-1.png
new file mode 100644
index 00000000..256b2b68
--- /dev/null
+++ b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-19-1.png
diff --git a/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-25-1.png b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-25-1.png
new file mode 100644
index 00000000..59011020
--- /dev/null
+++ b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-25-1.png
diff --git a/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-30-1.png b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-30-1.png
new file mode 100644
index 00000000..f427bc39
--- /dev/null
+++ b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-30-1.png
diff --git a/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-8-1.png b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-8-1.png
new file mode 100644
index 00000000..7c3b460b
--- /dev/null
+++ b/docs/articles/prebuilt/2023_mesotrione_parent_files/figure-html/unnamed-chunk-8-1.png
diff --git a/docs/authors.html b/docs/authors.html
index 313c86e9..8b91ee77 100644
--- a/docs/authors.html
+++ b/docs/authors.html
@@ -17,7 +17,7 @@
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
@@ -132,13 +132,13 @@
 
     <p>Ranke J (2023).
 <em>mkin: Kinetic Evaluation of Chemical Degradation Data</em>.
-R package version 1.2.4, <a href="https://pkgdown.jrwb.de/mkin/">https://pkgdown.jrwb.de/mkin/</a>. 
+R package version 1.2.5, <a href="https://pkgdown.jrwb.de/mkin/">https://pkgdown.jrwb.de/mkin/</a>. 
 </p>
     <pre>@Manual{,
   title = {mkin: Kinetic Evaluation of Chemical Degradation Data},
   author = {Johannes Ranke},
   year = {2023},
-  note = {R package version 1.2.4},
+  note = {R package version 1.2.5},
   url = {https://pkgdown.jrwb.de/mkin/},
 }</pre>
 
diff --git a/docs/index.html b/docs/index.html
index 06cdfdfb..79716087 100644
--- a/docs/index.html
+++ b/docs/index.html
@@ -44,7 +44,7 @@
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
diff --git a/docs/news/index.html b/docs/news/index.html
index 16fce355..7f7c0a6d 100644
--- a/docs/news/index.html
+++ b/docs/news/index.html
@@ -17,7 +17,7 @@
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="../index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
@@ -105,6 +105,12 @@
     </div>
 
     <div class="section level2">
+<h2 class="page-header" data-toc-text="1.2.5" id="mkin-125">mkin 1.2.5<a class="anchor" aria-label="anchor" href="#mkin-125"></a></h2>
+<ul><li><p>‘vignettes/mesotrione_parent_2023.rnw’: Prebuilt vignette showing how covariate modelling can be done for all relevant parent degradation models.</p></li>
+<li><p>‘inst/testdata/mesotrione_soil_efsa_2016}.xlsx’: Another example spreadsheets for use with ‘read_spreadsheet()’, featuring pH dependent degradation</p></li>
+<li><p>R/illparms.R: Fix the detection of ill-defined slope or error model parameters for the case that the estimate is negative</p></li>
+</ul></div>
+    <div class="section level2">
 <h2 class="page-header" data-toc-text="1.2.4" id="mkin-124">mkin 1.2.4<a class="anchor" aria-label="anchor" href="#mkin-124"></a></h2>
 <ul><li>R/endpoints.R: Fix the calculation of endpoints for user specified covariate values</li>
 </ul></div>
diff --git a/docs/pkgdown.yml b/docs/pkgdown.yml
index fdcd875b..65d77082 100644
--- a/docs/pkgdown.yml
+++ b/docs/pkgdown.yml
@@ -8,6 +8,7 @@ articles:
   2022_cyan_pathway: prebuilt/2022_cyan_pathway.html
   2022_dmta_parent: prebuilt/2022_dmta_parent.html
   2022_dmta_pathway: prebuilt/2022_dmta_pathway.html
+  2023_mesotrione_parent: prebuilt/2023_mesotrione_parent.html
   twa: twa.html
   FOCUS_Z: web_only/FOCUS_Z.html
   NAFTA_examples: web_only/NAFTA_examples.html
@@ -16,7 +17,7 @@ articles:
   dimethenamid_2018: web_only/dimethenamid_2018.html
   multistart: web_only/multistart.html
   saem_benchmarks: web_only/saem_benchmarks.html
-last_built: 2023-05-19T09:04Z
+last_built: 2023-08-09T15:42Z
 urls:
   reference: https://pkgdown.jrwb.de/mkin/reference
   article: https://pkgdown.jrwb.de/mkin/articles
diff --git a/docs/reference/D24_2014.html b/docs/reference/D24_2014.html
index fcc3cec3..1e35e864 100644
--- a/docs/reference/D24_2014.html
+++ b/docs/reference/D24_2014.html
@@ -22,7 +22,7 @@ constrained by data protection regulations."><!-- mathjax --><script src="https:
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="../index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
diff --git a/docs/reference/Rplot001.png b/docs/reference/Rplot001.png
index 5de2bdc7..fa2316b9 100644
--- a/docs/reference/Rplot001.png
+++ b/docs/reference/Rplot001.png
diff --git a/docs/reference/Rplot002.png b/docs/reference/Rplot002.png
index 556ca0a7..fa2316b9 100644
--- a/docs/reference/Rplot002.png
+++ b/docs/reference/Rplot002.png
diff --git a/docs/reference/Rplot003.png b/docs/reference/Rplot003.png
index 30cf38f3..fa2316b9 100644
--- a/docs/reference/Rplot003.png
+++ b/docs/reference/Rplot003.png
diff --git a/docs/reference/dimethenamid_2018.html b/docs/reference/dimethenamid_2018.html
index a1c800e8..961b28d2 100644
--- a/docs/reference/dimethenamid_2018.html
+++ b/docs/reference/dimethenamid_2018.html
@@ -22,7 +22,7 @@ constrained by data protection regulations."><!-- mathjax --><script src="https:
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="../index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
@@ -222,10 +222,10 @@ specific pieces of information in the comments.</p>
 <span class="r-in"><span><span class="co">#saemix::plot(f_dmta_saem_tc$so, plot.type = "convergence")</span></span></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_dmta_saem_tc</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.4 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> R version used for fitting:         4.3.0 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Date of fit:     Fri May 19 09:15:21 2023 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Date of summary: Fri May 19 09:15:21 2023 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> mkin version used for pre-fitting:  1.2.5 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> R version used for fitting:         4.3.1 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Date of fit:     Wed Aug  9 17:55:36 2023 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Date of summary: Wed Aug  9 17:55:36 2023 </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_DMTA/dt = - k_DMTA * DMTA</span>
@@ -238,7 +238,7 @@ specific pieces of information in the comments.</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 deSolve </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 282.941 s</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Fitted in 796.123 s</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> Using 300, 100 iterations and 9 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>
diff --git a/docs/reference/index.html b/docs/reference/index.html
index 4d59b3ab..a0d9e12b 100644
--- a/docs/reference/index.html
+++ b/docs/reference/index.html
@@ -17,7 +17,7 @@
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="../index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
diff --git a/docs/reference/read_spreadsheet.html b/docs/reference/read_spreadsheet.html
index 4d1aa706..7603cee2 100644
--- a/docs/reference/read_spreadsheet.html
+++ b/docs/reference/read_spreadsheet.html
@@ -22,7 +22,7 @@ factors can be given in columns named 'Temperature' and 'Moisture'."><!-- mathja
       </button>
       <span class="navbar-brand">
         <a class="navbar-link" href="../index.html">mkin</a>
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.4</span>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.2.5</span>
       </span>
     </div>
 
@@ -154,7 +154,7 @@ and moisture normalisation factors in the sheet 'Datasets'?</p></dd>
 The first row read after the header read in from this sheet is assumed
 to contain name and acronym of the parent compound.</p>
 <p>The dataset sheets should be named using the dataset numbers read in from
-the 'Datasets' sheet, i.e. '1', '2', ... . In each dataset sheet, name
+the 'Datasets' sheet, i.e. '1', '2', ... . In each dataset sheet, the name
 of the observed variable (e.g. the acronym of the parent compound or
 one of its transformation products) should be in the first column,
 the time values should be in the second colum, and the observed value
@@ -166,7 +166,7 @@ column and the column must have the same name as the second column in
 'Datasets'. Covariates will be read in from columns four and higher.
 Their names should preferably not contain special characters like spaces,
 so they can be easily used for specifying covariate models.</p>
-<p>An similar data structure is defined as the R6 class <a href="mkindsg.html">mkindsg</a>, but
+<p>A similar data structure is defined as the R6 class <a href="mkindsg.html">mkindsg</a>, but
 is probably more complicated to use.</p>
     </div>
 
diff --git a/docs/sitemap.xml b/docs/sitemap.xml
index f1b20197..9a84d4c0 100644
--- a/docs/sitemap.xml
+++ b/docs/sitemap.xml
@@ -25,6 +25,9 @@
     <loc>https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_pathway.html</loc>
   </url>
   <url>
+    <loc>https://pkgdown.jrwb.de/mkin/articles/prebuilt/2023_mesotrione_parent.html</loc>
+  </url>
+  <url>
     <loc>https://pkgdown.jrwb.de/mkin/articles/twa.html</loc>
   </url>
   <url>
author	Johannes Ranke <jranke@uni-bremen.de>	2023-08-09 17:59:42 +0200
committer	Johannes Ranke <jranke@uni-bremen.de>	2023-08-09 17:59:42 +0200
commit	9abab1e2d4385039b01ad3dc0d9c5966bbe94fee (patch)
tree	cc9f23c3abe70e289b40c2cb03a4cda4ed157122 /docs
parent	f21373a16d368c56234246417a32bdfb95589cf9 (diff)