diff options
Diffstat (limited to 'docs/articles/prebuilt')
72 files changed, 896 insertions, 1044 deletions
diff --git a/docs/articles/prebuilt/2022_cyan_pathway.html b/docs/articles/prebuilt/2022_cyan_pathway.html index c22c6735..3dd83fee 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway.html +++ b/docs/articles/prebuilt/2022_cyan_pathway.html @@ -4,142 +4,90 @@ <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"> +<meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"> <title>Testing hierarchical pathway kinetics with residue data on cyantraniliprole • 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 hierarchical pathway kinetics with residue data on cyantraniliprole"> -<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]--> +<script src="../../deps/jquery-3.6.0/jquery-3.6.0.min.js"></script><meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"> +<link href="../../deps/bootstrap-5.3.1/bootstrap.min.css" rel="stylesheet"> +<script src="../../deps/bootstrap-5.3.1/bootstrap.bundle.min.js"></script><link href="../../deps/font-awesome-6.5.2/css/all.min.css" rel="stylesheet"> +<link href="../../deps/font-awesome-6.5.2/css/v4-shims.min.css" rel="stylesheet"> +<script src="../../deps/headroom-0.11.0/headroom.min.js"></script><script src="../../deps/headroom-0.11.0/jQuery.headroom.min.js"></script><script src="../../deps/bootstrap-toc-1.0.1/bootstrap-toc.min.js"></script><script src="../../deps/clipboard.js-2.0.11/clipboard.min.js"></script><script src="../../deps/search-1.0.0/autocomplete.jquery.min.js"></script><script src="../../deps/search-1.0.0/fuse.min.js"></script><script src="../../deps/search-1.0.0/mark.min.js"></script><!-- pkgdown --><script src="../../pkgdown.js"></script><meta property="og:title" content="Testing hierarchical pathway kinetics with residue data on cyantraniliprole"> </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.4</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> +<body> + <a href="#main" class="visually-hidden-focusable">Skip to contents</a> + + + <nav class="navbar navbar-expand-lg fixed-top bg-light" data-bs-theme="default" aria-label="Site navigation"><div class="container"> + + <a class="navbar-brand me-2" href="../../index.html">mkin</a> + + <small class="nav-text text-muted me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="Released version">1.2.10</small> + + + <button class="navbar-toggler" type="button" data-bs-toggle="collapse" data-bs-target="#navbar" aria-controls="navbar" aria-expanded="false" aria-label="Toggle navigation"> + <span class="navbar-toggler-icon"></span> + </button> + + <div id="navbar" class="collapse navbar-collapse ms-3"> + <ul class="navbar-nav me-auto"> +<li class="nav-item"><a class="nav-link" href="../../reference/index.html">Reference</a></li> +<li class="active nav-item dropdown"> + <button class="nav-link dropdown-toggle" type="button" id="dropdown-articles" data-bs-toggle="dropdown" aria-expanded="false" aria-haspopup="true">Articles</button> + <ul class="dropdown-menu" aria-labelledby="dropdown-articles"> +<li><a class="dropdown-item" href="../../articles/mkin.html">Introduction to mkin</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Example evaluations with (generalised) nonlinear least squares</h6></li> + <li><a class="dropdown-item" href="../../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a></li> + <li><a class="dropdown-item" href="../../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Example evaluations with hierarchical models (nonlinear mixed-effects models)</h6></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_dmta_parent.html">Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_dmta_pathway.html">Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2023_mesotrione_parent.html">Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_cyan_pathway.html">Testing hierarchical pathway kinetics with residue data on cyantraniliprole</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/dimethenamid_2018.html">Comparison of saemix and nlme evaluations of dimethenamid data from 2018</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/multistart.html">Short demo of the multistart method</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Performance</h6></li> + <li><a class="dropdown-item" href="../../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/benchmarks.html">Benchmark timings for mkin</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/saem_benchmarks.html">Benchmark timings for saem.mmkin</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Miscellaneous</h6></li> + <li><a class="dropdown-item" href="../../articles/twa.html">Calculation of time weighted average concentrations with mkin</a></li> + <li><a class="dropdown-item" 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> +<li class="nav-item"><a class="nav-link" href="../../coverage/coverage.html">Test coverage</a></li> +<li class="nav-item"><a class="nav-link" 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 class="navbar-nav"> +<li class="nav-item"><form class="form-inline" role="search"> + <input class="form-control" type="search" name="search-input" id="search-input" autocomplete="off" aria-label="Search site" placeholder="Search for" data-search-index="../../search.json"> +</form></li> +<li class="nav-item"><a class="external-link nav-link" href="https://github.com/jranke/mkin/" aria-label="GitHub"><span class="fa fab fa-github fa-lg"></span></a></li> </ul> </div> -<!--/.nav-collapse --> + + </div> -<!--/.container --> -</div> -<!--/.navbar --> +</nav><div class="container template-article"> + - - </header><div class="row"> - <div class="col-md-9 contents"> - <div class="page-header toc-ignore"> - <h1 data-toc-skip>Testing hierarchical pathway kinetics with -residue data on cyantraniliprole</h1> + +<div class="row"> + <main id="main" class="col-md-9"><div class="page-header"> + + <h1>Testing hierarchical pathway kinetics with residue data on cyantraniliprole</h1> <h4 data-toc-skip class="author">Johannes Ranke</h4> - <h4 data-toc-skip class="date">Last change on 20 April 2023, -last compiled on 19 Mai 2023</h4> + <h4 data-toc-skip class="date">Last change on 13 February +2023, last compiled on 12 Mai 2025</h4> <small class="dont-index">Source: <a href="https://github.com/jranke/mkin/blob/HEAD/vignettes/prebuilt/2022_cyan_pathway.rmd" class="external-link"><code>vignettes/prebuilt/2022_cyan_pathway.rmd</code></a></small> - <div class="hidden name"><code>2022_cyan_pathway.rmd</code></div> - + <div class="d-none name"><code>2022_cyan_pathway.rmd</code></div> </div> @@ -155,7 +103,7 @@ be fitted with the mkin package.</p> 173340 (Application of nonlinear hierarchical models to the kinetic evaluation of chemical degradation data) of the German Environment Agency carried out in 2022 and 2023.</p> -<p>The mkin package is used in version 1.2.4 which is currently under +<p>The mkin package is used in version 1.2.10 which is currently under development. The newly introduced functionality that is used here is a simplification of excluding random effects for a set of fits based on a related set of fits with a reduced model, and the documentation of the @@ -1638,10 +1586,58 @@ sd(f_JCZ38_qlogis)</td> </tr> </tbody> </table> -<p>The model comparison below suggests that the pathway fits using DFOP +<p>The model comparisons below suggest that the pathway fits using DFOP or SFORB for the parent compound provide the best fit.</p> <div class="sourceCode" id="cb15"><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">|></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> +<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="st">"const"</span><span class="op">]</span><span class="op">)</span> <span class="op">|></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_path_1 const</td> +<td align="right">16</td> +<td align="right">2693.0</td> +<td align="right">2686.8</td> +<td align="right">-1330.5</td> +</tr> +<tr class="even"> +<td align="left">fomc_path_1 const</td> +<td align="right">18</td> +<td align="right">2427.9</td> +<td align="right">2420.9</td> +<td align="right">-1196.0</td> +</tr> +<tr class="odd"> +<td align="left">dfop_path_1 const</td> +<td align="right">20</td> +<td align="right">2403.2</td> +<td align="right">2395.4</td> +<td align="right">-1181.6</td> +</tr> +<tr class="even"> +<td align="left">sforb_path_1 const</td> +<td align="right">20</td> +<td align="right">2401.4</td> +<td align="right">2393.6</td> +<td align="right">-1180.7</td> +</tr> +<tr class="odd"> +<td align="left">hs_path_1 const</td> +<td align="right">20</td> +<td align="right">2427.2</td> +<td align="right">2419.4</td> +<td align="right">-1193.6</td> +</tr> +</tbody> +</table> +<div class="sourceCode" id="cb16"><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="fl">1</span><span class="op">:</span><span class="fl">4</span>, <span class="op">]</span><span class="op">)</span> <span class="op">|></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> @@ -1694,20 +1690,13 @@ or SFORB for the parent compound provide the best fit.</p> <td align="right">-1180.7</td> </tr> <tr class="odd"> -<td align="left">hs_path_1 const</td> -<td align="right">20</td> -<td align="right">2427.2</td> -<td align="right">2419.4</td> -<td align="right">-1193.6</td> -</tr> -<tr class="even"> <td align="left">dfop_path_1 tc</td> <td align="right">20</td> <td align="right">2398.0</td> <td align="right">2390.1</td> <td align="right">-1179.0</td> </tr> -<tr class="odd"> +<tr class="even"> <td align="left">sforb_path_1 tc</td> <td align="right">20</td> <td align="right">2399.9</td> @@ -1718,14 +1707,14 @@ or SFORB for the parent compound provide the best fit.</p> </table> <p>For these two parent model, successful fits are shown below. Plots of the fits with the other parent models are shown in the Appendix.</p> -<div class="sourceCode" id="cb16"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb17"><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">f_saem_1</span><span class="op">[[</span><span class="st">"dfop_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png" alt="DFOP pathway fit with two-component error" width="700"><p class="caption"> DFOP pathway fit with two-component error </p> </div> -<div class="sourceCode" id="cb17"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb18"><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">f_saem_1</span><span class="op">[[</span><span class="st">"sforb_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png" alt="SFORB pathway fit with two-component error" width="700"><p class="caption"> @@ -1735,7 +1724,7 @@ SFORB pathway fit with two-component error <p>A closer graphical analysis of these Figures shows that the residues of transformation product JCZ38 in the soils Tama and Nambsheim observed at later time points are strongly and systematically underestimated.</p> -<div class="sourceCode" id="cb18"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb19"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/parallel/makeCluster.html" class="external-link">stopCluster</a></span><span class="op">(</span><span class="va">cl_path_1</span><span class="op">)</span></span></code></pre></div> </div> <div class="section level3"> @@ -1749,7 +1738,7 @@ corresponding amide (Addendum 2014, p. 109). As FOMC provided the best fit for the parent, and the biexponential models DFOP and SFORB provided the best initial pathway fits, these three parent models are used in the alternative pathway fits.</p> -<div class="sourceCode" id="cb19"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb20"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">cyan_path_2</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> <span> fomc_path_2 <span class="op">=</span> <span class="fu"><a href="../../reference/mkinmod.html">mkinmod</a></span><span class="op">(</span></span> <span> cyan <span class="op">=</span> <span class="fu"><a href="../../reference/mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"FOMC"</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">"JCZ38"</span>, <span class="st">"J9Z38"</span><span class="op">)</span><span class="op">)</span>,</span> @@ -1827,7 +1816,7 @@ alternative pathway fits.</p> </table> <p>Using constant variance, separate fits converge with the exception of the fits to the Sassafras soil data.</p> -<div class="sourceCode" id="cb20"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb21"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">f_sep_2_tc</span> <span class="op"><-</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_2_const</span>, error_model <span class="op">=</span> <span class="st">"tc"</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_sep_2_tc</span><span class="op">)</span> <span class="op">|></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"> @@ -1869,11 +1858,11 @@ the fits to the Sassafras soil data.</p> <p>Using the two-component error model, all separate fits converge with the exception of the alternative pathway fit with DFOP used for the parent and the Sassafras dataset.</p> -<div class="sourceCode" id="cb21"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb22"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">f_saem_2</span> <span class="op"><-</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_2_const</span>, <span class="va">f_sep_2_tc</span><span class="op">)</span>,</span> <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">cyan_saem_full</span><span class="op">[</span><span class="fl">2</span><span class="op">:</span><span class="fl">4</span>, <span class="op">]</span><span class="op">)</span>,</span> <span> cluster <span class="op">=</span> <span class="va">cl_path_2</span><span class="op">)</span></span></code></pre></div> -<div class="sourceCode" id="cb22"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb23"><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_saem_2</span><span class="op">)</span> <span class="op">|></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"> @@ -1900,8 +1889,9 @@ parent and the Sassafras dataset.</p> </tbody> </table> <p>The hierarchical fits for the alternative pathway completed -successfully.</p> -<div class="sourceCode" id="cb23"><pre class="downlit sourceCode r"> +successfully, with the exception of the model using FOMC for the parent +compound and constant variance as the error model.</p> +<div class="sourceCode" id="cb24"><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">|></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> @@ -1932,11 +1922,12 @@ successfully.</p> </tr> </tbody> </table> -<p>In both fits, the random effects for the formation fractions for the -pathways from JCZ38 to JSE76, and for the reverse pathway from JSE76 to -JCZ38 are ill-defined.</p> -<div class="sourceCode" id="cb24"><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="op">|></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> +<p>In all biphasic fits (DFOP or SFORB for the parent compound), the +random effects for the formation fractions for the pathways from JCZ38 +to JSE76, and for the reverse pathway from JSE76 to JCZ38 are +ill-defined.</p> +<div class="sourceCode" id="cb25"><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">"tc"</span><span class="op">]</span><span class="op">)</span> <span class="op">|></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> @@ -1954,27 +1945,54 @@ JCZ38 are ill-defined.</p> <td align="right">-1103.5</td> </tr> <tr class="even"> +<td align="left">dfop_path_2 tc</td> +<td align="right">22</td> +<td align="right">2234.4</td> +<td align="right">2225.8</td> +<td align="right">-1095.2</td> +</tr> +<tr class="odd"> +<td align="left">sforb_path_2 tc</td> +<td align="right">22</td> +<td align="right">2239.7</td> +<td align="right">2231.1</td> +<td align="right">-1097.9</td> +</tr> +</tbody> +</table> +<div class="sourceCode" id="cb26"><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="fl">2</span><span class="op">:</span><span class="fl">3</span>,<span class="op">]</span><span class="op">)</span> <span class="op">|></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">dfop_path_2 const</td> <td align="right">22</td> <td align="right">2288.4</td> <td align="right">2279.8</td> <td align="right">-1122.2</td> </tr> -<tr class="odd"> +<tr class="even"> <td align="left">sforb_path_2 const</td> <td align="right">22</td> <td align="right">2283.3</td> <td align="right">2274.7</td> <td align="right">-1119.7</td> </tr> -<tr class="even"> +<tr class="odd"> <td align="left">dfop_path_2 tc</td> <td align="right">22</td> <td align="right">2234.4</td> <td align="right">2225.8</td> <td align="right">-1095.2</td> </tr> -<tr class="odd"> +<tr class="even"> <td align="left">sforb_path_2 tc</td> <td align="right">22</td> <td align="right">2239.7</td> @@ -1989,21 +2007,21 @@ and BIC values and are plotted below. Compared with the original pathway, the AIC and BIC values indicate a large improvement. This is confirmed by the plots, which show that the metabolite JCZ38 is fitted much better with this model.</p> -<div class="sourceCode" id="cb25"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb27"><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">f_saem_2</span><span class="op">[[</span><span class="st">"fomc_path_2"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png" alt="FOMC pathway fit with two-component error, alternative pathway" width="700"><p class="caption"> FOMC pathway fit with two-component error, alternative pathway </p> </div> -<div class="sourceCode" id="cb26"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb28"><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">f_saem_2</span><span class="op">[[</span><span class="st">"dfop_path_2"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png" alt="DFOP pathway fit with two-component error, alternative pathway" width="700"><p class="caption"> DFOP pathway fit with two-component error, alternative pathway </p> </div> -<div class="sourceCode" id="cb27"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb29"><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">f_saem_2</span><span class="op">[[</span><span class="st">"sforb_path_2"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png" alt="SFORB pathway fit with two-component error, alternative pathway" width="700"><p class="caption"> @@ -2021,7 +2039,7 @@ created below that can be indexed by row and column indices, and which contains the degradation parameter names for which random effects should be excluded for each of the hierarchical fits contained in <code>f_saem_2</code>.</p> -<div class="sourceCode" id="cb28"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb30"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">no_ranef</span> <span class="op"><-</span> <span class="fu"><a href="https://rdrr.io/r/base/matrix.html" class="external-link">matrix</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="op">)</span>, nrow <span class="op">=</span> <span class="fl">3</span>, ncol <span class="op">=</span> <span class="fl">2</span>, dimnames <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/dimnames.html" class="external-link">dimnames</a></span><span class="op">(</span><span class="va">f_saem_2</span><span class="op">)</span><span class="op">)</span></span> <span><span class="va">no_ranef</span><span class="op">[[</span><span class="st">"fomc_path_2"</span>, <span class="st">"const"</span><span class="op">]</span><span class="op">]</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">"log_beta"</span>, <span class="st">"f_JCZ38_qlogis"</span>, <span class="st">"f_JSE76_qlogis"</span><span class="op">)</span></span> <span><span class="va">no_ranef</span><span class="op">[[</span><span class="st">"fomc_path_2"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</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">"cyan_0"</span>, <span class="st">"f_JCZ38_qlogis"</span>, <span class="st">"f_JSE76_qlogis"</span><span class="op">)</span></span> @@ -2036,7 +2054,7 @@ be excluded for each of the hierarchical fits contained in <span><span class="va">f_saem_3</span> <span class="op"><-</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_2</span>,</span> <span> no_random_effect <span class="op">=</span> <span class="va">no_ranef</span>,</span> <span> cluster <span class="op">=</span> <span class="va">cl_path_2</span><span class="op">)</span></span></code></pre></div> -<div class="sourceCode" id="cb29"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb31"><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_saem_3</span><span class="op">)</span> <span class="op">|></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"> @@ -2066,7 +2084,7 @@ be excluded for each of the hierarchical fits contained in all updated fits completed successfully. However, the Fisher Information Matrix for the fixed effects (Fth) could not be inverted, so no confidence intervals for the optimised parameters are available.</p> -<div class="sourceCode" id="cb30"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb32"><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_3</span><span class="op">)</span> <span class="op">|></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"> @@ -2092,8 +2110,8 @@ confidence intervals for the optimised parameters are available.</p> </tr> </tbody> </table> -<div class="sourceCode" id="cb31"><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_3</span><span class="op">)</span> <span class="op">|></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> +<div class="sourceCode" id="cb33"><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_3</span><span class="op">[</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">)</span> <span class="op">|></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> @@ -2111,27 +2129,54 @@ confidence intervals for the optimised parameters are available.</p> <td align="right">-1105.5</td> </tr> <tr class="even"> +<td align="left">dfop_path_2 tc</td> +<td align="right">20</td> +<td align="right">2237.3</td> +<td align="right">2229.5</td> +<td align="right">-1098.6</td> +</tr> +<tr class="odd"> +<td align="left">sforb_path_2 tc</td> +<td align="right">20</td> +<td align="right">2241.3</td> +<td align="right">2233.5</td> +<td align="right">-1100.7</td> +</tr> +</tbody> +</table> +<div class="sourceCode" id="cb34"><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_3</span><span class="op">[</span><span class="fl">2</span><span class="op">:</span><span class="fl">3</span>,<span class="op">]</span><span class="op">)</span> <span class="op">|></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">dfop_path_2 const</td> <td align="right">20</td> <td align="right">2282.2</td> <td align="right">2274.4</td> <td align="right">-1121.1</td> </tr> -<tr class="odd"> +<tr class="even"> <td align="left">sforb_path_2 const</td> <td align="right">20</td> <td align="right">2279.7</td> <td align="right">2271.9</td> <td align="right">-1119.9</td> </tr> -<tr class="even"> +<tr class="odd"> <td align="left">dfop_path_2 tc</td> <td align="right">20</td> <td align="right">2237.3</td> <td align="right">2229.5</td> <td align="right">-1098.6</td> </tr> -<tr class="odd"> +<tr class="even"> <td align="left">sforb_path_2 tc</td> <td align="right">20</td> <td align="right">2241.3</td> @@ -2144,7 +2189,7 @@ confidence intervals for the optimised parameters are available.</p> two-component error) are lower than in the previous fits with the alternative pathway, the practical value of these refined evaluations is limited as no confidence intervals are obtained.</p> -<div class="sourceCode" id="cb32"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb35"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/parallel/makeCluster.html" class="external-link">stopCluster</a></span><span class="op">(</span><span class="va">cl_path_2</span><span class="op">)</span></span></code></pre></div> </div> </div> @@ -2172,21 +2217,21 @@ Agency are gratefully acknowledged.</p> <div class="section level3"> <h3 id="plots-of-fits-that-were-not-refined-further">Plots of fits that were not refined further<a class="anchor" aria-label="anchor" href="#plots-of-fits-that-were-not-refined-further"></a> </h3> -<div class="sourceCode" id="cb33"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb36"><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">f_saem_1</span><span class="op">[[</span><span class="st">"sfo_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png" alt="SFO pathway fit with two-component error" width="700"><p class="caption"> SFO pathway fit with two-component error </p> </div> -<div class="sourceCode" id="cb34"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb37"><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">f_saem_1</span><span class="op">[[</span><span class="st">"fomc_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png" alt="FOMC pathway fit with two-component error" width="700"><p class="caption"> FOMC pathway fit with two-component error </p> </div> -<div class="sourceCode" id="cb35"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb38"><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">f_saem_1</span><span class="op">[[</span><span class="st">"sforb_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png" alt="HS pathway fit with two-component error" width="700"><p class="caption"> @@ -2204,11 +2249,11 @@ HS pathway fit with two-component error Hierarchical SFO path 1 fit with constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:27:54 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:35:03 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - k_cyan * cyan @@ -2221,7 +2266,7 @@ Data: Model predictions using solution type deSolve -Fitted in 422.743 s +Fitted in 445.08 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -2332,11 +2377,11 @@ JSE76 143.91 478.04 Hierarchical SFO path 1 fit with two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:27:49 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:34:47 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - k_cyan * cyan @@ -2349,7 +2394,7 @@ Data: Model predictions using solution type deSolve -Fitted in 417.436 s +Fitted in 428.46 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -2462,11 +2507,11 @@ JSE76 114.20 379.35 Hierarchical FOMC path 1 fit with constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:28:29 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:35:21 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan @@ -2481,7 +2526,7 @@ Data: Model predictions using solution type deSolve -Fitted in 457.122 s +Fitted in 462.739 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -2607,11 +2652,11 @@ JSE76 318.86 1059.22 NA Hierarchical FOMC path 1 fit with two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:28:21 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:35:31 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan @@ -2626,7 +2671,7 @@ Data: Model predictions using solution type deSolve -Fitted in 449.531 s +Fitted in 472.067 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -2745,11 +2790,11 @@ JSE76 210.98 700.9 NA Hierarchical DFOP path 1 fit with constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:29:15 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:36:02 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -2768,7 +2813,7 @@ Data: Model predictions using solution type deSolve -Fitted in 503.737 s +Fitted in 503.671 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -2891,11 +2936,11 @@ JSE76 200.41 665.76 NA NA NA Hierarchical DFOP path 1 fit with two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:31:24 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:38:54 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -2914,7 +2959,7 @@ Data: Model predictions using solution type deSolve -Fitted in 632.55 s +Fitted in 674.81 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -3037,11 +3082,11 @@ JSE76 179.30 595.62 NA NA NA Hierarchical SFORB path 1 fit with constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:29:23 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:36:26 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * @@ -3059,7 +3104,7 @@ Data: Model predictions using solution type deSolve -Fitted in 511.715 s +Fitted in 527.639 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -3202,11 +3247,11 @@ JSE76 193.42 642.53 NA NA NA Hierarchical SFORB path 1 fit with two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:31:23 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:38:31 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * @@ -3224,7 +3269,7 @@ Data: Model predictions using solution type deSolve -Fitted in 630.627 s +Fitted in 651.67 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -3367,11 +3412,11 @@ JSE76 177.11 588.3 NA NA NA Hierarchical HS path 1 fit with constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:28:57 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:36:13 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ifelse(time <= tb, k1, k2) * cyan @@ -3386,7 +3431,7 @@ Data: Model predictions using solution type deSolve -Fitted in 485.304 s +Fitted in 514.353 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -3513,11 +3558,11 @@ JSE76 187.01 621.23 NA NA NA Hierarchical FOMC path 2 fit with two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:39:30 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:47:30 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan @@ -3532,7 +3577,7 @@ Data: Model predictions using solution type deSolve -Fitted in 474.942 s +Fitted in 505.533 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -3581,28 +3626,28 @@ Likelihood computed by importance sampling 2249 2241 -1104 Optimised parameters: - est. lower upper -cyan_0 101.55265 9.920e+01 103.90593 -log_k_JCZ38 -2.32302 -2.832e+00 -1.81416 -log_k_J9Z38 -5.13082 -5.942e+00 -4.31990 -log_k_JSE76 -3.01756 -4.262e+00 -1.77360 -f_cyan_ilr_1 0.70850 3.657e-01 1.05135 -f_cyan_ilr_2 0.95775 2.612e-01 1.65432 -f_JCZ38_qlogis 3.86105 9.248e-01 6.79733 -f_JSE76_qlogis 7.51583 -1.120e+02 127.03921 -log_alpha -0.15308 -4.508e-01 0.14462 -log_beta 2.99165 2.711e+00 3.27202 -a.1 2.04034 1.811e+00 2.26968 -b.1 0.06924 5.745e-02 0.08104 -SD.log_k_JCZ38 0.50818 1.390e-01 0.87736 -SD.log_k_J9Z38 0.86597 2.652e-01 1.46671 -SD.log_k_JSE76 1.38092 4.864e-01 2.27541 -SD.f_cyan_ilr_1 0.38204 1.354e-01 0.62864 -SD.f_cyan_ilr_2 0.55129 7.198e-02 1.03060 -SD.f_JCZ38_qlogis 1.88457 1.710e-02 3.75205 -SD.f_JSE76_qlogis 2.64018 -2.450e+03 2455.27887 -SD.log_alpha 0.31860 1.047e-01 0.53249 -SD.log_beta 0.24195 1.273e-02 0.47117 + est. lower upper +cyan_0 101.55265 9.920e+01 103.9059 +log_k_JCZ38 -2.32302 -2.832e+00 -1.8142 +log_k_J9Z38 -5.13082 -5.942e+00 -4.3199 +log_k_JSE76 -3.01756 -4.262e+00 -1.7736 +f_cyan_ilr_1 0.70850 3.657e-01 1.0513 +f_cyan_ilr_2 0.95775 2.612e-01 1.6543 +f_JCZ38_qlogis 3.86105 9.248e-01 6.7973 +f_JSE76_qlogis 7.51583 -1.120e+02 127.0392 +log_alpha -0.15308 -4.508e-01 0.1446 +log_beta 2.99165 2.711e+00 3.2720 +a.1 2.04034 1.843e+00 2.2382 +b.1 0.06924 5.749e-02 0.0810 +SD.log_k_JCZ38 0.50818 1.390e-01 0.8774 +SD.log_k_J9Z38 0.86597 2.652e-01 1.4667 +SD.log_k_JSE76 1.38092 4.864e-01 2.2754 +SD.f_cyan_ilr_1 0.38204 1.354e-01 0.6286 +SD.f_cyan_ilr_2 0.55129 7.198e-02 1.0306 +SD.f_JCZ38_qlogis 1.88457 1.711e-02 3.7520 +SD.f_JSE76_qlogis 2.64018 -2.450e+03 2454.9447 +SD.log_alpha 0.31860 1.047e-01 0.5325 +SD.log_beta 0.24195 1.273e-02 0.4712 Correlation: cyan_0 l__JCZ3 l__J9Z3 l__JSE7 f_cy__1 f_cy__2 f_JCZ38 f_JSE76 @@ -3633,15 +3678,15 @@ SD.log_k_J9Z38 0.8660 2.652e-01 1.4667 SD.log_k_JSE76 1.3809 4.864e-01 2.2754 SD.f_cyan_ilr_1 0.3820 1.354e-01 0.6286 SD.f_cyan_ilr_2 0.5513 7.198e-02 1.0306 -SD.f_JCZ38_qlogis 1.8846 1.710e-02 3.7520 -SD.f_JSE76_qlogis 2.6402 -2.450e+03 2455.2789 +SD.f_JCZ38_qlogis 1.8846 1.711e-02 3.7520 +SD.f_JSE76_qlogis 2.6402 -2.450e+03 2454.9447 SD.log_alpha 0.3186 1.047e-01 0.5325 SD.log_beta 0.2420 1.273e-02 0.4712 Variance model: - est. lower upper -a.1 2.04034 1.81101 2.26968 -b.1 0.06924 0.05745 0.08104 + est. lower upper +a.1 2.04034 1.84252 2.238 +b.1 0.06924 0.05749 0.081 Backtransformed parameters: est. lower upper @@ -3679,11 +3724,11 @@ JSE76 14.169 47.07 NA Hierarchical DFOP path 2 fit with constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:40:29 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:48:20 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -3703,7 +3748,7 @@ Data: Model predictions using solution type deSolve -Fitted in 533.901 s +Fitted in 555.413 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -3870,11 +3915,11 @@ JSE76 24.32 80.78 NA NA NA Hierarchical DFOP path 2 fit with two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:43:04 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:51:02 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -3894,7 +3939,7 @@ Data: Model predictions using solution type deSolve -Fitted in 688.913 s +Fitted in 717.675 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -3971,15 +4016,15 @@ f_JSE76_qlogis 4.54884 -7.25628 16.35395 log_k1 -2.38201 -2.51639 -2.24763 log_k2 -4.66741 -4.91865 -4.41617 g_qlogis -0.28446 -1.14192 0.57300 -a.1 2.05925 1.83267 2.28582 -b.1 0.06172 0.05076 0.07268 +a.1 2.05925 1.86481 2.25369 +b.1 0.06172 0.05062 0.07282 SD.log_k_JCZ38 0.81137 0.25296 1.36977 -SD.log_k_J9Z38 0.83542 0.25396 1.41689 +SD.log_k_J9Z38 0.83542 0.25395 1.41689 SD.log_k_JSE76 0.97903 0.30100 1.65707 SD.f_cyan_ilr_1 0.37878 0.13374 0.62382 SD.f_cyan_ilr_2 0.67274 0.10102 1.24446 -SD.f_JCZ38_qlogis 1.35327 -0.42361 3.13015 -SD.f_JSE76_qlogis 1.43956 -19.15140 22.03052 +SD.f_JCZ38_qlogis 1.35327 -0.42359 3.13012 +SD.f_JSE76_qlogis 1.43956 -19.14972 22.02884 SD.log_k2 0.25329 0.07521 0.43138 SD.g_qlogis 0.95167 0.35149 1.55184 @@ -4010,19 +4055,19 @@ g_qlogis -0.1656 -0.0928 Random effects: est. lower upper SD.log_k_JCZ38 0.8114 0.25296 1.3698 -SD.log_k_J9Z38 0.8354 0.25396 1.4169 +SD.log_k_J9Z38 0.8354 0.25395 1.4169 SD.log_k_JSE76 0.9790 0.30100 1.6571 SD.f_cyan_ilr_1 0.3788 0.13374 0.6238 SD.f_cyan_ilr_2 0.6727 0.10102 1.2445 -SD.f_JCZ38_qlogis 1.3533 -0.42361 3.1301 -SD.f_JSE76_qlogis 1.4396 -19.15140 22.0305 +SD.f_JCZ38_qlogis 1.3533 -0.42359 3.1301 +SD.f_JSE76_qlogis 1.4396 -19.14972 22.0288 SD.log_k2 0.2533 0.07521 0.4314 SD.g_qlogis 0.9517 0.35149 1.5518 Variance model: est. lower upper -a.1 2.05925 1.83267 2.28582 -b.1 0.06172 0.05076 0.07268 +a.1 2.05925 1.86481 2.25369 +b.1 0.06172 0.05062 0.07282 Backtransformed parameters: est. lower upper @@ -4061,11 +4106,11 @@ JSE76 17.96 59.66 NA NA NA Hierarchical SFORB path 2 fit with constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:40:32 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:48:15 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * @@ -4083,7 +4128,7 @@ Data: Model predictions using solution type deSolve -Fitted in 536.94 s +Fitted in 550.772 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -4257,11 +4302,11 @@ JSE76 14.287 47.46 NA NA NA Hierarchical SFORB path 2 fit with two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:42:47 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 20:50:59 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * @@ -4279,7 +4324,7 @@ Data: Model predictions using solution type deSolve -Fitted in 671.849 s +Fitted in 714.48 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -4358,17 +4403,17 @@ f_cyan_ilr_1 0.70640 3.562e-01 1.057e+00 f_cyan_ilr_2 1.42704 3.170e-01 2.537e+00 f_JCZ38_qlogis 2.84779 1.042e+00 4.654e+00 f_JSE76_qlogis 8.63674 -6.407e+02 6.580e+02 -a.1 2.07082 1.846e+00 2.296e+00 -b.1 0.06227 5.120e-02 7.334e-02 +a.1 2.07082 1.877e+00 2.265e+00 +b.1 0.06227 5.098e-02 7.355e-02 SD.log_k_cyan_free 0.49674 1.865e-01 8.069e-01 -SD.log_k_cyan_bound_free 0.28537 6.808e-02 5.027e-01 +SD.log_k_cyan_bound_free 0.28537 6.809e-02 5.027e-01 SD.log_k_JCZ38 0.74846 2.305e-01 1.266e+00 SD.log_k_J9Z38 0.86077 2.713e-01 1.450e+00 SD.log_k_JSE76 0.97613 3.030e-01 1.649e+00 SD.f_cyan_ilr_1 0.38994 1.382e-01 6.417e-01 SD.f_cyan_ilr_2 0.82869 3.917e-02 1.618e+00 -SD.f_JCZ38_qlogis 1.05000 -2.809e-02 2.128e+00 -SD.f_JSE76_qlogis 0.44681 -3.986e+05 3.986e+05 +SD.f_JCZ38_qlogis 1.05000 -2.808e-02 2.128e+00 +SD.f_JSE76_qlogis 0.44681 -3.985e+05 3.985e+05 Correlation: cyn_f_0 lg_k_c_ lg_k_cyn_f_ lg_k_cyn_b_ l__JCZ3 l__J9Z3 @@ -4397,19 +4442,19 @@ f_JSE76_qlogis 0.0064 0.0345 0.2015 -0.7058 Random effects: est. lower upper SD.log_k_cyan_free 0.4967 1.865e-01 8.069e-01 -SD.log_k_cyan_bound_free 0.2854 6.808e-02 5.027e-01 +SD.log_k_cyan_bound_free 0.2854 6.809e-02 5.027e-01 SD.log_k_JCZ38 0.7485 2.305e-01 1.266e+00 SD.log_k_J9Z38 0.8608 2.713e-01 1.450e+00 SD.log_k_JSE76 0.9761 3.030e-01 1.649e+00 SD.f_cyan_ilr_1 0.3899 1.382e-01 6.417e-01 SD.f_cyan_ilr_2 0.8287 3.917e-02 1.618e+00 -SD.f_JCZ38_qlogis 1.0500 -2.809e-02 2.128e+00 -SD.f_JSE76_qlogis 0.4468 -3.986e+05 3.986e+05 +SD.f_JCZ38_qlogis 1.0500 -2.808e-02 2.128e+00 +SD.f_JSE76_qlogis 0.4468 -3.985e+05 3.985e+05 Variance model: - est. lower upper -a.1 2.07082 1.8458 2.29588 -b.1 0.06227 0.0512 0.07334 + est. lower upper +a.1 2.07082 1.87680 2.26483 +b.1 0.06227 0.05098 0.07355 Backtransformed parameters: est. lower upper @@ -4458,11 +4503,11 @@ Hierarchical FOMC path 2 fit with reduced random effects, two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:55:35 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:04:10 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan @@ -4477,7 +4522,7 @@ Data: Model predictions using solution type deSolve -Fitted in 748.54 s +Fitted in 786.038 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -4526,26 +4571,26 @@ Likelihood computed by importance sampling 2249 2242 -1106 Optimised parameters: - est. lower upper -cyan_0 101.24524 NA NA -log_k_JCZ38 -2.85375 NA NA -log_k_J9Z38 -5.07729 NA NA -log_k_JSE76 -3.53511 NA NA -f_cyan_ilr_1 0.67478 NA NA -f_cyan_ilr_2 0.97152 NA NA -f_JCZ38_qlogis 213.48001 NA NA -f_JSE76_qlogis 2.02040 NA NA -log_alpha -0.11041 NA NA -log_beta 3.06575 NA NA -a.1 2.05279 1.82393 2.28166 -b.1 0.07116 0.05910 0.08322 -SD.log_k_JCZ38 1.21713 0.44160 1.99266 -SD.log_k_J9Z38 0.88268 0.27541 1.48995 -SD.log_k_JSE76 0.59452 0.15005 1.03898 -SD.f_cyan_ilr_1 0.35370 0.12409 0.58331 -SD.f_cyan_ilr_2 0.78186 0.18547 1.37824 -SD.log_alpha 0.27781 0.08168 0.47394 -SD.log_beta 0.32608 0.06490 0.58726 + est. lower upper +cyan_0 101.24524 NA NA +log_k_JCZ38 -2.85375 NA NA +log_k_J9Z38 -5.07729 NA NA +log_k_JSE76 -3.53511 NA NA +f_cyan_ilr_1 0.67478 NA NA +f_cyan_ilr_2 0.97152 NA NA +f_JCZ38_qlogis 213.48001 NA NA +f_JSE76_qlogis 2.02040 NA NA +log_alpha -0.11041 NA NA +log_beta 3.06575 NA NA +a.1 2.05279 1.85495 2.2506 +b.1 0.07116 0.05912 0.0832 +SD.log_k_JCZ38 1.21713 0.44160 1.9927 +SD.log_k_J9Z38 0.88268 0.27541 1.4900 +SD.log_k_JSE76 0.59452 0.15005 1.0390 +SD.f_cyan_ilr_1 0.35370 0.12409 0.5833 +SD.f_cyan_ilr_2 0.78186 0.18547 1.3782 +SD.log_alpha 0.27781 0.08168 0.4739 +SD.log_beta 0.32608 0.06490 0.5873 Correlation is not available @@ -4560,9 +4605,9 @@ SD.log_alpha 0.2778 0.08168 0.4739 SD.log_beta 0.3261 0.06490 0.5873 Variance model: - est. lower upper -a.1 2.05279 1.8239 2.28166 -b.1 0.07116 0.0591 0.08322 + est. lower upper +a.1 2.05279 1.85495 2.2506 +b.1 0.07116 0.05912 0.0832 Backtransformed parameters: est. lower upper @@ -4601,11 +4646,11 @@ Hierarchical DFOP path 2 fit with reduced random effects, constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:57:10 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:05:56 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -4625,7 +4670,7 @@ Data: Model predictions using solution type deSolve -Fitted in 843.793 s +Fitted in 892.139 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -4767,11 +4812,11 @@ Hierarchical DFOP path 2 fit with reduced random effects, two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:57:32 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:02 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -4791,7 +4836,7 @@ Data: Model predictions using solution type deSolve -Fitted in 865.636 s +Fitted in 898.534 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -4868,13 +4913,13 @@ f_JSE76_qlogis 414.80884 NA NA log_k1 -2.38601 NA NA log_k2 -4.63632 NA NA g_qlogis -0.33920 NA NA -a.1 2.10837 1.88051 2.33623 -b.1 0.06223 0.05108 0.07338 +a.1 2.10837 1.91261 2.30413 +b.1 0.06223 0.05085 0.07361 SD.log_k_JCZ38 1.30902 0.48128 2.13675 SD.log_k_J9Z38 0.83882 0.25790 1.41974 SD.log_k_JSE76 0.58104 0.14201 1.02008 SD.f_cyan_ilr_1 0.35421 0.12398 0.58443 -SD.f_cyan_ilr_2 0.79373 0.12007 1.46740 +SD.f_cyan_ilr_2 0.79373 0.12007 1.46739 SD.log_k2 0.27476 0.08557 0.46394 SD.g_qlogis 0.96170 0.35463 1.56878 @@ -4892,8 +4937,8 @@ SD.g_qlogis 0.9617 0.35463 1.5688 Variance model: est. lower upper -a.1 2.10837 1.88051 2.33623 -b.1 0.06223 0.05108 0.07338 +a.1 2.10837 1.91261 2.30413 +b.1 0.06223 0.05085 0.07361 Backtransformed parameters: est. lower upper @@ -4933,11 +4978,11 @@ Hierarchical SFORB path 2 fit with reduced random effects, constant variance </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:57:01 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:05:52 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * @@ -4955,7 +5000,7 @@ Data: Model predictions using solution type deSolve -Fitted in 834.906 s +Fitted in 888.333 s Using 300, 100 iterations and 10 chains Variance model: Constant variance @@ -5104,11 +5149,11 @@ Hierarchical SFORB path 2 fit with reduced random effects, two-component error </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.4 -R version used for fitting: 4.3.0 -Date of fit: Fri May 19 09:57:17 2023 -Date of summary: Fri May 19 09:57:33 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:21 2025 +Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * @@ -5126,7 +5171,7 @@ Data: Model predictions using solution type deSolve -Fitted in 850.751 s +Fitted in 916.619 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function @@ -5205,8 +5250,8 @@ f_cyan_ilr_1 0.72263 NA NA f_cyan_ilr_2 1.45352 NA NA f_JCZ38_qlogis 2.00778 NA NA f_JSE76_qlogis 941.58570 NA NA -a.1 2.11130 1.88299 2.33960 -b.1 0.06299 0.05176 0.07421 +a.1 2.11130 1.91479 2.30780 +b.1 0.06299 0.05152 0.07445 SD.log_k_cyan_free 0.50098 0.18805 0.81390 SD.log_k_cyan_bound_free 0.31671 0.08467 0.54875 SD.log_k_JCZ38 1.25865 0.45932 2.05798 @@ -5220,7 +5265,7 @@ Correlation is not available Random effects: est. lower upper SD.log_k_cyan_free 0.5010 0.18805 0.8139 -SD.log_k_cyan_bound_free 0.3167 0.08467 0.5488 +SD.log_k_cyan_bound_free 0.3167 0.08467 0.5487 SD.log_k_JCZ38 1.2587 0.45932 2.0580 SD.log_k_J9Z38 0.8683 0.27222 1.4644 SD.log_k_JSE76 0.5933 0.14711 1.0394 @@ -5229,8 +5274,8 @@ SD.f_cyan_ilr_2 0.8854 0.13797 1.6329 Variance model: est. lower upper -a.1 2.11130 1.88299 2.33960 -b.1 0.06299 0.05176 0.07421 +a.1 2.11130 1.91479 2.30780 +b.1 0.06299 0.05152 0.07445 Backtransformed parameters: est. lower upper @@ -5275,13 +5320,13 @@ JSE76 25.44 84.51 NA NA NA <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.0 Patched (2023-05-18 r84448) -Platform: x86_64-pc-linux-gnu (64-bit) +<pre><code>R version 4.5.0 (2025-04-11) +Platform: x86_64-pc-linux-gnu Running under: Debian GNU/Linux 12 (bookworm) Matrix products: default -BLAS: /home/jranke/svn/R/r-patched/build/lib/libRblas.so -LAPACK: /usr/lib/x86_64-linux-gnu/openblas-serial/liblapack.so.3; LAPACK version 3.11.0 +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 LAPACK version 3.11.0 locale: [1] LC_CTYPE=de_DE.UTF-8 LC_NUMERIC=C @@ -5299,63 +5344,52 @@ attached base packages: [8] base other attached packages: -[1] saemix_3.2 npde_3.3 knitr_1.42 mkin_1.2.4 +[1] rmarkdown_2.29 nvimcom_0.9-167 saemix_3.3 npde_3.5 +[5] knitr_1.49 mkin_1.2.10 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.21-8 digest_0.6.31 magrittr_2.0.3 evaluate_0.21 - [9] grid_4.3.0 fastmap_1.1.1 cellranger_1.1.0 rprojroot_2.0.3 -[13] jsonlite_1.8.4 processx_3.8.1 pkgbuild_1.4.0 deSolve_1.35 -[17] DBI_1.1.3 mclust_6.0.0 ps_1.7.5 gridExtra_2.3 -[21] purrr_1.0.1 fansi_1.0.4 scales_1.2.1 codetools_0.2-19 -[25] textshaping_0.3.6 jquerylib_0.1.4 cli_3.6.1 crayon_1.5.2 -[29] rlang_1.1.1 munsell_0.5.0 cachem_1.0.8 yaml_2.3.7 -[33] inline_0.3.19 tools_4.3.0 memoise_2.0.1 dplyr_1.1.2 -[37] colorspace_2.1-0 ggplot2_3.4.2 vctrs_0.6.2 R6_2.5.1 -[41] zoo_1.8-12 lifecycle_1.0.3 stringr_1.5.0 fs_1.6.2 -[45] ragg_1.2.5 callr_3.7.3 pkgconfig_2.0.3 desc_1.4.2 -[49] pkgdown_2.0.7 bslib_0.4.2 pillar_1.9.0 gtable_0.3.3 -[53] glue_1.6.2 systemfonts_1.0.4 highr_0.10 xfun_0.39 -[57] tibble_3.2.1 lmtest_0.9-40 tidyselect_1.2.0 htmltools_0.5.5 -[61] nlme_3.1-162 rmarkdown_2.21 compiler_4.3.0 prettyunits_1.1.1 -[65] readxl_1.4.2 </code></pre> + [1] sass_0.4.9 generics_0.1.3 lattice_0.22-6 digest_0.6.37 + [5] magrittr_2.0.3 evaluate_1.0.3 grid_4.5.0 fastmap_1.2.0 + [9] cellranger_1.1.0 jsonlite_1.9.0 processx_3.8.6 pkgbuild_1.4.6 +[13] deSolve_1.40 mclust_6.1.1 ps_1.9.0 gridExtra_2.3 +[17] scales_1.3.0 codetools_0.2-20 textshaping_1.0.0 jquerylib_0.1.4 +[21] cli_3.6.4 rlang_1.1.5 munsell_0.5.1 cachem_1.1.0 +[25] yaml_2.3.10 inline_0.3.21 tools_4.5.0 dplyr_1.1.4 +[29] colorspace_2.1-1 ggplot2_3.5.1 vctrs_0.6.5 R6_2.6.1 +[33] zoo_1.8-13 lifecycle_1.0.4 fs_1.6.5 htmlwidgets_1.6.4 +[37] MASS_7.3-65 ragg_1.3.3 pkgconfig_2.0.3 desc_1.4.3 +[41] callr_3.7.6 pkgdown_2.1.1 pillar_1.10.1 bslib_0.9.0 +[45] gtable_0.3.6 glue_1.8.0 systemfonts_1.2.1 xfun_0.51 +[49] tibble_3.2.1 lmtest_0.9-40 tidyselect_1.2.1 htmltools_0.5.8.1 +[53] nlme_3.1-168 compiler_4.5.0 readxl_1.4.4 </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: AMD Ryzen 9 7950X 16-Core Processor</code></pre> -<pre><code>MemTotal: 64925476 kB</code></pre> +<pre><code>MemTotal: 64927780 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> + </main><aside class="col-md-3"><nav id="toc" aria-label="Table of contents"><h2>On this page</h2> + </nav></aside> </div> -</div> - - <footer><div class="copyright"> - <p></p> -<p>Developed by Johannes Ranke.</p> + <footer><div class="pkgdown-footer-left"> + <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 class="pkgdown-footer-right"> + <p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.1.1.</p> </div> - </footer> + </footer> </div> - - + </body> </html> diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png Binary files differindex b3d6066b..d2201974 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png Binary files differindex 341d48b9..7380ba4c 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png Binary files differindex 92df5b6c..4de15105 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png Binary files differindex 1a17b3bc..d57badf1 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png Binary files differindex 15905501..eb629c4d 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png Binary files differindex cc7d0126..a2abb2f7 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png Binary files differindex 48bfc544..28ec82ba 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png Binary files differindex cc7d0126..a2abb2f7 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-11-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-11-1.pdf Binary files differnew file mode 100644 index 00000000..e38951ad --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-11-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-12-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-12-1.pdf Binary files differnew file mode 100644 index 00000000..3d43286e --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-12-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-13-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-13-1.pdf Binary files differnew file mode 100644 index 00000000..71c3a6a7 --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-13-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-14-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-14-1.pdf Binary files differnew file mode 100644 index 00000000..c7bc3667 --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-14-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-15-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-15-1.pdf Binary files differnew file mode 100644 index 00000000..8317cdd6 --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-15-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-17-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-17-1.pdf Binary files differnew file mode 100644 index 00000000..d270c72e --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-17-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-18-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-18-1.pdf Binary files differnew file mode 100644 index 00000000..75e8718c --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-18-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-19-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-19-1.pdf Binary files differnew file mode 100644 index 00000000..b6e1bf16 --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-19-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-20-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-20-1.pdf Binary files differnew file mode 100644 index 00000000..970cac53 --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-20-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-21-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-21-1.pdf Binary files differnew file mode 100644 index 00000000..9b9f59d0 --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-21-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-22-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-22-1.pdf Binary files differnew file mode 100644 index 00000000..df4c103f --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-22-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-6-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-6-1.pdf Binary files differnew file mode 100644 index 00000000..b0c7f51a --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-6-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-7-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-7-1.pdf Binary files differnew file mode 100644 index 00000000..01c707e3 --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-7-1.pdf diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-8-1.pdf b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-8-1.pdf Binary files differnew file mode 100644 index 00000000..e2f0820d --- /dev/null +++ b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-latex/unnamed-chunk-8-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent.html b/docs/articles/prebuilt/2022_dmta_parent.html index 2da41981..10035c99 100644 --- a/docs/articles/prebuilt/2022_dmta_parent.html +++ b/docs/articles/prebuilt/2022_dmta_parent.html @@ -4,142 +4,90 @@ <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"> +<meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"> <title>Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P • 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 hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P"> -<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]--> +<script src="../../deps/jquery-3.6.0/jquery-3.6.0.min.js"></script><meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"> +<link href="../../deps/bootstrap-5.3.1/bootstrap.min.css" rel="stylesheet"> +<script src="../../deps/bootstrap-5.3.1/bootstrap.bundle.min.js"></script><link href="../../deps/font-awesome-6.5.2/css/all.min.css" rel="stylesheet"> +<link href="../../deps/font-awesome-6.5.2/css/v4-shims.min.css" rel="stylesheet"> +<script src="../../deps/headroom-0.11.0/headroom.min.js"></script><script src="../../deps/headroom-0.11.0/jQuery.headroom.min.js"></script><script src="../../deps/bootstrap-toc-1.0.1/bootstrap-toc.min.js"></script><script src="../../deps/clipboard.js-2.0.11/clipboard.min.js"></script><script src="../../deps/search-1.0.0/autocomplete.jquery.min.js"></script><script src="../../deps/search-1.0.0/fuse.min.js"></script><script src="../../deps/search-1.0.0/mark.min.js"></script><!-- pkgdown --><script src="../../pkgdown.js"></script><meta property="og:title" content="Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P"> </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.4</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> +<body> + <a href="#main" class="visually-hidden-focusable">Skip to contents</a> + + + <nav class="navbar navbar-expand-lg fixed-top bg-light" data-bs-theme="default" aria-label="Site navigation"><div class="container"> + + <a class="navbar-brand me-2" href="../../index.html">mkin</a> + + <small class="nav-text text-muted me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="Released version">1.2.10</small> + + + <button class="navbar-toggler" type="button" data-bs-toggle="collapse" data-bs-target="#navbar" aria-controls="navbar" aria-expanded="false" aria-label="Toggle navigation"> + <span class="navbar-toggler-icon"></span> + </button> + + <div id="navbar" class="collapse navbar-collapse ms-3"> + <ul class="navbar-nav me-auto"> +<li class="nav-item"><a class="nav-link" href="../../reference/index.html">Reference</a></li> +<li class="active nav-item dropdown"> + <button class="nav-link dropdown-toggle" type="button" id="dropdown-articles" data-bs-toggle="dropdown" aria-expanded="false" aria-haspopup="true">Articles</button> + <ul class="dropdown-menu" aria-labelledby="dropdown-articles"> +<li><a class="dropdown-item" href="../../articles/mkin.html">Introduction to mkin</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Example evaluations with (generalised) nonlinear least squares</h6></li> + <li><a class="dropdown-item" href="../../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a></li> + <li><a class="dropdown-item" href="../../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Example evaluations with hierarchical models (nonlinear mixed-effects models)</h6></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_dmta_parent.html">Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_dmta_pathway.html">Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2023_mesotrione_parent.html">Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_cyan_pathway.html">Testing hierarchical pathway kinetics with residue data on cyantraniliprole</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/dimethenamid_2018.html">Comparison of saemix and nlme evaluations of dimethenamid data from 2018</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/multistart.html">Short demo of the multistart method</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Performance</h6></li> + <li><a class="dropdown-item" href="../../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/benchmarks.html">Benchmark timings for mkin</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/saem_benchmarks.html">Benchmark timings for saem.mmkin</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Miscellaneous</h6></li> + <li><a class="dropdown-item" href="../../articles/twa.html">Calculation of time weighted average concentrations with mkin</a></li> + <li><a class="dropdown-item" 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> +<li class="nav-item"><a class="nav-link" href="../../coverage/coverage.html">Test coverage</a></li> +<li class="nav-item"><a class="nav-link" 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 class="navbar-nav"> +<li class="nav-item"><form class="form-inline" role="search"> + <input class="form-control" type="search" name="search-input" id="search-input" autocomplete="off" aria-label="Search site" placeholder="Search for" data-search-index="../../search.json"> +</form></li> +<li class="nav-item"><a class="external-link nav-link" href="https://github.com/jranke/mkin/" aria-label="GitHub"><span class="fa fab fa-github fa-lg"></span></a></li> </ul> </div> -<!--/.nav-collapse --> + + </div> -<!--/.container --> -</div> -<!--/.navbar --> +</nav><div class="container template-article"> + + - - </header><div class="row"> - <div class="col-md-9 contents"> - <div class="page-header toc-ignore"> - <h1 data-toc-skip>Testing hierarchical parent degradation kinetics -with residue data on dimethenamid and dimethenamid-P</h1> +<div class="row"> + <main id="main" class="col-md-9"><div class="page-header"> + + <h1>Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P</h1> <h4 data-toc-skip class="author">Johannes Ranke</h4> <h4 data-toc-skip class="date">Last change on 5 January -2023, last compiled on 19 Mai 2023</h4> +2023, last compiled on 12 Mai 2025</h4> <small class="dont-index">Source: <a href="https://github.com/jranke/mkin/blob/HEAD/vignettes/prebuilt/2022_dmta_parent.rmd" class="external-link"><code>vignettes/prebuilt/2022_dmta_parent.rmd</code></a></small> - <div class="hidden name"><code>2022_dmta_parent.rmd</code></div> - + <div class="d-none name"><code>2022_dmta_parent.rmd</code></div> </div> @@ -154,7 +102,7 @@ FOMC, DFOP and HS can be fitted with the mkin package.</p> 173340 (Application of nonlinear hierarchical models to the kinetic evaluation of chemical degradation data) of the German Environment Agency carried out in 2022 and 2023.</p> -<p>The mkin package is used in version 1.2.4. It contains the test data +<p>The mkin package is used in version 1.2.10. It contains the test data and the functions used in 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> @@ -1005,7 +953,7 @@ updated assuming two-component error.</p> <td align="left">DFOP</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> @@ -1013,7 +961,7 @@ updated assuming two-component error.</p> <tr class="even"> <td align="left">HS</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> @@ -1111,9 +1059,9 @@ the best fits.</p> <tr class="even"> <td align="left">FOMC tc</td> <td align="right">8</td> -<td align="right">720.4</td> -<td align="right">718.8</td> -<td align="right">-352.2</td> +<td align="right">720.7</td> +<td align="right">719.1</td> +<td align="right">-352.4</td> </tr> <tr class="odd"> <td align="left">DFOP const</td> @@ -1132,9 +1080,9 @@ the best fits.</p> <tr class="odd"> <td align="left">DFOP tc</td> <td align="right">10</td> -<td align="right">665.5</td> -<td align="right">663.4</td> -<td align="right">-322.8</td> +<td align="right">665.7</td> +<td align="right">663.6</td> +<td align="right">-322.9</td> </tr> <tr class="even"> <td align="left">HS tc</td> @@ -1215,12 +1163,12 @@ achieved with the argument <code>test = TRUE</code> to the <span> <span class="fu"><a href="https://rdrr.io/pkg/knitr/man/kable.html" class="external-link">kable</a></span><span class="op">(</span>format.args <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>digits <span class="op">=</span> <span class="fl">4</span><span class="op">)</span><span class="op">)</span></span></code></pre></div> <table class="table"> <colgroup> -<col width="37%"> -<col width="6%"> +<col width="38%"> +<col width="7%"> <col width="8%"> <col width="8%"> <col width="9%"> -<col width="9%"> +<col width="8%"> <col width="4%"> <col width="15%"> </colgroup> @@ -1238,8 +1186,8 @@ achieved with the argument <code>test = TRUE</code> to the <tr class="odd"> <td align="left">f_saem_dfop_tc_no_ranef_k2</td> <td align="right">9</td> -<td align="right">663.8</td> -<td align="right">661.9</td> +<td align="right">663.7</td> +<td align="right">661.8</td> <td align="right">-322.9</td> <td align="right">NA</td> <td align="right">NA</td> @@ -1248,12 +1196,12 @@ achieved with the argument <code>test = TRUE</code> to the <tr class="even"> <td align="left">f_saem[[“DFOP”, “tc”]]</td> <td align="right">10</td> -<td align="right">665.5</td> -<td align="right">663.4</td> -<td align="right">-322.8</td> -<td align="right">0.2809</td> +<td align="right">665.7</td> +<td align="right">663.6</td> +<td align="right">-322.9</td> +<td align="right">0</td> +<td align="right">1</td> <td align="right">1</td> -<td align="right">0.5961</td> </tr> </tbody> </table> @@ -1285,11 +1233,11 @@ Plot of the final NLHM DFOP fit <p>Finally, a summary report of the fit is produced.</p> <div class="sourceCode" id="cb14"><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_dfop_tc_no_ranef_k2</span><span class="op">)</span></span></code></pre></div> -<pre><code>saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:09 2023 -Date of summary: Thu Apr 20 14:07:10 2023 +<pre><code>saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:46 2025 +Date of summary: Mon May 12 21:06:46 2025 Equations: d_DMTA/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -1301,21 +1249,21 @@ Data: Model predictions using solution type analytical -Fitted in 4.175 s +Fitted in 3.953 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function Starting values for degradation parameters: - DMTA_0 k1 k2 g -98.759266 0.087034 0.009933 0.930827 + DMTA_0 k1 k2 g +98.71186 0.08675 0.01374 0.93491 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k1 k2 g -DMTA_0 98.76 0 0 0 +DMTA_0 98.71 0 0 0 k1 0.00 1 0 0 k2 0.00 0 1 0 g 0.00 0 0 1 @@ -1328,40 +1276,40 @@ Results: Likelihood computed by importance sampling AIC BIC logLik - 663.8 661.9 -322.9 + 663.7 661.8 -322.9 Optimised parameters: est. lower upper -DMTA_0 98.228939 96.285869 100.17201 -k1 0.064063 0.033477 0.09465 -k2 0.008297 0.005824 0.01077 -g 0.953821 0.914328 0.99331 -a.1 1.068479 0.869538 1.26742 -b.1 0.029424 0.022406 0.03644 -SD.DMTA_0 2.030437 0.404824 3.65605 -SD.k1 0.594692 0.256660 0.93272 -SD.g 1.006754 0.361327 1.65218 +DMTA_0 98.256267 96.286112 100.22642 +k1 0.064037 0.033281 0.09479 +k2 0.008469 0.006002 0.01094 +g 0.954167 0.914460 0.99387 +a.1 1.061795 0.878608 1.24498 +b.1 0.029550 0.022593 0.03651 +SD.DMTA_0 2.068581 0.427178 3.70998 +SD.k1 0.598285 0.258235 0.93833 +SD.g 1.016689 0.360061 1.67332 Correlation: DMTA_0 k1 k2 -k1 0.0218 -k2 0.0556 0.0355 -g -0.0516 -0.0284 -0.2800 +k1 0.0213 +k2 0.0541 0.0344 +g -0.0521 -0.0286 -0.2744 Random effects: est. lower upper -SD.DMTA_0 2.0304 0.4048 3.6560 -SD.k1 0.5947 0.2567 0.9327 -SD.g 1.0068 0.3613 1.6522 +SD.DMTA_0 2.0686 0.4272 3.7100 +SD.k1 0.5983 0.2582 0.9383 +SD.g 1.0167 0.3601 1.6733 Variance model: est. lower upper -a.1 1.06848 0.86954 1.26742 -b.1 0.02942 0.02241 0.03644 +a.1 1.06180 0.87861 1.24498 +b.1 0.02955 0.02259 0.03651 Estimated disappearance times: - DT50 DT90 DT50back DT50_k1 DT50_k2 -DMTA 11.45 41.4 12.46 10.82 83.54</code></pre> + DT50 DT90 DT50back DT50_k1 DT50_k2 +DMTA 11.45 41.32 12.44 10.82 81.85</code></pre> </div> <div class="section level3"> <h3 id="alternative-check-of-parameter-identifiability">Alternative check of parameter identifiability<a class="anchor" aria-label="anchor" href="#alternative-check-of-parameter-identifiability"></a> @@ -1461,11 +1409,11 @@ Erythropoiesis.”</span> <em>BMC Bioinformatics</em> 22 (478). <a href="https:/ Hierarchical mkin fit of the SFO model with error model const </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:02 2023 -Date of summary: Thu Apr 20 14:08:16 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:38 2025 +Date of summary: Mon May 12 21:07:51 2025 Equations: d_DMTA/dt = - k_DMTA * DMTA @@ -1475,7 +1423,7 @@ Data: Model predictions using solution type analytical -Fitted in 0.982 s +Fitted in 0.793 s Using 300, 100 iterations and 9 chains Variance model: Constant variance @@ -1533,11 +1481,11 @@ DMTA 12.24 40.65 Hierarchical mkin fit of the SFO model with error model tc </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:03 2023 -Date of summary: Thu Apr 20 14:08:16 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:40 2025 +Date of summary: Mon May 12 21:07:51 2025 Equations: d_DMTA/dt = - k_DMTA * DMTA @@ -1547,7 +1495,7 @@ Data: Model predictions using solution type analytical -Fitted in 2.398 s +Fitted in 2.462 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function @@ -1575,13 +1523,13 @@ Likelihood computed by importance sampling 798.3 797.1 -393.2 Optimised parameters: - est. lower upper -DMTA_0 97.271822 95.703157 98.84049 -k_DMTA 0.056638 0.029110 0.08417 -a.1 2.660081 2.230398 3.08976 -b.1 0.001665 -0.006911 0.01024 -SD.DMTA_0 1.545520 0.145035 2.94601 -SD.k_DMTA 0.606422 0.262274 0.95057 + est. lower upper +DMTA_0 97.271822 95.70316 98.84049 +k_DMTA 0.056638 0.02911 0.08417 +a.1 2.660081 2.27492 3.04525 +b.1 0.001665 -0.14451 0.14784 +SD.DMTA_0 1.545520 0.14301 2.94803 +SD.k_DMTA 0.606422 0.26227 0.95057 Correlation: DMTA_0 @@ -1589,13 +1537,13 @@ k_DMTA 0.0169 Random effects: est. lower upper -SD.DMTA_0 1.5455 0.1450 2.9460 +SD.DMTA_0 1.5455 0.1430 2.9480 SD.k_DMTA 0.6064 0.2623 0.9506 Variance model: - est. lower upper -a.1 2.660081 2.230398 3.08976 -b.1 0.001665 -0.006911 0.01024 + est. lower upper +a.1 2.660081 2.2749 3.0452 +b.1 0.001665 -0.1445 0.1478 Estimated disappearance times: DT50 DT90 @@ -1607,11 +1555,11 @@ DMTA 12.24 40.65 Hierarchical mkin fit of the FOMC model with error model const </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:02 2023 -Date of summary: Thu Apr 20 14:08:16 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:39 2025 +Date of summary: Mon May 12 21:07:51 2025 Equations: d_DMTA/dt = - (alpha/beta) * 1/((time/beta) + 1) * DMTA @@ -1621,7 +1569,7 @@ Data: Model predictions using solution type analytical -Fitted in 1.398 s +Fitted in 1.408 s Using 300, 100 iterations and 9 chains Variance model: Constant variance @@ -1653,7 +1601,7 @@ Optimised parameters: est. lower upper DMTA_0 98.3435 96.9033 99.784 alpha 7.2007 2.5889 11.812 -beta 112.8746 34.8816 190.868 +beta 112.8745 34.8816 190.867 a.1 2.0459 1.8054 2.286 SD.DMTA_0 1.4795 0.2717 2.687 SD.alpha 0.6396 0.1509 1.128 @@ -1684,11 +1632,11 @@ DMTA 11.41 42.53 12.8 Hierarchical mkin fit of the FOMC model with error model tc </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:04 2023 -Date of summary: Thu Apr 20 14:08:16 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:40 2025 +Date of summary: Mon May 12 21:07:51 2025 Equations: d_DMTA/dt = - (alpha/beta) * 1/((time/beta) + 1) * DMTA @@ -1698,7 +1646,7 @@ Data: Model predictions using solution type analytical -Fitted in 3.044 s +Fitted in 2.735 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function @@ -1724,38 +1672,38 @@ Results: Likelihood computed by importance sampling AIC BIC logLik - 720.4 718.8 -352.2 + 720.7 719.1 -352.4 Optimised parameters: est. lower upper -DMTA_0 98.99136 97.26011 100.72261 -alpha 5.86312 2.57485 9.15138 -beta 88.55571 29.20889 147.90254 -a.1 1.51063 1.24384 1.77741 -b.1 0.02824 0.02040 0.03609 -SD.DMTA_0 1.57436 -0.04867 3.19739 -SD.alpha 0.59871 0.17132 1.02611 -SD.beta 0.72994 0.22849 1.23139 +DMTA_0 99.10577 97.33296 100.87859 +alpha 5.46260 2.52199 8.40321 +beta 81.66080 30.46664 132.85497 +a.1 1.50219 1.25801 1.74636 +b.1 0.02893 0.02048 0.03739 +SD.DMTA_0 1.61887 -0.03843 3.27618 +SD.alpha 0.58145 0.17364 0.98925 +SD.beta 0.68205 0.21108 1.15302 Correlation: DMTA_0 alpha -alpha -0.1363 -beta -0.1414 0.2542 +alpha -0.1321 +beta -0.1430 0.2467 Random effects: - est. lower upper -SD.DMTA_0 1.5744 -0.04867 3.197 -SD.alpha 0.5987 0.17132 1.026 -SD.beta 0.7299 0.22849 1.231 + est. lower upper +SD.DMTA_0 1.6189 -0.03843 3.2762 +SD.alpha 0.5814 0.17364 0.9892 +SD.beta 0.6821 0.21108 1.1530 Variance model: - est. lower upper -a.1 1.51063 1.2438 1.77741 -b.1 0.02824 0.0204 0.03609 + est. lower upper +a.1 1.50219 1.25801 1.74636 +b.1 0.02893 0.02048 0.03739 Estimated disappearance times: - DT50 DT90 DT50back -DMTA 11.11 42.6 12.82 + DT50 DT90 DT50back +DMTA 11.05 42.81 12.89 </code></pre> <p></p> @@ -1763,11 +1711,11 @@ DMTA 11.11 42.6 12.82 Hierarchical mkin fit of the DFOP model with error model const </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:02 2023 -Date of summary: Thu Apr 20 14:08:16 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:39 2025 +Date of summary: Mon May 12 21:07:51 2025 Equations: d_DMTA/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -1779,7 +1727,7 @@ Data: Model predictions using solution type analytical -Fitted in 1.838 s +Fitted in 1.769 s Using 300, 100 iterations and 9 chains Variance model: Constant variance @@ -1810,10 +1758,10 @@ Likelihood computed by importance sampling Optimised parameters: est. lower upper -DMTA_0 98.092481 96.573898 99.61106 +DMTA_0 98.092481 96.573899 99.61106 k1 0.062499 0.030336 0.09466 k2 0.009065 -0.005133 0.02326 -g 0.948967 0.862079 1.03586 +g 0.948967 0.862080 1.03586 a.1 1.821671 1.604774 2.03857 SD.DMTA_0 1.677785 0.472066 2.88350 SD.k1 0.634962 0.270788 0.99914 @@ -1847,11 +1795,11 @@ DMTA 11.79 42.8 12.88 11.09 76.46 Hierarchical mkin fit of the DFOP model with error model tc </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:04 2023 -Date of summary: Thu Apr 20 14:08:16 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:41 2025 +Date of summary: Mon May 12 21:07:51 2025 Equations: d_DMTA/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -1863,21 +1811,21 @@ Data: Model predictions using solution type analytical -Fitted in 3.297 s +Fitted in 3.069 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function Starting values for degradation parameters: - DMTA_0 k1 k2 g -98.759266 0.087034 0.009933 0.930827 + DMTA_0 k1 k2 g +98.71186 0.08675 0.01374 0.93491 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k1 k2 g -DMTA_0 98.76 0 0 0 +DMTA_0 98.71 0 0 0 k1 0.00 1 0 0 k2 0.00 0 1 0 g 0.00 0 0 1 @@ -1890,42 +1838,42 @@ Results: Likelihood computed by importance sampling AIC BIC logLik - 665.5 663.4 -322.8 + 665.7 663.6 -322.9 Optimised parameters: est. lower upper -DMTA_0 98.377019 96.447952 100.30609 -k1 0.064843 0.034607 0.09508 -k2 0.008895 0.006368 0.01142 -g 0.949696 0.903815 0.99558 -a.1 1.065241 0.865754 1.26473 -b.1 0.029340 0.022336 0.03634 -SD.DMTA_0 2.007754 0.387982 3.62753 -SD.k1 0.580473 0.250286 0.91066 -SD.k2 0.006105 -4.920337 4.93255 -SD.g 1.097149 0.412779 1.78152 +DMTA_0 98.347470 96.380815 100.31413 +k1 0.064524 0.034279 0.09477 +k2 0.008304 0.005843 0.01076 +g 0.952128 0.909578 0.99468 +a.1 1.068907 0.883665 1.25415 +b.1 0.029265 0.022318 0.03621 +SD.DMTA_0 2.065796 0.427951 3.70364 +SD.k1 0.583703 0.251796 0.91561 +SD.k2 0.004167 -7.831228 7.83956 +SD.g 1.064450 0.397479 1.73142 Correlation: DMTA_0 k1 k2 -k1 0.0235 -k2 0.0595 0.0424 -g -0.0470 -0.0278 -0.2731 +k1 0.0223 +k2 0.0568 0.0394 +g -0.0464 -0.0269 -0.2713 Random effects: est. lower upper -SD.DMTA_0 2.007754 0.3880 3.6275 -SD.k1 0.580473 0.2503 0.9107 -SD.k2 0.006105 -4.9203 4.9325 -SD.g 1.097149 0.4128 1.7815 +SD.DMTA_0 2.065796 0.4280 3.7036 +SD.k1 0.583703 0.2518 0.9156 +SD.k2 0.004167 -7.8312 7.8396 +SD.g 1.064450 0.3975 1.7314 Variance model: est. lower upper -a.1 1.06524 0.86575 1.26473 -b.1 0.02934 0.02234 0.03634 +a.1 1.06891 0.88367 1.25415 +b.1 0.02927 0.02232 0.03621 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 -DMTA 11.36 41.32 12.44 10.69 77.92 +DMTA 11.39 41.36 12.45 10.74 83.48 </code></pre> <p></p> @@ -1933,11 +1881,11 @@ DMTA 11.36 41.32 12.44 10.69 77.92 Hierarchical mkin fit of the HS model with error model const </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:03 2023 -Date of summary: Thu Apr 20 14:08:16 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:39 2025 +Date of summary: Mon May 12 21:07:51 2025 Equations: d_DMTA/dt = - ifelse(time <= tb, k1, k2) * DMTA @@ -1947,7 +1895,7 @@ Data: Model predictions using solution type analytical -Fitted in 1.972 s +Fitted in 1.996 s Using 300, 100 iterations and 9 chains Variance model: Constant variance @@ -2015,11 +1963,11 @@ DMTA 8.801 67.91 20.44 8.801 31.13 Hierarchical mkin fit of the HS model with error model tc </caption> <pre><code> -saemix version used for fitting: 3.2 -mkin version used for pre-fitting: 1.2.3 -R version used for fitting: 4.2.3 -Date of fit: Thu Apr 20 14:07:04 2023 -Date of summary: Thu Apr 20 14:08:16 2023 +saemix version used for fitting: 3.3 +mkin version used for pre-fitting: 1.2.10 +R version used for fitting: 4.5.0 +Date of fit: Mon May 12 21:06:41 2025 +Date of summary: Mon May 12 21:07:51 2025 Equations: d_DMTA/dt = - ifelse(time <= tb, k1, k2) * DMTA @@ -2029,7 +1977,7 @@ Data: Model predictions using solution type analytical -Fitted in 3.378 s +Fitted in 3.389 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function @@ -2060,16 +2008,16 @@ Likelihood computed by importance sampling Optimised parameters: est. lower upper -DMTA_0 97.76570 95.81350 99.71791 +DMTA_0 97.76571 95.81350 99.71791 k1 0.05855 0.03080 0.08630 k2 0.02337 0.01664 0.03010 tb 31.09638 29.38289 32.80987 -a.1 1.08835 0.88590 1.29080 -b.1 0.02964 0.02257 0.03671 -SD.DMTA_0 2.04877 0.42607 3.67147 +a.1 1.08835 0.90059 1.27611 +b.1 0.02964 0.02261 0.03667 +SD.DMTA_0 2.04877 0.42553 3.67200 SD.k1 0.59166 0.25621 0.92711 SD.k2 0.30698 0.09561 0.51835 -SD.tb 0.01274 -0.10914 0.13462 +SD.tb 0.01274 -0.10915 0.13464 Correlation: DMTA_0 k1 k2 @@ -2079,15 +2027,15 @@ tb -0.0668 -0.0103 -0.2013 Random effects: est. lower upper -SD.DMTA_0 2.04877 0.42607 3.6715 +SD.DMTA_0 2.04877 0.42553 3.6720 SD.k1 0.59166 0.25621 0.9271 SD.k2 0.30698 0.09561 0.5183 -SD.tb 0.01274 -0.10914 0.1346 +SD.tb 0.01274 -0.10915 0.1346 Variance model: est. lower upper -a.1 1.08835 0.88590 1.29080 -b.1 0.02964 0.02257 0.03671 +a.1 1.08835 0.90059 1.27611 +b.1 0.02964 0.02261 0.03667 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 @@ -2143,13 +2091,13 @@ Convergence plot for the NLHM HS fit with two-component error <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.2.3 (2023-03-15) -Platform: x86_64-pc-linux-gnu (64-bit) +<pre><code>R version 4.5.0 (2025-04-11) +Platform: x86_64-pc-linux-gnu Running under: Debian GNU/Linux 12 (bookworm) Matrix products: default -BLAS: /usr/lib/x86_64-linux-gnu/openblas-serial/libblas.so.3 -LAPACK: /usr/lib/x86_64-linux-gnu/openblas-serial/libopenblas-r0.3.21.so +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 LAPACK version 3.11.0 locale: [1] LC_CTYPE=de_DE.UTF-8 LC_NUMERIC=C @@ -2159,65 +2107,58 @@ locale: [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.42 mkin_1.2.3 +[1] rmarkdown_2.29 nvimcom_0.9-167 saemix_3.3 npde_3.5 +[5] knitr_1.49 mkin_1.2.10 loaded via a namespace (and not attached): - [1] highr_0.10 pillar_1.9.0 bslib_0.4.2 compiler_4.2.3 - [5] jquerylib_0.1.4 tools_4.2.3 mclust_6.0.0 digest_0.6.31 - [9] tibble_3.2.1 jsonlite_1.8.4 evaluate_0.20 memoise_2.0.1 -[13] lifecycle_1.0.3 nlme_3.1-162 gtable_0.3.3 lattice_0.21-8 -[17] pkgconfig_2.0.3 rlang_1.1.0 DBI_1.1.3 cli_3.6.1 -[21] yaml_2.3.7 pkgdown_2.0.7 xfun_0.38 fastmap_1.1.1 -[25] gridExtra_2.3 dplyr_1.1.1 stringr_1.5.0 generics_0.1.3 -[29] desc_1.4.2 fs_1.6.1 vctrs_0.6.1 sass_0.4.5 -[33] systemfonts_1.0.4 tidyselect_1.2.0 rprojroot_2.0.3 lmtest_0.9-40 -[37] grid_4.2.3 glue_1.6.2 R6_2.5.1 textshaping_0.3.6 -[41] fansi_1.0.4 rmarkdown_2.21 purrr_1.0.1 ggplot2_3.4.2 -[45] magrittr_2.0.3 codetools_0.2-19 scales_1.2.1 htmltools_0.5.5 -[49] colorspace_2.1-0 ragg_1.2.5 utf8_1.2.3 stringi_1.7.12 -[53] munsell_0.5.0 cachem_1.0.7 zoo_1.8-12 </code></pre> + [1] gtable_0.3.6 jsonlite_1.9.0 dplyr_1.1.4 compiler_4.5.0 + [5] tidyselect_1.2.1 gridExtra_2.3 jquerylib_0.1.4 systemfonts_1.2.1 + [9] scales_1.3.0 textshaping_1.0.0 yaml_2.3.10 fastmap_1.2.0 +[13] lattice_0.22-6 ggplot2_3.5.1 R6_2.6.1 generics_0.1.3 +[17] lmtest_0.9-40 MASS_7.3-65 htmlwidgets_1.6.4 tibble_3.2.1 +[21] desc_1.4.3 munsell_0.5.1 bslib_0.9.0 pillar_1.10.1 +[25] rlang_1.1.5 cachem_1.1.0 xfun_0.51 fs_1.6.5 +[29] sass_0.4.9 cli_3.6.4 pkgdown_2.1.1 magrittr_2.0.3 +[33] digest_0.6.37 grid_4.5.0 mclust_6.1.1 lifecycle_1.0.4 +[37] nlme_3.1-168 vctrs_0.6.5 evaluate_1.0.3 glue_1.8.0 +[41] codetools_0.2-20 ragg_1.3.3 zoo_1.8-13 colorspace_2.1-1 +[45] tools_4.5.0 pkgconfig_2.0.3 htmltools_0.5.8.1</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: AMD Ryzen 9 7950X 16-Core Processor</code></pre> -<pre><code>MemTotal: 64936316 kB</code></pre> +<pre><code>MemTotal: 64927780 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> - + </main><aside class="col-md-3"><nav id="toc" aria-label="Table of contents"><h2>On this page</h2> + </nav></aside> </div> - <footer><div class="copyright"> - <p></p> -<p>Developed by Johannes Ranke.</p> + <footer><div class="pkgdown-footer-left"> + <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 class="pkgdown-footer-right"> + <p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.1.1.</p> </div> - </footer> + </footer> </div> - - + </body> </html> diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-const-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-const-1.png Binary files differindex 3f145074..91c027a5 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-const-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-const-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-tc-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-tc-1.png Binary files differindex e5457fc9..678b76e2 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-tc-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-tc-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-tc-no-ranef-k2-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-tc-no-ranef-k2-1.png Binary files differindex 14707641..a3bc320e 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-tc-no-ranef-k2-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-dfop-tc-no-ranef-k2-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-fomc-const-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-fomc-const-1.png Binary files differindex c7ed69a3..da63e191 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-fomc-const-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-fomc-const-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-fomc-tc-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-fomc-tc-1.png Binary files differindex 1a48524c..00958654 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-fomc-tc-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-fomc-tc-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-hs-const-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-hs-const-1.png Binary files differindex 0f3b1184..72a996ac 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-hs-const-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-hs-const-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-hs-tc-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-hs-tc-1.png Binary files differindex 901a1579..46cfa3f8 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-hs-tc-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-hs-tc-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-sfo-const-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-sfo-const-1.png Binary files differindex a3e3a51f..6dc7a6e2 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-sfo-const-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-sfo-const-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-sfo-tc-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-sfo-tc-1.png Binary files differindex b85691eb..8df26d78 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-sfo-tc-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/convergence-saem-sfo-tc-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-full-par-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-full-par-1.png Binary files differindex a42950f0..f90b19e7 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-full-par-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-full-par-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-reduced-par-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-reduced-par-1.png Binary files differindex caebc768..f6f1e43d 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-reduced-par-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-reduced-par-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-reduced-par-llquant-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-reduced-par-llquant-1.png Binary files differindex 45ae57f1..ef6bc8bd 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-reduced-par-llquant-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/multistart-reduced-par-llquant-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/plot-saem-dfop-tc-no-ranef-k2-1.png b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/plot-saem-dfop-tc-no-ranef-k2-1.png Binary files differindex 1f8eb9f0..cd506f57 100644 --- a/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/plot-saem-dfop-tc-no-ranef-k2-1.png +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-html/plot-saem-dfop-tc-no-ranef-k2-1.png diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-const-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-const-1.pdf Binary files differnew file mode 100644 index 00000000..20f31bba --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-const-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-tc-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-tc-1.pdf Binary files differnew file mode 100644 index 00000000..079a4547 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-tc-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-tc-no-ranef-k2-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-tc-no-ranef-k2-1.pdf Binary files differnew file mode 100644 index 00000000..d732c9d6 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-dfop-tc-no-ranef-k2-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-fomc-const-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-fomc-const-1.pdf Binary files differnew file mode 100644 index 00000000..8a7cf103 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-fomc-const-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-fomc-tc-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-fomc-tc-1.pdf Binary files differnew file mode 100644 index 00000000..7556db36 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-fomc-tc-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-hs-const-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-hs-const-1.pdf Binary files differnew file mode 100644 index 00000000..6e8b9df1 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-hs-const-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-hs-tc-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-hs-tc-1.pdf Binary files differnew file mode 100644 index 00000000..08da9919 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-hs-tc-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-sfo-const-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-sfo-const-1.pdf Binary files differnew file mode 100644 index 00000000..ba3698f3 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-sfo-const-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-sfo-tc-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-sfo-tc-1.pdf Binary files differnew file mode 100644 index 00000000..7baff5bb --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/convergence-saem-sfo-tc-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-full-par-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-full-par-1.pdf Binary files differnew file mode 100644 index 00000000..3621b5a5 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-full-par-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-reduced-par-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-reduced-par-1.pdf Binary files differnew file mode 100644 index 00000000..93d185be --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-reduced-par-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-reduced-par-llquant-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-reduced-par-llquant-1.pdf Binary files differnew file mode 100644 index 00000000..31818088 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/multistart-reduced-par-llquant-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/plot-saem-dfop-tc-no-ranef-k2-1.pdf b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/plot-saem-dfop-tc-no-ranef-k2-1.pdf Binary files differnew file mode 100644 index 00000000..607a03ed --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_parent_files/figure-latex/plot-saem-dfop-tc-no-ranef-k2-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_pathway.html b/docs/articles/prebuilt/2022_dmta_pathway.html index ea4bd087..d2ef55fb 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway.html +++ b/docs/articles/prebuilt/2022_dmta_pathway.html @@ -4,142 +4,90 @@ <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"> +<meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"> <title>Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P • 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 hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P"> -<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]--> +<script src="../../deps/jquery-3.6.0/jquery-3.6.0.min.js"></script><meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"> +<link href="../../deps/bootstrap-5.3.1/bootstrap.min.css" rel="stylesheet"> +<script src="../../deps/bootstrap-5.3.1/bootstrap.bundle.min.js"></script><link href="../../deps/font-awesome-6.5.2/css/all.min.css" rel="stylesheet"> +<link href="../../deps/font-awesome-6.5.2/css/v4-shims.min.css" rel="stylesheet"> +<script src="../../deps/headroom-0.11.0/headroom.min.js"></script><script src="../../deps/headroom-0.11.0/jQuery.headroom.min.js"></script><script src="../../deps/bootstrap-toc-1.0.1/bootstrap-toc.min.js"></script><script src="../../deps/clipboard.js-2.0.11/clipboard.min.js"></script><script src="../../deps/search-1.0.0/autocomplete.jquery.min.js"></script><script src="../../deps/search-1.0.0/fuse.min.js"></script><script src="../../deps/search-1.0.0/mark.min.js"></script><!-- pkgdown --><script src="../../pkgdown.js"></script><meta property="og:title" content="Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P"> </head> -<body data-spy="scroll" data-target="#toc"> - +<body> + <a href="#main" class="visually-hidden-focusable">Skip to contents</a> - <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.4</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> + <nav class="navbar navbar-expand-lg fixed-top bg-light" data-bs-theme="default" aria-label="Site navigation"><div class="container"> + + <a class="navbar-brand me-2" href="../../index.html">mkin</a> + + <small class="nav-text text-muted me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="Released version">1.2.10</small> + + + <button class="navbar-toggler" type="button" data-bs-toggle="collapse" data-bs-target="#navbar" aria-controls="navbar" aria-expanded="false" aria-label="Toggle navigation"> + <span class="navbar-toggler-icon"></span> + </button> + + <div id="navbar" class="collapse navbar-collapse ms-3"> + <ul class="navbar-nav me-auto"> +<li class="nav-item"><a class="nav-link" href="../../reference/index.html">Reference</a></li> +<li class="active nav-item dropdown"> + <button class="nav-link dropdown-toggle" type="button" id="dropdown-articles" data-bs-toggle="dropdown" aria-expanded="false" aria-haspopup="true">Articles</button> + <ul class="dropdown-menu" aria-labelledby="dropdown-articles"> +<li><a class="dropdown-item" href="../../articles/mkin.html">Introduction to mkin</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Example evaluations with (generalised) nonlinear least squares</h6></li> + <li><a class="dropdown-item" href="../../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a></li> + <li><a class="dropdown-item" href="../../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Example evaluations with hierarchical models (nonlinear mixed-effects models)</h6></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_dmta_parent.html">Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_dmta_pathway.html">Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2023_mesotrione_parent.html">Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_cyan_pathway.html">Testing hierarchical pathway kinetics with residue data on cyantraniliprole</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/dimethenamid_2018.html">Comparison of saemix and nlme evaluations of dimethenamid data from 2018</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/multistart.html">Short demo of the multistart method</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Performance</h6></li> + <li><a class="dropdown-item" href="../../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/benchmarks.html">Benchmark timings for mkin</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/saem_benchmarks.html">Benchmark timings for saem.mmkin</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Miscellaneous</h6></li> + <li><a class="dropdown-item" href="../../articles/twa.html">Calculation of time weighted average concentrations with mkin</a></li> + <li><a class="dropdown-item" 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> +<li class="nav-item"><a class="nav-link" href="../../coverage/coverage.html">Test coverage</a></li> +<li class="nav-item"><a class="nav-link" 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 class="navbar-nav"> +<li class="nav-item"><form class="form-inline" role="search"> + <input class="form-control" type="search" name="search-input" id="search-input" autocomplete="off" aria-label="Search site" placeholder="Search for" data-search-index="../../search.json"> +</form></li> +<li class="nav-item"><a class="external-link nav-link" href="https://github.com/jranke/mkin/" aria-label="GitHub"><span class="fa fab fa-github fa-lg"></span></a></li> </ul> </div> -<!--/.nav-collapse --> + + </div> -<!--/.container --> -</div> -<!--/.navbar --> +</nav><div class="container template-article"> - - </header><div class="row"> - <div class="col-md-9 contents"> - <div class="page-header toc-ignore"> - <h1 data-toc-skip>Testing hierarchical pathway kinetics with -residue data on dimethenamid and dimethenamid-P</h1> + + +<div class="row"> + <main id="main" class="col-md-9"><div class="page-header"> + + <h1>Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</h1> <h4 data-toc-skip class="author">Johannes Ranke</h4> <h4 data-toc-skip class="date">Last change on 20 April 2023, -last compiled on 19 Mai 2023</h4> +last compiled on 12 Mai 2025</h4> <small class="dont-index">Source: <a href="https://github.com/jranke/mkin/blob/HEAD/vignettes/prebuilt/2022_dmta_pathway.rmd" class="external-link"><code>vignettes/prebuilt/2022_dmta_pathway.rmd</code></a></small> - <div class="hidden name"><code>2022_dmta_pathway.rmd</code></div> - + <div class="d-none name"><code>2022_dmta_pathway.rmd</code></div> </div> @@ -155,7 +103,7 @@ can be fitted with the mkin package.</p> 173340 (Application of nonlinear hierarchical models to the kinetic evaluation of chemical degradation data) of the German Environment Agency carried out in 2022 and 2023.</p> -<p>The mkin package is used in version 1.2.4, which is currently under +<p>The mkin package is used in version 1.2.10, which is currently under development. It contains the test data, and the functions used in 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 @@ -1847,7 +1795,9 @@ performed on the basis of the refined fit shown above.</p> <div class="sourceCode" id="cb13"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">saem_sforb_path_1_tc_reduced_multi</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/multistart.html">multistart</a></span><span class="op">(</span><span class="va">saem_sforb_path_1_tc_reduced</span>,</span> <span> n <span class="op">=</span> <span class="fl">32</span>, cores <span class="op">=</span> <span class="fl">10</span><span class="op">)</span></span></code></pre></div> -<div class="sourceCode" id="cb14"><pre class="downlit sourceCode r"> +<pre><code> + (subscript) logical subscript too long</code></pre> +<div class="sourceCode" id="cb15"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/print.html" class="external-link">print</a></span><span class="op">(</span><span class="va">saem_sforb_path_1_tc_reduced_multi</span><span class="op">)</span></span></code></pre></div> <pre><code><multistart> object with 32 fits: E OK @@ -1861,7 +1811,7 @@ the SAEM algorithm leads to parameter combinations for the degradation model that the numerical integration routine cannot cope with. Because of this variation of initial parameters, some of the model fits take up to two times more time than the original fit.</p> -<div class="sourceCode" id="cb16"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb17"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/graphics/par.html" class="external-link">par</a></span><span class="op">(</span>mar <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">12.1</span>, <span class="fl">4.1</span>, <span class="fl">2.1</span>, <span class="fl">2.1</span><span class="op">)</span><span class="op">)</span></span> <span><span class="fu"><a href="../../reference/parplot.html">parplot</a></span><span class="op">(</span><span class="va">saem_sforb_path_1_tc_reduced_multi</span>, ylim <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">0.5</span>, <span class="fl">2</span><span class="op">)</span>, las <span class="op">=</span> <span class="fl">2</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> @@ -1880,14 +1830,14 @@ ill-defined parameters.</p> </h2> <p>The SFORB pathway fits with full and reduced parameter distribution model are shown below.</p> -<div class="sourceCode" id="cb17"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb18"><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">saem_1</span><span class="op">[[</span><span class="st">"sforb_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png" alt="SFORB pathway fit with two-component error" width="700"><p class="caption"> SFORB pathway fit with two-component error </p> </div> -<div class="sourceCode" id="cb18"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb19"><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">saem_sforb_path_1_tc_reduced</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png" alt="SFORB pathway fit with two-component error, reduced parameter model" width="700"><p class="caption"> @@ -1896,7 +1846,7 @@ SFORB pathway fit with two-component error, reduced parameter model </div> <p>Plots of the remaining fits and listings for all successful fits are shown in the Appendix.</p> -<div class="sourceCode" id="cb19"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb20"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/parallel/makeCluster.html" class="external-link">stopCluster</a></span><span class="op">(</span><span class="va">cl</span><span class="op">)</span></span></code></pre></div> </div> <div class="section level2"> @@ -1936,21 +1886,21 @@ Level with Nonlinear Mixed-Effects Models.”</span> <em>Environments</em> <div class="section level3"> <h3 id="plots-of-hierarchical-fits-not-selected-for-refinement">Plots of hierarchical fits not selected for refinement<a class="anchor" aria-label="anchor" href="#plots-of-hierarchical-fits-not-selected-for-refinement"></a> </h3> -<div class="sourceCode" id="cb20"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb21"><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">saem_1</span><span class="op">[[</span><span class="st">"sfo_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png" alt="SFO pathway fit with two-component error" width="700"><p class="caption"> SFO pathway fit with two-component error </p> </div> -<div class="sourceCode" id="cb21"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb22"><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">saem_1</span><span class="op">[[</span><span class="st">"fomc_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png" alt="FOMC pathway fit with two-component error" width="700"><p class="caption"> FOMC pathway fit with two-component error </p> </div> -<div class="sourceCode" id="cb22"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb23"><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">saem_1</span><span class="op">[[</span><span class="st">"sforb_path_1"</span>, <span class="st">"tc"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <div class="figure" style="text-align: center"> <img src="2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png" alt="HS pathway fit with two-component error" width="700"><p class="caption"> @@ -1976,13 +1926,13 @@ error<a class="anchor" aria-label="anchor" href="#improved-fit-of-the-sforb-path <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.0 Patched (2023-05-18 r84448) -Platform: x86_64-pc-linux-gnu (64-bit) +<pre><code>R version 4.5.0 (2025-04-11) +Platform: x86_64-pc-linux-gnu Running under: Debian GNU/Linux 12 (bookworm) Matrix products: default -BLAS: /home/jranke/svn/R/r-patched/build/lib/libRblas.so -LAPACK: /usr/lib/x86_64-linux-gnu/openblas-serial/liblapack.so.3; LAPACK version 3.11.0 +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 LAPACK version 3.11.0 locale: [1] LC_CTYPE=de_DE.UTF-8 LC_NUMERIC=C @@ -2000,62 +1950,52 @@ attached base packages: [8] base other attached packages: -[1] saemix_3.2 npde_3.3 knitr_1.42 mkin_1.2.4 +[1] rmarkdown_2.29 nvimcom_0.9-167 saemix_3.3 npde_3.5 +[5] knitr_1.49 mkin_1.2.10 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.21-8 digest_0.6.31 magrittr_2.0.3 evaluate_0.21 - [9] grid_4.3.0 fastmap_1.1.1 rprojroot_2.0.3 jsonlite_1.8.4 -[13] processx_3.8.1 pkgbuild_1.4.0 deSolve_1.35 DBI_1.1.3 -[17] mclust_6.0.0 ps_1.7.5 gridExtra_2.3 purrr_1.0.1 -[21] fansi_1.0.4 scales_1.2.1 codetools_0.2-19 textshaping_0.3.6 -[25] jquerylib_0.1.4 cli_3.6.1 crayon_1.5.2 rlang_1.1.1 -[29] munsell_0.5.0 cachem_1.0.8 yaml_2.3.7 inline_0.3.19 -[33] tools_4.3.0 memoise_2.0.1 dplyr_1.1.2 colorspace_2.1-0 -[37] ggplot2_3.4.2 vctrs_0.6.2 R6_2.5.1 zoo_1.8-12 -[41] lifecycle_1.0.3 stringr_1.5.0 fs_1.6.2 ragg_1.2.5 -[45] callr_3.7.3 pkgconfig_2.0.3 desc_1.4.2 pkgdown_2.0.7 -[49] bslib_0.4.2 pillar_1.9.0 gtable_0.3.3 glue_1.6.2 -[53] systemfonts_1.0.4 highr_0.10 xfun_0.39 tibble_3.2.1 -[57] lmtest_0.9-40 tidyselect_1.2.0 htmltools_0.5.5 nlme_3.1-162 -[61] rmarkdown_2.21 compiler_4.3.0 prettyunits_1.1.1</code></pre> + [1] sass_0.4.9 generics_0.1.3 lattice_0.22-6 digest_0.6.37 + [5] magrittr_2.0.3 evaluate_1.0.3 grid_4.5.0 fastmap_1.2.0 + [9] jsonlite_1.9.0 processx_3.8.6 pkgbuild_1.4.6 deSolve_1.40 +[13] mclust_6.1.1 ps_1.9.0 gridExtra_2.3 scales_1.3.0 +[17] codetools_0.2-20 textshaping_1.0.0 jquerylib_0.1.4 cli_3.6.4 +[21] rlang_1.1.5 munsell_0.5.1 cachem_1.1.0 yaml_2.3.10 +[25] tools_4.5.0 inline_0.3.21 dplyr_1.1.4 colorspace_2.1-1 +[29] ggplot2_3.5.1 vctrs_0.6.5 R6_2.6.1 zoo_1.8-13 +[33] lifecycle_1.0.4 fs_1.6.5 htmlwidgets_1.6.4 MASS_7.3-65 +[37] ragg_1.3.3 pkgconfig_2.0.3 desc_1.4.3 callr_3.7.6 +[41] pkgdown_2.1.1 pillar_1.10.1 bslib_0.9.0 gtable_0.3.6 +[45] glue_1.8.0 systemfonts_1.2.1 xfun_0.51 tibble_3.2.1 +[49] lmtest_0.9-40 tidyselect_1.2.1 htmltools_0.5.8.1 nlme_3.1-168 +[53] compiler_4.5.0 </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: AMD Ryzen 9 7950X 16-Core Processor</code></pre> -<pre><code>MemTotal: 64925476 kB</code></pre> +<pre><code>MemTotal: 64927780 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> - + </main><aside class="col-md-3"><nav id="toc" aria-label="Table of contents"><h2>On this page</h2> + </nav></aside> </div> - <footer><div class="copyright"> - <p></p> -<p>Developed by Johannes Ranke.</p> + <footer><div class="pkgdown-footer-left"> + <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 class="pkgdown-footer-right"> + <p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.1.1.</p> </div> - </footer> + </footer> </div> - - + </body> </html> diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/saem-sforb-path-1-tc-reduced-convergence-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/saem-sforb-path-1-tc-reduced-convergence-1.png Binary files differindex a7187abc..0facfbe5 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/saem-sforb-path-1-tc-reduced-convergence-1.png +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/saem-sforb-path-1-tc-reduced-convergence-1.png diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-2-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-2-1.png Binary files differindex ba36beb2..4e363872 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-2-1.png +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-2-1.png diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png Binary files differindex 80376942..562d8d9d 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png Binary files differindex 15a30cfb..7bf4aa4a 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png Binary files differindex 5445ddb5..602c1138 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png Binary files differindex 067b144d..8db8e11c 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png Binary files differindex 80376942..562d8d9d 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/saem-sforb-path-1-tc-reduced-convergence-1.pdf b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/saem-sforb-path-1-tc-reduced-convergence-1.pdf Binary files differnew file mode 100644 index 00000000..ec73c393 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/saem-sforb-path-1-tc-reduced-convergence-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-2-1.pdf b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-2-1.pdf Binary files differnew file mode 100644 index 00000000..7bdec976 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-2-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-3-1.pdf b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-3-1.pdf Binary files differnew file mode 100644 index 00000000..1f8968f3 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-3-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-4-1.pdf b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-4-1.pdf Binary files differnew file mode 100644 index 00000000..1a65f828 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-4-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-5-1.pdf b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-5-1.pdf Binary files differnew file mode 100644 index 00000000..b8b103a8 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-5-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-6-1.pdf b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-6-1.pdf Binary files differnew file mode 100644 index 00000000..e4ef5ccb --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-6-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-7-1.pdf b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-7-1.pdf Binary files differnew file mode 100644 index 00000000..9b782707 --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-7-1.pdf diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-8-1.pdf b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-8-1.pdf Binary files differnew file mode 100644 index 00000000..96f2134a --- /dev/null +++ b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-latex/unnamed-chunk-8-1.pdf diff --git a/docs/articles/prebuilt/2023_mesotrione_parent.html b/docs/articles/prebuilt/2023_mesotrione_parent.html index 2b81cec4..92df6e51 100644 --- a/docs/articles/prebuilt/2023_mesotrione_parent.html +++ b/docs/articles/prebuilt/2023_mesotrione_parent.html @@ -4,145 +4,90 @@ <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"> +<meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"> <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]--> +<script src="../../deps/jquery-3.6.0/jquery-3.6.0.min.js"></script><meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"> +<link href="../../deps/bootstrap-5.3.1/bootstrap.min.css" rel="stylesheet"> +<script src="../../deps/bootstrap-5.3.1/bootstrap.bundle.min.js"></script><link href="../../deps/font-awesome-6.5.2/css/all.min.css" rel="stylesheet"> +<link href="../../deps/font-awesome-6.5.2/css/v4-shims.min.css" rel="stylesheet"> +<script src="../../deps/headroom-0.11.0/headroom.min.js"></script><script src="../../deps/headroom-0.11.0/jQuery.headroom.min.js"></script><script src="../../deps/bootstrap-toc-1.0.1/bootstrap-toc.min.js"></script><script src="../../deps/clipboard.js-2.0.11/clipboard.min.js"></script><script src="../../deps/search-1.0.0/autocomplete.jquery.min.js"></script><script src="../../deps/search-1.0.0/fuse.min.js"></script><script src="../../deps/search-1.0.0/mark.min.js"></script><!-- pkgdown --><script src="../../pkgdown.js"></script><meta property="og:title" content="Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione"> </head> -<body data-spy="scroll" data-target="#toc"> - +<body> + <a href="#main" class="visually-hidden-focusable">Skip to contents</a> - <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/2023_mesotrione_parent.html">Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione</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> + <nav class="navbar navbar-expand-lg fixed-top bg-light" data-bs-theme="default" aria-label="Site navigation"><div class="container"> + + <a class="navbar-brand me-2" href="../../index.html">mkin</a> + + <small class="nav-text text-muted me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="Released version">1.2.10</small> + + + <button class="navbar-toggler" type="button" data-bs-toggle="collapse" data-bs-target="#navbar" aria-controls="navbar" aria-expanded="false" aria-label="Toggle navigation"> + <span class="navbar-toggler-icon"></span> + </button> + + <div id="navbar" class="collapse navbar-collapse ms-3"> + <ul class="navbar-nav me-auto"> +<li class="nav-item"><a class="nav-link" href="../../reference/index.html">Reference</a></li> +<li class="active nav-item dropdown"> + <button class="nav-link dropdown-toggle" type="button" id="dropdown-articles" data-bs-toggle="dropdown" aria-expanded="false" aria-haspopup="true">Articles</button> + <ul class="dropdown-menu" aria-labelledby="dropdown-articles"> +<li><a class="dropdown-item" href="../../articles/mkin.html">Introduction to mkin</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Example evaluations with (generalised) nonlinear least squares</h6></li> + <li><a class="dropdown-item" href="../../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a></li> + <li><a class="dropdown-item" href="../../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Example evaluations with hierarchical models (nonlinear mixed-effects models)</h6></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_dmta_parent.html">Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_dmta_pathway.html">Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2023_mesotrione_parent.html">Testing covariate modelling in hierarchical parent degradation kinetics with residue data on mesotrione</a></li> + <li><a class="dropdown-item" href="../../articles/prebuilt/2022_cyan_pathway.html">Testing hierarchical pathway kinetics with residue data on cyantraniliprole</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/dimethenamid_2018.html">Comparison of saemix and nlme evaluations of dimethenamid data from 2018</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/multistart.html">Short demo of the multistart method</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Performance</h6></li> + <li><a class="dropdown-item" href="../../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/benchmarks.html">Benchmark timings for mkin</a></li> + <li><a class="dropdown-item" href="../../articles/web_only/saem_benchmarks.html">Benchmark timings for saem.mmkin</a></li> + <li><hr class="dropdown-divider"></li> + <li><h6 class="dropdown-header" data-toc-skip>Miscellaneous</h6></li> + <li><a class="dropdown-item" href="../../articles/twa.html">Calculation of time weighted average concentrations with mkin</a></li> + <li><a class="dropdown-item" 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> +<li class="nav-item"><a class="nav-link" href="../../coverage/coverage.html">Test coverage</a></li> +<li class="nav-item"><a class="nav-link" 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 class="navbar-nav"> +<li class="nav-item"><form class="form-inline" role="search"> + <input class="form-control" type="search" name="search-input" id="search-input" autocomplete="off" aria-label="Search site" placeholder="Search for" data-search-index="../../search.json"> +</form></li> +<li class="nav-item"><a class="external-link nav-link" href="https://github.com/jranke/mkin/" aria-label="GitHub"><span class="fa fab fa-github fa-lg"></span></a></li> </ul> </div> -<!--/.nav-collapse --> + + </div> -<!--/.container --> -</div> -<!--/.navbar --> +</nav><div class="container template-article"> - - </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> + + +<div class="row"> + <main id="main" class="col-md-9"><div class="page-header"> + + <h1>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 10 August 2023</h4> +last compiled on 12 Mai 2025</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 class="d-none name"><code>2023_mesotrione_parent.rmd</code></div> </div> @@ -158,7 +103,7 @@ 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 +<p>The mkin package is used in version 1.2.10, 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> @@ -1520,7 +1465,7 @@ according to the <code>illparms</code> function.</p> <tr class="odd"> <td align="left">HS</td> <td align="left"></td> -<td align="left"></td> +<td align="left">b.1</td> </tr> </tbody> </table> @@ -2480,13 +2425,13 @@ further refined to make them fully identifiable.</p> <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) +<pre><code>R version 4.5.0 (2025-04-11) +Platform: x86_64-pc-linux-gnu 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 +LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.11.0 LAPACK version 3.11.0 locale: [1] LC_CTYPE=de_DE.UTF-8 LC_NUMERIC=C @@ -2504,59 +2449,51 @@ attached base packages: [8] base other attached packages: -[1] saemix_3.2 npde_3.3 knitr_1.43 mkin_1.2.5 +[1] rmarkdown_2.29 nvimcom_0.9-167 saemix_3.3 npde_3.5 +[5] knitr_1.49 mkin_1.2.10 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 rprojroot_2.0.3 jsonlite_1.8.5 -[13] mclust_6.0.0 gridExtra_2.3 purrr_1.0.1 fansi_1.0.4 -[17] scales_1.2.1 textshaping_0.3.6 jquerylib_0.1.4 cli_3.6.1 -[21] rlang_1.1.1 munsell_0.5.0 cachem_1.0.8 yaml_2.3.7 -[25] tools_4.3.1 memoise_2.0.1 dplyr_1.1.2 colorspace_2.1-0 -[29] ggplot2_3.4.2 vctrs_0.6.2 R6_2.5.1 zoo_1.8-12 -[33] lifecycle_1.0.3 stringr_1.5.0 fs_1.6.2 ragg_1.2.5 -[37] pkgconfig_2.0.3 desc_1.4.2 pkgdown_2.0.7 bslib_0.4.2 -[41] pillar_1.9.0 gtable_0.3.3 glue_1.6.2 systemfonts_1.0.4 -[45] xfun_0.39 tibble_3.2.1 lmtest_0.9-40 tidyselect_1.2.0 -[49] htmltools_0.5.5 nlme_3.1-162 rmarkdown_2.22 compiler_4.3.1 </code></pre> + [1] gtable_0.3.6 jsonlite_1.9.0 dplyr_1.1.4 compiler_4.5.0 + [5] tidyselect_1.2.1 gridExtra_2.3 jquerylib_0.1.4 systemfonts_1.2.1 + [9] scales_1.3.0 textshaping_1.0.0 readxl_1.4.4 yaml_2.3.10 +[13] fastmap_1.2.0 lattice_0.22-6 ggplot2_3.5.1 R6_2.6.1 +[17] generics_0.1.3 lmtest_0.9-40 MASS_7.3-65 htmlwidgets_1.6.4 +[21] tibble_3.2.1 desc_1.4.3 munsell_0.5.1 bslib_0.9.0 +[25] pillar_1.10.1 rlang_1.1.5 cachem_1.1.0 xfun_0.51 +[29] fs_1.6.5 sass_0.4.9 cli_3.6.4 pkgdown_2.1.1 +[33] magrittr_2.0.3 digest_0.6.37 grid_4.5.0 mclust_6.1.1 +[37] lifecycle_1.0.4 nlme_3.1-168 vctrs_0.6.5 evaluate_1.0.3 +[41] glue_1.8.0 cellranger_1.1.0 codetools_0.2-20 ragg_1.3.3 +[45] zoo_1.8-13 colorspace_2.1-1 tools_4.5.0 pkgconfig_2.0.3 +[49] htmltools_0.5.8.1</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> +<pre><code>CPU model: AMD Ryzen 9 7950X 16-Core Processor</code></pre> +<pre><code>MemTotal: 64927780 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> - + </main><aside class="col-md-3"><nav id="toc" aria-label="Table of contents"><h2>On this page</h2> + </nav></aside> </div> - <footer><div class="copyright"> - <p></p> -<p>Developed by Johannes Ranke.</p> + <footer><div class="pkgdown-footer-left"> + <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 class="pkgdown-footer-right"> + <p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.1.1.</p> </div> - </footer> + </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 Binary files differindex 863a48bd..719bf0e9 100644 --- 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 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 Binary files differindex 256b2b68..6afcac48 100644 --- 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 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 Binary files differindex 59011020..50253fb1 100644 --- 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 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 Binary files differindex f427bc39..c8b77724 100644 --- 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 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 Binary files differindex 7c3b460b..8fa204a3 100644 --- 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 |
