nmode's Git Repositories - Rnaught/blob - man/idea.Rd

   1 % Generated by roxygen2: do not edit by hand
   2 % Please edit documentation in R/idea.R
   3 \name{idea}
   4 \alias{idea}
   5 \title{Incidence Decay and Exponential Adjustment (IDEA)}
   6 \usage{
   7 idea(cases, mu)
   8 }
   9 \arguments{
  10 \item{cases}{Vector of case counts. The vector must be of length at least two
  11 and only contain positive integers.}
  12
  13 \item{mu}{Mean of the serial distribution. This must be a positive number.
  14 The value should match the case counts in time units. For example, if case
  15 counts are weekly and the serial distribution has a mean of seven days,
  16 then \code{mu} should be set to \code{1}. If case counts are daily and the serial
  17 distribution has a mean of seven days, then \code{mu} should be set to \code{7}.}
  18 }
  19 \value{
  20 An estimate of the basic reproduction number (R0).
  21 }
  22 \description{
  23 This function implements a least squares estimation method of R0 due to
  24 Fisman et al. (PloS One, 2013). See details for implementation notes.
  25 }
  26 \details{
  27 This method is closely related to that implemented in \code{\link[=id]{id()}}. The method is
  28 based on an incidence decay model. The estimate of R0 is the value which
  29 minimizes the sum of squares between observed case counts and case counts
  30 expected under the model.
  31
  32 This method is based on an approximation of the SIR model, which is most
  33 valid at the beginning of an epidemic. The method assumes that the mean of
  34 the serial distribution (sometimes called the serial interval) is known. The
  35 final estimate can be quite sensitive to this value, so sensitivity testing
  36 is strongly recommended. Users should be careful about units of time (e.g.,
  37 are counts observed daily or weekly?) when implementing.
  38 }
  39 \examples{
  40 # Weekly data.
  41 cases <- c(1, 4, 10, 5, 3, 4, 19, 3, 3, 14, 4)
  42
  43 # Obtain R0 when the serial distribution has a mean of five days.
  44 idea(cases, mu = 5 / 7)
  45
  46 # Obtain R0 when the serial distribution has a mean of three days.
  47 idea(cases, mu = 3 / 7)
  48 }
  49 \references{
  50 \href{https://doi.org/10.1371/journal.pone.0083622}{Fisman et al. (PloS One, 2013)}
  51 }
  52 \seealso{
  53 \code{\link[=id]{id()}} for a similar method.
  54 }