Reproductive number of chickenpox among elderly in Egypt is 10.00 (Kitala et al, 2002)
Ro = 10.00
Generate sequence of numbers for fraction of population vaccinated from 0 to 1 with 0.1 interval
fraction_vaccinated_chickenpox = seq (0, 1, 0.1)
Print fraction of population vaccinated.
cat ("Fraction of population vaccinated: ", fraction_vaccinated_chickenpox)
## Fraction of population vaccinated: 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Compute effective reproductive number
Re = Ro * (1 - fraction_vaccinated_chickenpox)
Print effective reproductive number
cat ("Effective reproductive number:", Re)
## Effective reproductive number: 10 9 8 7 6 5 4 3 2 1 0
Compute herd immunity threshold
herd_immunity_threshold = 1 - (1/Ro)
Print herd immunity threshold
cat ("herd immunity threshold = ", herd_immunity_threshold)
## herd immunity threshold = 0.9
Plot fraction of population vaccinated (versus) effective reproductive number
subtitle = paste ("Ro = ", Ro, ", herd immunity threshold = ", round (herd_immunity_threshold, digits = 4), "; elderly (Egypt) - Kitala et al (2002)", sep="")
plot (fraction_vaccinated_chickenpox, Re, main = "Chicken pox" , sub = subtitle, xlab = "Fraction of Population Vaccinated \n", ylab = "Effective Reproductive Number (Re)")
Results and Discussion
Reproductive number of chickenpox among elderly in Egypt is 10.00 (Kitala et al, 2002). The graph illustrates that as fraction of population vaccinated increases, effective reproductive number decreases. Herd immunity threshold is 90%; that is, at this level of vaccination coverage, effective reproductive number is 1 (Re = 1). When vaccination coverage is above 90%, effective reproductive number is less than 1 (Re < 1); thereby, chicken pox will be eliminated at these higher levels of vaccination coverage
Public heath implications
Recommend chicken pox vaccination among elderly in Egypt at coverage levels of above 90%