plot(function(x) IORE.solution(x, 100, 0.2, 1.3), 0, 2, ylim = c(0, 100))
if (FALSE) {
fit.fomc <- mkinfit("FOMC", FOCUS_2006_C, quiet = TRUE)
fit.iore <- mkinfit("IORE", FOCUS_2006_C, quiet = TRUE)
fit.iore.deS <- mkinfit("IORE", FOCUS_2006_C, solution_type = "deSolve", quiet = TRUE)
print(data.frame(fit.fomc$par, fit.iore$par, fit.iore.deS$par,
row.names = paste("model par", 1:4)))
print(rbind(fomc = endpoints(fit.fomc)$distimes, iore = endpoints(fit.iore)$distimes,
iore.deS = endpoints(fit.iore)$distimes))
}
Run the code above in your browser using DataLab