From 9ae42bd20bc2543a94cf1581ba9820c2f9e3afbd Mon Sep 17 00:00:00 2001
From: Johannes Ranke
Date: Thu, 20 Apr 2023 19:53:28 +0200
Subject: Fix and rebuild documentation, see NEWS
I had to fix the two pathway vignettes, as they did not work with
the released version any more. So they and the multistart vignette
which got some small fixes as well were rebuilt.
Complete rebuild of the online docs with the released version. The
documentation of the 'hierarchial_kinetics' format had to be fixed
as well.
---
docs/reference/mkinpredict.html | 84 ++++++++++++++++++++++++-----------------
1 file changed, 49 insertions(+), 35 deletions(-)
(limited to 'docs/reference/mkinpredict.html')
diff --git a/docs/reference/mkinpredict.html b/docs/reference/mkinpredict.html
index 2aab0b50..20ace74a 100644
--- a/docs/reference/mkinpredict.html
+++ b/docs/reference/mkinpredict.html
@@ -19,13 +19,13 @@ kinetic parameters and initial values for the state variables.">
mkin
- 1.2.0
+ 1.2.3
News
@@ -107,10 +124,11 @@ kinetic parameters and initial values for the state variables.
outtimes = seq(0, 120, by = 0.1),
solution_type = "deSolve",
use_compiled = "auto",
+ use_symbols = FALSE,
method.ode = "lsoda",
atol = 1e-08,
rtol = 1e-10,
- maxsteps = 20000,
+ maxsteps = 20000L,
map_output = TRUE,
na_stop = TRUE,
...
@@ -166,8 +184,9 @@ solver is used.
The method that should be used for producing the
predictions. This should generally be "analytical" if there is only one
observed variable, and usually "deSolve" in the case of several observed
-variables. The third possibility "eigen" is faster but not applicable to
-some models e.g. using FOMC for the parent compound.
+variables. The third possibility "eigen" is fast in comparison to uncompiled
+ODE models, but not applicable to some models, e.g. using FOMC for the
+parent compound.
use_compiled
@@ -175,24 +194,27 @@ some models e.g. using FOMC for the parent compound.
mkinmod model is used, even if is present.
+use_symbols
+If set to TRUE
(default), symbol info present in
+the mkinmod object is used if available for accessing compiled code
+
+
method.ode
The solution method passed via mkinpredict to ode] in
-case the solution type is "deSolve". The default "lsoda" is performant, but
-sometimes fails to converge.
+case the solution type is "deSolve" and we are not using compiled code.
+When using compiled code, only lsoda is supported.
atol
-Absolute error tolerance, passed to ode. Default is 1e-8,
-lower than in lsoda.
+Absolute error tolerance, passed to the ode solver.
rtol
-Absolute error tolerance, passed to ode. Default is 1e-10,
-much lower than in lsoda.
+Absolute error tolerance, passed to the ode solver.
maxsteps
-Maximum number of steps, passed to ode.
+Maximum number of steps, passed to the ode solver.
map_output
@@ -325,15 +347,15 @@ as these always return mapped output.
#> time degradinol
#> 20.0000000 0.2478752
mkinpredict(SFO, c(k_degradinol = 0.3), c(degradinol = 100), 0:20,
- method = "lsoda")[21,]
+ method = "lsoda", use_compiled = FALSE)[21,]
#> time degradinol
#> 20.0000000 0.2478752
mkinpredict(SFO, c(k_degradinol = 0.3), c(degradinol = 100), 0:20,
- method = "ode45")[21,]
+ method = "ode45", use_compiled = FALSE)[21,]
#> time degradinol
#> 20.0000000 0.2478752
mkinpredict(SFO, c(k_degradinol = 0.3), c(degradinol = 100), 0:20,
- method = "rk4")[21,]
+ method = "rk4", use_compiled = FALSE)[21,]
#> time degradinol
#> 20.0000000 0.2480043
# rk4 is not as precise here
@@ -373,25 +395,17 @@ as these always return mapped output.
solution_type = "analytical", use_compiled = FALSE)[201,])
}
#> test relative elapsed
-#> 2 deSolve_compiled 1.0 0.004
-#> 4 analytical 1.0 0.004
-#> 1 eigen 5.5 0.022
-#> 3 deSolve 51.0 0.204
+#> 2 deSolve_compiled 1 0.002
+#> 4 analytical 1 0.002
+#> 1 eigen 4 0.008
+#> 3 deSolve 30 0.060
# \dontrun{
# Predict from a fitted model
f <- mkinfit(SFO_SFO, FOCUS_2006_C, quiet = TRUE)
f <- mkinfit(SFO_SFO, FOCUS_2006_C, quiet = TRUE, solution_type = "deSolve")
head(mkinpredict(f))
-#> DLSODA- At current T (=R1), MXSTEP (=I1) steps
-#> taken on this call before reaching TOUT
-#> In above message, I1 = 1
-#>
-#> In above message, R1 = 9.99904e-07
-#>
-#> Warning: an excessive amount of work (> maxsteps ) was done, but integration was not successful - increase maxsteps
-#> Warning: Returning early. Results are accurate, as far as they go
-#> Error in out[available, var]: (subscript) logical subscript too long
+#> Error in !is.null(x$symbols) & use_symbols: operations are possible only for numeric, logical or complex types
# }
@@ -408,7 +422,7 @@ as these always return mapped output.
--
cgit v1.2.1