diff options
Diffstat (limited to 'vignettes/FOCUS_L.html')
| -rw-r--r-- | vignettes/FOCUS_L.html | 357 |
1 files changed, 188 insertions, 169 deletions
diff --git a/vignettes/FOCUS_L.html b/vignettes/FOCUS_L.html index 96ea70ce..692caf93 100644 --- a/vignettes/FOCUS_L.html +++ b/vignettes/FOCUS_L.html @@ -214,13 +214,7 @@ hr { report, p. 284:</p> <pre><code class="r">library("mkin") -</code></pre> - -<pre><code>## Loading required package: minpack.lm -## Loading required package: rootSolve -</code></pre> - -<pre><code class="r">FOCUS_2006_L1 = data.frame( +FOCUS_2006_L1 = data.frame( t = rep(c(0, 1, 2, 3, 5, 7, 14, 21, 30), each = 2), parent = c(88.3, 91.4, 85.6, 84.5, 78.9, 77.6, 72.0, 71.9, 50.3, 59.4, 47.0, 45.1, @@ -242,17 +236,17 @@ given in the FOCUS report. </p> summary(m.L1.SFO) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:53 2015 -## Date of summary: Sat Feb 21 14:44:53 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:47:59 2015 +## Date of summary: Sun Jun 21 01:47:59 2015 ## ## Equations: ## d_parent = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 37 model solutions performed in 0.098 s +## Fitted with method Port using 37 model solutions performed in 0.093 s ## ## Weighting: none ## @@ -269,13 +263,10 @@ summary(m.L1.SFO) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) -## parent_0 92.470 1.36800 89.570 95.370 67.58 4.339e-21 -## log_k_parent_sink -2.347 0.04057 -2.433 -2.261 -57.86 5.155e-20 -## Pr(>t) -## parent_0 2.170e-21 -## log_k_parent_sink 2.577e-20 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 92.470 1.36800 89.570 95.370 +## log_k_parent_sink -2.347 0.04057 -2.433 -2.261 ## ## Parameter correlation: ## parent_0 log_k_parent_sink @@ -285,9 +276,12 @@ summary(m.L1.SFO) ## Residual standard error: 2.948 on 16 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 92.47000 89.57000 95.3700 -## k_parent_sink 0.09561 0.08773 0.1042 +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 92.47000 67.58 2.170e-21 89.57000 95.3700 +## k_parent_sink 0.09561 24.65 1.867e-14 0.08773 0.1042 ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -341,20 +335,31 @@ The residual plot can be easily obtained by</p> is checked.</p> <pre><code class="r">m.L1.FOMC <- mkinfit("FOMC", FOCUS_2006_L1_mkin, quiet=TRUE) -summary(m.L1.FOMC, data = FALSE) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:55 2015 -## Date of summary: Sat Feb 21 14:44:55 2015 +<pre><code>## Warning in mkinfit("FOMC", FOCUS_2006_L1_mkin, quiet = TRUE): Optimisation by method Port did not converge. +## Convergence code is 1 +</code></pre> + +<pre><code class="r">summary(m.L1.FOMC, data = FALSE) +</code></pre> + +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:00 2015 +## Date of summary: Sun Jun 21 01:48:00 2015 +## +## +## Warning: Optimisation by method Port did not converge. +## Convergence code is 1 +## ## ## Equations: -## d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent +## d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 611 model solutions performed in 1.509 s +## Fitted with method Port using 188 model solutions performed in 0.463 s ## ## Weighting: none ## @@ -373,29 +378,28 @@ summary(m.L1.FOMC, data = FALSE) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) -## parent_0 92.47 1.482 89.31 95.63 62.39000 1.546e-19 -## log_alpha 11.25 598.200 -1264.00 1286.00 0.01880 9.852e-01 -## log_beta 13.60 598.200 -1261.00 1289.00 0.02273 9.822e-01 -## Pr(>t) -## parent_0 7.730e-20 -## log_alpha 4.926e-01 -## log_beta 4.911e-01 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 92.47 1.422 89.44 95.50 +## log_alpha 15.43 15.080 -16.71 47.58 +## log_beta 17.78 15.090 -14.37 49.93 ## ## Parameter correlation: ## parent_0 log_alpha log_beta -## parent_0 1.0000 -0.3016 -0.3016 -## log_alpha -0.3016 1.0000 1.0000 -## log_beta -0.3016 1.0000 1.0000 +## parent_0 1.0000 0.1129 0.1112 +## log_alpha 0.1129 1.0000 1.0000 +## log_beta 0.1112 1.0000 1.0000 ## ## Residual standard error: 3.045 on 15 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 92.47 89.31 95.63 -## alpha 76830.00 0.00 Inf -## beta 803500.00 0.00 Inf +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 9.247e+01 65.150 4.044e-20 8.944e+01 9.550e+01 +## alpha 5.044e+06 1.271 1.115e-01 5.510e-08 4.618e+20 +## beta 5.276e+07 1.259 1.137e-01 5.732e-07 4.857e+21 ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -403,8 +407,8 @@ summary(m.L1.FOMC, data = FALSE) ## parent 3.619 3 6 ## ## Estimated disappearance times: -## DT50 DT90 DT50back -## parent 7.249 24.08 7.249 +## DT50 DT90 DT50back +## parent 7.25 24.08 7.25 </code></pre> <p>Due to the higher number of parameters, and the lower number of degrees of @@ -442,17 +446,17 @@ FOCUS_2006_L2_mkin <- mkin_wide_to_long(FOCUS_2006_L2) summary(m.L2.SFO) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:55 2015 -## Date of summary: Sat Feb 21 14:44:55 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:00 2015 +## Date of summary: Sun Jun 21 01:48:00 2015 ## ## Equations: ## d_parent = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 41 model solutions performed in 0.1 s +## Fitted with method Port using 41 model solutions performed in 0.097 s ## ## Weighting: none ## @@ -469,13 +473,10 @@ summary(m.L2.SFO) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) -## parent_0 91.4700 3.8070 82.9800 99.9500 24.030 3.545e-10 -## log_k_parent_sink -0.4112 0.1074 -0.6505 -0.1719 -3.828 3.329e-03 -## Pr(>t) -## parent_0 1.773e-10 -## log_k_parent_sink 1.664e-03 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 91.4700 3.8070 82.9800 99.9500 +## log_k_parent_sink -0.4112 0.1074 -0.6505 -0.1719 ## ## Parameter correlation: ## parent_0 log_k_parent_sink @@ -485,9 +486,12 @@ summary(m.L2.SFO) ## Residual standard error: 5.51 on 10 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 91.4700 82.9800 99.9500 -## k_parent_sink 0.6629 0.5218 0.8421 +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 91.4700 24.03 1.773e-10 82.9800 99.9500 +## k_parent_sink 0.6629 9.31 1.525e-06 0.5218 0.8421 ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -552,17 +556,17 @@ mkinresplot(m.L2.FOMC) <pre><code class="r">summary(m.L2.FOMC, data = FALSE) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:55 2015 -## Date of summary: Sat Feb 21 14:44:55 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:00 2015 +## Date of summary: Sun Jun 21 01:48:00 2015 ## ## Equations: -## d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent +## d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 81 model solutions performed in 0.201 s +## Fitted with method Port using 81 model solutions performed in 0.191 s ## ## Weighting: none ## @@ -581,11 +585,11 @@ mkinresplot(m.L2.FOMC) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) Pr(>t) -## parent_0 93.7700 1.8560 89.5700 97.9700 50.5100 2.345e-12 1.173e-12 -## log_alpha 0.3180 0.1867 -0.1044 0.7405 1.7030 1.227e-01 6.137e-02 -## log_beta 0.2102 0.2943 -0.4555 0.8759 0.7142 4.932e-01 2.466e-01 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 93.7700 1.8560 89.5700 97.9700 +## log_alpha 0.3180 0.1867 -0.1044 0.7405 +## log_beta 0.2102 0.2943 -0.4555 0.8759 ## ## Parameter correlation: ## parent_0 log_alpha log_beta @@ -596,10 +600,13 @@ mkinresplot(m.L2.FOMC) ## Residual standard error: 2.628 on 9 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 93.770 89.5700 97.970 -## alpha 1.374 0.9009 2.097 -## beta 1.234 0.6341 2.401 +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 93.770 50.510 1.173e-12 89.5700 97.970 +## alpha 1.374 5.355 2.296e-04 0.9009 2.097 +## beta 1.234 3.398 3.949e-03 0.6341 2.401 ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -638,10 +645,10 @@ plot(m.L2.DFOP) <pre><code class="r">summary(m.L2.DFOP, data = FALSE) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:57 2015 -## Date of summary: Sat Feb 21 14:44:57 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:02 2015 +## Date of summary: Sun Jun 21 01:48:02 2015 ## ## Equations: ## d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -650,7 +657,7 @@ plot(m.L2.DFOP) ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 336 model solutions performed in 0.856 s +## Fitted with method Port using 336 model solutions performed in 0.835 s ## ## Weighting: none ## @@ -671,12 +678,12 @@ plot(m.L2.DFOP) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) Pr(>t) -## parent_0 93.9500 NA NA NA NA NA NA -## log_k1 3.1210 NA NA NA NA NA NA -## log_k2 -1.0880 NA NA NA NA NA NA -## g_ilr -0.2821 NA NA NA NA NA NA +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 93.9500 NA NA NA +## log_k1 3.1210 NA NA NA +## log_k2 -1.0880 NA NA NA +## g_ilr -0.2821 NA NA NA ## ## Parameter correlation: ## Could not estimate covariance matrix; singular system: @@ -684,11 +691,14 @@ plot(m.L2.DFOP) ## Residual standard error: 1.732 on 8 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 93.9500 NA NA -## k1 22.6700 NA NA -## k2 0.3369 NA NA -## g 0.4016 NA NA +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 93.9500 NA NA NA NA +## k1 22.6700 NA NA NA NA +## k2 0.3369 NA NA NA NA +## g 0.4016 NA NA NA NA ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -727,17 +737,17 @@ plot(m.L3.SFO) <pre><code class="r">summary(m.L3.SFO) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:57 2015 -## Date of summary: Sat Feb 21 14:44:57 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:03 2015 +## Date of summary: Sun Jun 21 01:48:03 2015 ## ## Equations: ## d_parent = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 43 model solutions performed in 0.109 s +## Fitted with method Port using 43 model solutions performed in 0.104 s ## ## Weighting: none ## @@ -754,13 +764,10 @@ plot(m.L3.SFO) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) -## parent_0 74.870 8.4570 54.180 95.57 8.853 1.155e-04 -## log_k_parent_sink -3.678 0.3261 -4.476 -2.88 -11.280 2.903e-05 -## Pr(>t) -## parent_0 5.776e-05 -## log_k_parent_sink 1.451e-05 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 74.870 8.4570 54.180 95.57 +## log_k_parent_sink -3.678 0.3261 -4.476 -2.88 ## ## Parameter correlation: ## parent_0 log_k_parent_sink @@ -770,9 +777,12 @@ plot(m.L3.SFO) ## Residual standard error: 12.91 on 6 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 74.87000 54.18000 95.57000 -## k_parent_sink 0.02527 0.01138 0.05612 +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 74.87000 8.853 5.776e-05 54.18000 95.57000 +## k_parent_sink 0.02527 3.067 1.102e-02 0.01138 0.05612 ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -813,17 +823,17 @@ plot(m.L3.FOMC) <pre><code class="r">summary(m.L3.FOMC, data = FALSE) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:58 2015 -## Date of summary: Sat Feb 21 14:44:58 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:03 2015 +## Date of summary: Sun Jun 21 01:48:03 2015 ## ## Equations: -## d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent +## d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 83 model solutions performed in 0.203 s +## Fitted with method Port using 83 model solutions performed in 0.196 s ## ## Weighting: none ## @@ -842,11 +852,11 @@ plot(m.L3.FOMC) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) Pr(>t) -## parent_0 96.9700 4.5500 85.2800 108.7000 21.310 4.216e-06 2.108e-06 -## log_alpha -0.8619 0.1704 -1.3000 -0.4238 -5.057 3.911e-03 1.955e-03 -## log_beta 0.6193 0.4744 -0.6003 1.8390 1.305 2.486e-01 1.243e-01 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 96.9700 4.5500 85.2800 108.7000 +## log_alpha -0.8619 0.1704 -1.3000 -0.4238 +## log_beta 0.6193 0.4744 -0.6003 1.8390 ## ## Parameter correlation: ## parent_0 log_alpha log_beta @@ -857,10 +867,13 @@ plot(m.L3.FOMC) ## Residual standard error: 4.572 on 5 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 96.9700 85.2800 108.7000 -## alpha 0.4224 0.2725 0.6546 -## beta 1.8580 0.5487 6.2890 +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 96.9700 21.310 2.108e-06 85.2800 108.7000 +## alpha 0.4224 5.867 1.020e-03 0.2725 0.6546 +## beta 1.8580 2.108 4.444e-02 0.5487 6.2890 ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -886,10 +899,10 @@ plot(m.L3.DFOP) <pre><code class="r">summary(m.L3.DFOP, data = FALSE) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:58 2015 -## Date of summary: Sat Feb 21 14:44:58 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:03 2015 +## Date of summary: Sun Jun 21 01:48:03 2015 ## ## Equations: ## d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -898,7 +911,7 @@ plot(m.L3.DFOP) ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 137 model solutions performed in 0.346 s +## Fitted with method Port using 137 model solutions performed in 0.35 s ## ## Weighting: none ## @@ -919,12 +932,12 @@ plot(m.L3.DFOP) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) Pr(>t) -## parent_0 97.7500 1.43800 93.7500 101.70000 67.970 2.808e-07 1.404e-07 -## log_k1 -0.6612 0.13340 -1.0310 -0.29100 -4.958 7.715e-03 3.858e-03 -## log_k2 -4.2860 0.05902 -4.4500 -4.12200 -72.620 2.155e-07 1.077e-07 -## g_ilr -0.1229 0.05121 -0.2651 0.01925 -2.401 7.431e-02 3.716e-02 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 97.7500 1.43800 93.7500 101.70000 +## log_k1 -0.6612 0.13340 -1.0310 -0.29100 +## log_k2 -4.2860 0.05902 -4.4500 -4.12200 +## g_ilr -0.1229 0.05121 -0.2651 0.01925 ## ## Parameter correlation: ## parent_0 log_k1 log_k2 g_ilr @@ -936,11 +949,14 @@ plot(m.L3.DFOP) ## Residual standard error: 1.439 on 4 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 97.75000 93.75000 101.70000 -## k1 0.51620 0.35650 0.74750 -## k2 0.01376 0.01168 0.01621 -## g 0.45660 0.40730 0.50680 +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 97.75000 67.970 1.404e-07 93.75000 101.70000 +## k1 0.51620 7.499 8.460e-04 0.35650 0.74750 +## k2 0.01376 16.940 3.557e-05 0.01168 0.01621 +## g 0.45660 25.410 7.121e-06 0.40730 0.50680 ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -984,17 +1000,17 @@ plot(m.L4.SFO) <pre><code class="r">summary(m.L4.SFO, data = FALSE) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:58 2015 -## Date of summary: Sat Feb 21 14:44:58 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:04 2015 +## Date of summary: Sun Jun 21 01:48:04 2015 ## ## Equations: ## d_parent = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 46 model solutions performed in 0.109 s +## Fitted with method Port using 46 model solutions performed in 0.105 s ## ## Weighting: none ## @@ -1011,13 +1027,10 @@ plot(m.L4.SFO) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) -## parent_0 96.44 1.94900 91.670 101.200 49.49 4.566e-09 -## log_k_parent_sink -5.03 0.07999 -5.225 -4.834 -62.88 1.088e-09 -## Pr(>t) -## parent_0 2.283e-09 -## log_k_parent_sink 5.438e-10 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 96.44 1.94900 91.670 101.200 +## log_k_parent_sink -5.03 0.07999 -5.225 -4.834 ## ## Parameter correlation: ## parent_0 log_k_parent_sink @@ -1027,9 +1040,12 @@ plot(m.L4.SFO) ## Residual standard error: 3.651 on 6 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 96.440000 91.670000 1.012e+02 -## k_parent_sink 0.006541 0.005378 7.955e-03 +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 96.440000 49.49 2.283e-09 91.670000 1.012e+02 +## k_parent_sink 0.006541 12.50 8.008e-06 0.005378 7.955e-03 ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -1059,17 +1075,17 @@ plot(m.L4.FOMC) <pre><code class="r">summary(m.L4.FOMC, data = FALSE) </code></pre> -<pre><code>## mkin version: 0.9.35 -## R version: 3.1.2 -## Date of fit: Sat Feb 21 14:44:58 2015 -## Date of summary: Sat Feb 21 14:44:58 2015 +<pre><code>## mkin version: 0.9.36 +## R version: 3.2.1 +## Date of fit: Sun Jun 21 01:48:04 2015 +## Date of summary: Sun Jun 21 01:48:04 2015 ## ## Equations: -## d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent +## d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 66 model solutions performed in 0.161 s +## Fitted with method Port using 66 model solutions performed in 0.159 s ## ## Weighting: none ## @@ -1088,11 +1104,11 @@ plot(m.L4.FOMC) ## Fixed parameter values: ## None ## -## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) Pr(>t) -## parent_0 99.1400 1.6800 94.820 103.5000 59.0200 2.643e-08 1.322e-08 -## log_alpha -0.3506 0.3725 -1.308 0.6068 -0.9414 3.897e-01 1.949e-01 -## log_beta 4.1740 0.5635 2.726 5.6230 7.4070 7.059e-04 3.530e-04 +## Optimised, transformed parameters with symmetric confidence intervals: +## Estimate Std. Error Lower Upper +## parent_0 99.1400 1.6800 94.820 103.5000 +## log_alpha -0.3506 0.3725 -1.308 0.6068 +## log_beta 4.1740 0.5635 2.726 5.6230 ## ## Parameter correlation: ## parent_0 log_alpha log_beta @@ -1103,10 +1119,13 @@ plot(m.L4.FOMC) ## Residual standard error: 2.315 on 5 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 99.1400 94.8200 103.500 -## alpha 0.7042 0.2703 1.835 -## beta 64.9800 15.2600 276.600 +## Confidence intervals for internally transformed parameters are asymmetric. +## t-test (unrealistically) based on the assumption of normal distribution +## for estimators of untransformed parameters. +## Estimate t value Pr(>t) Lower Upper +## parent_0 99.1400 59.020 1.322e-08 94.8200 103.500 +## alpha 0.7042 2.685 2.178e-02 0.2703 1.835 +## beta 64.9800 1.775 6.807e-02 15.2600 276.600 ## ## Chi2 error levels in percent: ## err.min n.optim df |