From d89e3d22eb9dc383897b09e9c5aa1b57f65cdbf0 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Thu, 21 Feb 2019 14:34:45 +0100 Subject: Add the logistic model --- docs/articles/FOCUS_L.html | 58 +++++++++++++++++++++++----------------------- 1 file changed, 29 insertions(+), 29 deletions(-) (limited to 'docs/articles/FOCUS_L.html') diff --git a/docs/articles/FOCUS_L.html b/docs/articles/FOCUS_L.html index e3ca9aff..265dd96c 100644 --- a/docs/articles/FOCUS_L.html +++ b/docs/articles/FOCUS_L.html @@ -85,7 +85,7 @@

Example evaluation of FOCUS Laboratory Data L1 to L3

Johannes Ranke

-

2019-01-31

+

2019-02-21

@@ -109,17 +109,17 @@

Since mkin version 0.9-32 (July 2014), we can use shorthand notation like "SFO" for parent only degradation models. The following two lines fit the model and produce the summary report of the model fit. This covers the numerical analysis given in the FOCUS report.

m.L1.SFO <- mkinfit("SFO", FOCUS_2006_L1_mkin, quiet = TRUE)
 summary(m.L1.SFO)
-
## mkin version used for fitting:    0.9.47.5 
+
## mkin version used for fitting:    0.9.47.6 
 ## R version used for fitting:       3.5.2 
-## Date of fit:     Thu Jan 31 16:52:39 2019 
-## Date of summary: Thu Jan 31 16:52:39 2019 
+## Date of fit:     Thu Feb 21 14:32:06 2019 
+## Date of summary: Thu Feb 21 14:32:06 2019 
 ## 
 ## Equations:
 ## d_parent/dt = - k_parent_sink * parent
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 37 model solutions performed in 0.094 s
+## Fitted with method Port using 37 model solutions performed in 0.097 s
 ## 
 ## Weighting: none
 ## 
@@ -200,17 +200,17 @@
 plot(m.L1.FOMC, show_errmin = TRUE, main = "FOCUS L1 - FOMC")

 

 
summary(m.L1.FOMC, data = FALSE)

-
## mkin version used for fitting:    0.9.47.5 
+
## mkin version used for fitting:    0.9.47.6 
 ## R version used for fitting:       3.5.2 
-## Date of fit:     Thu Jan 31 16:52:41 2019 
-## Date of summary: Thu Jan 31 16:52:41 2019 
+## Date of fit:     Thu Feb 21 14:32:08 2019 
+## Date of summary: Thu Feb 21 14:32:08 2019 
 ## 
 ## Equations:
 ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 611 model solutions performed in 1.49 s
+## Fitted with method Port using 611 model solutions performed in 1.547 s
 ## 
 ## Weighting: none
 ## 
@@ -295,17 +295,17 @@
      main = "FOCUS L2 - FOMC")

 

 
summary(m.L2.FOMC, data = FALSE)

-
## mkin version used for fitting:    0.9.47.5 
+
## mkin version used for fitting:    0.9.47.6 
 ## R version used for fitting:       3.5.2 
-## Date of fit:     Thu Jan 31 16:52:42 2019 
-## Date of summary: Thu Jan 31 16:52:42 2019 
+## Date of fit:     Thu Feb 21 14:32:08 2019 
+## Date of summary: Thu Feb 21 14:32:08 2019 
 ## 
 ## Equations:
 ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 81 model solutions performed in 0.198 s
+## Fitted with method Port using 81 model solutions performed in 0.201 s
 ## 
 ## Weighting: none
 ## 
@@ -366,10 +366,10 @@
      main = "FOCUS L2 - DFOP")

 

 
summary(m.L2.DFOP, data = FALSE)

-
## mkin version used for fitting:    0.9.47.5 
+
## mkin version used for fitting:    0.9.47.6 
 ## R version used for fitting:       3.5.2 
-## Date of fit:     Thu Jan 31 16:52:43 2019 
-## Date of summary: Thu Jan 31 16:52:43 2019 
+## Date of fit:     Thu Feb 21 14:32:09 2019 
+## Date of summary: Thu Feb 21 14:32:09 2019 
 ## 
 ## Equations:
 ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) *
@@ -378,7 +378,7 @@
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 336 model solutions performed in 0.839 s
+## Fitted with method Port using 336 model solutions performed in 0.86 s
 ## 
 ## Weighting: none
 ## 
@@ -458,10 +458,10 @@
 
The objects returned by mmkin are arranged like a matrix, with models as a row index and datasets as a column index.

 
We can extract the summary and plot for e.g. the DFOP fit, using square brackets for indexing which will result in the use of the summary and plot functions working on mkinfit objects.

 
summary(mm.L3[["DFOP", 1]])

-
## mkin version used for fitting:    0.9.47.5 
+
## mkin version used for fitting:    0.9.47.6 
 ## R version used for fitting:       3.5.2 
-## Date of fit:     Thu Jan 31 16:52:44 2019 
-## Date of summary: Thu Jan 31 16:52:44 2019 
+## Date of fit:     Thu Feb 21 14:32:10 2019 
+## Date of summary: Thu Feb 21 14:32:11 2019 
 ## 
 ## Equations:
 ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) *
@@ -470,7 +470,7 @@
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 137 model solutions performed in 0.341 s
+## Fitted with method Port using 137 model solutions performed in 0.351 s
 ## 
 ## Weighting: none
 ## 
@@ -559,17 +559,17 @@
 

 
The \(\chi^2\) error level of 3.3% as well as the plot suggest that the SFO model fits very well. The error level at which the \(\chi^2\) test passes is slightly lower for the FOMC model. However, the difference appears negligible.

 
summary(mm.L4[["SFO", 1]], data = FALSE)

-
## mkin version used for fitting:    0.9.47.5 
+
## mkin version used for fitting:    0.9.47.6 
 ## R version used for fitting:       3.5.2 
-## Date of fit:     Thu Jan 31 16:52:44 2019 
-## Date of summary: Thu Jan 31 16:52:45 2019 
+## Date of fit:     Thu Feb 21 14:32:11 2019 
+## Date of summary: Thu Feb 21 14:32:11 2019 
 ## 
 ## Equations:
 ## d_parent/dt = - k_parent_sink * parent
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 46 model solutions performed in 0.11 s
+## Fitted with method Port using 46 model solutions performed in 0.113 s
 ## 
 ## Weighting: none
 ## 
@@ -619,17 +619,17 @@
 ##        DT50 DT90
 ## parent  106  352

 
summary(mm.L4[["FOMC", 1]], data = FALSE)

-
## mkin version used for fitting:    0.9.47.5 
+
## mkin version used for fitting:    0.9.47.6 
 ## R version used for fitting:       3.5.2 
-## Date of fit:     Thu Jan 31 16:52:45 2019 
-## Date of summary: Thu Jan 31 16:52:45 2019 
+## Date of fit:     Thu Feb 21 14:32:11 2019 
+## Date of summary: Thu Feb 21 14:32:11 2019 
 ## 
 ## Equations:
 ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 66 model solutions performed in 0.162 s
+## Fitted with method Port using 66 model solutions performed in 0.168 s
 ## 
 ## Weighting: none
 ## 
-- 
cgit v1.2.1