1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/data.R
\docType{data}
\name{mers_2014_15}
\alias{mers_2014_15}
\title{Data on Middle East Respiratory Syndrome (MERS) in Saudi Arabia.}
\format{A list of two elements:
\itemize{
\item \strong{incidence:} a dataframe containing 495 days of observations with dates
in the first column, and number of local (2nd column) and imported (3rd
column) cases of MERS,
\item \strong{si:} a list of estimates of the mean (mean_si) and standard deviation
(std_si) of the serial interval for MERS.
}}
\source{
The incidence data was extracted from the EMPRES I system from FAO
(Global Animal Disease Information System - Food and Agriculture Organization
of the United Nations, 2017). Note incidence on the first day was originally
made of one local case and zero imported cases; this has been modified to
zero local cases and one imported case in the dataset shown here so the
reproduction number can be estimated from the start using the function
\code{\link[=estimate_R]{estimate_R()}}. The serial interval parameters were those
estimated by Cauchemez et al. (2016).
}
\description{
This data set gives:
\enumerate{
\item the daily incidence of onset of symptoms of laboratory confirmed human
infections with MERS-CoV in Saudi Arabia between the beginning of July
2014 and the end of December 2015, and
\item estimates of the mean and standard deviation of the serial interval for
MERS.
}
}
\examples{
## load data
data("mers_2014_15")
## estimate the reproduction number (method "parametric_si")
bimonthly_R <- estimate_R(mers_2014_15$incidence[,c("local", "imported")],
method = "parametric_si",
config = make_config(
mean_si = mers_2014_15$si$mean_si,
std_si = mers_2014_15$si$std_si,
t_start = 2:(nrow(mers_2014_15$incidence)-8*7),
t_end = (2:(nrow(mers_2014_15$incidence)-8*7)) + 8*7))
plot(bimonthly_R, legend = FALSE, add_imported_cases = TRUE,
options_I = list(col = c("local" = "black",
"imported" = "red"),
interval = 7, # show weekly incidence
ylab = "Weekly incidence"),
options_R = list(ylab = "Bimonthly R"))
# The first plot shows the weekly incidence,
# with imported cases shown in red and local cases in black
}
\references{
Global Animal Disease Information System - Food and Agriculture Organization
of the United Nations, 2017
Cauchemez S, Nouvellet P, Cori A, Jombart T, Garske T, Clapham H, Moore S,
Linden Mills H, Salje H, Collins C, et al. 2016.
Unraveling the drivers of MERS-CoV transmission.
Proc Natl Acad Sci 113: 9081-9086.
}
|