aboutsummaryrefslogtreecommitdiff
path: root/R/mkinpredict.R
diff options
context:
space:
mode:
authorjranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb>2012-07-03 15:46:14 +0000
committerjranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb>2012-07-03 15:46:14 +0000
commit5b28712ad11470293e7ab6746718ab3edcfe966b (patch)
tree4d1b3329e04ca007048317536b5dd596bab29eec /R/mkinpredict.R
parente832b2255ad18d85e4ea9cc445ee87825a71c3a8 (diff)
- Add rtol argument in order to increase the precision of deSolve based
predictions for complex models - Improve default values for atol and set rtol default to 1e-10 git-svn-id: svn+ssh://svn.r-forge.r-project.org/svnroot/kinfit/pkg/mkin@49 edb9625f-4e0d-4859-8d74-9fd3b1da38cb
Diffstat (limited to 'R/mkinpredict.R')
-rw-r--r--R/mkinpredict.R3
1 files changed, 2 insertions, 1 deletions
diff --git a/R/mkinpredict.R b/R/mkinpredict.R
index 493c578..5170dc6 100644
--- a/R/mkinpredict.R
+++ b/R/mkinpredict.R
@@ -18,7 +18,7 @@
# You should have received a copy of the GNU General Public License along with
# this program. If not, see <http://www.gnu.org/licenses/>
-mkinpredict <- function(mkinmod, odeparms, odeini, outtimes, solution_type = "deSolve", map_output = TRUE, atol = 1e-6, ...) {
+mkinpredict <- function(mkinmod, odeparms, odeini, outtimes, solution_type = "deSolve", map_output = TRUE, atol = 1e-8, rtol = 1e-10, ...) {
# Get the names of the state variables in the model
mod_vars <- names(mkinmod$diffs)
@@ -92,6 +92,7 @@ mkinpredict <- function(mkinmod, odeparms, odeini, outtimes, solution_type = "de
func = mkindiff,
parms = odeparms,
atol = atol,
+ rtol = rtol,
...
)
}

Contact - Imprint