From 316f1360d1e12eefe491f86b0bbcd6dcf091c736 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Sat, 12 Jul 2014 08:32:08 +0200 Subject: First working fits with IORE model --- R/mkinfit.R | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) (limited to 'R/mkinfit.R') diff --git a/R/mkinfit.R b/R/mkinfit.R index b7ca1d74..26fc0e6b 100644 --- a/R/mkinfit.R +++ b/R/mkinfit.R @@ -90,13 +90,15 @@ mkinfit <- function(mkinmod, observed, defaultpar.names <- setdiff(mkinmod$parms, names(parms.ini)) for (parmname in defaultpar.names) { # Default values for rate constants, depending on the parameterisation - if (substr(parmname, 1, 2) == "k_") { + if (grepl("^k", parmname)) { parms.ini[parmname] = 0.1 + k_salt k_salt = k_salt + 1e-4 } # Default values for rate constants for reversible binding if (grepl("free_bound$", parmname)) parms.ini[parmname] = 0.1 if (grepl("bound_free$", parmname)) parms.ini[parmname] = 0.02 + # Default values for IORE exponents + if (grepl("^N", parmname)) parms.ini[parmname] = 1 # Default values for the FOMC, DFOP and HS models if (parmname == "alpha") parms.ini[parmname] = 1 if (parmname == "beta") parms.ini[parmname] = 10 -- cgit v1.2.1