diff options
Diffstat (limited to 'vignettes/FOCUS_L.html')
| -rw-r--r-- | vignettes/FOCUS_L.html | 129 |
1 files changed, 64 insertions, 65 deletions
diff --git a/vignettes/FOCUS_L.html b/vignettes/FOCUS_L.html index a06f360e..96ea70ce 100644 --- a/vignettes/FOCUS_L.html +++ b/vignettes/FOCUS_L.html @@ -214,7 +214,13 @@ hr { report, p. 284:</p> <pre><code class="r">library("mkin") -FOCUS_2006_L1 = data.frame( +</code></pre> + +<pre><code>## Loading required package: minpack.lm +## Loading required package: rootSolve +</code></pre> + +<pre><code class="r">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, @@ -237,16 +243,16 @@ summary(m.L1.SFO) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:55 2015 -## Date of summary: Fri May 15 13:29:55 2015 +## 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 ## ## Equations: ## d_parent = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 37 model solutions performed in 0.082 s +## Fitted with method Port using 37 model solutions performed in 0.098 s ## ## Weighting: none ## @@ -335,31 +341,20 @@ 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) -</code></pre> - -<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) +summary(m.L1.FOMC, data = FALSE) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:55 2015 -## Date of summary: Fri May 15 13:29:55 2015 -## -## -## Warning: Optimisation by method Port did not converge. -## Convergence code is 1 -## +## 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 ## ## Equations: ## d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 188 model solutions performed in 0.412 s +## Fitted with method Port using 611 model solutions performed in 1.509 s ## ## Weighting: none ## @@ -379,24 +374,28 @@ is checked.</p> ## None ## ## Optimised, transformed parameters: -## Estimate Std. Error Lower Upper t value Pr(>|t|) Pr(>t) -## parent_0 92.47 1.422 89.44 95.50 65.030 8.317e-20 4.158e-20 -## log_alpha 15.43 15.080 -16.71 47.58 1.023 3.224e-01 1.612e-01 -## log_beta 17.78 15.090 -14.37 49.93 1.179 2.569e-01 1.284e-01 +## 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 ## ## Parameter correlation: ## parent_0 log_alpha log_beta -## parent_0 1.0000 0.1129 0.1112 -## log_alpha 0.1129 1.0000 1.0000 -## log_beta 0.1112 1.0000 1.0000 +## parent_0 1.0000 -0.3016 -0.3016 +## log_alpha -0.3016 1.0000 1.0000 +## log_beta -0.3016 1.0000 1.0000 ## ## Residual standard error: 3.045 on 15 degrees of freedom ## ## Backtransformed parameters: -## Estimate Lower Upper -## parent_0 9.247e+01 8.944e+01 9.550e+01 -## alpha 5.044e+06 5.510e-08 4.618e+20 -## beta 5.276e+07 5.732e-07 4.857e+21 +## Estimate Lower Upper +## parent_0 92.47 89.31 95.63 +## alpha 76830.00 0.00 Inf +## beta 803500.00 0.00 Inf ## ## Chi2 error levels in percent: ## err.min n.optim df @@ -404,8 +403,8 @@ is checked.</p> ## parent 3.619 3 6 ## ## Estimated disappearance times: -## DT50 DT90 DT50back -## parent 7.25 24.08 7.25 +## DT50 DT90 DT50back +## parent 7.249 24.08 7.249 </code></pre> <p>Due to the higher number of parameters, and the lower number of degrees of @@ -444,16 +443,16 @@ summary(m.L2.SFO) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:56 2015 -## Date of summary: Fri May 15 13:29:56 2015 +## 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 ## ## Equations: ## d_parent = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 41 model solutions performed in 0.091 s +## Fitted with method Port using 41 model solutions performed in 0.1 s ## ## Weighting: none ## @@ -554,16 +553,16 @@ mkinresplot(m.L2.FOMC) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:56 2015 -## Date of summary: Fri May 15 13:29:56 2015 +## 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 ## ## Equations: ## d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 81 model solutions performed in 0.192 s +## Fitted with method Port using 81 model solutions performed in 0.201 s ## ## Weighting: none ## @@ -640,9 +639,9 @@ plot(m.L2.DFOP) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:58 2015 -## Date of summary: Fri May 15 13:29:58 2015 +## 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 ## ## Equations: ## d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -651,7 +650,7 @@ plot(m.L2.DFOP) ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 336 model solutions performed in 0.775 s +## Fitted with method Port using 336 model solutions performed in 0.856 s ## ## Weighting: none ## @@ -729,16 +728,16 @@ plot(m.L3.SFO) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:58 2015 -## Date of summary: Fri May 15 13:29:58 2015 +## 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 ## ## Equations: ## d_parent = - k_parent_sink * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 43 model solutions performed in 0.099 s +## Fitted with method Port using 43 model solutions performed in 0.109 s ## ## Weighting: none ## @@ -815,16 +814,16 @@ plot(m.L3.FOMC) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:58 2015 -## Date of summary: Fri May 15 13:29:58 2015 +## 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 ## ## Equations: ## d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 83 model solutions performed in 0.186 s +## Fitted with method Port using 83 model solutions performed in 0.203 s ## ## Weighting: none ## @@ -888,9 +887,9 @@ plot(m.L3.DFOP) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:59 2015 -## Date of summary: Fri May 15 13:29:59 2015 +## 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 ## ## Equations: ## d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -899,7 +898,7 @@ plot(m.L3.DFOP) ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 137 model solutions performed in 0.314 s +## Fitted with method Port using 137 model solutions performed in 0.346 s ## ## Weighting: none ## @@ -986,9 +985,9 @@ plot(m.L4.SFO) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:59 2015 -## Date of summary: Fri May 15 13:29:59 2015 +## 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 ## ## Equations: ## d_parent = - k_parent_sink * parent @@ -1061,16 +1060,16 @@ plot(m.L4.FOMC) </code></pre> <pre><code>## mkin version: 0.9.35 -## R version: 3.2.0 -## Date of fit: Fri May 15 13:29:59 2015 -## Date of summary: Fri May 15 13:29:59 2015 +## 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 ## ## Equations: ## d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent ## ## Model predictions using solution type analytical ## -## Fitted with method Port using 66 model solutions performed in 0.146 s +## Fitted with method Port using 66 model solutions performed in 0.161 s ## ## Weighting: none ## |