Extinction probability for negative binomial branching process

Extinction probability for negative binomial branching process

nbbp_ep(r, k, tol = .Machine$double.eps)

Arguments

r

effective reproduction number

k

dispersion parameter: when <1, overdispersed

tol

tolerance passed to numerical solver

Value

a list, with $prob the extinction prob and $error a measure of numerical error

Details

See Nishiura et al 2012 (10.1016/j.jtbi.2011.10.039) equation 4. Probability that a branching process goes extinct. Solved numerically unless R < 1.0 (in which case the extinction probability is 1).