chainbinomial (version 0.1.5)

estimate_sar: Estimate Secondary Attack Rate of the Chain Binomial Model

Description

Given data on the number of infected after a number of generation, initial number of susceptible, and initial number of infected, estimate the secondary attack rate (SAR) using maximum likelihood.

Usage

estimate_sar(infected, s0, i0 = 1, generations = Inf, se = TRUE)

Value

A list of class sar with the following components:

• sar_hat The point estimate of the secondary attack rate.

• se Standard error of the estimate (if se = TRUE).

• loglikelihood The log likelihood value at the point estimate.

• data The input data.

Arguments

infected

numeric.

s0

numeric.

i0

numeric.

generations

numeric.

se

logical. If TRUE (default), the standard error is computed.

See Also

confint.sar() for calculating confidence intervals.

Examples

set.seed(234)
mydata <- rchainbinom(n = 15, s0 = 5, sar = 0.2, i0 = 1, generations = Inf)
res <- estimate_sar(infected = mydata, s0 = 5, i0 = 1, generations = Inf)