aboutsummaryrefslogtreecommitdiff
path: root/docs/reference/mkinfit.html
diff options
context:
space:
mode:
Diffstat (limited to 'docs/reference/mkinfit.html')
-rw-r--r--docs/reference/mkinfit.html176
1 files changed, 104 insertions, 72 deletions
diff --git a/docs/reference/mkinfit.html b/docs/reference/mkinfit.html
index b75b7c95..09329a86 100644
--- a/docs/reference/mkinfit.html
+++ b/docs/reference/mkinfit.html
@@ -360,17 +360,17 @@ Per default, parameters in the kinetic models are internally transformed in
<h2 class="hasAnchor" id="examples"><a class="anchor" href="#examples"></a>Examples</h2>
<pre class="examples"><div class='input'><span class='co'># Use shorthand notation for parent only degradation</span>
<span class='no'>fit</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='st'>"FOMC"</span>, <span class='no'>FOCUS_2006_C</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
-<span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/summary'>summary</a></span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Could not calculate correlation; no covariance matrix</span></div><div class='output co'>#&gt; mkin version used for fitting: 0.9.49.4
+<span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/summary'>summary</a></span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.49.4
#&gt; R version used for fitting: 3.6.0
-#&gt; Date of fit: Tue May 7 08:08:23 2019
-#&gt; Date of summary: Tue May 7 08:08:23 2019
+#&gt; Date of fit: Tue May 7 08:36:16 2019
+#&gt; Date of summary: Tue May 7 08:36:16 2019
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent
#&gt;
#&gt; Model predictions using solution type analytical
#&gt;
-#&gt; Fitted using 66 model solutions performed in 0.139 s
+#&gt; Fitted using 222 model solutions performed in 0.89 s
#&gt;
#&gt; Error model:
#&gt; Constant variance
@@ -393,23 +393,28 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt; None
#&gt;
#&gt; Optimised, transformed parameters with symmetric confidence intervals:
-#&gt; Estimate Std. Error Lower Upper
-#&gt; parent_0 85.87000 NA NA NA
-#&gt; log_alpha 0.05192 NA NA NA
-#&gt; log_beta 0.65100 NA NA NA
-#&gt; sigma 1.85700 NA NA NA
+#&gt; Estimate Std. Error Lower Upper
+#&gt; parent_0 85.87000 1.8070 81.23000 90.5200
+#&gt; log_alpha 0.05192 0.1353 -0.29580 0.3996
+#&gt; log_beta 0.65100 0.2287 0.06315 1.2390
+#&gt; sigma 1.85700 0.4378 0.73200 2.9830
#&gt;
#&gt; Parameter correlation:
-#&gt; No covariance matrix
+#&gt; parent_0 log_alpha log_beta sigma
+#&gt; parent_0 1.000e+00 -1.565e-01 -3.142e-01 4.770e-08
+#&gt; log_alpha -1.565e-01 1.000e+00 9.564e-01 9.974e-08
+#&gt; log_beta -3.142e-01 9.564e-01 1.000e+00 8.468e-08
+#&gt; sigma 4.770e-08 9.974e-08 8.468e-08 1.000e+00
+#&gt;
#&gt; Backtransformed parameters:
#&gt; Confidence intervals for internally transformed parameters are asymmetric.
#&gt; t-test (unrealistically) based on the assumption of normal distribution
#&gt; for estimators of untransformed parameters.
-#&gt; Estimate t value Pr(&gt;t) Lower Upper
-#&gt; parent_0 85.870 NA NA NA NA
-#&gt; alpha 1.053 NA NA NA NA
-#&gt; beta 1.917 NA NA NA NA
-#&gt; sigma 1.857 NA NA NA NA
+#&gt; Estimate t value Pr(&gt;t) Lower Upper
+#&gt; parent_0 85.870 47.530 3.893e-08 81.2300 90.520
+#&gt; alpha 1.053 7.393 3.562e-04 0.7439 1.491
+#&gt; beta 1.917 4.373 3.601e-03 1.0650 3.451
+#&gt; sigma 1.857 4.243 4.074e-03 0.7320 2.983
#&gt;
#&gt; FOCUS Chi2 error levels in percent:
#&gt; err.min n.optim df
@@ -438,7 +443,7 @@ Per default, parameters in the kinetic models are internally transformed in
<span class='kw'>m1</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>))</div><div class='output co'>#&gt; <span class='message'>Successfully compiled differential equation model from auto-generated C code.</span></div><div class='input'><span class='co'># Fit the model to the FOCUS example dataset D using defaults</span>
<span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/print'>print</a></span>(<span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/system.time'>system.time</a></span>(<span class='no'>fit</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO</span>, <span class='no'>FOCUS_2006_D</span>,
<span class='kw'>solution_type</span> <span class='kw'>=</span> <span class='st'>"eigen"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)))</div><div class='output co'>#&gt; <span class='warning'>Warning: Observations with value of zero were removed from the data</span></div><div class='output co'>#&gt; User System verstrichen
-#&gt; 0.637 0.000 0.640 </div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/stats/topics/coef'>coef</a></span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; NULL</div><div class='input'><span class='fu'><a href='endpoints.html'>endpoints</a></span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; $ff
+#&gt; 2.251 0.000 2.253 </div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/stats/topics/coef'>coef</a></span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; NULL</div><div class='input'><span class='fu'><a href='endpoints.html'>endpoints</a></span>(<span class='no'>fit</span>)</div><div class='output co'>#&gt; $ff
#&gt; parent_sink parent_m1 m1_sink
#&gt; 0.485524 0.514476 1.000000
#&gt;
@@ -510,7 +515,7 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt; Sum of squared residuals at call 126: 371.2134
#&gt; Sum of squared residuals at call 135: 371.2134
#&gt; Negative log-likelihood at call 145: 97.22429</div><div class='output co'>#&gt; <span class='message'>Optimisation successfully terminated.</span></div><div class='output co'>#&gt; User System verstrichen
-#&gt; 0.544 0.000 0.550 </div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/stats/topics/coef'>coef</a></span>(<span class='no'>fit.deSolve</span>)</div><div class='output co'>#&gt; NULL</div><div class='input'><span class='fu'><a href='endpoints.html'>endpoints</a></span>(<span class='no'>fit.deSolve</span>)</div><div class='output co'>#&gt; $ff
+#&gt; 1.151 0.000 1.152 </div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/stats/topics/coef'>coef</a></span>(<span class='no'>fit.deSolve</span>)</div><div class='output co'>#&gt; NULL</div><div class='input'><span class='fu'><a href='endpoints.html'>endpoints</a></span>(<span class='no'>fit.deSolve</span>)</div><div class='output co'>#&gt; $ff
#&gt; parent_sink parent_m1 m1_sink
#&gt; 0.485524 0.514476 1.000000
#&gt;
@@ -540,10 +545,10 @@ Per default, parameters in the kinetic models are internally transformed in
<span class='no'>fit.SFORB_SFO</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFORB_SFO</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>parms.ini</span> <span class='kw'>=</span> <span class='no'>fit.SFORB</span>$<span class='no'>bparms.ode</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Observations with value of zero were removed from the data</span></div><div class='input'>
</div><div class='input'><span class='co'># Weighted fits, including IRLS</span>
<span class='no'>SFO_SFO.ff</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinmod.html'>mkinmod</a></span>(<span class='kw'>parent</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>, <span class='st'>"m1"</span>),
- <span class='kw'>m1</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>), <span class='kw'>use_of_ff</span> <span class='kw'>=</span> <span class='st'>"max"</span>)</div><div class='output co'>#&gt; <span class='message'>Successfully compiled differential equation model from auto-generated C code.</span></div><div class='input'><span class='no'>f.noweight</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Observations with value of zero were removed from the data</span></div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/summary'>summary</a></span>(<span class='no'>f.noweight</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Could not calculate correlation; no covariance matrix</span></div><div class='output co'>#&gt; mkin version used for fitting: 0.9.49.4
+ <span class='kw'>m1</span> <span class='kw'>=</span> <span class='fu'><a href='mkinsub.html'>mkinsub</a></span>(<span class='st'>"SFO"</span>), <span class='kw'>use_of_ff</span> <span class='kw'>=</span> <span class='st'>"max"</span>)</div><div class='output co'>#&gt; <span class='message'>Successfully compiled differential equation model from auto-generated C code.</span></div><div class='input'><span class='no'>f.noweight</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Observations with value of zero were removed from the data</span></div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/summary'>summary</a></span>(<span class='no'>f.noweight</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.49.4
#&gt; R version used for fitting: 3.6.0
-#&gt; Date of fit: Tue May 7 08:08:30 2019
-#&gt; Date of summary: Tue May 7 08:08:30 2019
+#&gt; Date of fit: Tue May 7 08:36:33 2019
+#&gt; Date of summary: Tue May 7 08:36:33 2019
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -551,7 +556,7 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted using 185 model solutions performed in 0.499 s
+#&gt; Fitted using 421 model solutions performed in 1.1 s
#&gt;
#&gt; Error model:
#&gt; Constant variance
@@ -577,25 +582,31 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt; m1_0 0 state
#&gt;
#&gt; Optimised, transformed parameters with symmetric confidence intervals:
-#&gt; Estimate Std. Error Lower Upper
-#&gt; parent_0 99.60000 NA NA NA
-#&gt; log_k_parent -2.31600 NA NA NA
-#&gt; log_k_m1 -5.24800 NA NA NA
-#&gt; f_parent_ilr_1 0.04096 NA NA NA
-#&gt; sigma 3.12600 NA NA NA
+#&gt; Estimate Std. Error Lower Upper
+#&gt; parent_0 99.60000 1.57000 96.40000 102.8000
+#&gt; log_k_parent -2.31600 0.04087 -2.39900 -2.2330
+#&gt; log_k_m1 -5.24800 0.13320 -5.51800 -4.9770
+#&gt; f_parent_ilr_1 0.04096 0.06312 -0.08746 0.1694
+#&gt; sigma 3.12600 0.35850 2.39600 3.8550
#&gt;
#&gt; Parameter correlation:
-#&gt; No covariance matrix
+#&gt; parent_0 log_k_parent log_k_m1 f_parent_ilr_1 sigma
+#&gt; parent_0 1.000e+00 5.174e-01 -1.688e-01 -5.471e-01 -2.265e-07
+#&gt; log_k_parent 5.174e-01 1.000e+00 -3.263e-01 -5.426e-01 3.785e-07
+#&gt; log_k_m1 -1.688e-01 -3.263e-01 1.000e+00 7.478e-01 -1.386e-07
+#&gt; f_parent_ilr_1 -5.471e-01 -5.426e-01 7.478e-01 1.000e+00 -3.641e-08
+#&gt; sigma -2.265e-07 3.785e-07 -1.386e-07 -3.641e-08 1.000e+00
+#&gt;
#&gt; Backtransformed parameters:
#&gt; Confidence intervals for internally transformed parameters are asymmetric.
#&gt; t-test (unrealistically) based on the assumption of normal distribution
#&gt; for estimators of untransformed parameters.
-#&gt; Estimate t value Pr(&gt;t) Lower Upper
-#&gt; parent_0 99.600000 NA NA NA NA
-#&gt; k_parent 0.098700 NA NA NA NA
-#&gt; k_m1 0.005261 NA NA NA NA
-#&gt; f_parent_to_m1 0.514500 NA NA NA NA
-#&gt; sigma 3.126000 NA NA NA NA
+#&gt; Estimate t value Pr(&gt;t) Lower Upper
+#&gt; parent_0 99.600000 63.430 2.298e-36 96.400000 1.028e+02
+#&gt; k_parent 0.098700 24.470 4.955e-23 0.090820 1.073e-01
+#&gt; k_m1 0.005261 7.510 6.165e-09 0.004012 6.898e-03
+#&gt; f_parent_to_m1 0.514500 23.070 3.104e-22 0.469100 5.596e-01
+#&gt; sigma 3.126000 8.718 2.235e-10 2.396000 3.855e+00
#&gt;
#&gt; FOCUS Chi2 error levels in percent:
#&gt; err.min n.optim df
@@ -652,10 +663,10 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt; 100 m1 31.04 31.98163 -9.416e-01
#&gt; 100 m1 33.13 31.98163 1.148e+00
#&gt; 120 m1 25.15 28.78984 -3.640e+00
-#&gt; 120 m1 33.31 28.78984 4.520e+00</div><div class='input'><span class='no'>f.obs</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>error_model</span> <span class='kw'>=</span> <span class='st'>"obs"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Observations with value of zero were removed from the data</span></div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/summary'>summary</a></span>(<span class='no'>f.obs</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Could not calculate correlation; no covariance matrix</span></div><div class='output co'>#&gt; mkin version used for fitting: 0.9.49.4
+#&gt; 120 m1 33.31 28.78984 4.520e+00</div><div class='input'><span class='no'>f.obs</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>error_model</span> <span class='kw'>=</span> <span class='st'>"obs"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Observations with value of zero were removed from the data</span></div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/summary'>summary</a></span>(<span class='no'>f.obs</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.49.4
#&gt; R version used for fitting: 3.6.0
-#&gt; Date of fit: Tue May 7 08:08:32 2019
-#&gt; Date of summary: Tue May 7 08:08:32 2019
+#&gt; Date of fit: Tue May 7 08:36:35 2019
+#&gt; Date of summary: Tue May 7 08:36:35 2019
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -663,7 +674,7 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted using 426 model solutions performed in 1.139 s
+#&gt; Fitted using 758 model solutions performed in 1.991 s
#&gt;
#&gt; Error model:
#&gt; Variance unique to each observed variable
@@ -691,27 +702,41 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt; m1_0 0 state
#&gt;
#&gt; Optimised, transformed parameters with symmetric confidence intervals:
-#&gt; Estimate Std. Error Lower Upper
-#&gt; parent_0 99.65000 NA NA NA
-#&gt; log_k_parent -2.31300 NA NA NA
-#&gt; log_k_m1 -5.25000 NA NA NA
-#&gt; f_parent_ilr_1 0.03861 NA NA NA
-#&gt; sigma_parent 3.40100 NA NA NA
-#&gt; sigma_m1 2.85500 NA NA NA
+#&gt; Estimate Std. Error Lower Upper
+#&gt; parent_0 99.65000 1.70200 96.19000 103.1000
+#&gt; log_k_parent -2.31300 0.04376 -2.40200 -2.2240
+#&gt; log_k_m1 -5.25000 0.12430 -5.50400 -4.9970
+#&gt; f_parent_ilr_1 0.03861 0.06171 -0.08708 0.1643
+#&gt; sigma_parent 3.40100 0.56820 2.24400 4.5590
+#&gt; sigma_m1 2.85500 0.45240 1.93400 3.7770
#&gt;
#&gt; Parameter correlation:
-#&gt; No covariance matrix
+#&gt; parent_0 log_k_parent log_k_m1 f_parent_ilr_1 sigma_parent
+#&gt; parent_0 1.00000 0.51078 -0.19133 -0.59997 0.035670
+#&gt; log_k_parent 0.51078 1.00000 -0.37458 -0.59239 0.069833
+#&gt; log_k_m1 -0.19133 -0.37458 1.00000 0.74398 -0.026158
+#&gt; f_parent_ilr_1 -0.59997 -0.59239 0.74398 1.00000 -0.041368
+#&gt; sigma_parent 0.03567 0.06983 -0.02616 -0.04137 1.000000
+#&gt; sigma_m1 -0.03385 -0.06627 0.02482 0.03925 -0.004628
+#&gt; sigma_m1
+#&gt; parent_0 -0.033847
+#&gt; log_k_parent -0.066265
+#&gt; log_k_m1 0.024821
+#&gt; f_parent_ilr_1 0.039255
+#&gt; sigma_parent -0.004628
+#&gt; sigma_m1 1.000000
+#&gt;
#&gt; Backtransformed parameters:
#&gt; Confidence intervals for internally transformed parameters are asymmetric.
#&gt; t-test (unrealistically) based on the assumption of normal distribution
#&gt; for estimators of untransformed parameters.
-#&gt; Estimate t value Pr(&gt;t) Lower Upper
-#&gt; parent_0 99.650000 NA NA NA NA
-#&gt; k_parent 0.098970 NA NA NA NA
-#&gt; k_m1 0.005245 NA NA NA NA
-#&gt; f_parent_to_m1 0.513600 NA NA NA NA
-#&gt; sigma_parent 3.401000 NA NA NA NA
-#&gt; sigma_m1 2.855000 NA NA NA NA
+#&gt; Estimate t value Pr(&gt;t) Lower Upper
+#&gt; parent_0 99.650000 58.560 2.004e-34 96.190000 1.031e+02
+#&gt; k_parent 0.098970 22.850 1.099e-21 0.090530 1.082e-01
+#&gt; k_m1 0.005245 8.046 1.732e-09 0.004072 6.756e-03
+#&gt; f_parent_to_m1 0.513600 23.560 4.352e-22 0.469300 5.578e-01
+#&gt; sigma_parent 3.401000 5.985 5.662e-07 2.244000 4.559e+00
+#&gt; sigma_m1 2.855000 6.311 2.215e-07 1.934000 3.777e+00
#&gt;
#&gt; FOCUS Chi2 error levels in percent:
#&gt; err.min n.optim df
@@ -768,10 +793,10 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt; 100 m1 31.04 31.98774 -9.477e-01
#&gt; 100 m1 33.13 31.98774 1.142e+00
#&gt; 120 m1 25.15 28.80430 -3.654e+00
-#&gt; 120 m1 33.31 28.80430 4.506e+00</div><div class='input'><span class='no'>f.tc</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>error_model</span> <span class='kw'>=</span> <span class='st'>"tc"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Observations with value of zero were removed from the data</span></div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/summary'>summary</a></span>(<span class='no'>f.tc</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Could not calculate correlation; no covariance matrix</span></div><div class='output co'>#&gt; mkin version used for fitting: 0.9.49.4
+#&gt; 120 m1 33.31 28.80430 4.506e+00</div><div class='input'><span class='no'>f.tc</span> <span class='kw'>&lt;-</span> <span class='fu'>mkinfit</span>(<span class='no'>SFO_SFO.ff</span>, <span class='no'>FOCUS_2006_D</span>, <span class='kw'>error_model</span> <span class='kw'>=</span> <span class='st'>"tc"</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#&gt; <span class='warning'>Warning: Observations with value of zero were removed from the data</span></div><div class='input'><span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/summary'>summary</a></span>(<span class='no'>f.tc</span>)</div><div class='output co'>#&gt; mkin version used for fitting: 0.9.49.4
#&gt; R version used for fitting: 3.6.0
-#&gt; Date of fit: Tue May 7 08:08:34 2019
-#&gt; Date of summary: Tue May 7 08:08:34 2019
+#&gt; Date of fit: Tue May 7 08:36:39 2019
+#&gt; Date of summary: Tue May 7 08:36:39 2019
#&gt;
#&gt; Equations:
#&gt; d_parent/dt = - k_parent * parent
@@ -779,7 +804,7 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt;
#&gt; Model predictions using solution type deSolve
#&gt;
-#&gt; Fitted using 489 model solutions performed in 2.013 s
+#&gt; Fitted using 821 model solutions performed in 3.304 s
#&gt;
#&gt; Error model:
#&gt; Two-component variance function
@@ -807,27 +832,34 @@ Per default, parameters in the kinetic models are internally transformed in
#&gt; m1_0 0 state
#&gt;
#&gt; Optimised, transformed parameters with symmetric confidence intervals:
-#&gt; Estimate Std. Error Lower Upper
-#&gt; parent_0 100.70000 NA NA NA
-#&gt; log_k_parent -2.29700 NA NA NA
-#&gt; log_k_m1 -5.26600 NA NA NA
-#&gt; f_parent_ilr_1 0.02374 NA NA NA
-#&gt; sigma_low 0.00305 NA NA NA
-#&gt; rsd_high 0.07928 NA NA NA
+#&gt; Estimate Std. Error Lower Upper
+#&gt; parent_0 100.70000 2.621000 95.400000 106.10000
+#&gt; log_k_parent -2.29700 0.008862 -2.315000 -2.27900
+#&gt; log_k_m1 -5.26600 0.091310 -5.452000 -5.08000
+#&gt; f_parent_ilr_1 0.02374 0.055300 -0.088900 0.13640
+#&gt; sigma_low 0.00305 0.004829 -0.006786 0.01289
+#&gt; rsd_high 0.07928 0.009418 0.060100 0.09847
#&gt;
#&gt; Parameter correlation:
-#&gt; No covariance matrix
+#&gt; parent_0 log_k_parent log_k_m1 f_parent_ilr_1 sigma_low rsd_high
+#&gt; parent_0 1.00000 0.67644 -0.10215 -0.76822 0.14294 -0.08783
+#&gt; log_k_parent 0.67644 1.00000 -0.15102 -0.59491 0.34611 -0.08125
+#&gt; log_k_m1 -0.10215 -0.15102 1.00000 0.51808 -0.05236 0.01240
+#&gt; f_parent_ilr_1 -0.76822 -0.59491 0.51808 1.00000 -0.13900 0.03248
+#&gt; sigma_low 0.14294 0.34611 -0.05236 -0.13900 1.00000 -0.16546
+#&gt; rsd_high -0.08783 -0.08125 0.01240 0.03248 -0.16546 1.00000
+#&gt;
#&gt; Backtransformed parameters:
#&gt; Confidence intervals for internally transformed parameters are asymmetric.
#&gt; t-test (unrealistically) based on the assumption of normal distribution
#&gt; for estimators of untransformed parameters.
-#&gt; Estimate t value Pr(&gt;t) Lower Upper
-#&gt; parent_0 1.007e+02 NA NA NA NA
-#&gt; k_parent 1.006e-01 NA NA NA NA
-#&gt; k_m1 5.167e-03 NA NA NA NA
-#&gt; f_parent_to_m1 5.084e-01 NA NA NA NA
-#&gt; sigma_low 3.050e-03 NA NA NA NA
-#&gt; rsd_high 7.928e-02 NA NA NA NA
+#&gt; Estimate t value Pr(&gt;t) Lower Upper
+#&gt; parent_0 1.007e+02 38.4300 1.180e-28 95.400000 1.061e+02
+#&gt; k_parent 1.006e-01 112.8000 1.718e-43 0.098760 1.024e-01
+#&gt; k_m1 5.167e-03 10.9500 1.172e-12 0.004290 6.223e-03
+#&gt; f_parent_to_m1 5.084e-01 26.0100 2.146e-23 0.468600 5.481e-01
+#&gt; sigma_low 3.050e-03 0.6314 2.661e-01 -0.006786 1.289e-02
+#&gt; rsd_high 7.928e-02 8.4170 6.418e-10 0.060100 9.847e-02
#&gt;
#&gt; FOCUS Chi2 error levels in percent:
#&gt; err.min n.optim df

Contact - Imprint