aboutsummaryrefslogtreecommitdiff
path: root/R/mkinfit.R
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2014-07-17 14:34:20 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2014-07-17 14:34:20 +0200
commit8def5006fc81c032c3fc99751e062cdb32a81cc1 (patch)
tree7e93b84c99c50fad930ef4edc642be3c108af91c /R/mkinfit.R
parentd2c1ab854491ff047135fa8377400a68499e72de (diff)
Return complete list of initial states after fitting
This is useful for specifying state.ini in a subsequent call to mkinfit
Diffstat (limited to 'R/mkinfit.R')
-rw-r--r--R/mkinfit.R5
1 files changed, 4 insertions, 1 deletions
diff --git a/R/mkinfit.R b/R/mkinfit.R
index d591c42a..a8fbfc78 100644
--- a/R/mkinfit.R
+++ b/R/mkinfit.R
@@ -417,8 +417,11 @@ mkinfit <- function(mkinmod, observed,
fit$bparms.optim <- bparms.optim
fit$bparms.fixed <- bparms.fixed
- # Return ode parameters for further fitting
+ # Return ode and state parameters for further fitting
fit$bparms.ode <- bparms.all[mkinmod$parms]
+ fit$bparms.state <- c(bparms.all[setdiff(names(bparms.all), names(fit$bparms.ode))],
+ state.ini.fixed)
+ names(fit$bparms.state) <- gsub("_0$", "", names(fit$bparms.state))
fit$date <- date()

Contact - Imprint