aboutsummaryrefslogtreecommitdiff
path: root/docs
diff options
context:
space:
mode:
Diffstat (limited to 'docs')
-rw-r--r--docs/articles/FOCUS_L.html56
-rw-r--r--docs/articles/web_only/compiled_models.html12
-rw-r--r--docs/news/index.html1
-rw-r--r--docs/reference/Extract.mmkin.html12
-rw-r--r--docs/reference/mccall81_245T.html12
-rw-r--r--docs/reference/mkinfit.html46
-rw-r--r--docs/reference/mkinmod.html2
-rw-r--r--docs/reference/mkinpredict.html76
-rw-r--r--docs/reference/mmkin.html4
-rw-r--r--docs/reference/summary.mkinfit.html4
-rw-r--r--docs/reference/transform_odeparms.html24
11 files changed, 145 insertions, 104 deletions
diff --git a/docs/articles/FOCUS_L.html b/docs/articles/FOCUS_L.html
index 8ef55b81..ea49ae10 100644
--- a/docs/articles/FOCUS_L.html
+++ b/docs/articles/FOCUS_L.html
@@ -108,17 +108,17 @@
<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" id="cb2"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb2-1" data-line-number="1">m.L1.SFO &lt;-<span class="st"> </span><span class="kw"><a href="../reference/mkinfit.html">mkinfit</a></span>(<span class="st">"SFO"</span>, FOCUS_<span class="dv">2006</span>_L1_mkin, <span class="dt">quiet =</span> <span class="ot">TRUE</span>)</a>
<a class="sourceLine" id="cb2-2" data-line-number="2"><span class="kw">summary</span>(m.L1.SFO)</a></code></pre></div>
-<pre><code>## mkin version used for fitting: 0.9.47.5
+<pre><code>## mkin version used for fitting: 0.9.47.6
## R version used for fitting: 3.5.1
-## Date of fit: Tue Sep 18 09:50:45 2018
-## Date of summary: Tue Sep 18 09:50:45 2018
+## Date of fit: Tue Sep 18 15:18:12 2018
+## Date of summary: Tue Sep 18 15:18:12 2018
##
## 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.081 s
+## Fitted with method Port using 37 model solutions performed in 0.084 s
##
## Weighting: none
##
@@ -199,17 +199,17 @@
<a class="sourceLine" id="cb6-2" data-line-number="2"><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>)</a></code></pre></div>
<p><img src="FOCUS_L_files/figure-html/unnamed-chunk-6-1.png" width="576"></p>
<div class="sourceCode" id="cb7"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb7-1" data-line-number="1"><span class="kw">summary</span>(m.L1.FOMC, <span class="dt">data =</span> <span class="ot">FALSE</span>)</a></code></pre></div>
-<pre><code>## mkin version used for fitting: 0.9.47.5
+<pre><code>## mkin version used for fitting: 0.9.47.6
## R version used for fitting: 3.5.1
-## Date of fit: Tue Sep 18 09:50:47 2018
-## Date of summary: Tue Sep 18 09:50:47 2018
+## Date of fit: Tue Sep 18 15:18:14 2018
+## Date of summary: Tue Sep 18 15:18:14 2018
##
## Equations:
## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent
##
## Model predictions using solution type analytical
##
-## Fitted with method Port using 611 model solutions performed in 1.294 s
+## Fitted with method Port using 611 model solutions performed in 1.321 s
##
## Weighting: none
##
@@ -294,17 +294,17 @@
<a class="sourceLine" id="cb11-3" data-line-number="3"> <span class="dt">main =</span> <span class="st">"FOCUS L2 - FOMC"</span>)</a></code></pre></div>
<p><img src="FOCUS_L_files/figure-html/unnamed-chunk-9-1.png" width="672"></p>
<div class="sourceCode" id="cb12"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb12-1" data-line-number="1"><span class="kw">summary</span>(m.L2.FOMC, <span class="dt">data =</span> <span class="ot">FALSE</span>)</a></code></pre></div>
-<pre><code>## mkin version used for fitting: 0.9.47.5
+<pre><code>## mkin version used for fitting: 0.9.47.6
## R version used for fitting: 3.5.1
-## Date of fit: Tue Sep 18 09:50:47 2018
-## Date of summary: Tue Sep 18 09:50:47 2018
+## Date of fit: Tue Sep 18 15:18:14 2018
+## Date of summary: Tue Sep 18 15:18:14 2018
##
## 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.166 s
+## Fitted with method Port using 81 model solutions performed in 0.169 s
##
## Weighting: none
##
@@ -365,10 +365,10 @@
<a class="sourceLine" id="cb14-3" data-line-number="3"> <span class="dt">main =</span> <span class="st">"FOCUS L2 - DFOP"</span>)</a></code></pre></div>
<p><img src="FOCUS_L_files/figure-html/unnamed-chunk-10-1.png" width="672"></p>
<div class="sourceCode" id="cb15"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb15-1" data-line-number="1"><span class="kw">summary</span>(m.L2.DFOP, <span class="dt">data =</span> <span class="ot">FALSE</span>)</a></code></pre></div>
-<pre><code>## mkin version used for fitting: 0.9.47.5
+<pre><code>## mkin version used for fitting: 0.9.47.6
## R version used for fitting: 3.5.1
-## Date of fit: Tue Sep 18 09:50:48 2018
-## Date of summary: Tue Sep 18 09:50:48 2018
+## Date of fit: Tue Sep 18 15:18:15 2018
+## Date of summary: Tue Sep 18 15:18:15 2018
##
## Equations:
## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) *
@@ -377,7 +377,7 @@
##
## Model predictions using solution type analytical
##
-## Fitted with method Port using 336 model solutions performed in 0.714 s
+## Fitted with method Port using 336 model solutions performed in 0.739 s
##
## Weighting: none
##
@@ -457,10 +457,10 @@
<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" id="cb21"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb21-1" data-line-number="1"><span class="kw">summary</span>(mm.L3[[<span class="st">"DFOP"</span>, <span class="dv">1</span>]])</a></code></pre></div>
-<pre><code>## mkin version used for fitting: 0.9.47.5
+<pre><code>## mkin version used for fitting: 0.9.47.6
## R version used for fitting: 3.5.1
-## Date of fit: Tue Sep 18 09:50:49 2018
-## Date of summary: Tue Sep 18 09:50:49 2018
+## Date of fit: Tue Sep 18 15:18:16 2018
+## Date of summary: Tue Sep 18 15:18:16 2018
##
## Equations:
## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) *
@@ -469,7 +469,7 @@
##
## Model predictions using solution type analytical
##
-## Fitted with method Port using 137 model solutions performed in 0.299 s
+## Fitted with method Port using 137 model solutions performed in 0.301 s
##
## Weighting: none
##
@@ -558,17 +558,17 @@
<p><img src="FOCUS_L_files/figure-html/unnamed-chunk-15-1.png" width="700"></p>
<p>The <span class="math inline">\(\chi^2\)</span> error level of 3.3% as well as the plot suggest that the SFO model fits very well. The error level at which the <span class="math inline">\(\chi^2\)</span> test passes is slightly lower for the FOMC model. However, the difference appears negligible.</p>
<div class="sourceCode" id="cb26"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb26-1" data-line-number="1"><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>)</a></code></pre></div>
-<pre><code>## mkin version used for fitting: 0.9.47.5
+<pre><code>## mkin version used for fitting: 0.9.47.6
## R version used for fitting: 3.5.1
-## Date of fit: Tue Sep 18 09:50:50 2018
-## Date of summary: Tue Sep 18 09:50:50 2018
+## Date of fit: Tue Sep 18 15:18:17 2018
+## Date of summary: Tue Sep 18 15:18:17 2018
##
## 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.093 s
+## Fitted with method Port using 46 model solutions performed in 0.107 s
##
## Weighting: none
##
@@ -618,17 +618,17 @@
## DT50 DT90
## parent 106 352</code></pre>
<div class="sourceCode" id="cb28"><pre class="sourceCode r"><code class="sourceCode r"><a class="sourceLine" id="cb28-1" data-line-number="1"><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>)</a></code></pre></div>
-<pre><code>## mkin version used for fitting: 0.9.47.5
+<pre><code>## mkin version used for fitting: 0.9.47.6
## R version used for fitting: 3.5.1
-## Date of fit: Tue Sep 18 09:50:50 2018
-## Date of summary: Tue Sep 18 09:50:50 2018
+## Date of fit: Tue Sep 18 15:18:17 2018
+## Date of summary: Tue Sep 18 15:18:17 2018
##
## 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.138 s
+## Fitted with method Port using 66 model solutions performed in 0.139 s
##
## Weighting: none
##
diff --git a/docs/articles/web_only/compiled_models.html b/docs/articles/web_only/compiled_models.html
index 3b3e9150..684f28cb 100644
--- a/docs/articles/web_only/compiled_models.html
+++ b/docs/articles/web_only/compiled_models.html
@@ -125,9 +125,9 @@
<a class="sourceLine" id="cb5-16" data-line-number="16">}</a></code></pre></div>
<pre><code>## Lade nötiges Paket: rbenchmark</code></pre>
<pre><code>## test replications elapsed relative user.self sys.self
-## 3 deSolve, compiled 3 2.001 1.000 1.999 0
-## 1 deSolve, not compiled 3 15.374 7.683 15.368 0
-## 2 Eigenvalue based 3 2.439 1.219 2.438 0
+## 3 deSolve, compiled 3 2.005 1.000 2.003 0
+## 1 deSolve, not compiled 3 15.480 7.721 15.474 0
+## 2 Eigenvalue based 3 2.460 1.227 2.459 0
## user.child sys.child
## 3 0 0
## 1 0 0
@@ -156,13 +156,13 @@
<a class="sourceLine" id="cb8-16" data-line-number="16">}</a></code></pre></div>
<pre><code>## Successfully compiled differential equation model from auto-generated C code.</code></pre>
<pre><code>## test replications elapsed relative user.self sys.self
-## 2 deSolve, compiled 3 3.608 1.000 3.607 0
-## 1 deSolve, not compiled 3 33.150 9.188 33.135 0
+## 2 deSolve, compiled 3 3.524 1.00 3.523 0.000
+## 1 deSolve, not compiled 3 33.021 9.37 33.001 0.004
## user.child sys.child
## 2 0 0
## 1 0 0</code></pre>
<p>Here we get a performance benefit of a factor of 9 using the version of the differential equation model compiled from C code!</p>
-<p>This vignette was built with mkin 0.9.47.5 on</p>
+<p>This vignette was built with mkin 0.9.47.6 on</p>
<pre><code>## R version 3.5.1 (2018-07-02)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Debian GNU/Linux 9 (stretch)</code></pre>
diff --git a/docs/news/index.html b/docs/news/index.html
index 2636b8b0..a316612a 100644
--- a/docs/news/index.html
+++ b/docs/news/index.html
@@ -124,6 +124,7 @@
<ul>
<li><p>‘add_err’: Respect the argument giving the number of replicates in the synthetic dataset</p></li>
<li><p>‘max_twa_parent’: Support maximum time weighted average concentration calculations for the hockey stick (HS) model</p></li>
+<li><p>‘mkinpredict’: Make the function generic and create a method for mkinfit objects</p></li>
</ul>
</div>
<div id="mkin-0-9-47-5-2018-09-14" class="section level1">
diff --git a/docs/reference/Extract.mmkin.html b/docs/reference/Extract.mmkin.html
index 39a8bb84..0447e2d9 100644
--- a/docs/reference/Extract.mmkin.html
+++ b/docs/reference/Extract.mmkin.html
@@ -278,7 +278,7 @@
#&gt;
#&gt; $time
#&gt; User System verstrichen
-#&gt; 0.161 0.000 0.162
+#&gt; 0.175 0.000 0.175
#&gt;
#&gt; $mkinmod
#&gt; &lt;mkinmod&gt; model generated with
@@ -467,8 +467,8 @@
#&gt; }
#&gt; return(mC)
#&gt; }
-#&gt; &lt;bytecode: 0x55555bd56848&gt;
-#&gt; &lt;environment: 0x55555bd53cb0&gt;
+#&gt; &lt;bytecode: 0x55555bd99510&gt;
+#&gt; &lt;environment: 0x55555bc57c88&gt;
#&gt;
#&gt; $cost_notrans
#&gt; function (P)
@@ -490,8 +490,8 @@
#&gt; scaleVar = scaleVar)
#&gt; return(mC)
#&gt; }
-#&gt; &lt;bytecode: 0x55555c9bc2c0&gt;
-#&gt; &lt;environment: 0x55555bd53cb0&gt;
+#&gt; &lt;bytecode: 0x55555cb90200&gt;
+#&gt; &lt;environment: 0x55555bc57c88&gt;
#&gt;
#&gt; $hessian_notrans
#&gt; parent_0 k_parent_sink
@@ -558,7 +558,7 @@
#&gt; 99.17407
#&gt;
#&gt; $date
-#&gt; [1] "Tue Sep 18 09:49:15 2018"
+#&gt; [1] "Tue Sep 18 15:16:36 2018"
#&gt;
#&gt; $version
#&gt; [1] "0.9.47.6"
diff --git a/docs/reference/mccall81_245T.html b/docs/reference/mccall81_245T.html
index 042073f2..be77df76 100644
--- a/docs/reference/mccall81_245T.html
+++ b/docs/reference/mccall81_245T.html
@@ -158,8 +158,8 @@
<span class='kw'>anisole</span> <span class='kw'>=</span> <span class='fu'>list</span>(<span class='kw'>type</span> <span class='kw'>=</span> <span class='st'>"SFO"</span>))</div><div class='output co'>#&gt; <span class='message'>Successfully compiled differential equation model from auto-generated C code.</span></div><div class='input'> </div><div class='input'> <span class='no'>fit.1</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinfit.html'>mkinfit</a></span>(<span class='no'>SFO_SFO_SFO</span>, <span class='fu'>subset</span>(<span class='no'>mccall81_245T</span>, <span class='no'>soil</span> <span class='kw'>==</span> <span class='st'>"Commerce"</span>), <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>fit.1</span>, <span class='kw'>data</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:26 2018
-#&gt; Date of summary: Tue Sep 18 09:49:26 2018
+#&gt; Date of fit: Tue Sep 18 15:16:48 2018
+#&gt; Date of summary: Tue Sep 18 15:16:48 2018
#&gt;
#&gt; Equations:
#&gt; d_T245/dt = - k_T245_sink * T245 - k_T245_phenol * T245
@@ -169,7 +169,7 @@
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 574 model solutions performed in 3.213 s
+#&gt; Fitted with method Port using 574 model solutions performed in 3.602 s
#&gt;
#&gt; Weighting: none
#&gt;
@@ -247,8 +247,8 @@
<span class='kw'>fixed_parms</span> <span class='kw'>=</span> <span class='st'>"k_phenol_sink"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>fit.2</span>, <span class='kw'>data</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:27 2018
-#&gt; Date of summary: Tue Sep 18 09:49:27 2018
+#&gt; Date of fit: Tue Sep 18 15:16:50 2018
+#&gt; Date of summary: Tue Sep 18 15:16:50 2018
#&gt;
#&gt; Equations:
#&gt; d_T245/dt = - k_T245_sink * T245 - k_T245_phenol * T245
@@ -258,7 +258,7 @@
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 246 model solutions performed in 1.363 s
+#&gt; Fitted with method Port using 246 model solutions performed in 1.513 s
#&gt;
#&gt; Weighting: none
#&gt;
diff --git a/docs/reference/mkinfit.html b/docs/reference/mkinfit.html
index 1ae53f8a..3c8442b6 100644
--- a/docs/reference/mkinfit.html
+++ b/docs/reference/mkinfit.html
@@ -434,15 +434,15 @@
<span class='no'>fit</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='st'>"FOMC"</span>, <span class='no'>FOCUS_2006_C</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:30 2018
-#&gt; Date of summary: Tue Sep 18 09:49:30 2018
+#&gt; Date of fit: Tue Sep 18 15:16:52 2018
+#&gt; Date of summary: Tue Sep 18 15:16:52 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent
#&gt;
#&gt; Model predictions using solution type analytical
#&gt;
-#&gt; Fitted with method Port using 64 model solutions performed in 0.14 s
+#&gt; Fitted with method Port using 64 model solutions performed in 0.158 s
#&gt;
#&gt; Weighting: none
#&gt;
@@ -511,7 +511,7 @@
<span class='kw'>m1</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>))</div><div class='output co'>#&gt; <span class='message'>Successfully compiled differential equation model from auto-generated C code.</span></div><div class='input'><span class='co'># Fit the model to the FOCUS example dataset D using defaults</span>
<span class='fu'>print</span>(<span class='fu'>system.time</span>(<span class='no'>fit</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO</span>, <span class='no'>FOCUS_2006_D</span>,
<span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"eigen"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)))</div><div class='output co'>#&gt; User System verstrichen
-#&gt; 0.872 0.001 0.874 </div><div class='input'><span class='fu'>coef</span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; parent_0 log_k_parent_sink log_k_parent_m1 log_k_m1_sink
+#&gt; 1.011 0.000 1.012 </div><div class='input'><span class='fu'>coef</span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; parent_0 log_k_parent_sink log_k_parent_m1 log_k_m1_sink
#&gt; 99.59848 -3.03822 -2.98030 -5.24750 </div><div class='input'><span class='fu'><a href='endpoints.html'>endpoints</a></span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; $ff
#&gt; parent_sink parent_m1 m1_sink
#&gt; 0.485524 0.514476 1.000000
@@ -586,7 +586,7 @@
#&gt; Model cost at call 146 : 371.2134
#&gt; Optimisation by method Port successfully terminated.
#&gt; User System verstrichen
-#&gt; 0.690 0.000 0.689 </div><div class='input'><span class='fu'>coef</span>(<span class='no'>fit.deSolve</span>)</div><div class='output co'>#&gt; parent_0 log_k_parent_sink log_k_parent_m1 log_k_m1_sink
+#&gt; 0.755 0.000 0.755 </div><div class='input'><span class='fu'>coef</span>(<span class='no'>fit.deSolve</span>)</div><div class='output co'>#&gt; parent_0 log_k_parent_sink log_k_parent_m1 log_k_m1_sink
#&gt; 99.59848 -3.03822 -2.98030 -5.24750 </div><div class='input'><span class='fu'><a href='endpoints.html'>endpoints</a></span>(<span class='no'>fit.deSolve</span>)</div><div class='output co'>#&gt; $ff
#&gt; parent_sink parent_m1 m1_sink
#&gt; 0.485524 0.514476 1.000000
@@ -624,8 +624,8 @@
<span class='kw'>m1</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>), <span class='kw'>use_of_ff</span> <span class='kw'>=</span> <span class='st'>"max"</span>)</div><div class='output co'>#&gt; <span class='message'>Successfully compiled differential equation model from auto-generated C code.</span></div><div class='input'><span class='no'>f.noweight</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>f.noweight</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:40 2018
-#&gt; Date of summary: Tue Sep 18 09:49:40 2018
+#&gt; Date of fit: Tue Sep 18 15:17:04 2018
+#&gt; Date of summary: Tue Sep 18 15:17:04 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -633,7 +633,7 @@
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 186 model solutions performed in 0.725 s
+#&gt; Fitted with method Port using 186 model solutions performed in 0.823 s
#&gt;
#&gt; Weighting: none
#&gt;
@@ -741,8 +741,8 @@
#&gt; 120 m1 33.31 28.78984 4.520e+00</div><div class='input'><span class='no'>f.irls</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>reweight.method</span> <span class='kw'>=</span> <span class='st'>"obs"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>f.irls</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:42 2018
-#&gt; Date of summary: Tue Sep 18 09:49:42 2018
+#&gt; Date of fit: Tue Sep 18 15:17:06 2018
+#&gt; Date of summary: Tue Sep 18 15:17:06 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -750,7 +750,7 @@
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 551 model solutions performed in 2.147 s
+#&gt; Fitted with method Port using 551 model solutions performed in 2.374 s
#&gt;
#&gt; Weighting: none
#&gt;
@@ -863,8 +863,8 @@
#&gt; 120 m1 33.31 28.80898 4.501e+00 2.722</div><div class='input'><span class='no'>f.w.mean</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>weight</span> <span class='kw'>=</span> <span class='st'>"mean"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>f.w.mean</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:43 2018
-#&gt; Date of summary: Tue Sep 18 09:49:43 2018
+#&gt; Date of fit: Tue Sep 18 15:17:07 2018
+#&gt; Date of summary: Tue Sep 18 15:17:07 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -872,7 +872,7 @@
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 155 model solutions performed in 0.603 s
+#&gt; Fitted with method Port using 155 model solutions performed in 0.636 s
#&gt;
#&gt; Weighting: mean
#&gt;
@@ -981,8 +981,8 @@
<span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>f.w.value</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:44 2018
-#&gt; Date of summary: Tue Sep 18 09:49:44 2018
+#&gt; Date of fit: Tue Sep 18 15:17:08 2018
+#&gt; Date of summary: Tue Sep 18 15:17:08 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -990,7 +990,7 @@
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 174 model solutions performed in 0.69 s
+#&gt; Fitted with method Port using 174 model solutions performed in 0.717 s
#&gt;
#&gt; Weighting: manual
#&gt;
@@ -1101,8 +1101,8 @@
<span class='no'>f.w.man</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>dw</span>, <span class='kw'>err</span> <span class='kw'>=</span> <span class='st'>"err.man"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>f.w.man</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:45 2018
-#&gt; Date of summary: Tue Sep 18 09:49:45 2018
+#&gt; Date of fit: Tue Sep 18 15:17:09 2018
+#&gt; Date of summary: Tue Sep 18 15:17:09 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -1110,7 +1110,7 @@
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 270 model solutions performed in 1.133 s
+#&gt; Fitted with method Port using 270 model solutions performed in 1.181 s
#&gt;
#&gt; Weighting: manual
#&gt;
@@ -1219,8 +1219,8 @@
<span class='kw'>reweight.method</span> <span class='kw'>=</span> <span class='st'>"obs"</span>)
<span class='fu'>summary</span>(<span class='no'>f.w.man.irls</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:49:48 2018
-#&gt; Date of summary: Tue Sep 18 09:49:48 2018
+#&gt; Date of fit: Tue Sep 18 15:17:13 2018
+#&gt; Date of summary: Tue Sep 18 15:17:13 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -1228,7 +1228,7 @@
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 692 model solutions performed in 2.77 s
+#&gt; Fitted with method Port using 692 model solutions performed in 3.077 s
#&gt;
#&gt; Weighting: manual
#&gt;
diff --git a/docs/reference/mkinmod.html b/docs/reference/mkinmod.html
index c1e7b840..af364403 100644
--- a/docs/reference/mkinmod.html
+++ b/docs/reference/mkinmod.html
@@ -229,7 +229,7 @@ For the definition of model types and their parameters, the equations given
<span class='no'>SFO_SFO</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinmod</span>(
<span class='kw'>parent</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>, <span class='st'>"m1"</span>),
<span class='kw'>m1</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>), <span class='kw'>verbose</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#&gt; Compilation argument:
-#&gt; /usr/lib/R/bin/R CMD SHLIB file73da22034c44.c 2&gt; file73da22034c44.c.err.txt
+#&gt; /usr/lib/R/bin/R CMD SHLIB file20267955fc75.c 2&gt; file20267955fc75.c.err.txt
#&gt; Program source:
#&gt; 1: #include &lt;R.h&gt;
#&gt; 2:
diff --git a/docs/reference/mkinpredict.html b/docs/reference/mkinpredict.html
index 5fad64de..91e9289d 100644
--- a/docs/reference/mkinpredict.html
+++ b/docs/reference/mkinpredict.html
@@ -131,16 +131,18 @@
</div>
- <pre class="usage"><span class='fu'>mkinpredict</span>(<span class='no'>mkinmod</span>, <span class='no'>odeparms</span>, <span class='no'>odeini</span>, <span class='no'>outtimes</span>, <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"deSolve"</span>,
- <span class='kw'>use_compiled</span> <span class='kw'>=</span> <span class='st'>"auto"</span>, <span class='kw'>method.ode</span> <span class='kw'>=</span> <span class='st'>"lsoda"</span>, <span class='kw'>atol</span> <span class='kw'>=</span> <span class='fl'>1e-08</span>, <span class='kw'>rtol</span> <span class='kw'>=</span> <span class='fl'>1e-10</span>,
- <span class='kw'>map_output</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>, <span class='no'>...</span>)</pre>
+ <pre class="usage"><span class='fu'>mkinpredict</span>(<span class='no'>x</span>, <span class='no'>odeparms</span>, <span class='no'>odeini</span>, <span class='kw'>outtimes</span> <span class='kw'>=</span> <span class='fu'>seq</span>(<span class='fl'>0</span>, <span class='fl'>120</span>, <span class='kw'>by</span> <span class='kw'>=</span> <span class='fl'>0.1</span>),
+ <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"deSolve"</span>, <span class='kw'>use_compiled</span> <span class='kw'>=</span> <span class='st'>"auto"</span>, <span class='kw'>method.ode</span> <span class='kw'>=</span> <span class='st'>"lsoda"</span>,
+ <span class='kw'>atol</span> <span class='kw'>=</span> <span class='fl'>1e-08</span>, <span class='kw'>rtol</span> <span class='kw'>=</span> <span class='fl'>1e-10</span>, <span class='kw'>map_output</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>, <span class='no'>...</span>)</pre>
<h2 class="hasAnchor" id="arguments"><a class="anchor" href="#arguments"></a>Arguments</h2>
<table class="ref-arguments">
<colgroup><col class="name" /><col class="desc" /></colgroup>
<tr>
- <th>mkinmod</th>
- <td><p>A kinetic model as produced by <code><a href='mkinmod.html'>mkinmod</a></code>.</p></td>
+ <th>x</th>
+ <td><p>A kinetic model as produced by <code><a href='mkinmod.html'>mkinmod</a></code>, or a kinetic
+ fit as fitted by <code><a href='mkinfit.html'>mkinfit</a></code>. In the latter case, the fitted
+ parameters are used for the prediction.</p></td>
</tr>
<tr>
<th>odeparms</th>
@@ -204,10 +206,10 @@
<h2 class="hasAnchor" id="examples"><a class="anchor" href="#examples"></a>Examples</h2>
- <pre class="examples"><div class='input'> <span class='no'>SFO</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinmod.html'>mkinmod</a></span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fu'>list</span>(<span class='kw'>type</span> <span class='kw'>=</span> <span class='st'>"SFO"</span>))
+ <pre class="examples"><div class='input'> <span class='no'>SFO</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinmod.html'>mkinmod</a></span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>))
<span class='co'># Compare solution types</span>
<span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>), <span class='fl'>0</span>:<span class='fl'>20</span>,
- <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"analytical"</span>)</div><div class='output co'>#&gt; time degradinol
+ <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"analytical"</span>)</div><div class='output co'>#&gt; time degradinol
#&gt; 1 0 100.0000000
#&gt; 2 1 74.0818221
#&gt; 3 2 54.8811636
@@ -229,7 +231,7 @@
#&gt; 19 18 0.4516581
#&gt; 20 19 0.3345965
#&gt; 21 20 0.2478752</div><div class='input'> <span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>), <span class='fl'>0</span>:<span class='fl'>20</span>,
- <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"deSolve"</span>)</div><div class='output co'>#&gt; time degradinol
+ <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"deSolve"</span>)</div><div class='output co'>#&gt; time degradinol
#&gt; 1 0 100.0000000
#&gt; 2 1 74.0818221
#&gt; 3 2 54.8811636
@@ -251,7 +253,7 @@
#&gt; 19 18 0.4516581
#&gt; 20 19 0.3345965
#&gt; 21 20 0.2478752</div><div class='input'> <span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>), <span class='fl'>0</span>:<span class='fl'>20</span>,
- <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"deSolve"</span>, <span class='kw'>use_compiled</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>)</div><div class='output co'>#&gt; time degradinol
+ <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"deSolve"</span>, <span class='kw'>use_compiled</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>)</div><div class='output co'>#&gt; time degradinol
#&gt; 1 0 100.0000000
#&gt; 2 1 74.0818221
#&gt; 3 2 54.8811636
@@ -273,7 +275,7 @@
#&gt; 19 18 0.4516581
#&gt; 20 19 0.3345965
#&gt; 21 20 0.2478752</div><div class='input'> <span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>), <span class='fl'>0</span>:<span class='fl'>20</span>,
- <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"eigen"</span>)</div><div class='output co'>#&gt; time degradinol
+ <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"eigen"</span>)</div><div class='output co'>#&gt; time degradinol
#&gt; 1 0 100.0000000
#&gt; 2 1 74.0818221
#&gt; 3 2 54.8811636
@@ -298,21 +300,21 @@
<span class='co'># Compare integration methods to analytical solution</span>
<span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>), <span class='fl'>0</span>:<span class='fl'>20</span>,
- <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"analytical"</span>)[<span class='fl'>21</span>,]</div><div class='output co'>#&gt; time degradinol
+ <span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"analytical"</span>)[<span class='fl'>21</span>,]</div><div class='output co'>#&gt; time degradinol
#&gt; 21 20 0.2478752</div><div class='input'> <span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>), <span class='fl'>0</span>:<span class='fl'>20</span>,
- <span class='kw'>method</span> <span class='kw'>=</span> <span class='st'>"lsoda"</span>)[<span class='fl'>21</span>,]</div><div class='output co'>#&gt; time degradinol
+ <span class='kw'>method</span> <span class='kw'>=</span> <span class='st'>"lsoda"</span>)[<span class='fl'>21</span>,]</div><div class='output co'>#&gt; time degradinol
#&gt; 21 20 0.2478752</div><div class='input'> <span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>), <span class='fl'>0</span>:<span class='fl'>20</span>,
- <span class='kw'>method</span> <span class='kw'>=</span> <span class='st'>"ode45"</span>)[<span class='fl'>21</span>,]</div><div class='output co'>#&gt; time degradinol
+ <span class='kw'>method</span> <span class='kw'>=</span> <span class='st'>"ode45"</span>)[<span class='fl'>21</span>,]</div><div class='output co'>#&gt; time degradinol
#&gt; 21 20 0.2478752</div><div class='input'> <span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>), <span class='fl'>0</span>:<span class='fl'>20</span>,
- <span class='kw'>method</span> <span class='kw'>=</span> <span class='st'>"rk4"</span>)[<span class='fl'>21</span>,]</div><div class='output co'>#&gt; time degradinol
+ <span class='kw'>method</span> <span class='kw'>=</span> <span class='st'>"rk4"</span>)[<span class='fl'>21</span>,]</div><div class='output co'>#&gt; time degradinol
#&gt; 21 20 0.2480043</div><div class='input'> <span class='co'># rk4 is not as precise here</span>
<span class='co'># The number of output times used to make a lot of difference until the</span>
<span class='co'># default for atol was adjusted</span>
<span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>),
- <span class='fu'>seq</span>(<span class='fl'>0</span>, <span class='fl'>20</span>, <span class='kw'>by</span> <span class='kw'>=</span> <span class='fl'>0.1</span>))[<span class='fl'>201</span>,]</div><div class='output co'>#&gt; time degradinol
+ <span class='fu'>seq</span>(<span class='fl'>0</span>, <span class='fl'>20</span>, <span class='kw'>by</span> <span class='kw'>=</span> <span class='fl'>0.1</span>))[<span class='fl'>201</span>,]</div><div class='output co'>#&gt; time degradinol
#&gt; 201 20 0.2478752</div><div class='input'> <span class='fu'>mkinpredict</span>(<span class='no'>SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_degradinol_sink</span> <span class='kw'>=</span> <span class='fl'>0.3</span>), <span class='fu'>c</span>(<span class='kw'>degradinol</span> <span class='kw'>=</span> <span class='fl'>100</span>),
- <span class='fu'>seq</span>(<span class='fl'>0</span>, <span class='fl'>20</span>, <span class='kw'>by</span> <span class='kw'>=</span> <span class='fl'>0.01</span>))[<span class='fl'>2001</span>,]</div><div class='output co'>#&gt; time degradinol
+ <span class='fu'>seq</span>(<span class='fl'>0</span>, <span class='fl'>20</span>, <span class='kw'>by</span> <span class='kw'>=</span> <span class='fl'>0.01</span>))[<span class='fl'>2001</span>,]</div><div class='output co'>#&gt; time degradinol
#&gt; 2001 20 0.2478752</div><div class='input'>
<span class='co'># Check compiled model versions - they are faster than the eigenvalue based solutions!</span>
<span class='no'>SFO_SFO</span> <span class='kw'>=</span> <span class='fu'><a href='mkinmod.html'>mkinmod</a></span>(<span class='kw'>parent</span> <span class='kw'>=</span> <span class='fu'>list</span>(<span class='kw'>type</span> <span class='kw'>=</span> <span class='st'>"SFO"</span>, <span class='kw'>to</span> <span class='kw'>=</span> <span class='st'>"m1"</span>),
@@ -321,7 +323,7 @@
<span class='fu'>c</span>(<span class='kw'>parent</span> <span class='kw'>=</span> <span class='fl'>100</span>, <span class='kw'>m1</span> <span class='kw'>=</span> <span class='fl'>0</span>), <span class='fu'>seq</span>(<span class='fl'>0</span>, <span class='fl'>20</span>, <span class='kw'>by</span> <span class='kw'>=</span> <span class='fl'>0.1</span>),
<span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"eigen"</span>)[<span class='fl'>201</span>,]))</div><div class='output co'>#&gt; time parent m1
#&gt; 201 20 4.978707 27.46227</div><div class='output co'>#&gt; User System verstrichen
-#&gt; 0.003 0.000 0.003 </div><div class='input'> <span class='fu'>system.time</span>(
+#&gt; 0.003 0.000 0.004 </div><div class='input'> <span class='fu'>system.time</span>(
<span class='fu'>print</span>(<span class='fu'>mkinpredict</span>(<span class='no'>SFO_SFO</span>, <span class='fu'>c</span>(<span class='kw'>k_parent_m1</span> <span class='kw'>=</span> <span class='fl'>0.05</span>, <span class='kw'>k_parent_sink</span> <span class='kw'>=</span> <span class='fl'>0.1</span>, <span class='kw'>k_m1_sink</span> <span class='kw'>=</span> <span class='fl'>0.01</span>),
<span class='fu'>c</span>(<span class='kw'>parent</span> <span class='kw'>=</span> <span class='fl'>100</span>, <span class='kw'>m1</span> <span class='kw'>=</span> <span class='fl'>0</span>), <span class='fu'>seq</span>(<span class='fl'>0</span>, <span class='fl'>20</span>, <span class='kw'>by</span> <span class='kw'>=</span> <span class='fl'>0.1</span>),
<span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"deSolve"</span>)[<span class='fl'>201</span>,]))</div><div class='output co'>#&gt; time parent m1
@@ -331,7 +333,45 @@
<span class='fu'>c</span>(<span class='kw'>parent</span> <span class='kw'>=</span> <span class='fl'>100</span>, <span class='kw'>m1</span> <span class='kw'>=</span> <span class='fl'>0</span>), <span class='fu'>seq</span>(<span class='fl'>0</span>, <span class='fl'>20</span>, <span class='kw'>by</span> <span class='kw'>=</span> <span class='fl'>0.1</span>),
<span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"deSolve"</span>, <span class='kw'>use_compiled</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>)[<span class='fl'>201</span>,]))</div><div class='output co'>#&gt; time parent m1
#&gt; 201 20 4.978707 27.46227</div><div class='output co'>#&gt; User System verstrichen
-#&gt; 0.033 0.000 0.034 </div></pre>
+#&gt; 0.034 0.000 0.034 </div><div class='input'>
+ </div><div class='input'> <span class='co'># Predict from a fitted model</span>
+ <span class='no'>f</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinfit.html'>mkinfit</a></span>(<span class='no'>SFO_SFO</span>, <span class='no'>FOCUS_2006_C</span>)</div><div class='output co'>#&gt; Model cost at call 1 : 552.5739
+#&gt; Model cost at call 3 : 552.5739
+#&gt; Model cost at call 4 : 552.5739
+#&gt; Model cost at call 6 : 279.9345
+#&gt; Model cost at call 7 : 279.9344
+#&gt; Model cost at call 8 : 279.9294
+#&gt; Model cost at call 9 : 279.9294
+#&gt; Model cost at call 12 : 200.3629
+#&gt; Model cost at call 13 : 200.3629
+#&gt; Model cost at call 18 : 197.9039
+#&gt; Model cost at call 23 : 197.9039
+#&gt; Model cost at call 25 : 196.6754
+#&gt; Model cost at call 27 : 196.6754
+#&gt; Model cost at call 32 : 196.5742
+#&gt; Model cost at call 33 : 196.5742
+#&gt; Model cost at call 34 : 196.5742
+#&gt; Model cost at call 38 : 196.5361
+#&gt; Model cost at call 40 : 196.5361
+#&gt; Model cost at call 44 : 196.5336
+#&gt; Model cost at call 45 : 196.5336
+#&gt; Model cost at call 50 : 196.5334
+#&gt; Model cost at call 51 : 196.5334
+#&gt; Model cost at call 52 : 196.5334
+#&gt; Model cost at call 56 : 196.5334
+#&gt; Model cost at call 58 : 196.5334
+#&gt; Model cost at call 59 : 196.5334
+#&gt; Model cost at call 65 : 196.5334
+#&gt; Model cost at call 73 : 196.5334
+#&gt; Model cost at call 78 : 196.5334
+#&gt; Model cost at call 80 : 196.5334
+#&gt; Optimisation by method Port successfully terminated.</div><div class='input'> <span class='fu'>head</span>(<span class='fu'>mkinpredict</span>(<span class='no'>f</span>))</div><div class='output co'>#&gt; time parent m1
+#&gt; 1 0.0 82.49216 0.000000
+#&gt; 2 0.1 80.00563 1.179955
+#&gt; 3 0.2 77.59404 2.312580
+#&gt; 4 0.3 75.25515 3.399419
+#&gt; 5 0.4 72.98675 4.441969
+#&gt; 6 0.5 70.78673 5.441679</div><div class='input'> </div></pre>
</div>
<div class="col-md-3 hidden-xs hidden-sm" id="sidebar">
<h2>Contents</h2>
diff --git a/docs/reference/mmkin.html b/docs/reference/mmkin.html
index 13cdedff..ad011b75 100644
--- a/docs/reference/mmkin.html
+++ b/docs/reference/mmkin.html
@@ -189,8 +189,8 @@
<span class='no'>time_1</span> <span class='kw'>&lt;-</span> <span class='fu'>system.time</span>(<span class='no'>fits.4</span> <span class='kw'>&lt;-</span> <span class='fu'>mmkin</span>(<span class='no'>models</span>, <span class='no'>datasets</span>, <span class='kw'>cores</span> <span class='kw'>=</span> <span class='fl'>1</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>))
<span class='no'>time_default</span></div><div class='output co'>#&gt; User System verstrichen
-#&gt; 0.043 0.040 6.475 </div><div class='input'><span class='no'>time_1</span></div><div class='output co'>#&gt; User System verstrichen
-#&gt; 19.165 0.000 19.179 </div><div class='input'>
+#&gt; 0.054 0.027 6.836 </div><div class='input'><span class='no'>time_1</span></div><div class='output co'>#&gt; User System verstrichen
+#&gt; 20.444 0.000 20.459 </div><div class='input'>
<span class='fu'><a href='endpoints.html'>endpoints</a></span>(<span class='no'>fits.0</span><span class='kw'>[[</span><span class='st'>"SFO_lin"</span>, <span class='fl'>2</span>]])</div><div class='output co'>#&gt; $ff
#&gt; parent_M1 parent_sink M1_M2 M1_sink
#&gt; 0.7340480 0.2659520 0.7505686 0.2494314
diff --git a/docs/reference/summary.mkinfit.html b/docs/reference/summary.mkinfit.html
index 0af65ed8..f8e09b45 100644
--- a/docs/reference/summary.mkinfit.html
+++ b/docs/reference/summary.mkinfit.html
@@ -206,8 +206,8 @@
<h2 class="hasAnchor" id="examples"><a class="anchor" href="#examples"></a>Examples</h2>
<pre class="examples"><div class='input'> <span class='fu'>summary</span>(<span class='fu'><a href='mkinfit.html'>mkinfit</a></span>(<span class='fu'><a href='mkinmod.html'>mkinmod</a></span>(<span class='kw'>parent</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>)), <span class='no'>FOCUS_2006_A</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>))</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:50:32 2018
-#&gt; Date of summary: Tue Sep 18 09:50:32 2018
+#&gt; Date of fit: Tue Sep 18 15:17:59 2018
+#&gt; Date of summary: Tue Sep 18 15:17:59 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent_sink * parent
diff --git a/docs/reference/transform_odeparms.html b/docs/reference/transform_odeparms.html
index ce6f6c8d..f2602792 100644
--- a/docs/reference/transform_odeparms.html
+++ b/docs/reference/transform_odeparms.html
@@ -200,8 +200,8 @@ The transformation of sets of formation fractions is fragile, as it supposes
<span class='no'>fit</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinfit.html'>mkinfit</a></span>(<span class='no'>SFO_SFO</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>fit</span>, <span class='kw'>data</span><span class='kw'>=</span><span class='fl'>FALSE</span>) <span class='co'># See transformed and backtransformed parameters</span></div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:50:39 2018
-#&gt; Date of summary: Tue Sep 18 09:50:39 2018
+#&gt; Date of fit: Tue Sep 18 15:18:06 2018
+#&gt; Date of summary: Tue Sep 18 15:18:06 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent_sink * parent - k_parent_m1 * parent
@@ -209,7 +209,7 @@ The transformation of sets of formation fractions is fragile, as it supposes
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 153 model solutions performed in 0.586 s
+#&gt; Fitted with method Port using 153 model solutions performed in 0.662 s
#&gt;
#&gt; Weighting: none
#&gt;
@@ -276,8 +276,8 @@ The transformation of sets of formation fractions is fragile, as it supposes
</div><div class='input'><span class='no'>fit.2</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinfit.html'>mkinfit</a></span>(<span class='no'>SFO_SFO</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>transform_rates</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>fit.2</span>, <span class='kw'>data</span><span class='kw'>=</span><span class='fl'>FALSE</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:50:41 2018
-#&gt; Date of summary: Tue Sep 18 09:50:41 2018
+#&gt; Date of fit: Tue Sep 18 15:18:08 2018
+#&gt; Date of summary: Tue Sep 18 15:18:08 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent_sink * parent - k_parent_m1 * parent
@@ -285,7 +285,7 @@ The transformation of sets of formation fractions is fragile, as it supposes
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 350 model solutions performed in 1.341 s
+#&gt; Fitted with method Port using 350 model solutions performed in 1.342 s
#&gt;
#&gt; Weighting: none
#&gt;
@@ -364,8 +364,8 @@ The transformation of sets of formation fractions is fragile, as it supposes
<span class='no'>fit.ff</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinfit.html'>mkinfit</a></span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>fit.ff</span>, <span class='kw'>data</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:50:42 2018
-#&gt; Date of summary: Tue Sep 18 09:50:42 2018
+#&gt; Date of fit: Tue Sep 18 15:18:09 2018
+#&gt; Date of summary: Tue Sep 18 15:18:09 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -373,7 +373,7 @@ The transformation of sets of formation fractions is fragile, as it supposes
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 186 model solutions performed in 0.73 s
+#&gt; Fitted with method Port using 186 model solutions performed in 0.737 s
#&gt;
#&gt; Weighting: none
#&gt;
@@ -448,8 +448,8 @@ The transformation of sets of formation fractions is fragile, as it supposes
<span class='no'>fit.ff.2</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinfit.html'>mkinfit</a></span>(<span class='no'>SFO_SFO.ff.2</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
<span class='fu'>summary</span>(<span class='no'>fit.ff.2</span>, <span class='kw'>data</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.47.6
#&gt; R version used for fitting: 3.5.1
-#&gt; Date of fit: Tue Sep 18 09:50:42 2018
-#&gt; Date of summary: Tue Sep 18 09:50:42 2018
+#&gt; Date of fit: Tue Sep 18 15:18:10 2018
+#&gt; Date of summary: Tue Sep 18 15:18:10 2018
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -457,7 +457,7 @@ The transformation of sets of formation fractions is fragile, as it supposes
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted with method Port using 104 model solutions performed in 0.401 s
+#&gt; Fitted with method Port using 104 model solutions performed in 0.414 s
#&gt;
#&gt; Weighting: none
#&gt;

Contact - Imprint