From cb318ea1067bf18d8ff896cef8d1e847a49e3aa7 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Tue, 28 Jun 2016 01:43:13 +0200 Subject: Static documentation rebuilt by staticdocs::build_site() --- vignettes/FOCUS_D.html | 6 +++--- vignettes/FOCUS_L.html | 42 ++++++++++++++++++++--------------------- vignettes/FOCUS_Z.pdf | Bin 237947 -> 238788 bytes vignettes/compiled_models.html | 32 +++++++++++++++---------------- 4 files changed, 40 insertions(+), 40 deletions(-) (limited to 'vignettes') diff --git a/vignettes/FOCUS_D.html b/vignettes/FOCUS_D.html index 48ae4654..7f564ad8 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:22:20 2016 
-## Date of summary: Tue Jun 28 01:22:21 2016 
+## Date of fit:     Tue Jun 28 01:37:23 2016 
+## Date of summary: Tue Jun 28 01:37:23 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.759 s +## Fitted with method Port using 153 model solutions performed in 1.704 s ## ## Weighting: none ## diff --git a/vignettes/FOCUS_L.html b/vignettes/FOCUS_L.html index 7026e3ce..154ac896 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:22:22 2016 
-## Date of summary: Tue Jun 28 01:22:22 2016 
+## Date of fit:     Tue Jun 28 01:37:25 2016 
+## Date of summary: Tue Jun 28 01:37:25 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.253 s
+## Fitted with method Port using 37 model solutions performed in 0.244 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:22:24 2016 
-## Date of summary: Tue Jun 28 01:22:24 2016 
+## Date of fit:     Tue Jun 28 01:37:27 2016 
+## Date of summary: Tue Jun 28 01:37:27 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.241 s +## Fitted with method Port using 188 model solutions performed in 1.236 s ## ## Weighting: none ## @@ -425,15 +425,15 @@ 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:22:27 2016 
-## Date of summary: Tue Jun 28 01:22:27 2016 
+## Date of fit:     Tue Jun 28 01:37:29 2016 
+## Date of summary: Tue Jun 28 01:37:29 2016 
 ## 
 ## Equations:
 ## 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.541 s
+## Fitted with method Port using 81 model solutions performed in 0.534 s
 ## 
 ## Weighting: none
 ## 
@@ -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:22:30 2016 
-## Date of summary: Tue Jun 28 01:22:30 2016 
+## Date of fit:     Tue Jun 28 01:37:32 2016 
+## Date of summary: Tue Jun 28 01:37:32 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.29 s
+## Fitted with method Port using 336 model solutions performed in 2.323 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:22:33 2016 
-## Date of summary: Tue Jun 28 01:22:33 2016 
+## Date of fit:     Tue Jun 28 01:37:35 2016 
+## Date of summary: Tue Jun 28 01:37:35 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.959 s
+## Fitted with method Port using 137 model solutions performed in 0.909 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:22:34 2016 
-## Date of summary: Tue Jun 28 01:22:35 2016 
+## Date of fit:     Tue Jun 28 01:37:36 2016 
+## Date of summary: Tue Jun 28 01:37:37 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.302 s
+## Fitted with method Port using 46 model solutions performed in 0.312 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:22:35 2016 
-## Date of summary: Tue Jun 28 01:22:35 2016 
+## Date of fit:     Tue Jun 28 01:37:37 2016 
+## Date of summary: Tue Jun 28 01:37:37 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.438 s
+## Fitted with method Port using 66 model solutions performed in 0.417 s
 ## 
 ## Weighting: none
 ## 
diff --git a/vignettes/FOCUS_Z.pdf b/vignettes/FOCUS_Z.pdf
index d21f7cb9..ac3af921 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 3a81c6c3..9e2dff36 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.634645 25.717260 25.925956 25.799875 26.071612
-##       Eigenvalue based  2.235689  2.237675  2.246796  2.239661  2.252349
-##      deSolve, compiled  1.869556  1.881345  1.887495  1.893133  1.896465
+##  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
 ##        max neval cld
-##  26.343349     3   b
-##   2.265038     3  a 
-##   1.899797     3  a
+## 25.736998 3 b +## 2.262771 3 a +## 1.880606 3 a 
autoplot(mb.1)
-

-

We see that using the compiled model is by a factor of 13.6 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.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:

rownames(smb.1) <- smb.1$expr
 smb.1["median"]/smb.1["deSolve, compiled", "median"]
##                          median
-## deSolve, not compiled 13.628133
-## Eigenvalue based       1.183044
+## deSolve, not compiled 13.750973
+## Eigenvalue based       1.207962
 ## 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 53.993961 54.012379 54.225054 54.030797 54.340600
-##      deSolve, compiled  3.450737  3.535551  3.583034  3.620364  3.649183
+##                   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
 ##        max neval cld
-##  54.650403     3   b
-##   3.678001     3  a
+## 55.105678 3 b +## 3.722071 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