diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2018-07-04 08:50:43 +0200 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2018-07-04 08:50:43 +0200 |
| commit | 94b9a6980dce02c7c32c3bc3b61aad61a3ce3151 (patch) | |
| tree | 6c1ddd37ea4f6e7e9c893a9e66bb005c66f5adff /R | |
| parent | 0b27c116b44d09460f3151b87762193d67d3b09e (diff) | |
Improve IORE fitting, use 1.1 as N starting value
Diffstat (limited to 'R')
| -rw-r--r-- | R/mkinfit.R | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/R/mkinfit.R b/R/mkinfit.R index c471b01c..840faf53 100644 --- a/R/mkinfit.R +++ b/R/mkinfit.R @@ -145,7 +145,7 @@ mkinfit <- function(mkinmod, observed, 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
+ if (grepl("^N", parmname)) parms.ini[parmname] = 1.1
# Default values for the FOMC, DFOP and HS models
if (parmname == "alpha") parms.ini[parmname] = 1
if (parmname == "beta") parms.ini[parmname] = 10
@@ -389,6 +389,9 @@ mkinfit <- function(mkinmod, observed, upper[other_fraction_parms] <- 1
}
+ # Show parameter names if tracing is requested
+ if(trace_parms) cat(names(c(state.ini.optim, transparms.optim)), "\n")
+
# Do the fit and take the time
fit_time <- system.time({
fit <- modFit(cost, c(state.ini.optim, transparms.optim),
|