diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2014-07-12 08:32:08 +0200 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2014-07-12 08:32:08 +0200 |
| commit | 316f1360d1e12eefe491f86b0bbcd6dcf091c736 (patch) | |
| tree | f155d7bcf82738b5623b4bc3a96aa431b007fceb /R/mkinfit.R | |
| parent | 65987c9b1ff9206859105c0f94be8346d457262e (diff) | |
First working fits with IORE model
Diffstat (limited to 'R/mkinfit.R')
| -rw-r--r-- | R/mkinfit.R | 4 |
1 files changed, 3 insertions, 1 deletions
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
|