As the COVID-19 pandemic continues worsening in the US, it is of critical importance todevelop a health information system that provides timely risk evaluation and prediction of theCOVID-19 infection in communities. We propose a spatiotemporal epidemiological forecastmodel that combines a spatial cellular automata (CA) with a temporal extended Susceptible-Antibody-Infectious-Removed (eSAIR) model under time-varying state-specific control mea-sures. This new toolbox enables the projection of the county-level COVID-19 prevalence over3,109 counties in the continental US, includingt-day ahead risk forecast and the risk relatedto a travel route. In comparison to the existing temporal risk prediction models, the proposedCA-eSAIR model informs the projected county-level risk to governments and residents of thelocal coronavirus spread patterns and the associated personal risks at specific geolocations. Suchhigh-resolution risk projection is useful for decision-making on business reopening and resourceallocation for COVID-19 tests.
To install and use this R package from Github, you will need to first install the R package devtools. Please uncomment the codes to install them. eSIR depends on three other packages, rjags (an interface to the JAGS library), chron and gtools, which could be installed with CA-eSAIR if not yet.
An error may occur if you have not yet installed JAGS-4.x.y.exe (for any x >= 0, y >=0). Windows users may download and install JAGS from here. Mac users may follow steps at casallas/8411082.
# install.packages("devtools") # library(devtools) # install_github("leyaozh/CA-eSAIR") library(CA-eSAIR)
Our data are collected daily from dxy.com. Alternatively, we notice some convenient access to COVID-19 data from GuangchuangYu/nCov2019.
For data outside China, we use JHU CSSE GitHub data. Another package coronavirus has its GitHub version udpated daily, which is also quite useful.
We also pre-calculate some parameters in
data("confirmed") # From JHU CSSE data("death") # From JHU CSSE data("recovered") # partly from 1Point3Acres data("USA_state_N") #population in each state
knitr::opts_chunk$set(echo = TRUE)
This is an R Markdown format used for publishing markdown documents to GitHub. When you click the Knit button all R code chunks are run and a markdown file (.md) suitable for publishing to GitHub is generated.
You can include R code in the document as follows:
summary(cars)
You can also embed plots, for example:
plot(pressure)
Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.