diff options
Diffstat (limited to 'docs/articles/FOCUS_L.html')
| -rw-r--r-- | docs/articles/FOCUS_L.html | 59 |
1 files changed, 31 insertions, 28 deletions
diff --git a/docs/articles/FOCUS_L.html b/docs/articles/FOCUS_L.html index feceab6f..a6ec5532 100644 --- a/docs/articles/FOCUS_L.html +++ b/docs/articles/FOCUS_L.html @@ -37,6 +37,9 @@ <a href="../articles/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a> </li> </ul></li> +<li> + <a href="../news/index.html">News</a> +</li> </ul><ul class="nav navbar-nav navbar-right"><li> <a href="http://github.com/jranke/mkin"> <span class="fa fa-github fa-lg"></span> @@ -74,17 +77,17 @@ FOCUS_2006_L1_mkin <-<span class="st"> </span><span class="kw"><a href="../re <p>Since mkin version 0.9-32 (July 2014), we can use shorthand notation like <code>"SFO"</code> for parent only degradation models. The following two lines fit the model and produce the summary report of the model fit. This covers the numerical analysis given in the FOCUS report.</p> <div class="sourceCode"><pre class="sourceCode r"><code class="sourceCode r">m.L1.SFO <-<span class="st"> </span><span class="kw"><a href="../reference/mkinfit.html">mkinfit</a></span>(<span class="st">"SFO"</span>, FOCUS_2006_L1_mkin, <span class="dt">quiet =</span> <span class="ot">TRUE</span>) <span class="kw">summary</span>(m.L1.SFO)</code></pre></div> -<pre><code>## mkin version: 0.9.44.9000 +<pre><code>## mkin version: 0.9.45 ## R version: 3.3.2 -## Date of fit: Thu Dec 8 08:18:47 2016 -## Date of summary: Thu Dec 8 08:18:47 2016 +## Date of fit: Thu Dec 8 09:23:13 2016 +## Date of summary: Thu Dec 8 09:23:13 2016 ## ## Equations: ## d_parent/dt = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 37 model solutions performed in 0.084 s +## Fitted with method Port using 37 model solutions performed in 0.088 s ## ## Weighting: none ## @@ -167,10 +170,10 @@ FOCUS_2006_L1_mkin <-<span class="st"> </span><span class="kw"><a href="../re <div class="sourceCode"><pre class="sourceCode r"><code class="sourceCode r"><span class="kw">plot</span>(m.L1.FOMC, <span class="dt">show_errmin =</span> <span class="ot">TRUE</span>, <span class="dt">main =</span> <span class="st">"FOCUS L1 - FOMC"</span>)</code></pre></div> <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-6-1.png" width="576"></p> <div class="sourceCode"><pre class="sourceCode r"><code class="sourceCode r"><span class="kw">summary</span>(m.L1.FOMC, <span class="dt">data =</span> <span class="ot">FALSE</span>)</code></pre></div> -<pre><code>## mkin version: 0.9.44.9000 +<pre><code>## mkin version: 0.9.45 ## R version: 3.3.2 -## Date of fit: Thu Dec 8 08:18:47 2016 -## Date of summary: Thu Dec 8 08:18:47 2016 +## Date of fit: Thu Dec 8 09:23:14 2016 +## Date of summary: Thu Dec 8 09:23:14 2016 ## ## ## Warning: Optimisation by method Port did not converge. @@ -182,7 +185,7 @@ FOCUS_2006_L1_mkin <-<span class="st"> </span><span class="kw"><a href="../re ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 188 model solutions performed in 0.462 s +## Fitted with method Port using 188 model solutions performed in 0.427 s ## ## Weighting: none ## @@ -264,17 +267,17 @@ FOCUS_2006_L2_mkin <-<span class="st"> </span><span class="kw"><a href="../re <span class="dt">main =</span> <span class="st">"FOCUS L2 - FOMC"</span>)</code></pre></div> <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-9-1.png" width="672"></p> <div class="sourceCode"><pre class="sourceCode r"><code class="sourceCode r"><span class="kw">summary</span>(m.L2.FOMC, <span class="dt">data =</span> <span class="ot">FALSE</span>)</code></pre></div> -<pre><code>## mkin version: 0.9.44.9000 +<pre><code>## mkin version: 0.9.45 ## R version: 3.3.2 -## Date of fit: Thu Dec 8 08:18:48 2016 -## Date of summary: Thu Dec 8 08:18:48 2016 +## Date of fit: Thu Dec 8 09:23:14 2016 +## Date of summary: Thu Dec 8 09:23:14 2016 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 81 model solutions performed in 0.201 s +## Fitted with method Port using 81 model solutions performed in 0.185 s ## ## Weighting: none ## @@ -334,10 +337,10 @@ FOCUS_2006_L2_mkin <-<span class="st"> </span><span class="kw"><a href="../re <span class="dt">main =</span> <span class="st">"FOCUS L2 - DFOP"</span>)</code></pre></div> <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-10-1.png" width="672"></p> <div class="sourceCode"><pre class="sourceCode r"><code class="sourceCode r"><span class="kw">summary</span>(m.L2.DFOP, <span class="dt">data =</span> <span class="ot">FALSE</span>)</code></pre></div> -<pre><code>## mkin version: 0.9.44.9000 +<pre><code>## mkin version: 0.9.45 ## R version: 3.3.2 -## Date of fit: Thu Dec 8 08:18:49 2016 -## Date of summary: Thu Dec 8 08:18:49 2016 +## Date of fit: Thu Dec 8 09:23:15 2016 +## Date of summary: Thu Dec 8 09:23:15 2016 ## ## Equations: ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * @@ -346,7 +349,7 @@ FOCUS_2006_L2_mkin <-<span class="st"> </span><span class="kw"><a href="../re ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 336 model solutions performed in 0.767 s +## Fitted with method Port using 336 model solutions performed in 0.765 s ## ## Weighting: none ## @@ -423,10 +426,10 @@ mm.L3 <-<span class="st"> </span><span class="kw"><a href="../reference/mmkin <p>The objects returned by mmkin are arranged like a matrix, with models as a row index and datasets as a column index.</p> <p>We can extract the summary and plot for <em>e.g.</em> the DFOP fit, using square brackets for indexing which will result in the use of the summary and plot functions working on mkinfit objects.</p> <div class="sourceCode"><pre class="sourceCode r"><code class="sourceCode r"><span class="kw">summary</span>(mm.L3[[<span class="st">"DFOP"</span>, <span class="dv">1</span>]])</code></pre></div> -<pre><code>## mkin version: 0.9.44.9000 +<pre><code>## mkin version: 0.9.45 ## R version: 3.3.2 -## Date of fit: Thu Dec 8 08:18:50 2016 -## Date of summary: Thu Dec 8 08:18:51 2016 +## Date of fit: Thu Dec 8 09:23:17 2016 +## Date of summary: Thu Dec 8 09:23:17 2016 ## ## Equations: ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * @@ -435,7 +438,7 @@ mm.L3 <-<span class="st"> </span><span class="kw"><a href="../reference/mmkin ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 137 model solutions performed in 0.311 s +## Fitted with method Port using 137 model solutions performed in 0.309 s ## ## Weighting: none ## @@ -523,17 +526,17 @@ mm.L4 <-<span class="st"> </span><span class="kw"><a href="../reference/mmkin <p><img src="FOCUS_L_files/figure-html/unnamed-chunk-15-1.png" width="672"></p> <p>The <span class="math inline">\(\chi^2\)</span> error level of 3.3% as well as the plot suggest that the SFO model fits very well. The error level at which the <span class="math inline">\(\chi^2\)</span> test passes is slightly lower for the FOMC model. However, the difference appears negligible.</p> <div class="sourceCode"><pre class="sourceCode r"><code class="sourceCode r"><span class="kw">summary</span>(mm.L4[[<span class="st">"SFO"</span>, <span class="dv">1</span>]], <span class="dt">data =</span> <span class="ot">FALSE</span>)</code></pre></div> -<pre><code>## mkin version: 0.9.44.9000 +<pre><code>## mkin version: 0.9.45 ## R version: 3.3.2 -## Date of fit: Thu Dec 8 08:18:51 2016 -## Date of summary: Thu Dec 8 08:18:51 2016 +## Date of fit: Thu Dec 8 09:23:17 2016 +## Date of summary: Thu Dec 8 09:23:17 2016 ## ## Equations: ## d_parent/dt = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 46 model solutions performed in 0.118 s +## Fitted with method Port using 46 model solutions performed in 0.106 s ## ## Weighting: none ## @@ -583,17 +586,17 @@ mm.L4 <-<span class="st"> </span><span class="kw"><a href="../reference/mmkin ## DT50 DT90 ## parent 106 352</code></pre> <div class="sourceCode"><pre class="sourceCode r"><code class="sourceCode r"><span class="kw">summary</span>(mm.L4[[<span class="st">"FOMC"</span>, <span class="dv">1</span>]], <span class="dt">data =</span> <span class="ot">FALSE</span>)</code></pre></div> -<pre><code>## mkin version: 0.9.44.9000 +<pre><code>## mkin version: 0.9.45 ## R version: 3.3.2 -## Date of fit: Thu Dec 8 08:18:51 2016 -## Date of summary: Thu Dec 8 08:18:51 2016 +## Date of fit: Thu Dec 8 09:23:17 2016 +## Date of summary: Thu Dec 8 09:23:17 2016 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 66 model solutions performed in 0.151 s +## Fitted with method Port using 66 model solutions performed in 0.148 s ## ## Weighting: none ## |