diff options
Diffstat (limited to 'docs/dev/reference/summary.saem.mmkin.html')
| -rw-r--r-- | docs/dev/reference/summary.saem.mmkin.html | 677 | 
1 files changed, 0 insertions, 677 deletions
| diff --git a/docs/dev/reference/summary.saem.mmkin.html b/docs/dev/reference/summary.saem.mmkin.html deleted file mode 100644 index e434ad8d..00000000 --- a/docs/dev/reference/summary.saem.mmkin.html +++ /dev/null @@ -1,677 +0,0 @@ -<!DOCTYPE html> -<!-- Generated by pkgdown: do not edit by hand --><html lang="en"><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><meta charset="utf-8"><meta http-equiv="X-UA-Compatible" content="IE=edge"><meta name="viewport" content="width=device-width, initial-scale=1.0"><title>Summary method for class "saem.mmkin" — summary.saem.mmkin • mkin</title><!-- jquery --><script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.4.1/jquery.min.js" integrity="sha256-CSXorXvZcTkaix6Yvo6HppcZGetbYMGWSFlBw8HfCJo=" crossorigin="anonymous"></script><!-- Bootstrap --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/css/bootstrap.min.css" integrity="sha256-bZLfwXAP04zRMK2BjiO8iu9pf4FbLqX6zitd+tIvLhE=" crossorigin="anonymous"><script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/js/bootstrap.min.js" integrity="sha256-nuL8/2cJ5NDSSwnKD8VqreErSWHtnEP9E7AySL+1ev4=" crossorigin="anonymous"></script><!-- bootstrap-toc --><link rel="stylesheet" href="../bootstrap-toc.css"><script src="../bootstrap-toc.js"></script><!-- Font Awesome icons --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/all.min.css" integrity="sha256-mmgLkCYLUQbXn0B1SRqzHar6dCnv9oZFPEC1g1cwlkk=" crossorigin="anonymous"><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous"><!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- pkgdown --><link href="../pkgdown.css" rel="stylesheet"><script src="../pkgdown.js"></script><meta property="og:title" content="Summary method for class " saem.mmkin summary.saem.mmkin><meta property="og:description" content="Lists model equations, initial parameter values, optimised parameters -for fixed effects (population), random effects (deviations from the -population mean) and residual error model, as well as the resulting -endpoints such as formation fractions and DT50 values. Optionally -(default is FALSE), the data are listed in full."><meta name="robots" content="noindex"><!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]> -<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script> -<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script> -<![endif]--></head><body data-spy="scroll" data-target="#toc"> -     - -    <div class="container template-reference-topic"> -      <header><div class="navbar navbar-default navbar-fixed-top" role="navigation"> -  <div class="container"> -    <div class="navbar-header"> -      <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar" aria-expanded="false"> -        <span class="sr-only">Toggle navigation</span> -        <span class="icon-bar"></span> -        <span class="icon-bar"></span> -        <span class="icon-bar"></span> -      </button> -      <span class="navbar-brand"> -        <a class="navbar-link" href="../index.html">mkin</a> -        <span class="version label label-info" data-toggle="tooltip" data-placement="bottom" title="In-development version">1.2.3</span> -      </span> -    </div> - -    <div id="navbar" class="navbar-collapse collapse"> -      <ul class="nav navbar-nav"><li> -  <a href="../reference/index.html">Reference</a> -</li> -<li class="dropdown"> -  <a href="#" class="dropdown-toggle" data-toggle="dropdown" role="button" data-bs-toggle="dropdown" aria-expanded="false"> -    Articles -      -    <span class="caret"></span> -  </a> -  <ul class="dropdown-menu" role="menu"><li> -      <a href="../articles/mkin.html">Introduction to mkin</a> -    </li> -    <li class="divider"> -    <li class="dropdown-header">Example evaluations with (generalised) nonlinear least squares</li> -    <li> -      <a href="../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a> -    </li> -    <li> -      <a href="../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a> -    </li> -    <li> -      <a href="../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a> -    </li> -    <li class="divider"> -    <li class="dropdown-header">Example evaluations with hierarchical models (nonlinear mixed-effects models)</li> -    <li> -      <a href="../articles/prebuilt/2022_dmta_parent.html">Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P</a> -    </li> -    <li> -      <a href="../articles/prebuilt/2022_dmta_pathway.html">Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P</a> -    </li> -    <li> -      <a href="../articles/prebuilt/2022_cyan_pathway.html">Testing hierarchical pathway kinetics with residue data on cyantraniliprole</a> -    </li> -    <li> -      <a href="../articles/web_only/dimethenamid_2018.html">Comparison of saemix and nlme evaluations of dimethenamid data from 2018</a> -    </li> -    <li> -      <a href="../articles/web_only/multistart.html">Short demo of the multistart method</a> -    </li> -    <li class="divider"> -    <li class="dropdown-header">Performance</li> -    <li> -      <a href="../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a> -    </li> -    <li> -      <a href="../articles/web_only/benchmarks.html">Benchmark timings for mkin</a> -    </li> -    <li> -      <a href="../articles/web_only/saem_benchmarks.html">Benchmark timings for saem.mmkin</a> -    </li> -    <li class="divider"> -    <li class="dropdown-header">Miscellaneous</li> -    <li> -      <a href="../articles/twa.html">Calculation of time weighted average concentrations with mkin</a> -    </li> -    <li> -      <a href="../articles/web_only/NAFTA_examples.html">Example evaluation of NAFTA SOP Attachment examples</a> -    </li> -  </ul></li> -<li> -  <a href="../news/index.html">News</a> -</li> -      </ul><ul class="nav navbar-nav navbar-right"><li> -  <a href="https://github.com/jranke/mkin/" class="external-link"> -    <span class="fab fa-github fa-lg"></span> -      -  </a> -</li> -      </ul></div><!--/.nav-collapse --> -  </div><!--/.container --> -</div><!--/.navbar --> - -       - -      </header><div class="row"> -  <div class="col-md-9 contents"> -    <div class="page-header"> -    <h1>Summary method for class "saem.mmkin"</h1> -    <small class="dont-index">Source: <a href="https://github.com/jranke/mkin/blob/HEAD/R/summary.saem.mmkin.R" class="external-link"><code>R/summary.saem.mmkin.R</code></a></small> -    <div class="hidden name"><code>summary.saem.mmkin.Rd</code></div> -    </div> - -    <div class="ref-description"> -    <p>Lists model equations, initial parameter values, optimised parameters -for fixed effects (population), random effects (deviations from the -population mean) and residual error model, as well as the resulting -endpoints such as formation fractions and DT50 values. Optionally -(default is FALSE), the data are listed in full.</p> -    </div> - -    <div id="ref-usage"> -    <div class="sourceCode"><pre class="sourceCode r"><code><span><span class="co"># S3 method for saem.mmkin</span></span> -<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> -<span>  <span class="va">object</span>,</span> -<span>  data <span class="op">=</span> <span class="cn">FALSE</span>,</span> -<span>  verbose <span class="op">=</span> <span class="cn">FALSE</span>,</span> -<span>  covariates <span class="op">=</span> <span class="cn">NULL</span>,</span> -<span>  covariate_quantile <span class="op">=</span> <span class="fl">0.5</span>,</span> -<span>  distimes <span class="op">=</span> <span class="cn">TRUE</span>,</span> -<span>  <span class="va">...</span></span> -<span><span class="op">)</span></span> -<span></span> -<span><span class="co"># S3 method for summary.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>, verbose <span class="op">=</span> <span class="va">x</span><span class="op">$</span><span class="va">verbose</span>, <span class="va">...</span><span class="op">)</span></span></code></pre></div> -    </div> - -    <div id="arguments"> -    <h2>Arguments</h2> -    <dl><dt>object</dt> -<dd><p>an object of class <a href="saem.html">saem.mmkin</a></p></dd> - - -<dt>data</dt> -<dd><p>logical, indicating whether the full data should be included in -the summary.</p></dd> - - -<dt>verbose</dt> -<dd><p>Should the summary be verbose?</p></dd> - - -<dt>covariates</dt> -<dd><p>Numeric vector with covariate values for all variables in -any covariate models in the object. If given, it overrides 'covariate_quantile'.</p></dd> - - -<dt>covariate_quantile</dt> -<dd><p>This argument only has an effect if the fitted -object has covariate models. If so, the default is to show endpoints -for the median of the covariate values (50th percentile).</p></dd> - - -<dt>distimes</dt> -<dd><p>logical, indicating whether DT50 and DT90 values should be -included.</p></dd> - - -<dt>...</dt> -<dd><p>optional arguments passed to methods like <code>print</code>.</p></dd> - - -<dt>x</dt> -<dd><p>an object of class summary.saem.mmkin</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>The summary function returns a list based on the <a href="https://rdrr.io/pkg/saemix/man/SaemixObject-class.html" class="external-link">saemix::SaemixObject</a></p> - - -<p>obtained in the fit, with at least the following additional components</p> -<dl><dt>saemixversion, mkinversion, Rversion</dt> -<dd><p>The saemix, mkin and R versions used</p></dd> - -<dt>date.fit, date.summary</dt> -<dd><p>The dates where the fit and the summary were -produced</p></dd> - -<dt>diffs</dt> -<dd><p>The differential equations used in the degradation model</p></dd> - -<dt>use_of_ff</dt> -<dd><p>Was maximum or minimum use made of formation fractions</p></dd> - -<dt>data</dt> -<dd><p>The data</p></dd> - -<dt>confint_trans</dt> -<dd><p>Transformed parameters as used in the optimisation, with confidence intervals</p></dd> - -<dt>confint_back</dt> -<dd><p>Backtransformed parameters, with confidence intervals if available</p></dd> - -<dt>confint_errmod</dt> -<dd><p>Error model parameters with confidence intervals</p></dd> - -<dt>ff</dt> -<dd><p>The estimated formation fractions derived from the fitted -model.</p></dd> - -<dt>distimes</dt> -<dd><p>The DT50 and DT90 values for each observed variable.</p></dd> - -<dt>SFORB</dt> -<dd><p>If applicable, eigenvalues of SFORB components of the model.</p></dd> - -</dl><p>The print method is called for its side effect, i.e. printing the summary.</p> -    </div> -    <div id="author"> -    <h2>Author</h2> -    <p>Johannes Ranke for the mkin specific parts -saemix authors for the parts inherited from saemix.</p> -    </div> - -    <div id="ref-examples"> -    <h2>Examples</h2> -    <div class="sourceCode"><pre class="sourceCode r"><code><span class="r-in"><span><span class="co"># Generate five datasets following DFOP-SFO kinetics</span></span></span> -<span class="r-in"><span><span class="va">sampling_times</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">0</span>, <span class="fl">1</span>, <span class="fl">3</span>, <span class="fl">7</span>, <span class="fl">14</span>, <span class="fl">28</span>, <span class="fl">60</span>, <span class="fl">90</span>, <span class="fl">120</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="va">dfop_sfo</span> <span class="op"><-</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">"m1"</span><span class="op">)</span>,</span></span> -<span class="r-in"><span> m1 <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>, quiet <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/Random.html" class="external-link">set.seed</a></span><span class="op">(</span><span class="fl">1234</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="va">k1_in</span> <span class="op"><-</span> <span class="fu"><a href="https://rdrr.io/r/stats/Lognormal.html" class="external-link">rlnorm</a></span><span class="op">(</span><span class="fl">5</span>, <span class="fu"><a href="https://rdrr.io/r/base/Log.html" class="external-link">log</a></span><span class="op">(</span><span class="fl">0.1</span><span class="op">)</span>, <span class="fl">0.3</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="va">k2_in</span> <span class="op"><-</span> <span class="fu"><a href="https://rdrr.io/r/stats/Lognormal.html" class="external-link">rlnorm</a></span><span class="op">(</span><span class="fl">5</span>, <span class="fu"><a href="https://rdrr.io/r/base/Log.html" class="external-link">log</a></span><span class="op">(</span><span class="fl">0.02</span><span class="op">)</span>, <span class="fl">0.3</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="va">g_in</span> <span class="op"><-</span> <span class="fu"><a href="https://rdrr.io/r/stats/Logistic.html" class="external-link">plogis</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/stats/Normal.html" class="external-link">rnorm</a></span><span class="op">(</span><span class="fl">5</span>, <span class="fu"><a href="https://rdrr.io/r/stats/Logistic.html" class="external-link">qlogis</a></span><span class="op">(</span><span class="fl">0.5</span><span class="op">)</span>, <span class="fl">0.3</span><span class="op">)</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="va">f_parent_to_m1_in</span> <span class="op"><-</span> <span class="fu"><a href="https://rdrr.io/r/stats/Logistic.html" class="external-link">plogis</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/stats/Normal.html" class="external-link">rnorm</a></span><span class="op">(</span><span class="fl">5</span>, <span class="fu"><a href="https://rdrr.io/r/stats/Logistic.html" class="external-link">qlogis</a></span><span class="op">(</span><span class="fl">0.3</span><span class="op">)</span>, <span class="fl">0.3</span><span class="op">)</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="va">k_m1_in</span> <span class="op"><-</span> <span class="fu"><a href="https://rdrr.io/r/stats/Lognormal.html" class="external-link">rlnorm</a></span><span class="op">(</span><span class="fl">5</span>, <span class="fu"><a href="https://rdrr.io/r/base/Log.html" class="external-link">log</a></span><span class="op">(</span><span class="fl">0.02</span><span class="op">)</span>, <span class="fl">0.3</span><span class="op">)</span></span></span> -<span class="r-in"><span></span></span> -<span class="r-in"><span><span class="va">pred_dfop_sfo</span> <span class="op"><-</span> <span class="kw">function</span><span class="op">(</span><span class="va">k1</span>, <span class="va">k2</span>, <span class="va">g</span>, <span class="va">f_parent_to_m1</span>, <span class="va">k_m1</span><span class="op">)</span> <span class="op">{</span></span></span> -<span class="r-in"><span>  <span class="fu"><a href="mkinpredict.html">mkinpredict</a></span><span class="op">(</span><span class="va">dfop_sfo</span>,</span></span> -<span class="r-in"><span>    <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span>k1 <span class="op">=</span> <span class="va">k1</span>, k2 <span class="op">=</span> <span class="va">k2</span>, g <span class="op">=</span> <span class="va">g</span>, f_parent_to_m1 <span class="op">=</span> <span class="va">f_parent_to_m1</span>, k_m1 <span class="op">=</span> <span class="va">k_m1</span><span class="op">)</span>,</span></span> -<span class="r-in"><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>, m1 <span class="op">=</span> <span class="fl">0</span><span class="op">)</span>,</span></span> -<span class="r-in"><span>    <span class="va">sampling_times</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="op">}</span></span></span> -<span class="r-in"><span></span></span> -<span class="r-in"><span><span class="va">ds_mean_dfop_sfo</span> <span class="op"><-</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="fl">1</span><span class="op">:</span><span class="fl">5</span>, <span class="kw">function</span><span class="op">(</span><span class="va">i</span><span class="op">)</span> <span class="op">{</span></span></span> -<span class="r-in"><span>  <span class="fu"><a href="mkinpredict.html">mkinpredict</a></span><span class="op">(</span><span class="va">dfop_sfo</span>,</span></span> -<span class="r-in"><span>    <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span>k1 <span class="op">=</span> <span class="va">k1_in</span><span class="op">[</span><span class="va">i</span><span class="op">]</span>, k2 <span class="op">=</span> <span class="va">k2_in</span><span class="op">[</span><span class="va">i</span><span class="op">]</span>, g <span class="op">=</span> <span class="va">g_in</span><span class="op">[</span><span class="va">i</span><span class="op">]</span>,</span></span> -<span class="r-in"><span>      f_parent_to_m1 <span class="op">=</span> <span class="va">f_parent_to_m1_in</span><span class="op">[</span><span class="va">i</span><span class="op">]</span>, k_m1 <span class="op">=</span> <span class="va">k_m1_in</span><span class="op">[</span><span class="va">i</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/c.html" class="external-link">c</a></span><span class="op">(</span>parent <span class="op">=</span> <span class="fl">100</span>, m1 <span class="op">=</span> <span class="fl">0</span><span class="op">)</span>,</span></span> -<span class="r-in"><span>    <span class="va">sampling_times</span><span class="op">)</span></span></span> -<span class="r-in"><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_mean_dfop_sfo</span><span class="op">)</span> <span class="op"><-</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">"ds"</span>, <span class="fl">1</span><span class="op">:</span><span class="fl">5</span><span class="op">)</span></span></span> -<span class="r-in"><span></span></span> -<span class="r-in"><span><span class="va">ds_syn_dfop_sfo</span> <span class="op"><-</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">ds_mean_dfop_sfo</span>, <span class="kw">function</span><span class="op">(</span><span class="va">ds</span><span class="op">)</span> <span class="op">{</span></span></span> -<span class="r-in"><span>  <span class="fu"><a href="add_err.html">add_err</a></span><span class="op">(</span><span class="va">ds</span>,</span></span> -<span class="r-in"><span>    sdfunc <span class="op">=</span> <span class="kw">function</span><span class="op">(</span><span class="va">value</span><span class="op">)</span> <span class="fu"><a href="https://rdrr.io/r/base/MathFun.html" class="external-link">sqrt</a></span><span class="op">(</span><span class="fl">1</span><span class="op">^</span><span class="fl">2</span> <span class="op">+</span> <span class="va">value</span><span class="op">^</span><span class="fl">2</span> <span class="op">*</span> <span class="fl">0.07</span><span class="op">^</span><span class="fl">2</span><span class="op">)</span>,</span></span> -<span class="r-in"><span>    n <span class="op">=</span> <span class="fl">1</span><span class="op">)</span><span class="op">[[</span><span class="fl">1</span><span class="op">]</span><span class="op">]</span></span></span> -<span class="r-in"><span><span class="op">}</span><span class="op">)</span></span></span> -<span class="r-in"><span></span></span> -<span class="r-in"><span><span class="co"># \dontrun{</span></span></span> -<span class="r-in"><span><span class="co"># Evaluate using mmkin and saem</span></span></span> -<span class="r-in"><span><span class="va">f_mmkin_dfop_sfo</span> <span class="op"><-</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 class="va">dfop_sfo</span><span class="op">)</span>, <span class="va">ds_syn_dfop_sfo</span>,</span></span> -<span class="r-in"><span>  quiet <span class="op">=</span> <span class="cn">TRUE</span>, error_model <span class="op">=</span> <span class="st">"tc"</span>, cores <span class="op">=</span> <span class="fl">5</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="va">f_saem_dfop_sfo</span> <span class="op"><-</span> <span class="fu"><a href="saem.html">saem</a></span><span class="op">(</span><span class="va">f_mmkin_dfop_sfo</span><span class="op">)</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">#></span> Kinetic nonlinear mixed-effects model fit by SAEM</span> -<span class="r-out co"><span class="r-pr">#></span> Structural model:</span> -<span class="r-out co"><span class="r-pr">#></span> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *</span> -<span class="r-out co"><span class="r-pr">#></span>            time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time)))</span> -<span class="r-out co"><span class="r-pr">#></span>            * parent</span> -<span class="r-out co"><span class="r-pr">#></span> d_m1/dt = + f_parent_to_m1 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g)</span> -<span class="r-out co"><span class="r-pr">#></span>            * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) *</span> -<span class="r-out co"><span class="r-pr">#></span>            exp(-k2 * time))) * parent - k_m1 * m1</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Data:</span> -<span class="r-out co"><span class="r-pr">#></span> 171 observations of 2 variable(s) grouped in 5 datasets</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Likelihood computed by importance sampling</span> -<span class="r-out co"><span class="r-pr">#></span>     AIC   BIC logLik</span> -<span class="r-out co"><span class="r-pr">#></span>   810.8 805.4 -391.4</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Fitted parameters:</span> -<span class="r-out co"><span class="r-pr">#></span>                     estimate     lower     upper</span> -<span class="r-out co"><span class="r-pr">#></span> parent_0           100.86947  97.81542 103.92353</span> -<span class="r-out co"><span class="r-pr">#></span> log_k_m1            -4.06947  -4.16944  -3.96950</span> -<span class="r-out co"><span class="r-pr">#></span> f_parent_qlogis     -0.93256  -1.34200  -0.52312</span> -<span class="r-out co"><span class="r-pr">#></span> log_k1              -2.37017  -2.72660  -2.01375</span> -<span class="r-out co"><span class="r-pr">#></span> log_k2              -4.06264  -4.21344  -3.91184</span> -<span class="r-out co"><span class="r-pr">#></span> g_qlogis            -0.02174  -0.45898   0.41549</span> -<span class="r-out co"><span class="r-pr">#></span> a.1                  0.87598   0.67275   1.07922</span> -<span class="r-out co"><span class="r-pr">#></span> b.1                  0.07949   0.06389   0.09509</span> -<span class="r-out co"><span class="r-pr">#></span> SD.parent_0          0.19170 -30.36286  30.74626</span> -<span class="r-out co"><span class="r-pr">#></span> SD.log_k_m1          0.01883  -0.28736   0.32502</span> -<span class="r-out co"><span class="r-pr">#></span> SD.f_parent_qlogis   0.44300   0.16391   0.72209</span> -<span class="r-out co"><span class="r-pr">#></span> SD.log_k1            0.35320   0.09661   0.60978</span> -<span class="r-out co"><span class="r-pr">#></span> SD.log_k2            0.13707   0.02359   0.25056</span> -<span class="r-out co"><span class="r-pr">#></span> SD.g_qlogis          0.37478   0.04490   0.70467</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_sfo</span><span class="op">)</span></span></span> -<span class="r-out co"><span class="r-pr">#></span> [1] "sd(parent_0)" "sd(log_k_m1)"</span> -<span class="r-in"><span><span class="va">f_saem_dfop_sfo_2</span> <span class="op"><-</span> <span class="fu"><a href="https://rdrr.io/r/stats/update.html" class="external-link">update</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>,</span></span> -<span class="r-in"><span>  no_random_effect <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">"parent_0"</span>, <span class="st">"log_k_m1"</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_sfo_2</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/pkg/nlme/man/intervals.html" class="external-link">intervals</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo_2</span><span class="op">)</span></span></span> -<span class="r-out co"><span class="r-pr">#></span> Approximate 95% confidence intervals</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span>  Fixed effects:</span> -<span class="r-out co"><span class="r-pr">#></span>                      lower         est.        upper</span> -<span class="r-out co"><span class="r-pr">#></span> parent_0       98.36731429 101.42508066 104.48284703</span> -<span class="r-out co"><span class="r-pr">#></span> k_m1            0.01513234   0.01670094   0.01843214</span> -<span class="r-out co"><span class="r-pr">#></span> f_parent_to_m1  0.20221431   0.27608850   0.36461630</span> -<span class="r-out co"><span class="r-pr">#></span> k1              0.06915073   0.09759718   0.13774560</span> -<span class="r-out co"><span class="r-pr">#></span> k2              0.01487068   0.01740389   0.02036863</span> -<span class="r-out co"><span class="r-pr">#></span> g               0.37365671   0.48384821   0.59563299</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span>  Random effects:</span> -<span class="r-out co"><span class="r-pr">#></span>                          lower      est.     upper</span> -<span class="r-out co"><span class="r-pr">#></span> sd(f_parent_qlogis) 0.16439770 0.4427585 0.7211193</span> -<span class="r-out co"><span class="r-pr">#></span> sd(log_k1)          0.08304243 0.3345213 0.5860002</span> -<span class="r-out co"><span class="r-pr">#></span> sd(log_k2)          0.03146410 0.1490210 0.2665779</span> -<span class="r-out co"><span class="r-pr">#></span> sd(g_qlogis)        0.06216385 0.4023430 0.7425221</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span>  </span> -<span class="r-out co"><span class="r-pr">#></span>          lower       est.      upper</span> -<span class="r-out co"><span class="r-pr">#></span> a.1 0.67696663 0.87777355 1.07858048</span> -<span class="r-out co"><span class="r-pr">#></span> b.1 0.06363957 0.07878001 0.09392044</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_2</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">#></span> saemix version used for fitting:      3.2 </span> -<span class="r-out co"><span class="r-pr">#></span> mkin version used for pre-fitting:  1.2.3 </span> -<span class="r-out co"><span class="r-pr">#></span> R version used for fitting:         4.2.3 </span> -<span class="r-out co"><span class="r-pr">#></span> Date of fit:     Sun Apr 16 08:34:58 2023 </span> -<span class="r-out co"><span class="r-pr">#></span> Date of summary: Sun Apr 16 08:34:58 2023 </span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Equations:</span> -<span class="r-out co"><span class="r-pr">#></span> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *</span> -<span class="r-out co"><span class="r-pr">#></span>            time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time)))</span> -<span class="r-out co"><span class="r-pr">#></span>            * parent</span> -<span class="r-out co"><span class="r-pr">#></span> d_m1/dt = + f_parent_to_m1 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g)</span> -<span class="r-out co"><span class="r-pr">#></span>            * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) *</span> -<span class="r-out co"><span class="r-pr">#></span>            exp(-k2 * time))) * parent - k_m1 * m1</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Data:</span> -<span class="r-out co"><span class="r-pr">#></span> 171 observations of 2 variable(s) grouped in 5 datasets</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Model predictions using solution type analytical </span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Fitted in 9.384 s</span> -<span class="r-out co"><span class="r-pr">#></span> Using 300, 100 iterations and 10 chains</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Variance model: Two-component variance function </span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Starting values for degradation parameters:</span> -<span class="r-out co"><span class="r-pr">#></span>        parent_0        log_k_m1 f_parent_qlogis          log_k1          log_k2 </span> -<span class="r-out co"><span class="r-pr">#></span>       101.65645        -4.05368        -0.94311        -2.35943        -4.07006 </span> -<span class="r-out co"><span class="r-pr">#></span>        g_qlogis </span> -<span class="r-out co"><span class="r-pr">#></span>        -0.01132 </span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Fixed degradation parameter values:</span> -<span class="r-out co"><span class="r-pr">#></span> None</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Starting values for random effects (square root of initial entries in omega):</span> -<span class="r-out co"><span class="r-pr">#></span>                 parent_0 log_k_m1 f_parent_qlogis log_k1 log_k2 g_qlogis</span> -<span class="r-out co"><span class="r-pr">#></span> parent_0           6.742   0.0000          0.0000 0.0000 0.0000    0.000</span> -<span class="r-out co"><span class="r-pr">#></span> log_k_m1           0.000   0.2236          0.0000 0.0000 0.0000    0.000</span> -<span class="r-out co"><span class="r-pr">#></span> f_parent_qlogis    0.000   0.0000          0.5572 0.0000 0.0000    0.000</span> -<span class="r-out co"><span class="r-pr">#></span> log_k1             0.000   0.0000          0.0000 0.8031 0.0000    0.000</span> -<span class="r-out co"><span class="r-pr">#></span> log_k2             0.000   0.0000          0.0000 0.0000 0.2931    0.000</span> -<span class="r-out co"><span class="r-pr">#></span> g_qlogis           0.000   0.0000          0.0000 0.0000 0.0000    0.807</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Starting values for error model parameters:</span> -<span class="r-out co"><span class="r-pr">#></span> a.1 b.1 </span> -<span class="r-out co"><span class="r-pr">#></span>   1   1 </span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Results:</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Likelihood computed by importance sampling</span> -<span class="r-out co"><span class="r-pr">#></span>   AIC   BIC logLik</span> -<span class="r-out co"><span class="r-pr">#></span>   807 802.3 -391.5</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Optimised parameters:</span> -<span class="r-out co"><span class="r-pr">#></span>                         est.    lower     upper</span> -<span class="r-out co"><span class="r-pr">#></span> parent_0           101.42508 98.36731 104.48285</span> -<span class="r-out co"><span class="r-pr">#></span> log_k_m1            -4.09229 -4.19092  -3.99366</span> -<span class="r-out co"><span class="r-pr">#></span> f_parent_qlogis     -0.96395 -1.37251  -0.55538</span> -<span class="r-out co"><span class="r-pr">#></span> log_k1              -2.32691 -2.67147  -1.98235</span> -<span class="r-out co"><span class="r-pr">#></span> log_k2              -4.05106 -4.20836  -3.89376</span> -<span class="r-out co"><span class="r-pr">#></span> g_qlogis            -0.06463 -0.51656   0.38730</span> -<span class="r-out co"><span class="r-pr">#></span> a.1                  0.87777  0.67697   1.07858</span> -<span class="r-out co"><span class="r-pr">#></span> b.1                  0.07878  0.06364   0.09392</span> -<span class="r-out co"><span class="r-pr">#></span> SD.f_parent_qlogis   0.44276  0.16440   0.72112</span> -<span class="r-out co"><span class="r-pr">#></span> SD.log_k1            0.33452  0.08304   0.58600</span> -<span class="r-out co"><span class="r-pr">#></span> SD.log_k2            0.14902  0.03146   0.26658</span> -<span class="r-out co"><span class="r-pr">#></span> SD.g_qlogis          0.40234  0.06216   0.74252</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Correlation: </span> -<span class="r-out co"><span class="r-pr">#></span>                 parnt_0 lg_k_m1 f_prnt_ log_k1  log_k2 </span> -<span class="r-out co"><span class="r-pr">#></span> log_k_m1        -0.4693                                </span> -<span class="r-out co"><span class="r-pr">#></span> f_parent_qlogis -0.2378  0.2595                        </span> -<span class="r-out co"><span class="r-pr">#></span> log_k1           0.1720 -0.1593 -0.0669                </span> -<span class="r-out co"><span class="r-pr">#></span> log_k2           0.0179  0.0594  0.0035  0.1995        </span> -<span class="r-out co"><span class="r-pr">#></span> g_qlogis         0.1073 -0.1060 -0.0322 -0.2299 -0.3168</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Random effects:</span> -<span class="r-out co"><span class="r-pr">#></span>                      est.   lower  upper</span> -<span class="r-out co"><span class="r-pr">#></span> SD.f_parent_qlogis 0.4428 0.16440 0.7211</span> -<span class="r-out co"><span class="r-pr">#></span> SD.log_k1          0.3345 0.08304 0.5860</span> -<span class="r-out co"><span class="r-pr">#></span> SD.log_k2          0.1490 0.03146 0.2666</span> -<span class="r-out co"><span class="r-pr">#></span> SD.g_qlogis        0.4023 0.06216 0.7425</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Variance model:</span> -<span class="r-out co"><span class="r-pr">#></span>        est.   lower   upper</span> -<span class="r-out co"><span class="r-pr">#></span> a.1 0.87777 0.67697 1.07858</span> -<span class="r-out co"><span class="r-pr">#></span> b.1 0.07878 0.06364 0.09392</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Backtransformed parameters:</span> -<span class="r-out co"><span class="r-pr">#></span>                    est.    lower     upper</span> -<span class="r-out co"><span class="r-pr">#></span> parent_0       101.4251 98.36731 104.48285</span> -<span class="r-out co"><span class="r-pr">#></span> k_m1             0.0167  0.01513   0.01843</span> -<span class="r-out co"><span class="r-pr">#></span> f_parent_to_m1   0.2761  0.20221   0.36462</span> -<span class="r-out co"><span class="r-pr">#></span> k1               0.0976  0.06915   0.13775</span> -<span class="r-out co"><span class="r-pr">#></span> k2               0.0174  0.01487   0.02037</span> -<span class="r-out co"><span class="r-pr">#></span> g                0.4838  0.37366   0.59563</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Resulting formation fractions:</span> -<span class="r-out co"><span class="r-pr">#></span>                 ff</span> -<span class="r-out co"><span class="r-pr">#></span> parent_m1   0.2761</span> -<span class="r-out co"><span class="r-pr">#></span> parent_sink 0.7239</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Estimated disappearance times:</span> -<span class="r-out co"><span class="r-pr">#></span>         DT50   DT90 DT50back DT50_k1 DT50_k2</span> -<span class="r-out co"><span class="r-pr">#></span> parent 15.54  94.33     28.4   7.102   39.83</span> -<span class="r-out co"><span class="r-pr">#></span> m1     41.50 137.87       NA      NA      NA</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Data:</span> -<span class="r-out co"><span class="r-pr">#></span>    ds   name time observed predicted  residual    std standardized</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent    0     89.8 1.014e+02 -11.62508 8.0383     -1.44620</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent    0    104.1 1.014e+02   2.67492 8.0383      0.33277</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent    1     88.7 9.650e+01  -7.80311 7.6530     -1.01961</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent    1     95.5 9.650e+01  -1.00311 7.6530     -0.13107</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent    3     81.8 8.753e+01  -5.72638 6.9510     -0.82382</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent    3     94.5 8.753e+01   6.97362 6.9510      1.00326</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent    7     71.5 7.254e+01  -1.04133 5.7818     -0.18010</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent    7     70.3 7.254e+01  -2.24133 5.7818     -0.38765</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent   14     54.2 5.349e+01   0.71029 4.3044      0.16502</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent   14     49.6 5.349e+01  -3.88971 4.3044     -0.90366</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent   28     31.5 3.167e+01  -0.16616 2.6446     -0.06283</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent   28     28.8 3.167e+01  -2.86616 2.6446     -1.08379</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent   60     12.1 1.279e+01  -0.69287 1.3365     -0.51843</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent   60     13.6 1.279e+01   0.80713 1.3365      0.60392</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent   90      6.2 6.397e+00  -0.19718 1.0122     -0.19481</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent   90      8.3 6.397e+00   1.90282 1.0122      1.87996</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent  120      2.2 3.323e+00  -1.12320 0.9160     -1.22623</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1 parent  120      2.4 3.323e+00  -0.92320 0.9160     -1.00788</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1    1      0.3 1.179e+00  -0.87919 0.8827     -0.99605</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1    1      0.2 1.179e+00  -0.97919 0.8827     -1.10935</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1    3      2.2 3.273e+00  -1.07272 0.9149     -1.17256</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1    3      3.0 3.273e+00  -0.27272 0.9149     -0.29811</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1    7      6.5 6.559e+00  -0.05872 1.0186     -0.05765</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1    7      5.0 6.559e+00  -1.55872 1.0186     -1.53032</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1   14     10.2 1.016e+01   0.03787 1.1880      0.03188</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1   14      9.5 1.016e+01  -0.66213 1.1880     -0.55734</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1   28     12.2 1.268e+01  -0.47913 1.3297     -0.36032</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1   28     13.4 1.268e+01   0.72087 1.3297      0.54211</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1   60     11.8 1.078e+01   1.02493 1.2211      0.83936</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1   60     13.2 1.078e+01   2.42493 1.2211      1.98588</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1   90      6.6 7.705e+00  -1.10464 1.0672     -1.03509</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1   90      9.3 7.705e+00   1.59536 1.0672      1.49491</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1  120      3.5 5.236e+00  -1.73617 0.9699     -1.79010</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 1     m1  120      5.4 5.236e+00   0.16383 0.9699      0.16892</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent    0    118.0 1.014e+02  16.57492 8.0383      2.06198</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent    0     99.8 1.014e+02  -1.62508 8.0383     -0.20217</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent    1     90.2 9.599e+01  -5.79045 7.6129     -0.76061</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent    1     94.6 9.599e+01  -1.39045 7.6129     -0.18264</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent    3     96.1 8.652e+01   9.57931 6.8724      1.39388</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent    3     78.4 8.652e+01  -8.12069 6.8724     -1.18164</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent    7     77.9 7.197e+01   5.93429 5.7370      1.03439</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent    7     77.7 7.197e+01   5.73429 5.7370      0.99953</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent   14     56.0 5.555e+01   0.44657 4.4637      0.10005</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent   14     54.7 5.555e+01  -0.85343 4.4637     -0.19120</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent   28     36.6 3.853e+01  -1.93170 3.1599     -0.61132</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent   28     36.8 3.853e+01  -1.73170 3.1599     -0.54803</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent   60     22.1 2.110e+01   1.00360 1.8795      0.53396</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent   60     24.7 2.110e+01   3.60360 1.8795      1.91728</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent   90     12.4 1.250e+01  -0.09712 1.3190     -0.07363</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent   90     10.8 1.250e+01  -1.69712 1.3190     -1.28667</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent  120      6.8 7.419e+00  -0.61913 1.0546     -0.58709</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2 parent  120      7.9 7.419e+00   0.48087 1.0546      0.45599</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1    1      1.3 1.422e+00  -0.12194 0.8849     -0.13781</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1    3      3.7 3.831e+00  -0.13149 0.9282     -0.14166</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1    3      4.7 3.831e+00   0.86851 0.9282      0.93567</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1    7      8.1 7.292e+00   0.80812 1.0490      0.77034</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1    7      7.9 7.292e+00   0.60812 1.0490      0.57969</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1   14     10.1 1.055e+01  -0.45332 1.2090     -0.37495</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1   14     10.3 1.055e+01  -0.25332 1.2090     -0.20953</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1   28     10.7 1.230e+01  -1.59960 1.3074     -1.22347</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1   28     12.2 1.230e+01  -0.09960 1.3074     -0.07618</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1   60     10.7 1.065e+01   0.05342 1.2141      0.04400</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1   60     12.5 1.065e+01   1.85342 1.2141      1.52661</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1   90      9.1 8.196e+00   0.90368 1.0897      0.82930</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1   90      7.4 8.196e+00  -0.79632 1.0897     -0.73078</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1  120      6.1 5.997e+00   0.10252 0.9969      0.10284</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 2     m1  120      4.5 5.997e+00  -1.49748 0.9969     -1.50220</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent    0    106.2 1.014e+02   4.77492 8.0383      0.59402</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent    0    106.9 1.014e+02   5.47492 8.0383      0.68110</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent    1    107.4 9.390e+01  13.49935 7.4494      1.81214</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent    1     96.1 9.390e+01   2.19935 7.4494      0.29524</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent    3     79.4 8.152e+01  -2.12307 6.4821     -0.32753</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent    3     82.6 8.152e+01   1.07693 6.4821      0.16614</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent    7     63.9 6.446e+01  -0.55834 5.1533     -0.10834</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent    7     62.4 6.446e+01  -2.05834 5.1533     -0.39942</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent   14     51.0 4.826e+01   2.74073 3.9019      0.70241</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent   14     47.1 4.826e+01  -1.15927 3.9019     -0.29711</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent   28     36.1 3.424e+01   1.86399 2.8364      0.65718</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent   28     36.6 3.424e+01   2.36399 2.8364      0.83346</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent   60     20.1 1.968e+01   0.42172 1.7815      0.23672</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent   60     19.8 1.968e+01   0.12172 1.7815      0.06833</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent   90     11.3 1.195e+01  -0.64633 1.2869     -0.50222</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent   90     10.7 1.195e+01  -1.24633 1.2869     -0.96844</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent  120      8.2 7.255e+00   0.94532 1.0474      0.90251</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3 parent  120      7.3 7.255e+00   0.04532 1.0474      0.04327</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1    0      0.8 2.956e-11   0.80000 0.8778      0.91140</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1    1      1.8 1.758e+00   0.04187 0.8886      0.04712</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1    1      2.3 1.758e+00   0.54187 0.8886      0.60978</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1    3      4.2 4.567e+00  -0.36697 0.9486     -0.38683</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1    3      4.1 4.567e+00  -0.46697 0.9486     -0.49224</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1    7      6.8 8.151e+00  -1.35124 1.0876     -1.24242</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1    7     10.1 8.151e+00   1.94876 1.0876      1.79182</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1   14     11.4 1.083e+01   0.57098 1.2240      0.46647</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1   14     12.8 1.083e+01   1.97098 1.2240      1.61022</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1   28     11.5 1.147e+01   0.03175 1.2597      0.02520</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1   28     10.6 1.147e+01  -0.86825 1.2597     -0.68928</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1   60      7.5 9.298e+00  -1.79834 1.1433     -1.57298</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1   60      8.6 9.298e+00  -0.69834 1.1433     -0.61083</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1   90      7.3 7.038e+00   0.26249 1.0382      0.25283</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1   90      8.1 7.038e+00   1.06249 1.0382      1.02340</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1  120      5.3 5.116e+00   0.18417 0.9659      0.19068</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 3     m1  120      3.8 5.116e+00  -1.31583 0.9659     -1.36232</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent    0    104.7 1.014e+02   3.27492 8.0383      0.40741</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent    0     88.3 1.014e+02 -13.12508 8.0383     -1.63281</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent    1     94.2 9.781e+01  -3.61183 7.7555     -0.46572</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent    1     94.6 9.781e+01  -3.21183 7.7555     -0.41414</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent    3     78.1 9.110e+01 -13.00467 7.2307     -1.79853</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent    3     96.5 9.110e+01   5.39533 7.2307      0.74617</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent    7     76.2 7.951e+01  -3.30511 6.3246     -0.52258</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent    7     77.8 7.951e+01  -1.70511 6.3246     -0.26960</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent   14     70.8 6.376e+01   7.03783 5.0993      1.38016</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent   14     67.3 6.376e+01   3.53783 5.0993      0.69379</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent   28     43.1 4.340e+01  -0.30456 3.5303     -0.08627</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent   28     45.1 4.340e+01   1.69544 3.5303      0.48026</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent   60     21.3 2.142e+01  -0.12077 1.9022     -0.06349</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent   60     23.5 2.142e+01   2.07923 1.9022      1.09308</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent   90     11.8 1.207e+01  -0.26813 1.2940     -0.20721</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent   90     12.1 1.207e+01   0.03187 1.2940      0.02463</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent  120      7.0 6.954e+00   0.04554 1.0347      0.04402</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4 parent  120      6.2 6.954e+00  -0.75446 1.0347     -0.72914</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1    0      1.6 1.990e-13   1.60000 0.8778      1.82279</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1    1      0.9 7.305e-01   0.16949 0.8797      0.19267</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1    3      3.7 2.051e+00   1.64896 0.8925      1.84753</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1    3      2.0 2.051e+00  -0.05104 0.8925     -0.05719</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1    7      3.6 4.204e+00  -0.60375 0.9382     -0.64354</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1    7      3.8 4.204e+00  -0.40375 0.9382     -0.43036</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1   14      7.1 6.760e+00   0.34021 1.0267      0.33137</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1   14      6.6 6.760e+00  -0.15979 1.0267     -0.15563</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1   28      9.5 9.011e+00   0.48856 1.1289      0.43277</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1   28      9.3 9.011e+00   0.28856 1.1289      0.25561</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1   60      8.3 8.611e+00  -0.31077 1.1093     -0.28014</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1   60      9.0 8.611e+00   0.38923 1.1093      0.35086</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1   90      6.6 6.678e+00  -0.07753 1.0233     -0.07576</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1   90      7.7 6.678e+00   1.02247 1.0233      0.99915</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1  120      3.7 4.847e+00  -1.14679 0.9572     -1.19804</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 4     m1  120      3.5 4.847e+00  -1.34679 0.9572     -1.40698</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent    0    110.4 1.014e+02   8.97492 8.0383      1.11651</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent    0    112.1 1.014e+02  10.67492 8.0383      1.32800</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent    1     93.5 9.466e+01  -1.16118 7.5089     -0.15464</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent    1     91.0 9.466e+01  -3.66118 7.5089     -0.48758</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent    3     71.0 8.302e+01 -12.01844 6.5988     -1.82130</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent    3     89.7 8.302e+01   6.68156 6.5988      1.01254</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent    7     60.4 6.563e+01  -5.22574 5.2440     -0.99652</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent    7     59.1 6.563e+01  -6.52574 5.2440     -1.24442</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent   14     56.5 4.727e+01   9.22621 3.8263      2.41128</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent   14     47.0 4.727e+01  -0.27379 3.8263     -0.07156</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent   28     30.2 3.103e+01  -0.83405 2.5977     -0.32108</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent   28     23.9 3.103e+01  -7.13405 2.5977     -2.74634</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent   60     17.0 1.800e+01  -0.99696 1.6675     -0.59787</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent   60     18.7 1.800e+01   0.70304 1.6675      0.42161</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent   90     11.3 1.167e+01  -0.36809 1.2710     -0.28961</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent   90     11.9 1.167e+01   0.23191 1.2710      0.18246</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent  120      9.0 7.595e+00   1.40496 1.0623      1.32256</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5 parent  120      8.1 7.595e+00   0.50496 1.0623      0.47535</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1    0      0.7 0.000e+00   0.70000 0.8778      0.79747</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1    1      3.0 3.158e+00  -0.15799 0.9123     -0.17317</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1    1      2.6 3.158e+00  -0.55799 0.9123     -0.61160</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1    3      5.1 8.443e+00  -3.34286 1.1013     -3.03535</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1    3      7.5 8.443e+00  -0.94286 1.1013     -0.85613</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1    7     16.5 1.580e+01   0.69781 1.5232      0.45811</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1    7     19.0 1.580e+01   3.19781 1.5232      2.09935</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1   14     22.9 2.216e+01   0.73604 1.9543      0.37663</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1   14     23.2 2.216e+01   1.03604 1.9543      0.53014</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1   28     22.2 2.423e+01  -2.03128 2.1011     -0.96678</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1   28     24.4 2.423e+01   0.16872 2.1011      0.08030</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1   60     15.5 1.876e+01  -3.25610 1.7187     -1.89455</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1   60     19.8 1.876e+01   1.04390 1.7187      0.60739</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1   90     14.9 1.366e+01   1.23585 1.3890      0.88976</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1   90     14.2 1.366e+01   0.53585 1.3890      0.38579</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1  120     10.9 9.761e+00   1.13911 1.1670      0.97613</span> -<span class="r-out co"><span class="r-pr">#></span>  ds 5     m1  120     10.4 9.761e+00   0.63911 1.1670      0.54767</span> -<span class="r-in"><span><span class="co"># Add a correlation between random effects of g and k2</span></span></span> -<span class="r-in"><span><span class="va">cov_model_3</span> <span class="op"><-</span> <span class="va">f_saem_dfop_sfo_2</span><span class="op">$</span><span class="va">so</span><span class="op">@</span><span class="va">model</span><span class="op">@</span><span class="va">covariance.model</span></span></span> -<span class="r-in"><span><span class="va">cov_model_3</span><span class="op">[</span><span class="st">"log_k2"</span>, <span class="st">"g_qlogis"</span><span class="op">]</span> <span class="op"><-</span> <span class="fl">1</span></span></span> -<span class="r-in"><span><span class="va">cov_model_3</span><span class="op">[</span><span class="st">"g_qlogis"</span>, <span class="st">"log_k2"</span><span class="op">]</span> <span class="op"><-</span> <span class="fl">1</span></span></span> -<span class="r-in"><span><span class="va">f_saem_dfop_sfo_3</span> <span class="op"><-</span> <span class="fu"><a href="https://rdrr.io/r/stats/update.html" class="external-link">update</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>,</span></span> -<span class="r-in"><span>  covariance.model <span class="op">=</span> <span class="va">cov_model_3</span><span class="op">)</span></span></span> -<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/pkg/nlme/man/intervals.html" class="external-link">intervals</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo_3</span><span class="op">)</span></span></span> -<span class="r-out co"><span class="r-pr">#></span> Approximate 95% confidence intervals</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span>  Fixed effects:</span> -<span class="r-out co"><span class="r-pr">#></span>                      lower         est.        upper</span> -<span class="r-out co"><span class="r-pr">#></span> parent_0       98.39888363 101.48951337 104.58014311</span> -<span class="r-out co"><span class="r-pr">#></span> k_m1            0.01508704   0.01665986   0.01839665</span> -<span class="r-out co"><span class="r-pr">#></span> f_parent_to_m1  0.20141557   0.27540583   0.36418131</span> -<span class="r-out co"><span class="r-pr">#></span> k1              0.07708759   0.10430866   0.14114200</span> -<span class="r-out co"><span class="r-pr">#></span> k2              0.01476621   0.01786384   0.02161129</span> -<span class="r-out co"><span class="r-pr">#></span> g               0.33679867   0.45083525   0.57028162</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span>  Random effects:</span> -<span class="r-out co"><span class="r-pr">#></span>                             lower       est.      upper</span> -<span class="r-out co"><span class="r-pr">#></span> sd(f_parent_qlogis)    0.38085375  0.4441841  0.5075145</span> -<span class="r-out co"><span class="r-pr">#></span> sd(log_k1)             0.04774819  0.2660384  0.4843286</span> -<span class="r-out co"><span class="r-pr">#></span> sd(log_k2)            -0.63842736  0.1977024  1.0338321</span> -<span class="r-out co"><span class="r-pr">#></span> sd(g_qlogis)           0.22711289  0.4502227  0.6733326</span> -<span class="r-out co"><span class="r-pr">#></span> corr(log_k2,g_qlogis) -0.83271473 -0.6176939 -0.4026730</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span>  </span> -<span class="r-out co"><span class="r-pr">#></span>          lower       est.      upper</span> -<span class="r-out co"><span class="r-pr">#></span> a.1 0.67347568 0.87437392 1.07527216</span> -<span class="r-out co"><span class="r-pr">#></span> b.1 0.06393032 0.07912417 0.09431802</span> -<span class="r-in"><span><span class="co"># The correlation does not improve the fit judged by AIC and BIC, although</span></span></span> -<span class="r-in"><span><span class="co"># the likelihood is higher with the additional parameter</span></span></span> -<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/stats/anova.html" class="external-link">anova</a></span><span class="op">(</span><span class="va">f_saem_dfop_sfo</span>, <span class="va">f_saem_dfop_sfo_2</span>, <span class="va">f_saem_dfop_sfo_3</span><span class="op">)</span></span></span> -<span class="r-out co"><span class="r-pr">#></span> Data: 171 observations of 2 variable(s) grouped in 5 datasets</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span>                   npar    AIC    BIC     Lik</span> -<span class="r-out co"><span class="r-pr">#></span> f_saem_dfop_sfo_2   12 806.96 802.27 -391.48</span> -<span class="r-out co"><span class="r-pr">#></span> f_saem_dfop_sfo_3   13 807.99 802.91 -391.00</span> -<span class="r-out co"><span class="r-pr">#></span> f_saem_dfop_sfo     14 810.83 805.36 -391.42</span> -<span class="r-in"><span><span class="co"># }</span></span></span> -<span class="r-in"><span></span></span> -</code></pre></div> -    </div> -  </div> -  <div class="col-md-3 hidden-xs hidden-sm" id="pkgdown-sidebar"> -    <nav id="toc" data-toggle="toc" class="sticky-top"><h2 data-toc-skip>Contents</h2> -    </nav></div> -</div> - - -      <footer><div class="copyright"> -  <p></p><p>Developed by Johannes Ranke.</p> -</div> - -<div class="pkgdown"> -  <p></p><p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.0.7.</p> -</div> - -      </footer></div> - -   - - -   - -  </body></html> - | 
