390 likes | 484 Vues
Learn the basic concepts, models, dynamics, and more. Dive into outbreak thresholds, immunity, and complex strategies with real-world examples and simulation tools. Understand the impact of vaccination coverage and effectiveness in controlling epidemics.
E N D
Modelling as a tool for planning vaccination policies Kari Auranen Department of Vaccines National Public Health Institute, Finland Department of Mathematics and Statistics University of Helsinki, Finland
Outline • Basic concepts and models • dynamics of transmission • herd immunity threshold • basic reproduction number • herd immunity and critical coverage of vaccination • mass action principle
Outline continues • Heterogeneously mixing populations • more complex models and new survey data to learn about routes of transmission • Use of models in decision making • example: varicella vaccinations in Finland • Ude of models in planning contaiment strategies • example: a simulation tool for pandemic influenza
A simple epidemic model (Hamer, 1906) • Consider an infection that • involves three “compartments” of infection: • proceeds in discrete generations (of infection) • is transmitted in a homogeneously mixing (“everyone meets everyone”) population of size N Susceptible Case Immune
Dynamics of transmission • Numbers of cases and susceptibles at generation t+1 • C = R * C * S / N • S = S - C + B t + 1 0 t t t+1 t t+1 t S = number of susceptibles at time t (i.e. generation t) C = number of cases (infectious individuals) at time t B = number of new susceptibles (by birth) t t t
Dynamics of transmission Epidemic threshold : S = N/R e 0
Epidemic threshold S e • S - S = - C + B • the number of susceptibles increases when C < B • decreases when C > B • the number of susceptibles cycles around the epidemic threshold S = N / R • this pattern is sustained as long as transmission is possible t+1 t t+1 t t+1 t t t+1 e 0
Epidemic threshold • C / C = R x S / N = S / S • the number of cases increases when S > S • decreases when S < S • the number of cases cycles around B (influx of new susceptibles) t+1 t 0 t t e e e t
Herd immunity threshold • incidence of infection decreases as long as the proportion of immunes exceeds the herd immunity threshold • H = 1- S / N • a complementary concept to the epidemic threshold • implies a critical vaccination coverage e
Basic reproduction number (R ) 0 • the average number of secondary cases that an infected individual produces in a totally susceptible population during his/her infectious period • in the Hamer model : R = R x 1 x N / N = R • herd immunity threshold H = 1 - 1 / R • in the endemic equilibrium: S = N / R , i.e., 0 0 0 0 e 0 R x S / N = 1 e e 0 0
Basic reproduction number (2) R = 3 0
Basic reproduction number (3) R = 3 endemic equilibrium 0 R x S / N = 1 0 e
Herd immunity threshold and R 0 H = 1-1/R 0 (Assumes homogeneous mixing)
Effect of vaccination Hamer model under vaccination S = S - C + B (1- VCxVE) Vaccine effectiveness (VE) x Vaccine coverage (VC) = 80% t t+1 t+1 Epidemic threshold sustained: S = N / R e 0
Mass action principle • all epidemic/transmission models are variations of the use of the mass action principle which • captures the effect of contacts between individuals • uses an analogy to modelling the rate of chemical reactions • is responsible for indirect effects of vaccination • assumes homogenous mixing • in the whole population • in appropriate subpopulations (defined by usually by age categories)
The SIR epidemic model • a continuous time model: overlapping generations • permanent immunity after infection • the system descibes the flow of individuals between the epidemiological compartments • uses a set of differential equations Susceptiple Infectious Removed
The SIR model equations = birth rate = rate of clearing infection = rate of infectious contacts by one individual = force of infection
Endemic equilibrium (SIR) N = 10,000 = 300/10000 (per time unit) = 10 (per time unit) = 1 (per time unit) 0
The basic reproduction number • Under the SIR model, Ro given by the ratio of two rates: • R = = rate of infectious contacts x • mean duration of infection • R not directly observable • need to derive relations to observable quantities 0 0
Force of infection • the number of infective contacts in the population per susceptible per time unit: • (t) = x I(t) / N • incidence rate of infection: (t) x S(t) • endemic force of infection = x (R - 1) 0
Estimation of R 0 Relation between the average age at infection and R (SIR model) 0 = 1/75 (per year)
A simple alternative formula • Assume everyone is infected at age A • everyone dies at age L (rectangular age distribution) Proportion 100 % Susceptibles Immunes Proportion of susceptibles: S / N = A / L R = N / S = L / A e A L e 0 Age (years)
Estimation of and Ro from seroprevalence data 1) Assume equilibrium 2) Parameterise force of infection 3) Estimate 4) Calculate Ro Ex. constant Proportion not yet infected: 1 - exp(- a) , estimate = 0.1 per year gives reasonable fit to the data
Estimates of R 0 * * * * Anderson and May: Infectious Diseases of Humans, 1991
Critical vaccination coverage to obtain herd immunity • Immunise a proportion p of newborns with a vaccine that offers complete protection against infection • R = (1-p) x R • If the proportion of vaccinated exceeds the herd immunity threshold, i.e., if p > H = 1-1/R , infection cannot persist in the population (herd immunity) vacc 0 0
Critical vaccination coverage as a function of R0 p = 1 – 1/R0
Indirect effects of vaccination • If p < H = 1-1/R , in the new endemic equilibrium: • S = N/R , = (R -1) • proportion of susceptibles remains untouched • force of infection decreases 0 e 0 vacc vacc e
Effect of vaccination on average age A’ at infection (SIR) • Life length L; proportion p vaccinated at birth, complete protection • every susceptible infected at age A Susceptibles S / N = (1-p) A’/L S / N = A/ L => A’ = A/(1-p) i.e., increase in the average age of infection Proportion e 1 e Immunes p A ’ L Age (years)
Vaccination at age V > 0 (SIR) • Assume proportion p vaccinated at age V • Every susceptible infected at age A • How big should p be to obtain herd immunity threshold H Proportion H = 1 - 1/R = 1 - A/L H = p (L-V)/L => p = (L-A)/(L-V) 1 Susceptibles p i.e., p bigger than when vaccination at birth Immunes V A L Age (years)
Modelling transmission in a heterogeneously mixing population
More complex mixing patterns • So far we have assumes (so called) homogeneous mixing • “everyone meets everyone” • More realistic models incorporate some form of heterogeneity in mixing (“who meets whom”) • e.g. individuals of the same age meet more often each other than individual from other age classes (assortative mixing)
Example: WAIFW matrix • structure of the Who Acquires Infection From Whom matrix for varicella , five age groups (0-4, 5-9, 10-14, 15-19, 20-75 years) table entry = rate of transmission between an infective and a susceptible of respective age groups e.g., force of infection in age group 0-4: a*I1 + a*I2 + c*I3 + d*I4 + e*I5 I1 = equilibrium number of infectives in age group 0-4, etc.
POLYMOD contact survey • Records the number of daily conversations in study participants in 7 European countries • Use the number of contacts between individuals from different age categories as a proxy for chances of transmission • Is currently being used to aid in modelling the impact of varicella vaccination in Finland
Use of models in policy making • Large-scale vaccinations usually bring along indirect effects • the mean age at disease increases • population immunity changes • Population-level experiments are impossible • Need for mathematicl modelling • to predict indirect effects of vaccination • to summarise the epidemiology of the infection • to identify missing data or knowledge about the natural history of the infection
References • Fine P.E.M, "Herd immunity: History, Theory, Practice", Epidemiologic Reviews, 15, 265-302,1993 • Fine P.E.M., "The contribution of modelling to vaccination policy, Vaccination and World Health, Eds. F.T. Cutts and P.G. Smith, Wiley and Sons, 1994. 3 Nokes D.J., Anderson R.M., "The use of mathematical models in the epidemiological study of infectious diseases and in the desing of mass immunization programmes", Epidemiology and Infection, 101, 1-20, 1988 • Anderson R.M. and May R.M., ”Infectious Diseases of Humans”; Oxford University Press, 1992. • Mossong J et al, Social contacts and mixing patterns relevant to the spread of infectious diseases: a multi-country population-based survey, Plos Medicine, in press • Duerr et al, Influenza pandemic intervention planning using InfluSim, BMC Infect Dis, 2007