aboutsummaryrefslogtreecommitdiff
path: root/docs/reference/saem.html
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2022-08-10 13:21:34 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2022-08-10 13:21:34 +0200
commit9ee3d9f025ec7f5effddb0bcf9cf6e054c99794b (patch)
tree3ea151f52dbc17b85839100b5639d64301460d5e /docs/reference/saem.html
parent9e346fabe99de71b21ef085be102027cfa774910 (diff)
Update static docs
Diffstat (limited to 'docs/reference/saem.html')
-rw-r--r--docs/reference/saem.html295
1 files changed, 178 insertions, 117 deletions
diff --git a/docs/reference/saem.html b/docs/reference/saem.html
index e9898a73..1bc6308e 100644
--- a/docs/reference/saem.html
+++ b/docs/reference/saem.html
@@ -19,7 +19,7 @@ Expectation Maximisation algorithm (SAEM)."><!-- mathjax --><script src="https:/
</button>
<span class="navbar-brand">
<a class="navbar-link" href="../index.html">mkin</a>
- <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.1.0</span>
+ <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.1.2</span>
</span>
</div>
@@ -91,40 +91,40 @@ Expectation Maximisation algorithm (SAEM).</p>
</div>
<div id="ref-usage">
- <div class="sourceCode"><pre class="sourceCode r"><code><span class="fu">saem</span><span class="op">(</span><span class="va">object</span>, <span class="va">...</span><span class="op">)</span>
-
-<span class="co"># S3 method for mmkin</span>
-<span class="fu">saem</span><span class="op">(</span>
- <span class="va">object</span>,
- transformations <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"mkin"</span>, <span class="st">"saemix"</span><span class="op">)</span>,
- degparms_start <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/numeric.html" class="external-link">numeric</a></span><span class="op">(</span><span class="op">)</span>,
- test_log_parms <span class="op">=</span> <span class="cn">TRUE</span>,
- conf.level <span class="op">=</span> <span class="fl">0.6</span>,
- solution_type <span class="op">=</span> <span class="st">"auto"</span>,
- nbiter.saemix <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">300</span>, <span class="fl">100</span><span class="op">)</span>,
- control <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span>displayProgress <span class="op">=</span> <span class="cn">FALSE</span>, print <span class="op">=</span> <span class="cn">FALSE</span>, nbiter.saemix <span class="op">=</span> <span class="va">nbiter.saemix</span>,
- save <span class="op">=</span> <span class="cn">FALSE</span>, save.graphs <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span>,
- fail_with_errors <span class="op">=</span> <span class="cn">TRUE</span>,
- verbose <span class="op">=</span> <span class="cn">FALSE</span>,
- quiet <span class="op">=</span> <span class="cn">FALSE</span>,
- <span class="va">...</span>
-<span class="op">)</span>
-
-<span class="co"># S3 method for saem.mmkin</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">x</span>, digits <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/Extremes.html" class="external-link">max</a></span><span class="op">(</span><span class="fl">3</span>, <span class="fu"><a href="https://rdrr.io/r/base/options.html" class="external-link">getOption</a></span><span class="op">(</span><span class="st">"digits"</span><span class="op">)</span> <span class="op">-</span> <span class="fl">3</span><span class="op">)</span>, <span class="va">...</span><span class="op">)</span>
-
-<span class="fu">saemix_model</span><span class="op">(</span>
- <span class="va">object</span>,
- solution_type <span class="op">=</span> <span class="st">"auto"</span>,
- transformations <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"mkin"</span>, <span class="st">"saemix"</span><span class="op">)</span>,
- degparms_start <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/numeric.html" class="external-link">numeric</a></span><span class="op">(</span><span class="op">)</span>,
- test_log_parms <span class="op">=</span> <span class="cn">FALSE</span>,
- conf.level <span class="op">=</span> <span class="fl">0.6</span>,
- verbose <span class="op">=</span> <span class="cn">FALSE</span>,
- <span class="va">...</span>
-<span class="op">)</span>
-
-<span class="fu">saemix_data</span><span class="op">(</span><span class="va">object</span>, verbose <span class="op">=</span> <span class="cn">FALSE</span>, <span class="va">...</span><span class="op">)</span></code></pre></div>
+ <div class="sourceCode"><pre class="sourceCode r"><code><span><span class="fu">saem</span><span class="op">(</span><span class="va">object</span>, <span class="va">...</span><span class="op">)</span></span>
+<span></span>
+<span><span class="co"># S3 method for mmkin</span></span>
+<span><span class="fu">saem</span><span class="op">(</span></span>
+<span> <span class="va">object</span>,</span>
+<span> transformations <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"mkin"</span>, <span class="st">"saemix"</span><span class="op">)</span>,</span>
+<span> degparms_start <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/numeric.html" class="external-link">numeric</a></span><span class="op">(</span><span class="op">)</span>,</span>
+<span> test_log_parms <span class="op">=</span> <span class="cn">TRUE</span>,</span>
+<span> conf.level <span class="op">=</span> <span class="fl">0.6</span>,</span>
+<span> solution_type <span class="op">=</span> <span class="st">"auto"</span>,</span>
+<span> nbiter.saemix <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">300</span>, <span class="fl">100</span><span class="op">)</span>,</span>
+<span> control <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span>displayProgress <span class="op">=</span> <span class="cn">FALSE</span>, print <span class="op">=</span> <span class="cn">FALSE</span>, nbiter.saemix <span class="op">=</span> <span class="va">nbiter.saemix</span>,</span>
+<span> save <span class="op">=</span> <span class="cn">FALSE</span>, save.graphs <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span>,</span>
+<span> fail_with_errors <span class="op">=</span> <span class="cn">TRUE</span>,</span>
+<span> verbose <span class="op">=</span> <span class="cn">FALSE</span>,</span>
+<span> quiet <span class="op">=</span> <span class="cn">FALSE</span>,</span>
+<span> <span class="va">...</span></span>
+<span><span class="op">)</span></span>
+<span></span>
+<span><span class="co"># S3 method for saem.mmkin</span></span>
+<span><span class="fu"><a href="https://rdrr.io/r/base/print.html" class="external-link">print</a></span><span class="op">(</span><span class="va">x</span>, digits <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/Extremes.html" class="external-link">max</a></span><span class="op">(</span><span class="fl">3</span>, <span class="fu"><a href="https://rdrr.io/r/base/options.html" class="external-link">getOption</a></span><span class="op">(</span><span class="st">"digits"</span><span class="op">)</span> <span class="op">-</span> <span class="fl">3</span><span class="op">)</span>, <span class="va">...</span><span class="op">)</span></span>
+<span></span>
+<span><span class="fu">saemix_model</span><span class="op">(</span></span>
+<span> <span class="va">object</span>,</span>
+<span> solution_type <span class="op">=</span> <span class="st">"auto"</span>,</span>
+<span> transformations <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"mkin"</span>, <span class="st">"saemix"</span><span class="op">)</span>,</span>
+<span> degparms_start <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/numeric.html" class="external-link">numeric</a></span><span class="op">(</span><span class="op">)</span>,</span>
+<span> test_log_parms <span class="op">=</span> <span class="cn">FALSE</span>,</span>
+<span> conf.level <span class="op">=</span> <span class="fl">0.6</span>,</span>
+<span> verbose <span class="op">=</span> <span class="cn">FALSE</span>,</span>
+<span> <span class="va">...</span></span>
+<span><span class="op">)</span></span>
+<span></span>
+<span><span class="fu">saemix_data</span><span class="op">(</span><span class="va">object</span>, verbose <span class="op">=</span> <span class="cn">FALSE</span>, <span class="va">...</span><span class="op">)</span></span></code></pre></div>
</div>
<div id="arguments">
@@ -132,54 +132,87 @@ Expectation Maximisation algorithm (SAEM).</p>
<dl><dt>object</dt>
<dd><p>An <a href="mmkin.html">mmkin</a> row object containing several fits of the same
<a href="mkinmod.html">mkinmod</a> model to different datasets</p></dd>
+
+
<dt>...</dt>
<dd><p>Further parameters passed to <a href="https://rdrr.io/pkg/saemix/man/saemixModel.html" class="external-link">saemix::saemixModel</a>.</p></dd>
+
+
<dt>transformations</dt>
<dd><p>Per default, all parameter transformations are done
in mkin. If this argument is set to 'saemix', parameter transformations
-are done in 'saemix' for the supported cases. Currently this is only
-supported in cases where the initial concentration of the parent is not fixed,
-SFO or DFOP is used for the parent and there is either no metabolite or one.</p></dd>
+are done in 'saemix' for the supported cases, i.e. (as of version 1.1.2)
+SFO, FOMC, DFOP and HS without fixing <code>parent_0</code>, and SFO or DFOP with
+one SFO metabolite.</p></dd>
+
+
<dt>degparms_start</dt>
<dd><p>Parameter values given as a named numeric vector will
be used to override the starting values obtained from the 'mmkin' object.</p></dd>
+
+
<dt>test_log_parms</dt>
<dd><p>If TRUE, an attempt is made to use more robust starting
values for population parameters fitted as log parameters in mkin (like
rate constants) by only considering rate constants that pass the t-test
when calculating mean degradation parameters using <a href="mean_degparms.html">mean_degparms</a>.</p></dd>
+
+
<dt>conf.level</dt>
<dd><p>Possibility to adjust the required confidence level
for parameter that are tested if requested by 'test_log_parms'.</p></dd>
+
+
<dt>solution_type</dt>
<dd><p>Possibility to specify the solution type in case the
automatic choice is not desired</p></dd>
+
+
<dt>nbiter.saemix</dt>
<dd><p>Convenience option to increase the number of
iterations</p></dd>
+
+
<dt>control</dt>
<dd><p>Passed to <a href="https://rdrr.io/pkg/saemix/man/saemix.html" class="external-link">saemix::saemix</a>.</p></dd>
+
+
<dt>fail_with_errors</dt>
<dd><p>Should a failure to compute standard errors
from the inverse of the Fisher Information Matrix be a failure?</p></dd>
+
+
<dt>verbose</dt>
<dd><p>Should we print information about created objects of
type <a href="https://rdrr.io/pkg/saemix/man/SaemixModel-class.html" class="external-link">saemix::SaemixModel</a> and <a href="https://rdrr.io/pkg/saemix/man/SaemixData-class.html" class="external-link">saemix::SaemixData</a>?</p></dd>
+
+
<dt>quiet</dt>
<dd><p>Should we suppress the messages saemix prints at the beginning
and the end of the optimisation process?</p></dd>
+
+
<dt>x</dt>
<dd><p>An saem.mmkin object to print</p></dd>
+
+
<dt>digits</dt>
<dd><p>Number of digits to use for printing</p></dd>
+
</dl></div>
<div id="value">
<h2>Value</h2>
- <p>An S3 object of class 'saem.mmkin', containing the fitted
+
+
+<p>An S3 object of class 'saem.mmkin', containing the fitted
<a href="https://rdrr.io/pkg/saemix/man/SaemixObject-class.html" class="external-link">saemix::SaemixObject</a> as a list component named 'so'. The
-object also inherits from 'mixed.mmkin'.
-An <a href="https://rdrr.io/pkg/saemix/man/SaemixModel-class.html" class="external-link">saemix::SaemixModel</a> object.
-An <a href="https://rdrr.io/pkg/saemix/man/SaemixData-class.html" class="external-link">saemix::SaemixData</a> object.</p>
+object also inherits from 'mixed.mmkin'.</p>
+
+
+<p>An <a href="https://rdrr.io/pkg/saemix/man/SaemixModel-class.html" class="external-link">saemix::SaemixModel</a> object.</p>
+
+
+<p>An <a href="https://rdrr.io/pkg/saemix/man/SaemixData-class.html" class="external-link">saemix::SaemixData</a> object.</p>
</div>
<div id="details">
<h2>Details</h2>
@@ -196,77 +229,105 @@ using <a href="mmkin.html">mmkin</a>.</p>
<div id="ref-examples">
<h2>Examples</h2>
- <div class="sourceCode"><pre class="sourceCode r"><code><span class="r-in"><span class="co"># \dontrun{</span></span>
-<span class="r-in"><span class="va">ds</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/lapply.html" class="external-link">lapply</a></span><span class="op">(</span><span class="va">experimental_data_for_UBA_2019</span><span class="op">[</span><span class="fl">6</span><span class="op">:</span><span class="fl">10</span><span class="op">]</span>,</span>
-<span class="r-in"> <span class="kw">function</span><span class="op">(</span><span class="va">x</span><span class="op">)</span> <span class="fu"><a href="https://rdrr.io/r/base/subset.html" class="external-link">subset</a></span><span class="op">(</span><span class="va">x</span><span class="op">$</span><span class="va">data</span><span class="op">[</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"name"</span>, <span class="st">"time"</span>, <span class="st">"value"</span><span class="op">)</span><span class="op">]</span><span class="op">)</span><span class="op">)</span></span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/r/base/names.html" class="external-link">names</a></span><span class="op">(</span><span class="va">ds</span><span class="op">)</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/paste.html" class="external-link">paste</a></span><span class="op">(</span><span class="st">"Dataset"</span>, <span class="fl">6</span><span class="op">:</span><span class="fl">10</span><span class="op">)</span></span>
-<span class="r-in"><span class="va">f_mmkin_parent_p0_fixed</span> <span class="op">&lt;-</span> <span class="fu"><a href="mmkin.html">mmkin</a></span><span class="op">(</span><span class="st">"FOMC"</span>, <span class="va">ds</span>,</span>
-<span class="r-in"> state.ini <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fl">100</span><span class="op">)</span>, fixed_initials <span class="op">=</span> <span class="st">"parent"</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span>
-<span class="r-in"><span class="va">f_saem_p0_fixed</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent_p0_fixed</span><span class="op">)</span></span>
-<span class="r-in"></span>
-<span class="r-in"><span class="va">f_mmkin_parent</span> <span class="op">&lt;-</span> <span class="fu"><a href="mmkin.html">mmkin</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"SFO"</span>, <span class="st">"FOMC"</span>, <span class="st">"DFOP"</span><span class="op">)</span>, <span class="va">ds</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span>
-<span class="r-in"><span class="va">f_saem_sfo</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent</span><span class="op">[</span><span class="st">"SFO"</span>, <span class="op">]</span><span class="op">)</span></span>
-<span class="r-in"><span class="va">f_saem_fomc</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent</span><span class="op">[</span><span class="st">"FOMC"</span>, <span class="op">]</span><span class="op">)</span></span>
-<span class="r-in"><span class="va">f_saem_dfop</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent</span><span class="op">[</span><span class="st">"DFOP"</span>, <span class="op">]</span><span class="op">)</span></span>
-<span class="r-in"></span>
-<span class="r-in"><span class="co"># The returned saem.mmkin object contains an SaemixObject, therefore we can use</span></span>
-<span class="r-in"><span class="co"># functions from saemix</span></span>
-<span class="r-in"><span class="kw"><a href="https://rdrr.io/r/base/library.html" class="external-link">library</a></span><span class="op">(</span><span class="va">saemix</span><span class="op">)</span></span>
+ <div class="sourceCode"><pre class="sourceCode r"><code><span class="r-in"><span><span class="co"># \dontrun{</span></span></span>
+<span class="r-in"><span><span class="va">ds</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/lapply.html" class="external-link">lapply</a></span><span class="op">(</span><span class="va">experimental_data_for_UBA_2019</span><span class="op">[</span><span class="fl">6</span><span class="op">:</span><span class="fl">10</span><span class="op">]</span>,</span></span>
+<span class="r-in"><span> <span class="kw">function</span><span class="op">(</span><span class="va">x</span><span class="op">)</span> <span class="fu"><a href="https://rdrr.io/r/base/subset.html" class="external-link">subset</a></span><span class="op">(</span><span class="va">x</span><span class="op">$</span><span class="va">data</span><span class="op">[</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"name"</span>, <span class="st">"time"</span>, <span class="st">"value"</span><span class="op">)</span><span class="op">]</span><span class="op">)</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/names.html" class="external-link">names</a></span><span class="op">(</span><span class="va">ds</span><span class="op">)</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/paste.html" class="external-link">paste</a></span><span class="op">(</span><span class="st">"Dataset"</span>, <span class="fl">6</span><span class="op">:</span><span class="fl">10</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="va">f_mmkin_parent_p0_fixed</span> <span class="op">&lt;-</span> <span class="fu"><a href="mmkin.html">mmkin</a></span><span class="op">(</span><span class="st">"FOMC"</span>, <span class="va">ds</span>,</span></span>
+<span class="r-in"><span> state.ini <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fl">100</span><span class="op">)</span>, fixed_initials <span class="op">=</span> <span class="st">"parent"</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="va">f_saem_p0_fixed</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent_p0_fixed</span><span class="op">)</span></span></span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="va">f_mmkin_parent</span> <span class="op">&lt;-</span> <span class="fu"><a href="mmkin.html">mmkin</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"SFO"</span>, <span class="st">"FOMC"</span>, <span class="st">"DFOP"</span><span class="op">)</span>, <span class="va">ds</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="va">f_saem_sfo</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent</span><span class="op">[</span><span class="st">"SFO"</span>, <span class="op">]</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="va">f_saem_fomc</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent</span><span class="op">[</span><span class="st">"FOMC"</span>, <span class="op">]</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="va">f_saem_dfop</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent</span><span class="op">[</span><span class="st">"DFOP"</span>, <span class="op">]</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="fu"><a href="illparms.html">illparms</a></span><span class="op">(</span><span class="va">f_saem_dfop</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [1] "sd(g_qlogis)"</span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/update.html" class="external-link">update</a></span><span class="op">(</span><span class="va">f_saem_dfop</span>, covariance.model <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/diag.html" class="external-link">diag</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">1</span>, <span class="fl">1</span>, <span class="fl">1</span>, <span class="fl">0</span><span class="op">)</span><span class="op">)</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Kinetic nonlinear mixed-effects model fit by SAEM</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Structural model:</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time)))</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> * parent</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Data:</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 90 observations of 1 variable(s) grouped in 5 datasets</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Likelihood computed by importance sampling</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> AIC BIC logLik</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 490.6 487.5 -237.3</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Fitted parameters:</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> estimate lower upper</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> parent_0 93.902 91.3695 96.4339</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 -2.936 -3.9950 -1.8762</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 -3.091 -4.9290 -1.2523</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis -0.366 -0.6484 -0.0836</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 2.385 2.0033 2.7664</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.parent_0 2.476 0.3890 4.5623</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 1.195 0.4381 1.9517</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 2.092 0.7906 3.3932</span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/AIC.html" class="external-link">AIC</a></span><span class="op">(</span><span class="va">f_saem_dfop</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [1] 493.9811</span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="co"># The returned saem.mmkin object contains an SaemixObject, therefore we can use</span></span></span>
+<span class="r-in"><span><span class="co"># functions from saemix</span></span></span>
+<span class="r-in"><span><span class="kw"><a href="https://rdrr.io/r/base/library.html" class="external-link">library</a></span><span class="op">(</span><span class="va">saemix</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Loading required package: npde</span>
-<span class="r-msg co"><span class="r-pr">#&gt;</span> Package saemix, version 3.0</span>
+<span class="r-msg co"><span class="r-pr">#&gt;</span> Package saemix, version 3.1</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> please direct bugs, questions and feedback to emmanuelle.comets@inserm.fr</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> </span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Attaching package: ‘saemix’</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> The following objects are masked from ‘package:npde’:</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> </span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> kurtosis, skewness</span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/compare.saemix.html" class="external-link">compare.saemix</a></span><span class="op">(</span><span class="va">f_saem_sfo</span><span class="op">$</span><span class="va">so</span>, <span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, <span class="va">f_saem_dfop</span><span class="op">$</span><span class="va">so</span><span class="op">)</span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/compare.saemix.html" class="external-link">compare.saemix</a></span><span class="op">(</span><span class="va">f_saem_sfo</span><span class="op">$</span><span class="va">so</span>, <span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, <span class="va">f_saem_dfop</span><span class="op">$</span><span class="va">so</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Likelihoods calculated by importance sampling</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> AIC BIC</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> 1 624.2598 622.3070</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> 2 467.8664 465.1324</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> 3 493.9811 490.4660</span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/plot-SaemixObject-ANY-method.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, plot.type <span class="op">=</span> <span class="st">"convergence"</span><span class="op">)</span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, plot.type <span class="op">=</span> <span class="st">"convergence"</span><span class="op">)</span></span></span>
<span class="r-plt img"><img src="saem-1.png" alt="" width="700" height="433"></span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/plot-SaemixObject-ANY-method.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, plot.type <span class="op">=</span> <span class="st">"individual.fit"</span><span class="op">)</span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, plot.type <span class="op">=</span> <span class="st">"individual.fit"</span><span class="op">)</span></span></span>
<span class="r-plt img"><img src="saem-2.png" alt="" width="700" height="433"></span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/plot-SaemixObject-ANY-method.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, plot.type <span class="op">=</span> <span class="st">"npde"</span><span class="op">)</span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, plot.type <span class="op">=</span> <span class="st">"npde"</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Simulating data using nsim = 1000 simulated datasets</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Computing WRES and npde .</span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Please use npdeSaemix to obtain VPC and npde</span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/plot-SaemixObject-ANY-method.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, plot.type <span class="op">=</span> <span class="st">"vpc"</span><span class="op">)</span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, plot.type <span class="op">=</span> <span class="st">"vpc"</span><span class="op">)</span></span></span>
<span class="r-plt img"><img src="saem-3.png" alt="" width="700" height="433"></span>
-<span class="r-in"></span>
-<span class="r-in"><span class="va">f_mmkin_parent_tc</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/stats/update.html" class="external-link">update</a></span><span class="op">(</span><span class="va">f_mmkin_parent</span>, error_model <span class="op">=</span> <span class="st">"tc"</span><span class="op">)</span></span>
-<span class="r-in"><span class="va">f_saem_fomc_tc</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent_tc</span><span class="op">[</span><span class="st">"FOMC"</span>, <span class="op">]</span><span class="op">)</span></span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/compare.saemix.html" class="external-link">compare.saemix</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, <span class="va">f_saem_fomc_tc</span><span class="op">$</span><span class="va">so</span><span class="op">)</span></span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="va">f_mmkin_parent_tc</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/stats/update.html" class="external-link">update</a></span><span class="op">(</span><span class="va">f_mmkin_parent</span>, error_model <span class="op">=</span> <span class="st">"tc"</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="va">f_saem_fomc_tc</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin_parent_tc</span><span class="op">[</span><span class="st">"FOMC"</span>, <span class="op">]</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/compare.saemix.html" class="external-link">compare.saemix</a></span><span class="op">(</span><span class="va">f_saem_fomc</span><span class="op">$</span><span class="va">so</span>, <span class="va">f_saem_fomc_tc</span><span class="op">$</span><span class="va">so</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Likelihoods calculated by importance sampling</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> AIC BIC</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> 1 467.8664 465.1324</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> 2 469.9096 466.7851</span>
-<span class="r-in"></span>
-<span class="r-in"><span class="va">sfo_sfo</span> <span class="op">&lt;-</span> <span class="fu"><a href="mkinmod.html">mkinmod</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"SFO"</span>, <span class="st">"A1"</span><span class="op">)</span>,</span>
-<span class="r-in"> A1 <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"SFO"</span><span class="op">)</span><span class="op">)</span></span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="va">sfo_sfo</span> <span class="op">&lt;-</span> <span class="fu"><a href="mkinmod.html">mkinmod</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"SFO"</span>, <span class="st">"A1"</span><span class="op">)</span>,</span></span>
+<span class="r-in"><span> A1 <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"SFO"</span><span class="op">)</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Temporary DLL for differentials generated and loaded</span>
-<span class="r-in"><span class="va">fomc_sfo</span> <span class="op">&lt;-</span> <span class="fu"><a href="mkinmod.html">mkinmod</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"FOMC"</span>, <span class="st">"A1"</span><span class="op">)</span>,</span>
-<span class="r-in"> A1 <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"SFO"</span><span class="op">)</span><span class="op">)</span></span>
+<span class="r-in"><span><span class="va">fomc_sfo</span> <span class="op">&lt;-</span> <span class="fu"><a href="mkinmod.html">mkinmod</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"FOMC"</span>, <span class="st">"A1"</span><span class="op">)</span>,</span></span>
+<span class="r-in"><span> A1 <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"SFO"</span><span class="op">)</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Temporary DLL for differentials generated and loaded</span>
-<span class="r-in"><span class="va">dfop_sfo</span> <span class="op">&lt;-</span> <span class="fu"><a href="mkinmod.html">mkinmod</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"DFOP"</span>, <span class="st">"A1"</span><span class="op">)</span>,</span>
-<span class="r-in"> A1 <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"SFO"</span><span class="op">)</span><span class="op">)</span></span>
+<span class="r-in"><span><span class="va">dfop_sfo</span> <span class="op">&lt;-</span> <span class="fu"><a href="mkinmod.html">mkinmod</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"DFOP"</span>, <span class="st">"A1"</span><span class="op">)</span>,</span></span>
+<span class="r-in"><span> A1 <span class="op">=</span> <span class="fu"><a href="mkinmod.html">mkinsub</a></span><span class="op">(</span><span class="st">"SFO"</span><span class="op">)</span><span class="op">)</span></span></span>
<span class="r-msg co"><span class="r-pr">#&gt;</span> Temporary DLL for differentials generated and loaded</span>
-<span class="r-in"><span class="co"># The following fit uses analytical solutions for SFO-SFO and DFOP-SFO,</span></span>
-<span class="r-in"><span class="co"># and compiled ODEs for FOMC that are much slower</span></span>
-<span class="r-in"><span class="va">f_mmkin</span> <span class="op">&lt;-</span> <span class="fu"><a href="mmkin.html">mmkin</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span></span>
-<span class="r-in"> <span class="st">"SFO-SFO"</span> <span class="op">=</span> <span class="va">sfo_sfo</span>, <span class="st">"FOMC-SFO"</span> <span class="op">=</span> <span class="va">fomc_sfo</span>, <span class="st">"DFOP-SFO"</span> <span class="op">=</span> <span class="va">dfop_sfo</span><span class="op">)</span>,</span>
-<span class="r-in"> <span class="va">ds</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span>
-<span class="r-in"><span class="co"># saem fits of SFO-SFO and DFOP-SFO to these data take about five seconds</span></span>
-<span class="r-in"><span class="co"># each on this system, as we use analytical solutions written for saemix.</span></span>
-<span class="r-in"><span class="co"># When using the analytical solutions written for mkin this took around</span></span>
-<span class="r-in"><span class="co"># four minutes</span></span>
-<span class="r-in"><span class="va">f_saem_sfo_sfo</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin</span><span class="op">[</span><span class="st">"SFO-SFO"</span>, <span class="op">]</span><span class="op">)</span></span>
-<span class="r-in"><span class="va">f_saem_dfop_sfo</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin</span><span class="op">[</span><span class="st">"DFOP-SFO"</span>, <span class="op">]</span><span class="op">)</span></span>
-<span class="r-in"><span class="co"># We can use print, plot and summary methods to check the results</span></span>
-<span class="r-in"><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">f_saem_dfop_sfo</span><span class="op">)</span></span>
+<span class="r-in"><span><span class="co"># The following fit uses analytical solutions for SFO-SFO and DFOP-SFO,</span></span></span>
+<span class="r-in"><span><span class="co"># and compiled ODEs for FOMC that are much slower</span></span></span>
+<span class="r-in"><span><span class="va">f_mmkin</span> <span class="op">&lt;-</span> <span class="fu"><a href="mmkin.html">mmkin</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/list.html" class="external-link">list</a></span><span class="op">(</span></span></span>
+<span class="r-in"><span> <span class="st">"SFO-SFO"</span> <span class="op">=</span> <span class="va">sfo_sfo</span>, <span class="st">"FOMC-SFO"</span> <span class="op">=</span> <span class="va">fomc_sfo</span>, <span class="st">"DFOP-SFO"</span> <span class="op">=</span> <span class="va">dfop_sfo</span><span class="op">)</span>,</span></span>
+<span class="r-in"><span> <span class="va">ds</span>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="co"># saem fits of SFO-SFO and DFOP-SFO to these data take about five seconds</span></span></span>
+<span class="r-in"><span><span class="co"># each on this system, as we use analytical solutions written for saemix.</span></span></span>
+<span class="r-in"><span><span class="co"># When using the analytical solutions written for mkin this took around</span></span></span>
+<span class="r-in"><span><span class="co"># four minutes</span></span></span>
+<span class="r-in"><span><span class="va">f_saem_sfo_sfo</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin</span><span class="op">[</span><span class="st">"SFO-SFO"</span>, <span class="op">]</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="va">f_saem_dfop_sfo</span> <span class="op">&lt;-</span> <span class="fu">saem</span><span class="op">(</span><span class="va">f_mmkin</span><span class="op">[</span><span class="st">"DFOP-SFO"</span>, <span class="op">]</span><span class="op">)</span></span></span>
+<span class="r-in"><span><span class="co"># We can use print, plot and summary methods to check the results</span></span></span>
+<span class="r-in"><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">f_saem_dfop_sfo</span><span class="op">)</span></span></span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Kinetic nonlinear mixed-effects model fit by SAEM</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Structural model:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *</span>
@@ -291,21 +352,21 @@ using <a href="mmkin.html">mmkin</a>.</p>
<span class="r-out co"><span class="r-pr">#&gt;</span> log_k1 -2.5841 -3.6876 -1.4805</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> log_k2 -3.5228 -5.3254 -1.7203</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> g_qlogis -0.1027 -0.8719 0.6665</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 1.8856 1.6770 2.0943</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.parent_0 2.7682 0.7669 4.7695</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 1.8856 1.6676 2.1037</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.parent_0 2.7682 0.7668 4.7695</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k_A1 1.7447 0.4047 3.0848</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> SD.f_parent_qlogis 0.4525 0.1620 0.7431</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 1.2423 0.4560 2.0285</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 2.0390 0.7601 3.3180</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.g_qlogis 0.4439 -0.3068 1.1946</span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/plot-SaemixObject-ANY-method.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span><span class="op">)</span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.g_qlogis 0.4439 -0.3069 1.1947</span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/plot.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span><span class="op">)</span></span></span>
<span class="r-plt img"><img src="saem-4.png" alt="" width="700" height="433"></span>
-<span class="r-in"><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/summary-methods.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>, data <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> saemix version used for fitting: 3.0 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> mkin version used for pre-fitting: 1.1.0 </span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/pkg/saemix/man/summary-methods.html" class="external-link">summary</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>, data <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> saemix version used for fitting: 3.1 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> mkin version used for pre-fitting: 1.1.2 </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> R version used for fitting: 4.2.1 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Date of fit: Thu Jun 30 10:23:53 2022 </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Date of summary: Thu Jun 30 10:23:53 2022 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Date of fit: Wed Aug 10 13:15:20 2022 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Date of summary: Wed Aug 10 13:15:20 2022 </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Equations:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *</span>
@@ -320,7 +381,7 @@ using <a href="mmkin.html">mmkin</a>.</p>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Model predictions using solution type analytical </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Fitted in 8.85 s</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Fitted in 9.258 s</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Using 300, 100 iterations and 10 chains</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Variance model: Constant variance </span>
@@ -359,16 +420,16 @@ using <a href="mmkin.html">mmkin</a>.</p>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Random effects:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.parent_0 2.7682 0.7669 4.7695</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.parent_0 2.7682 0.7668 4.7695</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k_A1 1.7447 0.4047 3.0848</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> SD.f_parent_qlogis 0.4525 0.1620 0.7431</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k1 1.2423 0.4560 2.0285</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> SD.log_k2 2.0390 0.7601 3.3180</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> SD.g_qlogis 0.4439 -0.3068 1.1946</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> SD.g_qlogis 0.4439 -0.3069 1.1947</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Variance model:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 1.886 1.677 2.094</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> a.1 1.886 1.668 2.104</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> </span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Backtransformed parameters:</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> est. lower upper</span>
@@ -554,27 +615,27 @@ using <a href="mmkin.html">mmkin</a>.</p>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 21 6.4 7.59693 -1.19693 1.886 -0.634761</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 21 7.2 7.59693 -0.39693 1.886 -0.210502</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 41 9.1 9.86436 -0.76436 1.886 -0.405361</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 41 8.5 9.86436 -1.36436 1.886 -0.723556</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 41 8.5 9.86436 -1.36436 1.886 -0.723555</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 63 11.7 10.99397 0.70603 1.886 0.374425</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 63 12.0 10.99397 1.00603 1.886 0.533522</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 91 13.3 11.91274 1.38726 1.886 0.735696</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 91 13.2 11.91274 1.28726 1.886 0.682663</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 120 14.3 12.66519 1.63481 1.886 0.866981</span>
<span class="r-out co"><span class="r-pr">#&gt;</span> Dataset 10 A1 120 12.1 12.66519 -0.56519 1.886 -0.299733</span>
-<span class="r-in"></span>
-<span class="r-in"><span class="co"># The following takes about 6 minutes</span></span>
-<span class="r-in"><span class="co">#f_saem_dfop_sfo_deSolve &lt;- saem(f_mmkin["DFOP-SFO", ], solution_type = "deSolve",</span></span>
-<span class="r-in"><span class="co"># control = list(nbiter.saemix = c(200, 80), nbdisplay = 10))</span></span>
-<span class="r-in"></span>
-<span class="r-in"><span class="co">#saemix::compare.saemix(list(</span></span>
-<span class="r-in"><span class="co"># f_saem_dfop_sfo$so,</span></span>
-<span class="r-in"><span class="co"># f_saem_dfop_sfo_deSolve$so))</span></span>
-<span class="r-in"></span>
-<span class="r-in"><span class="co"># If the model supports it, we can also use eigenvalue based solutions, which</span></span>
-<span class="r-in"><span class="co"># take a similar amount of time</span></span>
-<span class="r-in"><span class="co">#f_saem_sfo_sfo_eigen &lt;- saem(f_mmkin["SFO-SFO", ], solution_type = "eigen",</span></span>
-<span class="r-in"><span class="co"># control = list(nbiter.saemix = c(200, 80), nbdisplay = 10))</span></span>
-<span class="r-in"><span class="co"># }</span></span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="co"># The following takes about 6 minutes</span></span></span>
+<span class="r-in"><span><span class="co">#f_saem_dfop_sfo_deSolve &lt;- saem(f_mmkin["DFOP-SFO", ], solution_type = "deSolve",</span></span></span>
+<span class="r-in"><span><span class="co"># control = list(nbiter.saemix = c(200, 80), nbdisplay = 10))</span></span></span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="co">#saemix::compare.saemix(list(</span></span></span>
+<span class="r-in"><span><span class="co"># f_saem_dfop_sfo$so,</span></span></span>
+<span class="r-in"><span><span class="co"># f_saem_dfop_sfo_deSolve$so))</span></span></span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="co"># If the model supports it, we can also use eigenvalue based solutions, which</span></span></span>
+<span class="r-in"><span><span class="co"># take a similar amount of time</span></span></span>
+<span class="r-in"><span><span class="co">#f_saem_sfo_sfo_eigen &lt;- saem(f_mmkin["SFO-SFO", ], solution_type = "eigen",</span></span></span>
+<span class="r-in"><span><span class="co"># control = list(nbiter.saemix = c(200, 80), nbdisplay = 10))</span></span></span>
+<span class="r-in"><span><span class="co"># }</span></span></span>
</code></pre></div>
</div>
</div>
@@ -589,7 +650,7 @@ using <a href="mmkin.html">mmkin</a>.</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.3.</p>
+ <p></p><p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.0.6.</p>
</div>
</footer></div>

Contact - Imprint