From 1b2bdb8c02b1090f6c20f3ce9f7b461aab26e21c Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Tue, 28 Jun 2016 05:32:21 +0200 Subject: Static documentation rebuilt by staticdocs::build_site() --- vignettes/FOCUS_D.html | 6 +++--- vignettes/FOCUS_L.html | 40 ++++++++++++++++++++-------------------- vignettes/FOCUS_Z.pdf | Bin 238788 -> 238788 bytes vignettes/compiled_models.html | 32 ++++++++++++++++---------------- 4 files changed, 39 insertions(+), 39 deletions(-) (limited to 'vignettes') diff --git a/vignettes/FOCUS_D.html b/vignettes/FOCUS_D.html index 7f564ad8..e01c5869 100644 --- a/vignettes/FOCUS_D.html +++ b/vignettes/FOCUS_D.html @@ -192,8 +192,8 @@ print(FOCUS_2006_D) 
summary(fit)
## mkin version:    0.9.43 
 ## R version:       3.3.1 
-## Date of fit:     Tue Jun 28 01:37:23 2016 
-## Date of summary: Tue Jun 28 01:37:23 2016 
+## Date of fit:     Tue Jun 28 01:59:43 2016 
+## Date of summary: Tue Jun 28 01:59:44 2016 
 ## 
 ## Equations:
 ## d_parent = - k_parent_sink * parent - k_parent_m1 * parent
@@ -201,7 +201,7 @@ print(FOCUS_2006_D)
## ## Model predictions using solution type deSolve ## -## Fitted with method Port using 153 model solutions performed in 1.704 s +## Fitted with method Port using 153 model solutions performed in 1.698 s ## ## Weighting: none ## diff --git a/vignettes/FOCUS_L.html b/vignettes/FOCUS_L.html index 154ac896..828d9def 100644 --- a/vignettes/FOCUS_L.html +++ b/vignettes/FOCUS_L.html @@ -235,15 +235,15 @@ FOCUS_2006_L1_mkin <- mkin_wide_to_long(FOCUS_2006_L1) summary(m.L1.SFO) 
## mkin version:    0.9.43 
 ## R version:       3.3.1 
-## Date of fit:     Tue Jun 28 01:37:25 2016 
-## Date of summary: Tue Jun 28 01:37:25 2016 
+## Date of fit:     Tue Jun 28 01:59:45 2016 
+## Date of summary: Tue Jun 28 01:59:45 2016 
 ## 
 ## Equations:
 ## d_parent = - k_parent_sink * parent
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 37 model solutions performed in 0.244 s
+## Fitted with method Port using 37 model solutions performed in 0.243 s
 ## 
 ## Weighting: none
 ## 
@@ -328,8 +328,8 @@ summary(m.L1.SFO)
summary(m.L1.FOMC, data = FALSE)
## mkin version:    0.9.43 
 ## R version:       3.3.1 
-## Date of fit:     Tue Jun 28 01:37:27 2016 
-## Date of summary: Tue Jun 28 01:37:27 2016 
+## Date of fit:     Tue Jun 28 01:59:47 2016 
+## Date of summary: Tue Jun 28 01:59:47 2016 
 ## 
 ## 
 ## Warning: Optimisation by method Port did not converge.
@@ -341,7 +341,7 @@ summary(m.L1.SFO)
## ## Model predictions using solution type analytical ## -## Fitted with method Port using 188 model solutions performed in 1.236 s +## Fitted with method Port using 188 model solutions performed in 1.241 s ## ## Weighting: none ## @@ -425,8 +425,8 @@ plot(m.L2.FOMC, show_residuals = TRUE, 
summary(m.L2.FOMC, data = FALSE)
## mkin version:    0.9.43 
 ## R version:       3.3.1 
-## Date of fit:     Tue Jun 28 01:37:29 2016 
-## Date of summary: Tue Jun 28 01:37:29 2016 
+## Date of fit:     Tue Jun 28 01:59:49 2016 
+## Date of summary: Tue Jun 28 01:59:49 2016 
 ## 
 ## Equations:
 ## d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent
@@ -495,8 +495,8 @@ plot(m.L2.DFOP, show_residuals = TRUE, show_errmin = TRUE,
 
summary(m.L2.DFOP, data = FALSE)

 
## mkin version:    0.9.43 
 ## R version:       3.3.1 
-## Date of fit:     Tue Jun 28 01:37:32 2016 
-## Date of summary: Tue Jun 28 01:37:32 2016 
+## Date of fit:     Tue Jun 28 01:59:52 2016 
+## Date of summary: Tue Jun 28 01:59:52 2016 
 ## 
 ## Equations:
 ## d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -505,7 +505,7 @@ plot(m.L2.DFOP, show_residuals = TRUE, show_errmin = TRUE,
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 336 model solutions performed in 2.323 s
+## Fitted with method Port using 336 model solutions performed in 2.286 s
 ## 
 ## Weighting: none
 ## 
@@ -584,8 +584,8 @@ plot(mm.L3)

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

 
## mkin version:    0.9.43 
 ## R version:       3.3.1 
-## Date of fit:     Tue Jun 28 01:37:35 2016 
-## Date of summary: Tue Jun 28 01:37:35 2016 
+## Date of fit:     Tue Jun 28 01:59:55 2016 
+## Date of summary: Tue Jun 28 01:59:56 2016 
 ## 
 ## Equations:
 ## d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -594,7 +594,7 @@ plot(mm.L3)

 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 137 model solutions performed in 0.909 s
+## Fitted with method Port using 137 model solutions performed in 0.907 s
 ## 
 ## Weighting: none
 ## 
@@ -684,15 +684,15 @@ plot(mm.L4)
summary(mm.L4[["SFO", 1]], data = FALSE)
## mkin version:    0.9.43 
 ## R version:       3.3.1 
-## Date of fit:     Tue Jun 28 01:37:36 2016 
-## Date of summary: Tue Jun 28 01:37:37 2016 
+## Date of fit:     Tue Jun 28 01:59:57 2016 
+## Date of summary: Tue Jun 28 01:59:58 2016 
 ## 
 ## Equations:
 ## d_parent = - k_parent_sink * parent
 ## 
 ## Model predictions using solution type analytical 
 ## 
-## Fitted with method Port using 46 model solutions performed in 0.312 s
+## Fitted with method Port using 46 model solutions performed in 0.306 s
 ## 
 ## Weighting: none
 ## 
@@ -744,15 +744,15 @@ plot(mm.L4)
summary(mm.L4[["FOMC", 1]], data = FALSE)
## mkin version:    0.9.43 
 ## R version:       3.3.1 
-## Date of fit:     Tue Jun 28 01:37:37 2016 
-## Date of summary: Tue Jun 28 01:37:37 2016 
+## Date of fit:     Tue Jun 28 01:59:57 2016 
+## Date of summary: Tue Jun 28 01:59:58 2016 
 ## 
 ## Equations:
 ## 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.417 s
+## Fitted with method Port using 66 model solutions performed in 0.42 s
 ## 
 ## Weighting: none
 ## 
diff --git a/vignettes/FOCUS_Z.pdf b/vignettes/FOCUS_Z.pdf
index ac3af921..efda950f 100644
Binary files a/vignettes/FOCUS_Z.pdf and b/vignettes/FOCUS_Z.pdf differ
diff --git a/vignettes/compiled_models.html b/vignettes/compiled_models.html
index 9e2dff36..5e426a7f 100644
--- a/vignettes/compiled_models.html
+++ b/vignettes/compiled_models.html
@@ -250,21 +250,21 @@ mb.1 <- microbenchmark(
 print(mb.1)
## Unit: seconds
 ##                   expr       min        lq      mean    median        uq
-##  deSolve, not compiled 25.160343 25.384579 25.502052 25.608814 25.672906
-##       Eigenvalue based  2.219737  2.234679  2.244043  2.249621  2.256196
-##      deSolve, compiled  1.825299  1.843813  1.856078  1.862327  1.871467
+##  deSolve, not compiled 25.042204 25.078629 25.467550 25.115054 25.680223
+##       Eigenvalue based  2.273059  2.277424  2.285719  2.281790  2.292049
+##      deSolve, compiled  1.878785  1.883750  1.891594  1.888716  1.897998
 ##        max neval cld
-##  25.736998     3   b
-##   2.262771     3  a 
-##   1.880606     3  a
+## 26.245391 3 b +## 2.302308 3 a +## 1.907281 3 a 
autoplot(mb.1)
-

-

We see that using the compiled model is by a factor of 13.8 faster than using the R version with the default ode solver, and it is even faster than the Eigenvalue based solution implemented in R which does not need iterative solution of the ODEs:

+

+

We see that using the compiled model is by a factor of 13.3 faster than using the R version with the default ode solver, and it is even faster than the Eigenvalue based solution implemented in R which does not need iterative solution of the ODEs:

rownames(smb.1) <- smb.1$expr
 smb.1["median"]/smb.1["deSolve, compiled", "median"]
##                          median
-## deSolve, not compiled 13.750973
-## Eigenvalue based       1.207962
+## deSolve, not compiled 13.297425
+## Eigenvalue based       1.208117
 ## deSolve, compiled      1.000000
@@ -285,18 +285,18 @@ smb.1["median"]/smb.1["deSolve, compiled", "median" 
smb.2 <- summary(mb.2)
 print(mb.2)
## Unit: seconds
-##                   expr       min        lq      mean   median        uq
-##  deSolve, not compiled 54.725198 54.787875 54.893809 54.85055 54.978114
-##      deSolve, compiled  3.618315  3.644838  3.670582  3.67136  3.696716
+##                   expr      min        lq      mean    median        uq
+##  deSolve, not compiled 53.69252 53.938844 54.137601 54.185167 54.360141
+##      deSolve, compiled  3.42508  3.526298  3.588392  3.627516  3.670048
 ##        max neval cld
-##  55.105678     3   b
-##   3.722071     3  a
+## 54.535116 3 b +## 3.712579 3 a 
smb.2["median"]/smb.2["deSolve, compiled", "median"]
##   median
 ## 1     NA
 ## 2     NA
autoplot(mb.2)
-

+

Here we get a performance benefit of a factor of 14.9 using the version of the differential equation model compiled from C code!

This vignette was built with mkin 0.9.43 on

## R version 3.3.1 (2016-06-21)
-- 
cgit v1.2.1