From 94b9a6980dce02c7c32c3bc3b61aad61a3ce3151 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Wed, 4 Jul 2018 08:50:43 +0200 Subject: Improve IORE fitting, use 1.1 as N starting value --- R/mkinfit.R | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) (limited to 'R') 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), -- cgit v1.2.1