Compute the likelihood of a transmission pair given the serial interval distribution

Two distributions can be specified for the serial interval distribution: Normal and Gamma. A Normal distribution with mean \(\mu\) and standard deviation \(\sigma\) is assumed by default. If a Gamma distribution is specified, the mean \(\mu\) and standard deviation \(\sigma\) are converted to the shape and rate parameters for a Gamma distribution as: $$ \beta = \frac{\mu}{\sigma^2} \alpha = \left(\frac{\mu}{\beta}\right)^2 $$

Usage

get_likelihood(serial_intervals, mu, sigma, distn = "normal")

Arguments

serial_intervals

numeric vector of serial intervals

mu

numeric; mean of serial interval distribution

sigma

numeric; standard deviation of serial interval distribution

distn

string; distribution to be assumed for serial interval. Defaults to "normal", but also accepts "gamma".

Value

Likelihood of the transmission pair.