secondary_attack_rate {SIRmcmc}

Compute the secondary attack rate

Description

Using the value of a call to household_transmission, computes the secondary attack rate for households.

Usage

secondary_attack_rate(household_size,SIRmcmc,probs=c(0.5,0.025,0.975))

Arguments

household_size

A numeric vector containing the number of people in a household.

SIRmcmc

The value of a call to household_transmission.

probs

A numeric vector of the quantiles of the posterior distribution. The default is the median and the central 95% credible region.

Details

Computes the posterior probability distribution of the secondary attack rate from the Metropolis algorithm. Returns quantiles of the distribution specified in the probs argument.

Value

An array with secondary attack rates.

The first dimension of the array is the household size.

The second dimension of the array is the quantiles of the posterior distribution in the probs argument.

The third dimension of the array is the value of epsilon in SIRmcmc.