diff options
Diffstat (limited to 'docs/reference/saem.html')
| -rw-r--r-- | docs/reference/saem.html | 39 |
1 files changed, 11 insertions, 28 deletions
diff --git a/docs/reference/saem.html b/docs/reference/saem.html index 159bc00f..4e81388b 100644 --- a/docs/reference/saem.html +++ b/docs/reference/saem.html @@ -11,7 +11,7 @@ Expectation Maximisation algorithm (SAEM)."></head><body> <a class="navbar-brand me-2" href="../index.html">mkin</a> - <small class="nav-text text-muted me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="Released version">1.2.11</small> + <small class="nav-text text-muted me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="Released version">1.2.10</small> <button class="navbar-toggler" type="button" data-bs-toggle="collapse" data-bs-target="#navbar" aria-controls="navbar" aria-expanded="false" aria-label="Toggle navigation"> @@ -402,10 +402,10 @@ using <a href="mmkin.html">mmkin</a>.</p> <span class="r-plt img"><img src="saem-4.png" alt="" width="700" height="433"></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">#></span> saemix version used for fitting: 3.4 </span> -<span class="r-out co"><span class="r-pr">#></span> mkin version used for pre-fitting: 1.2.11 </span> -<span class="r-out co"><span class="r-pr">#></span> R version used for fitting: 4.5.1 </span> -<span class="r-out co"><span class="r-pr">#></span> Date of fit: Fri Sep 12 20:34:55 2025 </span> -<span class="r-out co"><span class="r-pr">#></span> Date of summary: Fri Sep 12 20:34:55 2025 </span> +<span class="r-out co"><span class="r-pr">#></span> mkin version used for pre-fitting: 1.2.10 </span> +<span class="r-out co"><span class="r-pr">#></span> R version used for fitting: 4.5.2 </span> +<span class="r-out co"><span class="r-pr">#></span> Date of fit: Fri Nov 28 09:18:42 2025 </span> +<span class="r-out co"><span class="r-pr">#></span> Date of summary: Fri Nov 28 09:18:42 2025 </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> @@ -420,7 +420,7 @@ using <a href="mmkin.html">mmkin</a>.</p> <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 4.017 s</span> +<span class="r-out co"><span class="r-pr">#></span> Fitted in 3.833 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: Constant variance </span> @@ -682,28 +682,11 @@ using <a href="mmkin.html">mmkin</a>.</p> <span class="r-out co"><span class="r-pr">#></span> Dataset 10 A1 120 14.3 12.79238 1.50762 1.882 0.80111</span> <span class="r-out co"><span class="r-pr">#></span> Dataset 10 A1 120 12.1 12.79238 -0.69238 1.882 -0.36791</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="va">f_saem_dfop_sfo_deSolve</span> <span class="op"><-</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>, solution_type <span class="op">=</span> <span class="st">"deSolve"</span>,</span></span> -<span class="r-in"><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">200</span>, <span class="fl">80</span><span class="op">)</span><span class="op">)</span></span></span> -<span class="r-out co"><span class="r-pr">#></span> DINTDY- T (=R1) illegal </span> -<span class="r-out co"><span class="r-pr">#></span> In above message, R1 = 70</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> T not in interval TCUR - HU (= R1) to TCUR (=R2) </span> -<span class="r-out co"><span class="r-pr">#></span> In above message, R1 = 53.1042, R2 = 56.6326</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> DINTDY- T (=R1) illegal </span> -<span class="r-out co"><span class="r-pr">#></span> In above message, R1 = 91</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> T not in interval TCUR - HU (= R1) to TCUR (=R2) </span> -<span class="r-out co"><span class="r-pr">#></span> In above message, R1 = 53.1042, R2 = 56.6326</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> DLSODA- Trouble in DINTDY. ITASK = I1, TOUT = R1</span> -<span class="r-out co"><span class="r-pr">#></span> In above message, I1 = 1</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> In above message, R1 = 91</span> -<span class="r-out co"><span class="r-pr">#></span> </span> -<span class="r-out co"><span class="r-pr">#></span> Error in deSolve::lsoda(y = odeini, times = outtimes_deSolve, func = lsoda_func, : </span> -<span class="r-out co"><span class="r-pr">#></span> illegal input detected before taking any integration steps - see written message</span> +<span class="r-in"><span><span class="co"># The following took about 6 minutes and currently fails with</span></span></span> +<span class="r-in"><span><span class="co"># "Error in deSolve::lsoda(y = odeini, times = outtimes_deSolve, func = lsoda_func, : </span></span></span> +<span class="r-in"><span><span class="co"># illegal input detected before taking any integration steps - see written message"</span></span></span> +<span class="r-in"><span><span class="co">#f_saem_dfop_sfo_deSolve <- saem(f_mmkin["DFOP-SFO", ], solution_type = "deSolve",</span></span></span> +<span class="r-in"><span><span class="co"># nbiter.saemix = c(200, 80))</span></span></span> <span class="r-in"><span></span></span> <span class="r-in"><span><span class="co">#anova(</span></span></span> <span class="r-in"><span><span class="co"># f_saem_dfop_sfo,</span></span></span> |