ontario/ont_single_var.R

In bbolker/McMasterPandemic: Pandemic Model

library(McMasterPandemic)
library(splines)
library(dplyr)
library(parallel)

L <- load(".ont_keep.RData")
print(unique(ont_all_sub$var))
ont_noICU <- dplyr::filter(ont_all_sub, var != "ICU")
params <- fix_pars(read_params("ICU1.csv"))
dat <- ont_noICU
datH <- dat %>% filter(var == "H") %>% mutate(t_vec = as.numeric(date-min(date)))
datreport <- dat %>% filter(var == "report") %>% mutate(t_vec = as.numeric(date-min(date)))



Honly <- do.call(calibrate_comb,c(nlist(params, data=datH, use_spline=TRUE, use_phenomhet=TRUE,use_DEoptim=FALSE)))
reportonly <- do.call(calibrate_comb,c(nlist(params, data=datreport, use_spline=TRUE, use_phenomhet=TRUE,use_DEoptim=FALSE)))

plot(Honly)
plot(reportonly)

#rdsave(reportonly,Honly)