From 97b4b0a109d55080c781d9bbecec59b56f7037b5 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Mon, 9 Nov 2015 10:16:17 +0100 Subject: Update staticdocs --- inst/web/Extract.mmkin.html | 24 ++++++------- inst/web/mccall81_245T.html | 6 ++-- inst/web/mkinfit.html | 8 ++--- inst/web/mkinmod.html | 61 +++++++++----------------------- inst/web/mkinpredict.html | 6 ++-- inst/web/summary.mkinfit.html | 6 ++-- inst/web/transform_odeparms.html | 6 ++-- inst/web/vignettes/FOCUS_Z.pdf | Bin 224802 -> 224775 bytes inst/web/vignettes/compiled_models.html | 34 +++++++++--------- inst/web/vignettes/mkin.pdf | Bin 160268 -> 160268 bytes 10 files changed, 62 insertions(+), 89 deletions(-) (limited to 'inst') diff --git a/inst/web/Extract.mmkin.html b/inst/web/Extract.mmkin.html index c6ae3ef7..6a175a02 100644 --- a/inst/web/Extract.mmkin.html +++ b/inst/web/Extract.mmkin.html @@ -181,7 +181,7 @@ $calls $time user system elapsed - 0.260 0.000 0.262 + 0.268 0.000 0.269 $mkinmod <mkinmod> model generated with @@ -367,7 +367,7 @@ function (P) } return(mC) } -<environment: 0x2a98780> +<environment: 0x2d96b00> $cost_notrans function (P) @@ -389,7 +389,7 @@ function (P) scaleVar = scaleVar) return(mC) } -<environment: 0x2a98780> +<environment: 0x2d96b00> $hessian_notrans parent_0 alpha beta @@ -455,7 +455,7 @@ $bparms.state 99.66619 $date -[1] "Mon Nov 9 09:05:53 2015" +[1] "Mon Nov 9 10:11:49 2015" attr(,"class") [1] "mkinfit" "modFit" @@ -540,7 +540,7 @@ $calls $time user system elapsed - 0.088 0.000 0.086 + 0.112 0.000 0.112 $mkinmod <mkinmod> model generated with @@ -727,7 +727,7 @@ function (P) } return(mC) } -<environment: 0x2ca6db8> +<environment: 0x2fa3d28> $cost_notrans function (P) @@ -749,7 +749,7 @@ function (P) scaleVar = scaleVar) return(mC) } -<environment: 0x2ca6db8> +<environment: 0x2fa3d28> $hessian_notrans parent_0 k_parent_sink @@ -812,7 +812,7 @@ $bparms.state 99.17407 $date -[1] "Mon Nov 9 09:05:53 2015" +[1] "Mon Nov 9 10:11:48 2015" attr(,"class") [1] "mkinfit" "modFit" @@ -890,7 +890,7 @@ $calls $time user system elapsed - 0.088 0.000 0.086 + 0.112 0.000 0.112 $mkinmod <mkinmod> model generated with @@ -1077,7 +1077,7 @@ function (P) } return(mC) } -<environment: 0x2ca6db8> +<environment: 0x2fa3d28> $cost_notrans function (P) @@ -1099,7 +1099,7 @@ function (P) scaleVar = scaleVar) return(mC) } -<environment: 0x2ca6db8> +<environment: 0x2fa3d28> $hessian_notrans parent_0 k_parent_sink @@ -1162,7 +1162,7 @@ $bparms.state 99.17407 $date -[1] "Mon Nov 9 09:05:53 2015" +[1] "Mon Nov 9 10:11:48 2015" attr(,"class") [1] "mkinfit" "modFit" diff --git a/inst/web/mccall81_245T.html b/inst/web/mccall81_245T.html index f7a454ca..cd8d4d73 100644 --- a/inst/web/mccall81_245T.html +++ b/inst/web/mccall81_245T.html @@ -114,8 +114,8 @@
mkin version: 0.9.41 R version: 3.2.2 -Date of fit: Mon Nov 9 09:06:00 2015 -Date of summary: Mon Nov 9 09:06:00 2015 +Date of fit: Mon Nov 9 10:11:56 2015 +Date of summary: Mon Nov 9 10:11:56 2015 Equations: d_T245 = - k_T245_sink * T245 - k_T245_phenol * T245 @@ -124,7 +124,7 @@ d_anisole = + k_phenol_anisole * phenol - k_anisole_sink * anisole Model predictions using solution type deSolve -Fitted with method Port using 246 model solutions performed in 1.372 s +Fitted with method Port using 246 model solutions performed in 1.39 s Weighting: none diff --git a/inst/web/mkinfit.html b/inst/web/mkinfit.html index 96bdf95f..0550ee65 100644 --- a/inst/web/mkinfit.html +++ b/inst/web/mkinfit.html @@ -317,15 +317,15 @@ summary(fit)
mkin version: 0.9.41 R version: 3.2.2 -Date of fit: Mon Nov 9 09:06:03 2015 -Date of summary: Mon Nov 9 09:06:03 2015 +Date of fit: Mon Nov 9 10:11:58 2015 +Date of summary: Mon Nov 9 10:11:58 2015 Equations: d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent Model predictions using solution type analytical -Fitted with method Port using 64 model solutions performed in 0.18 s +Fitted with method Port using 64 model solutions performed in 0.194 s Weighting: none @@ -401,7 +401,7 @@ print(system.time(fit <- mkinfit(SFO_SFO, FOCUS_2006_D, solution_type = "eigen", quiet = TRUE)))
user system elapsed - 1.180 1.196 0.889 + 1.188 1.228 0.910
coef(fit)
diff --git a/inst/web/mkinmod.html b/inst/web/mkinmod.html index 797e96d5..7b43beb7 100644 --- a/inst/web/mkinmod.html +++ b/inst/web/mkinmod.html @@ -170,51 +170,24 @@ SFO_SFO <- mkinmod( Successfully compiled differential equation model from auto-generated C code.
-# The above model used to be specified like this, before the advent of mkinsub() -SFO_SFO <- mkinmod( - parent = list(type = "SFO", to = "m1"), - m1 = list(type = "SFO")) -
-Successfully compiled differential equation model from auto-generated C code. -
-# Show details of creating the C function -SFO_SFO <- mkinmod( - parent = mkinsub("SFO", "m1"), - m1 = mkinsub("SFO"), verbose = TRUE) -
-
Compilation argument: - /usr/lib/R/bin/R CMD SHLIB file16d352645eb9.c 2> file16d352645eb9.c.err.txt -Program source: - 1: #include <R.h> - 2: - 3: - 4: static double parms [3]; - 5: #define k_parent_sink parms[0] - 6: #define k_parent_m1 parms[1] - 7: #define k_m1_sink parms[2] - 8: - 9: void initpar(void (* odeparms)(int *, double *)) { - 10: int N = 3; - 11: odeparms(&N, parms); - 12: } - 13: - 14: - 15: void func ( int * n, double * t, double * y, double * f, double * rpar, int * ipar ) { - 16: - 17: f[0] = - k_parent_sink * y[0] - k_parent_m1 * y[0]; - 18: f[1] = + k_parent_m1 * y[0] - k_m1_sink * y[1]; - 19: } -
-Successfully compiled differential equation model from auto-generated C code. -
-# If we have several parallel metabolites -# (compare tests/testthat/test_synthetic_data_for_UBA_2014.R) -m_synth_DFOP_par <- mkinmod(parent = mkinsub("DFOP", c("M1", "M2")), - M1 = mkinsub("SFO"), - M2 = mkinsub("SFO"), - use_of_ff = "max", quiet = TRUE) - ## Not run: +# # The above model used to be specified like this, before the advent of mkinsub() +# SFO_SFO <- mkinmod( +# parent = list(type = "SFO", to = "m1"), +# m1 = list(type = "SFO")) +# +# # Show details of creating the C function +# SFO_SFO <- mkinmod( +# parent = mkinsub("SFO", "m1"), +# m1 = mkinsub("SFO"), verbose = TRUE) +# +# # If we have several parallel metabolites +# # (compare tests/testthat/test_synthetic_data_for_UBA_2014.R) +# m_synth_DFOP_par <- mkinmod(parent = mkinsub("DFOP", c("M1", "M2")), +# M1 = mkinsub("SFO"), +# M2 = mkinsub("SFO"), +# use_of_ff = "max", quiet = TRUE) +# # fit_DFOP_par_c <- mkinfit(m_synth_DFOP_par, # synthetic_data_for_UBA_2014[[12]]$data, # quiet = TRUE) diff --git a/inst/web/mkinpredict.html b/inst/web/mkinpredict.html index 3a6300d0..900ee26f 100644 --- a/inst/web/mkinpredict.html +++ b/inst/web/mkinpredict.html @@ -304,7 +304,7 @@ 201 20 4.978707 27.46227
user system elapsed - 0.008 0.024 0.004 + 0.004 0.028 0.005
system.time( print(mkinpredict(SFO_SFO, c(k_parent_m1 = 0.05, k_parent_sink = 0.1, k_m1_sink = 0.01), @@ -315,7 +315,7 @@ 201 20 4.978707 27.46227
user system elapsed - 0.004 0.020 0.002 + 0.020 0.000 0.003
system.time( print(mkinpredict(SFO_SFO, c(k_parent_m1 = 0.05, k_parent_sink = 0.1, k_m1_sink = 0.01), @@ -326,7 +326,7 @@ 201 20 4.978707 27.46227
user system elapsed - 0.064 0.000 0.066 + 0.092 0.000 0.092
diff --git a/inst/web/summary.mkinfit.html b/inst/web/summary.mkinfit.html index 02469bb6..2958b020 100644 --- a/inst/web/summary.mkinfit.html +++ b/inst/web/summary.mkinfit.html @@ -159,15 +159,15 @@
mkin version: 0.9.41 R version: 3.2.2 -Date of fit: Mon Nov 9 09:06:13 2015 -Date of summary: Mon Nov 9 09:06:13 2015 +Date of fit: Mon Nov 9 10:12:08 2015 +Date of summary: Mon Nov 9 10:12:08 2015 Equations: d_parent = - k_parent_sink * parent Model predictions using solution type analytical -Fitted with method Port using 35 model solutions performed in 0.1 s +Fitted with method Port using 35 model solutions performed in 0.105 s Weighting: none diff --git a/inst/web/transform_odeparms.html b/inst/web/transform_odeparms.html index 25dcaba9..36dfecf8 100644 --- a/inst/web/transform_odeparms.html +++ b/inst/web/transform_odeparms.html @@ -135,8 +135,8 @@ summary(fit, data=FALSE) # See transformed and backtransformed parameters
mkin version: 0.9.41 R version: 3.2.2 -Date of fit: Mon Nov 9 09:06:14 2015 -Date of summary: Mon Nov 9 09:06:14 2015 +Date of fit: Mon Nov 9 10:12:09 2015 +Date of summary: Mon Nov 9 10:12:09 2015 Equations: d_parent = - k_parent_sink * parent - k_parent_m1 * parent @@ -144,7 +144,7 @@ d_m1 = + k_parent_m1 * parent - k_m1_sink * m1 Model predictions using solution type deSolve -Fitted with method Port using 153 model solutions performed in 0.619 s +Fitted with method Port using 153 model solutions performed in 0.632 s Weighting: none diff --git a/inst/web/vignettes/FOCUS_Z.pdf b/inst/web/vignettes/FOCUS_Z.pdf index 1234de20..31480019 100644 Binary files a/inst/web/vignettes/FOCUS_Z.pdf and b/inst/web/vignettes/FOCUS_Z.pdf differ diff --git a/inst/web/vignettes/compiled_models.html b/inst/web/vignettes/compiled_models.html index 92919da9..7722d95a 100644 --- a/inst/web/vignettes/compiled_models.html +++ b/inst/web/vignettes/compiled_models.html @@ -99,18 +99,18 @@ smb.1 <- summary(mb.1)[-1] rownames(smb.1) <- c("deSolve, not compiled", "Eigenvalue based", "deSolve, compiled") print(smb.1) 
##                             min        lq      mean    median        uq
-## deSolve, not compiled 9307.3194 9319.9546 9332.8171 9332.5899 9345.5659
-## Eigenvalue based       855.3608  855.8081  869.4725  856.2555  876.5283
-## deSolve, compiled      686.6143  687.9256  698.0279  689.2369  703.7346
+## deSolve, not compiled 9442.5119 9447.2060 9458.3420 9451.9001 9466.2571
+## Eigenvalue based       868.6312  872.4552  895.3422  876.2792  908.6977
+## deSolve, compiled      691.9663  697.5653  701.1004  703.1643  705.6674
 ##                             max neval cld
-## deSolve, not compiled 9358.5420     3   c
-## Eigenvalue based       896.8012     3  b 
-## deSolve, compiled      718.2324     3 a
-

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

+## deSolve, not compiled 9480.6141 3 c +## Eigenvalue based 941.1163 3 b +## deSolve, compiled 708.1706 3 a +

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

smb.1["median"]/smb.1["deSolve, compiled", "median"]
##                          median
-## deSolve, not compiled 13.540468
-## Eigenvalue based       1.242324
+## deSolve, not compiled 13.441951
+## Eigenvalue based       1.246194
 ## deSolve, compiled      1.000000
@@ -127,17 +127,17 @@ print(smb.1) smb.2 <- summary(mb.2)[-1] rownames(smb.2) <- c("deSolve, not compiled", "deSolve, compiled") print(smb.2) -
##                             min        lq      mean    median        uq
-## deSolve, not compiled 20.130709 20.147223 20.180429 20.163737 20.205289
-## deSolve, compiled      1.235864  1.255748  1.267458  1.275632  1.283255
-##                             max neval cld
-## deSolve, not compiled 20.246841     3   b
-## deSolve, compiled      1.290878     3  a
+
##                             min       lq      mean    median        uq
+## deSolve, not compiled 20.444632 20.48824 20.557595 20.531857 20.614077
+## deSolve, compiled      1.251733  1.25179  1.275227  1.251846  1.286973
+##                           max neval cld
+## deSolve, not compiled 20.6963     3   b
+## deSolve, compiled      1.3221     3  a
smb.2["median"]/smb.2["deSolve, compiled", "median"]
##                         median
-## deSolve, not compiled 15.80686
+## deSolve, not compiled 16.40126
 ## deSolve, compiled      1.00000
-

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

+

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

This vignette was built with mkin 0.9.41 on

## R version 3.2.2 (2015-08-14)
 ## Platform: x86_64-pc-linux-gnu (64-bit)
diff --git a/inst/web/vignettes/mkin.pdf b/inst/web/vignettes/mkin.pdf
index 9c1a03f3..9f9f9af0 100644
Binary files a/inst/web/vignettes/mkin.pdf and b/inst/web/vignettes/mkin.pdf differ
-- 
cgit v1.2.1