diff options
Diffstat (limited to 'docs/dev/articles/web_only/NAFTA_examples.html')
| -rw-r--r-- | docs/dev/articles/web_only/NAFTA_examples.html | 375 |
1 files changed, 153 insertions, 222 deletions
diff --git a/docs/dev/articles/web_only/NAFTA_examples.html b/docs/dev/articles/web_only/NAFTA_examples.html index e4e457fc..54c45b3e 100644 --- a/docs/dev/articles/web_only/NAFTA_examples.html +++ b/docs/dev/articles/web_only/NAFTA_examples.html @@ -4,143 +4,91 @@ <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>Evaluation of example datasets from Attachment 1 to the US EPA SOP for the NAFTA guidance • 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="Evaluation of example datasets from Attachment 1 to the US EPA SOP for the NAFTA guidance"> -<meta property="og:description" content="mkin"> +<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="Evaluation of example datasets from Attachment 1 to the US EPA SOP for the NAFTA guidance"> <meta name="robots" content="noindex"> -<!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]> -<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script> -<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script> -<![endif]--> </head> -<body data-spy="scroll" data-target="#toc"> - +<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-info" data-toggle="tooltip" data-placement="bottom" title="In-development version">1.2.5</span> - </span> - </div> - <div id="navbar" class="navbar-collapse collapse"> - <ul class="nav navbar-nav"> -<li> - <a href="../../reference/index.html">Reference</a> -</li> -<li class="dropdown"> - <a href="#" class="dropdown-toggle" data-toggle="dropdown" role="button" data-bs-toggle="dropdown" aria-expanded="false"> - Articles - - <span class="caret"></span> - </a> - <ul class="dropdown-menu" role="menu"> -<li> - <a href="../../articles/mkin.html">Introduction to mkin</a> - </li> - <li class="divider"> - </li> -<li class="dropdown-header">Example evaluations with (generalised) nonlinear least squares</li> - <li> - <a href="../../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a> - </li> - <li> - <a href="../../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a> - </li> - <li> - <a href="../../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a> - </li> - <li class="divider"> - </li> -<li class="dropdown-header">Example evaluations with hierarchical models (nonlinear mixed-effects models)</li> - <li> - <a href="../../articles/prebuilt/2022_dmta_parent.html">Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P</a> - </li> - <li> - <a href="../../articles/prebuilt/2022_dmta_pathway.html">Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</a> - </li> - <li> - <a href="../../articles/prebuilt/2022_cyan_pathway.html">Testing hierarchical pathway kinetics with residue data on cyantraniliprole</a> - </li> - <li> - <a href="../../articles/web_only/dimethenamid_2018.html">Comparison of saemix and nlme evaluations of dimethenamid data from 2018</a> - </li> - <li> - <a href="../../articles/web_only/multistart.html">Short demo of the multistart method</a> - </li> - <li class="divider"> - </li> -<li class="dropdown-header">Performance</li> - <li> - <a href="../../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a> - </li> - <li> - <a href="../../articles/web_only/benchmarks.html">Benchmark timings for mkin</a> - </li> - <li> - <a href="../../articles/web_only/saem_benchmarks.html">Benchmark timings for saem.mmkin</a> - </li> - <li class="divider"> - </li> -<li class="dropdown-header">Miscellaneous</li> - <li> - <a href="../../articles/twa.html">Calculation of time weighted average concentrations with mkin</a> - </li> - <li> - <a href="../../articles/web_only/NAFTA_examples.html">Example evaluation of NAFTA SOP Attachment examples</a> - </li> + <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-info me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="In-development 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>Evaluation of example datasets from Attachment 1 -to the US EPA SOP for the NAFTA guidance</h1> + + +<div class="row"> + <main id="main" class="col-md-9"><div class="page-header"> + + <h1>Evaluation of example datasets from Attachment 1 to the US EPA SOP for the NAFTA guidance</h1> <h4 data-toc-skip class="author">Johannes Ranke</h4> <h4 data-toc-skip class="date">26 February 2019 (rebuilt -2023-05-19)</h4> +2025-02-16)</h4> <small class="dont-index">Source: <a href="https://github.com/jranke/mkin/blob/HEAD/vignettes/web_only/NAFTA_examples.rmd" class="external-link"><code>vignettes/web_only/NAFTA_examples.rmd</code></a></small> - <div class="hidden name"><code>NAFTA_examples.rmd</code></div> - + <div class="d-none name"><code>NAFTA_examples.rmd</code></div> </div> @@ -201,7 +149,7 @@ same.</p> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> <span><span class="co">## parent_0 9.99e+01 1.41e-26 98.8116 101.0810</span></span> <span><span class="co">## k1 2.67e-02 5.05e-06 0.0243 0.0295</span></span> -<span><span class="co">## k2 2.26e-12 5.00e-01 0.0000 Inf</span></span> +<span><span class="co">## k2 3.41e-12 5.00e-01 0.0000 Inf</span></span> <span><span class="co">## g 6.47e-01 3.67e-06 0.6248 0.6677</span></span> <span><span class="co">## sigma 1.27e+00 8.91e-06 0.8395 1.6929</span></span> <span><span class="co">## </span></span> @@ -210,7 +158,7 @@ same.</p> <span><span class="co">## DT50 DT90 DT50_rep</span></span> <span><span class="co">## SFO 67.7 2.25e+02 6.77e+01</span></span> <span><span class="co">## IORE 58.2 1.07e+03 3.22e+02</span></span> -<span><span class="co">## DFOP 55.5 5.59e+11 3.07e+11</span></span> +<span><span class="co">## DFOP 55.5 3.70e+11 2.03e+11</span></span> <span><span class="co">## </span></span> <span><span class="co">## Representative half-life:</span></span> <span><span class="co">## [1] 321.51</span></span></code></pre> @@ -252,7 +200,7 @@ same.</p> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> <span><span class="co">## parent_0 9.84e+01 1.24e-27 97.8078 98.9187</span></span> <span><span class="co">## k1 1.55e-02 4.10e-04 0.0143 0.0167</span></span> -<span><span class="co">## k2 8.63e-12 5.00e-01 0.0000 Inf</span></span> +<span><span class="co">## k2 9.07e-12 5.00e-01 0.0000 Inf</span></span> <span><span class="co">## g 6.89e-01 2.92e-03 0.6626 0.7142</span></span> <span><span class="co">## sigma 6.48e-01 2.38e-05 0.4147 0.8813</span></span> <span><span class="co">## </span></span> @@ -261,7 +209,7 @@ same.</p> <span><span class="co">## DT50 DT90 DT50_rep</span></span> <span><span class="co">## SFO 86.6 2.88e+02 8.66e+01</span></span> <span><span class="co">## IORE 85.5 7.17e+02 2.16e+02</span></span> -<span><span class="co">## DFOP 83.6 1.32e+11 8.04e+10</span></span> +<span><span class="co">## DFOP 83.6 1.25e+11 7.64e+10</span></span> <span><span class="co">## </span></span> <span><span class="co">## Representative half-life:</span></span> <span><span class="co">## [1] 215.87</span></span></code></pre> @@ -303,7 +251,7 @@ same.</p> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> <span><span class="co">## parent_0 9.66e+01 1.57e-25 95.3476 97.8979</span></span> <span><span class="co">## k1 2.55e-02 7.33e-06 0.0233 0.0278</span></span> -<span><span class="co">## k2 3.22e-11 5.00e-01 0.0000 Inf</span></span> +<span><span class="co">## k2 3.84e-11 5.00e-01 0.0000 Inf</span></span> <span><span class="co">## g 8.61e-01 7.55e-06 0.8314 0.8867</span></span> <span><span class="co">## sigma 1.46e+00 6.93e-06 0.9661 1.9483</span></span> <span><span class="co">## </span></span> @@ -312,7 +260,7 @@ same.</p> <span><span class="co">## DT50 DT90 DT50_rep</span></span> <span><span class="co">## SFO 38.6 1.28e+02 3.86e+01</span></span> <span><span class="co">## IORE 34.0 1.77e+02 5.32e+01</span></span> -<span><span class="co">## DFOP 34.1 1.01e+10 2.15e+10</span></span> +<span><span class="co">## DFOP 34.1 8.50e+09 1.80e+10</span></span> <span><span class="co">## </span></span> <span><span class="co">## Representative half-life:</span></span> <span><span class="co">## [1] 53.17</span></span></code></pre> @@ -354,7 +302,7 @@ same.</p> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> <span><span class="co">## parent_0 9.89e+01 9.44e-49 95.4640 102.2573</span></span> <span><span class="co">## k1 1.81e-02 1.75e-01 0.0116 0.0281</span></span> -<span><span class="co">## k2 3.63e-10 5.00e-01 0.0000 Inf</span></span> +<span><span class="co">## k2 3.62e-10 5.00e-01 0.0000 Inf</span></span> <span><span class="co">## g 6.06e-01 2.19e-01 0.4826 0.7178</span></span> <span><span class="co">## sigma 7.40e+00 2.97e-15 6.0201 8.7754</span></span> <span><span class="co">## </span></span> @@ -363,7 +311,7 @@ same.</p> <span><span class="co">## DT50 DT90 DT50_rep</span></span> <span><span class="co">## SFO 94.3 3.13e+02 9.43e+01</span></span> <span><span class="co">## IORE 96.7 1.51e+03 4.55e+02</span></span> -<span><span class="co">## DFOP 96.4 3.77e+09 1.91e+09</span></span> +<span><span class="co">## DFOP 96.4 3.79e+09 1.92e+09</span></span> <span><span class="co">## </span></span> <span><span class="co">## Representative half-life:</span></span> <span><span class="co">## [1] 454.55</span></span></code></pre> @@ -477,7 +425,7 @@ lower value for the rate constant is used here.</p> <span><span class="co">## DT50 DT90 DT50_rep</span></span> <span><span class="co">## SFO 16.9 5.63e+01 1.69e+01</span></span> <span><span class="co">## IORE 11.6 3.37e+02 1.01e+02</span></span> -<span><span class="co">## DFOP 10.5 1.38e+12 7.69e+11</span></span> +<span><span class="co">## DFOP 10.5 1.38e+12 7.68e+11</span></span> <span><span class="co">## </span></span> <span><span class="co">## Representative half-life:</span></span> <span><span class="co">## [1] 101.43</span></span></code></pre> @@ -490,17 +438,12 @@ suggest a simple exponential decline.</p> </h3> <div class="sourceCode" id="cb37"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p9b</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p9b"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> -<pre><code><span><span class="co">## Warning in sqrt(diag(covar)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in sqrt(diag(covar_notrans)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in sqrt(1/diag(V)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(.) had 0 or NA entries; non-finite result</span></span> -<span><span class="co">## is doubtful</span></span></code></pre> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb44"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb40"><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">p9b</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p9b-1.png" width="700"></p> -<div class="sourceCode" id="cb45"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb41"><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">p9b</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -527,8 +470,8 @@ suggest a simple exponential decline.</p> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> <span><span class="co">## parent_0 94.7123 1.61e-16 93.1355 96.2891</span></span> <span><span class="co">## k1 0.0389 1.08e-04 0.0266 0.0569</span></span> -<span><span class="co">## k2 0.0389 2.23e-04 0.0255 0.0592</span></span> -<span><span class="co">## g 0.5256 NaN NA NA</span></span> +<span><span class="co">## k2 0.0389 2.24e-04 0.0255 0.0592</span></span> +<span><span class="co">## g 0.5256 5.00e-01 0.0000 1.0000</span></span> <span><span class="co">## sigma 1.5957 2.50e-04 0.9135 2.2779</span></span> <span><span class="co">## </span></span> <span><span class="co">## </span></span> @@ -550,18 +493,15 @@ in PestDF and g in mkin. In mkin, it is restricted to the interval from <div class="section level3"> <h3 id="example-on-page-10">Example on page 10<a class="anchor" aria-label="anchor" href="#example-on-page-10"></a> </h3> -<div class="sourceCode" id="cb47"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb43"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p10</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p10"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> -<pre><code><span><span class="co">## Warning in sqrt(diag(covar)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in sqrt(1/diag(V)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(.) had 0 or NA entries; non-finite result</span></span> -<span><span class="co">## is doubtful</span></span></code></pre> +<pre><code><span><span class="co">## Warning in sqrt(diag(covar_notrans)): NaNs produced</span></span></code></pre> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb53"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb47"><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">p10</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p10-1.png" width="700"></p> -<div class="sourceCode" id="cb54"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb48"><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">p10</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -585,12 +525,12 @@ in PestDF and g in mkin. In mkin, it is restricted to the interval from <span><span class="co">## sigma 4.90 1.77e-04 2.837 6.968</span></span> <span><span class="co">## </span></span> <span><span class="co">## $DFOP</span></span> -<span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> -<span><span class="co">## parent_0 101.7315 1.41e-09 91.6534 111.8097</span></span> -<span><span class="co">## k1 0.0495 6.58e-03 0.0303 0.0809</span></span> -<span><span class="co">## k2 0.0495 2.60e-03 0.0410 0.0598</span></span> -<span><span class="co">## g 0.4487 5.00e-01 NA NA</span></span> -<span><span class="co">## sigma 8.0152 2.50e-04 4.5886 11.4418</span></span> +<span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> +<span><span class="co">## parent_0 101.7315 1.41e-09 91.6534 111.810</span></span> +<span><span class="co">## k1 0.0495 3.04e-03 0.0188 0.131</span></span> +<span><span class="co">## k2 0.0495 4.92e-04 0.0197 0.124</span></span> +<span><span class="co">## g 0.4487 NaN 0.0000 1.000</span></span> +<span><span class="co">## sigma 8.0152 2.50e-04 4.5886 11.442</span></span> <span><span class="co">## </span></span> <span><span class="co">## </span></span> <span><span class="co">## DTx values:</span></span> @@ -614,14 +554,14 @@ difference in IORE model parameters between PestDF and mkin.</p> <div class="section level3"> <h3 id="example-on-page-11">Example on page 11<a class="anchor" aria-label="anchor" href="#example-on-page-11"></a> </h3> -<div class="sourceCode" id="cb56"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb50"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p11</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p11"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb59"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb53"><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">p11</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p11-1.png" width="700"></p> -<div class="sourceCode" id="cb60"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb54"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="fu"><a href="https://rdrr.io/r/base/print.html" class="external-link">print</a></span><span class="op">(</span><span class="va">p11</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -660,7 +600,7 @@ difference in IORE model parameters between PestDF and mkin.</p> <span><span class="co">## DFOP 3.07e+11 1.93e+12 6.98e+11</span></span> <span><span class="co">## </span></span> <span><span class="co">## Representative half-life:</span></span> -<span><span class="co">## [1] 41148170</span></span></code></pre> +<span><span class="co">## [1] 41148169</span></span></code></pre> <p>In this case, the DFOP fit reported for PestDF resulted in a negative value for the slower rate constant, which is not possible in mkin. The other results are in agreement.</p> @@ -677,21 +617,17 @@ overparameterisation.</p> <div class="section level3"> <h3 id="example-on-page-12-upper-panel">Example on page 12, upper panel<a class="anchor" aria-label="anchor" href="#example-on-page-12-upper-panel"></a> </h3> -<div class="sourceCode" id="cb62"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb56"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p12a</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p12a"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Warning in summary.mkinfit(x): Could not calculate correlation; no covariance</span></span> <span><span class="co">## matrix</span></span></code></pre> -<pre><code><span><span class="co">## Warning in sqrt(diag(covar)): NaNs produced</span></span></code></pre> <pre><code><span><span class="co">## Warning in sqrt(diag(covar_notrans)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in sqrt(1/diag(V)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(.) had 0 or NA entries; non-finite result</span></span> -<span><span class="co">## is doubtful</span></span></code></pre> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb70"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb61"><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">p12a</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p12a-1.png" width="700"></p> -<div class="sourceCode" id="cb71"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb62"><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">p12a</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -719,7 +655,7 @@ overparameterisation.</p> <span><span class="co">## parent_0 100.521 2.74e-10 92.2366 108.805</span></span> <span><span class="co">## k1 0.124 2.53e-05 0.0908 0.170</span></span> <span><span class="co">## k2 0.124 2.52e-02 0.0456 0.339</span></span> -<span><span class="co">## g 0.793 NaN NA NA</span></span> +<span><span class="co">## g 0.793 NaN 0.0000 1.000</span></span> <span><span class="co">## sigma 7.048 2.50e-04 4.0349 10.061</span></span> <span><span class="co">## </span></span> <span><span class="co">## </span></span> @@ -735,18 +671,20 @@ overparameterisation.</p> <div class="section level3"> <h3 id="example-on-page-12-lower-panel">Example on page 12, lower panel<a class="anchor" aria-label="anchor" href="#example-on-page-12-lower-panel"></a> </h3> -<div class="sourceCode" id="cb73"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb64"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p12b</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p12b"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> +<pre><code><span><span class="co">## Warning in sqrt(diag(covar)): NaNs produced</span></span></code></pre> <pre><code><span><span class="co">## Warning in qt(alpha/2, rdf): NaNs produced</span></span></code></pre> <pre><code><span><span class="co">## Warning in qt(1 - alpha/2, rdf): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in sqrt(diag(covar_notrans)): NaNs produced</span></span></code></pre> <pre><code><span><span class="co">## Warning in pt(abs(tval), rdf, lower.tail = FALSE): NaNs produced</span></span></code></pre> +<pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(V) had non-positive or NA entries; the</span></span> +<span><span class="co">## non-finite result may be dubious</span></span></code></pre> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb80"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb72"><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">p12b</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p12b-1.png" width="700"></p> -<div class="sourceCode" id="cb81"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb73"><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">p12b</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -790,14 +728,14 @@ overparameterisation.</p> <div class="section level3"> <h3 id="example-on-page-13">Example on page 13<a class="anchor" aria-label="anchor" href="#example-on-page-13"></a> </h3> -<div class="sourceCode" id="cb83"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb75"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p13</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p13"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb86"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb78"><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">p13</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p13-1.png" width="700"></p> -<div class="sourceCode" id="cb87"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb79"><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">p13</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -823,8 +761,8 @@ overparameterisation.</p> <span><span class="co">## $DFOP</span></span> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> <span><span class="co">## parent_0 92.73500 NA 8.95e+01 95.92118</span></span> -<span><span class="co">## k1 0.00258 NA 4.14e-04 0.01611</span></span> -<span><span class="co">## k2 0.00258 NA 1.74e-03 0.00383</span></span> +<span><span class="co">## k1 0.00258 NA 4.18e-04 0.01592</span></span> +<span><span class="co">## k2 0.00258 NA 1.75e-03 0.00381</span></span> <span><span class="co">## g 0.16452 NA 0.00e+00 1.00000</span></span> <span><span class="co">## sigma 3.41172 NA 2.02e+00 4.79960</span></span> <span><span class="co">## </span></span> @@ -842,18 +780,17 @@ overparameterisation.</p> <div class="section level2"> <h2 id="dt50-not-observed-in-the-study-and-dfop-problems-in-pestdf">DT50 not observed in the study and DFOP problems in PestDF<a class="anchor" aria-label="anchor" href="#dt50-not-observed-in-the-study-and-dfop-problems-in-pestdf"></a> </h2> -<div class="sourceCode" id="cb89"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb81"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p14</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p14"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Warning in sqrt(diag(covar)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in sqrt(1/diag(V)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(.) had 0 or NA entries; non-finite result</span></span> -<span><span class="co">## is doubtful</span></span></code></pre> +<pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(V) had non-positive or NA entries; the</span></span> +<span><span class="co">## non-finite result may be dubious</span></span></code></pre> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb95"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb86"><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">p14</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p14-1.png" width="700"></p> -<div class="sourceCode" id="cb96"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb87"><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">p14</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -880,7 +817,7 @@ overparameterisation.</p> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> <span><span class="co">## parent_0 1.00e+02 2.96e-28 99.40280 101.2768</span></span> <span><span class="co">## k1 9.53e-03 1.20e-01 0.00638 0.0143</span></span> -<span><span class="co">## k2 6.08e-12 5.00e-01 0.00000 Inf</span></span> +<span><span class="co">## k2 5.21e-12 5.00e-01 0.00000 Inf</span></span> <span><span class="co">## g 3.98e-01 2.19e-01 0.30481 0.4998</span></span> <span><span class="co">## sigma 1.17e+00 7.68e-06 0.77406 1.5610</span></span> <span><span class="co">## </span></span> @@ -889,7 +826,7 @@ overparameterisation.</p> <span><span class="co">## DT50 DT90 DT50_rep</span></span> <span><span class="co">## SFO 2.48e+02 8.25e+02 2.48e+02</span></span> <span><span class="co">## IORE 4.34e+02 2.22e+04 6.70e+03</span></span> -<span><span class="co">## DFOP 3.05e+10 2.95e+11 1.14e+11</span></span> +<span><span class="co">## DFOP 3.55e+10 3.44e+11 1.33e+11</span></span> <span><span class="co">## </span></span> <span><span class="co">## Representative half-life:</span></span> <span><span class="co">## [1] 6697.44</span></span></code></pre> @@ -900,14 +837,17 @@ same results in mkin and PestDF.</p> <div class="section level2"> <h2 id="n-is-less-than-1-and-dfop-fraction-parameter-is-below-zero">N is less than 1 and DFOP fraction parameter is below zero<a class="anchor" aria-label="anchor" href="#n-is-less-than-1-and-dfop-fraction-parameter-is-below-zero"></a> </h2> -<div class="sourceCode" id="cb98"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb89"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p15a</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p15a"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> +<pre><code><span><span class="co">## Warning in sqrt(diag(covar)): NaNs produced</span></span></code></pre> +<pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(V) had non-positive or NA entries; the</span></span> +<span><span class="co">## non-finite result may be dubious</span></span></code></pre> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb101"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb94"><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">p15a</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p15a-1.png" width="700"></p> -<div class="sourceCode" id="cb102"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb95"><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">p15a</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -933,9 +873,9 @@ same results in mkin and PestDF.</p> <span><span class="co">## $DFOP</span></span> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> <span><span class="co">## parent_0 97.96751 2.85e-13 94.21913 101.7159</span></span> -<span><span class="co">## k1 0.00952 6.28e-02 0.00250 0.0363</span></span> -<span><span class="co">## k2 0.00952 1.27e-04 0.00646 0.0140</span></span> -<span><span class="co">## g 0.21241 5.00e-01 0.00000 1.0000</span></span> +<span><span class="co">## k1 0.00952 6.28e-02 0.00260 0.0349</span></span> +<span><span class="co">## k2 0.00952 1.27e-04 0.00652 0.0139</span></span> +<span><span class="co">## g 0.21241 5.00e-01 NA NA</span></span> <span><span class="co">## sigma 4.18778 2.50e-04 2.39747 5.9781</span></span> <span><span class="co">## </span></span> <span><span class="co">## </span></span> @@ -947,18 +887,16 @@ same results in mkin and PestDF.</p> <span><span class="co">## </span></span> <span><span class="co">## Representative half-life:</span></span> <span><span class="co">## [1] 41.33</span></span></code></pre> -<div class="sourceCode" id="cb104"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb97"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p15b</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p15b"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> -<pre><code><span><span class="co">## Warning in sqrt(diag(covar)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in sqrt(1/diag(V)): NaNs produced</span></span></code></pre> -<pre><code><span><span class="co">## Warning in cov2cor(ans$covar): diag(.) had 0 or NA entries; non-finite result</span></span> -<span><span class="co">## is doubtful</span></span></code></pre> +<pre><code><span><span class="co">## Warning in summary.mkinfit(x): Could not calculate correlation; no covariance</span></span> +<span><span class="co">## matrix</span></span></code></pre> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The half-life obtained from the IORE model may be used</span></span></code></pre> -<div class="sourceCode" id="cb110"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb101"><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">p15b</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p15b-1.png" width="700"></p> -<div class="sourceCode" id="cb111"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb102"><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">p15b</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -976,18 +914,18 @@ same results in mkin and PestDF.</p> <span><span class="co">## </span></span> <span><span class="co">## $IORE</span></span> <span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> -<span><span class="co">## parent_0 99.83 1.81e-16 97.51349 102.14</span></span> +<span><span class="co">## parent_0 99.83 1.81e-16 97.51348 102.14</span></span> <span><span class="co">## k__iore_parent 0.38 3.22e-01 0.00352 41.05</span></span> <span><span class="co">## N_parent 0.00 5.00e-01 -1.07696 1.08</span></span> <span><span class="co">## sigma 2.21 2.57e-04 1.23245 3.19</span></span> <span><span class="co">## </span></span> <span><span class="co">## $DFOP</span></span> -<span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> -<span><span class="co">## parent_0 1.01e+02 NA 9.82e+01 1.04e+02</span></span> -<span><span class="co">## k1 4.86e-03 NA 8.63e-04 2.73e-02</span></span> -<span><span class="co">## k2 4.86e-03 NA 3.21e-03 7.35e-03</span></span> -<span><span class="co">## g 1.88e-01 NA NA NA</span></span> -<span><span class="co">## sigma 2.76e+00 NA 1.58e+00 3.94e+00</span></span> +<span><span class="co">## Estimate Pr(>t) Lower Upper</span></span> +<span><span class="co">## parent_0 1.01e+02 NA NA NA</span></span> +<span><span class="co">## k1 4.86e-03 NA NA NA</span></span> +<span><span class="co">## k2 4.86e-03 NA NA NA</span></span> +<span><span class="co">## g 1.88e-01 NA NA NA</span></span> +<span><span class="co">## sigma 2.76e+00 NA NA NA</span></span> <span><span class="co">## </span></span> <span><span class="co">## </span></span> <span><span class="co">## DTx values:</span></span> @@ -1006,16 +944,16 @@ mkin and PestDF.</p> <div class="section level2"> <h2 id="the-dfop-fraction-parameter-is-greater-than-1">The DFOP fraction parameter is greater than 1<a class="anchor" aria-label="anchor" href="#the-dfop-fraction-parameter-is-greater-than-1"></a> </h2> -<div class="sourceCode" id="cb113"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb104"><pre class="downlit sourceCode r"> <code class="sourceCode R"><span><span class="va">p16</span> <span class="op"><-</span> <span class="fu"><a href="../../reference/nafta.html">nafta</a></span><span class="op">(</span><span class="va">NAFTA_SOP_Attachment</span><span class="op">[[</span><span class="st">"p16"</span><span class="op">]</span><span class="op">]</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## The SFO model is rejected as S_SFO is equal or higher than the critical value S_c</span></span></code></pre> <pre><code><span><span class="co">## The representative half-life of the IORE model is longer than the one corresponding</span></span></code></pre> <pre><code><span><span class="co">## to the terminal degradation rate found with the DFOP model.</span></span></code></pre> <pre><code><span><span class="co">## The representative half-life obtained from the DFOP model may be used</span></span></code></pre> -<div class="sourceCode" id="cb118"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb109"><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">p16</span><span class="op">)</span></span></code></pre></div> <p><img src="NAFTA_examples_files/figure-html/p16-1.png" width="700"></p> -<div class="sourceCode" id="cb119"><pre class="downlit sourceCode r"> +<div class="sourceCode" id="cb110"><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">p16</span><span class="op">)</span></span></code></pre></div> <pre><code><span><span class="co">## Sums of squares:</span></span> <span><span class="co">## SFO IORE DFOP </span></span> @@ -1056,10 +994,12 @@ mkin and PestDF.</p> <span><span class="co">## Representative half-life:</span></span> <span><span class="co">## [1] 8.93</span></span></code></pre> <p>In PestDF, the DFOP fit seems to have stuck in a local minimum, as -mkin finds a solution with a much lower <span class="math inline">\(\chi^2\)</span> error level. As the half-life from -the slower rate constant of the DFOP model is larger than the IORE -derived half-life, the NAFTA recommendation obtained with mkin is to use -the DFOP representative half-life of 8.9 days.</p> +mkin finds a solution with a much lower +<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><msup><mi>χ</mi><mn>2</mn></msup><annotation encoding="application/x-tex">\chi^2</annotation></semantics></math> +error level. As the half-life from the slower rate constant of the DFOP +model is larger than the IORE derived half-life, the NAFTA +recommendation obtained with mkin is to use the DFOP representative +half-life of 8.9 days.</p> </div> <div class="section level2"> <h2 id="conclusions">Conclusions<a class="anchor" aria-label="anchor" href="#conclusions"></a> @@ -1084,35 +1024,26 @@ Pesticide Degradation.”</span> <a href="https://www.epa.gov/pesticide-science- </div> </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> |