Critical vaccination threshold for herd immunity

Formula: p_c = 1 - 1/R0; with vaccine efficacy e the coverage needed is p_c/e

Usage

Vaccthresh(R0, efficacy = 1)

Arguments

R0: See Usage.
efficacy: See Usage.

Value

List with threshold, coverage, feasible, R0, efficacy.

References

Anderson and May (1991), Infectious Diseases of Humans: Dynamics and Control, Oxford University Press. P_c = 1 - 1/R0 is the standard published result and is stated in the same form in every open source consulted.

Examples

V <- c(1, 2, 3, 4, 5, 6, 7, 8)
Vaccthresh(V)
#> $threshold
#> [1] 0.0000000 0.5000000 0.6666667 0.7500000 0.8000000 0.8333333 0.8571429
#> [8] 0.8750000
#> 
#> $coverage
#> [1] 0.0000000 0.5000000 0.6666667 0.7500000 0.8000000 0.8333333 0.8571429
#> [8] 0.8750000
#> 
#> $feasible
#> [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
#> 
#> $R0
#> [1] 1 2 3 4 5 6 7 8
#> 
#> $efficacy
#> [1] 1
#> 
#> $method
#> [1] "Critical vaccination threshold"
#> 
#> attr(,"class")
#> [1] "morie_rich_result" "list"