scripts/job_submissions/runs_18112017/rio_forecasts.R

In jameshay218/zikaProj: zikaInfer

datFile <- "~/Documents/Zika/Data/brazil/microceph_reports_2016.csv"
incFile <- "~/Documents/Zika/Data/brazil/zika_inc_reports.csv"

local <- "riograndedonorte"
microDat <- read.csv(datFile,
                     stringsAsFactors=FALSE)
incDat <- read.csv(incFile,
                   stringsAsFactors=FALSE)
microDat <- microDat[microDat$local == local,]
incDat <- incDat[incDat$local == local,]
parTab[parTab$names %in% c("L_H","N_H"),"values"] <- as.numeric(incDat[1,c("L_H","N_H")])
parTab[parTab$names == "L_H","values"] <- parTab[parTab$names == "L_H","values"]*365
parTab[parTab$names == "N_H","values"] <- as.integer(parTab[parTab$names == "N_H","values"])
parTab$local <- local
combos <- expand.grid(runName = "riograndedonorte_forecast",chainNo=1:maxChain,stringsAsFactors=FALSE)

jobs_rio_forecasts <- queuer::enqueue_bulk(obj1,combos,"model_forecasting",
                                           microDat=microDat,
                                           incDat=incDat,
                                           parTab=parTab,
                                           mcmcPars1=mcmcPars1,
                                           mcmcPars2=mcmcPars2,
                                           do_call=TRUE,timeout=0)