Skip to main content
  • Research article
  • Open access
  • Published:

Evaluating human papillomavirus vaccination programs in Canada: should provincial healthcare pay for voluntary adult vaccination?



Recently, provincial health programs in Canada and elsewhere have begun rolling out vaccination against human papillomavirus for girls aged 9–13. While vaccination is voluntary, the cost of vaccination is waived, to encourage parents to have their daughters vaccinated. Adult women who are eligible for the vaccine may still receive it, but at a cost of approximately CAN$400. Given the high efficacy and immunogenicity of the vaccine, the possibility of eradicating targeted types of the virus may be feasible, assuming the vaccination programs are undertaken strategically.


We develop a mathematical model to describe the epidemiology of vaccination against human papillomavirus, accounting for a widespread childhood vaccination program that may be supplemented by voluntary adult vaccination. A stability analysis is performed to determine the stability of the disease-free equilibrium. The critical vaccine efficacy and immunogenicity thresholds are derived, and the minimum level of adult vaccination required for eradication of targeted types is determined.


We demonstrate that eradication of targeted types is indeed feasible, although the burden of coverage for a childhood-only vaccination program may be high. However, if a small, but non-negligible, proportion of eligible adults can be vaccinated, then the possibility of eradication of targeted types becomes much more favourable. We provide a threshold for eradication in general communities and illustrate the results with numerical simulations. We also investigate the effects of suboptimal efficacy and immunogenicity and show that there is a critical efficacy below which eradication of targeted types is not possible. If eradication is possible, then there is a critical immunogenicity such that even 100% childhood vaccination will not eradicate the targeted types of the virus and must be supplemented with voluntary adult vaccination. However, the level of adult vaccination coverage required is modest and may be achieved simply by removing the cost burden to vaccination.


We recommend that provincial healthcare programs should pay for voluntary adult vaccination for women aged 14–26. However, it should be noted that our model results are preliminary, in that we have made a number of simplifying assumptions, including a lack of age-dependency in sexual partner rates, a lack of sexual activity outside of the vaccine age-range among females and a uniform age of sexual debut; thus, further work is desired to enhance the external generalisability of our results.

Peer Review reports


Human papillomavirus (HPV) produces epithelial tumors of the skin and mucous membranes [1]. There are over 100 types, many of which are relatively benign. However, some types have emerged as high risk because they produce lesions that may lead to carcinomas [2]. Resulting disease includes genital warts, respiratory papillomatosis, and cancer of the cervix, vulva, vagina, anus and the penis [2], as well as cancers of the head and neck [3]. Prevalence of HPV in Canada has been estimated at 24% in female university students [4]. In Ontario, 500 women are diagnosed with cervical cancer annually, leading to 140 deaths [5]. Between 30 and 40 types are transmitted through sexual contact. Without condoms, risk of transmission, given contact with an infected partner, is close to 90%; this risk is still high (40%) when condoms are used [6]. No antivirals have been developed for HPV and detection has largely relied on the recommended yearly pap smear, which locates cellular abnormalities that indicate that HPV may be present [7].

Cervical cancer is the second most common cause of death in women (after breast cancer) and it accounts for 10% of all cancers in women [8]. Progression to malignancy after acquisition of HPV usually takes at least 10 years [8]. Types 16 and 18 account for approximately 70% of these cervical cancers [810]. Merck and GlaxoSmithKline (GSK) have developed commercial vaccines which target types 16 and 18 [8]. Merck's vaccine also protects against types 6 and 11, which are responsible for 90% of external genital warts [8]. The combination of a successful vaccine and vaccination strategy, in combination with the yearly pap smear, seems to be the best approach towards preventing cervical cancer. The vaccine has been approved for women aged 9–26 [11]. Current studies are showing 90–100% efficacy as well as over 98% immunogenicity rates [8, 12] for existing vaccines, with no loss of immunity for at least 5 years [12]. Current vaccination programs against HPV [1319] consider two distinct groups: girls who have not yet begun to be sexually active and women under 26 who are sexually active.

Vaccination programs in Canada have begun in several provinces (Ontario, Nova Scotia, Prince Edward Island) in Fall 2007, targeting school children aged 9–13 [20]. The vaccine is covered by Canadian healthcare for girls in this age group, but is voluntary [5]. Sociological research suggests that about 77% of parents would immunize their children with the HPV vaccine [21]. Vaccination of women aged 14–26 is available, but not covered by provincial healthcare, costing approximately CAN$400 [5] for a full course (three doses) [8]. State-funded voluntary HPV vaccination for adults has been underway in other countries, such as France and Australia, since early 2007 [22]. Possible limitations of a vaccination program include: i) the vaccine may only be delivered to a proportion p of the population, ii) the vaccine may only confer immunogenicity in a proportion ε of the vaccinated population (ie the vaccine may not always take when administered to the patient), iii) the vaccine may have incomplete efficacy ψ (ie the vaccine doesn't always protect against infection during sexual intercourse), iv) the vaccine may wane over time and v) the vaccine does not target all HPV types.

In this paper, we develop a mathematical model to evaluate the effectiveness of supplementing a childhood vaccination program with voluntary adult vaccination. We address the following research questions: 1. Can a childhood-only vaccination program eradicate targeted types of HPV? 2. Should an adult vaccination program supplement childhood vaccination? 3. Is eradication of targeted types possible for vaccines with suboptimal efficacy or immunogenicity?

Results and Discussion

The mathematical model considers female children who may be vaccinated or unvaccinated. Female children grow into sexually active adult women, whereupon unvaccinated adult women may become vaccinated, or they may be infected with HPV. Vaccinated women may still be infected, although they will have a lower probability of infection, due to the vaccine. Men are also included in the model and may be susceptible or infected, but cannot be vaccinated. See Methods for details. The model is illustrated in Figure 1.

Figure 1
figure 1

The model. Schematic representation of the model representing flow of individuals from pre- to post-vaccination. School-aged girls may be vaccinated (C V ) or unvaccinated (C U ), depending on the proportion vaccinated (p) and the immunogenicity of the vaccine (ε). Children progress to sexually active adults at rate α. Unvaccinated adult women (A U ) can either be vaccinated (f) or become infected (I U ), with probability β N , when they meet an infected man N . Vaccinated adult women (A V ) can also become infected (I V ), but with reduced transmissibility, due to the efficacy of the vaccine ψ. Unvaccinated men (M) become infected upon contact with infected women, with transmission probability β M . The mortality rate of children is μ C and the leaving rate of adults is μ.

There is a threshold of eradication of targeted types given by

ε p = 1 ψ μ [ μ + f ( ε ¯ p ¯ ) ( 1 ψ ) μ 4 ( μ + f ( ε ¯ p ¯ ) ) ( α + μ C ) β M β N π M π W α ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqyTduMaemiCaaNaeyypa0tcfa4aaSaaaeaacqaIXaqmaeaacqaHipqEcqaH8oqBaaGcdaWadaqaaiabeY7aTjabgUcaRiabdAgaMjabcIcaOiqbew7aLzaaraGafmiCaaNbaebacqGGPaqkcqGGOaakcqaIXaqmcqGHsislcqaHipqEcqGGPaqkcqGHsisljuaGdaWcaaqaaiabeY7aTnaaCaaabeqaaiabisda0aaacqGGOaakcqaH8oqBcqGHRaWkcqWGMbGzcqGGOaakcuaH1oqzgaqeaiqbdchaWzaaraGaeiykaKIaeiykaKIaeiikaGIaeqySdeMaey4kaSIaeqiVd02aaSbaaeaacqWGdbWqaeqaaiabcMcaPaqaaiabek7aInaaBaaabaGaemyta0eabeaacqaHYoGydaWgaaqaaiabd6eaobqabaGaeqiWda3aaSbaaeaacqWGnbqtaeqaaiabec8aWnaaBaaabaGaem4vaCfabeaacqaHXoqyaaaakiaawUfacaGLDbaacqGGSaalaaa@69F1@

with parameters as given in Table 1. See Methods for mathematical details. Note that, if the vaccination is not at all efficacious (ψ = 0), or if the vaccine confers no immunogenicity to children (ε = 0), then it is not possible to eradicate targeted types of the disease, regardless of how many children or adults are vaccinated.

Table 1 List of symbols

We chose sample parameters to illustrate this threshold, using the following key assumptions: females enter the model as children at age 13, the mean rate of progression to sexually active adults is 3 years, women are in the sexually active pool for 10 years (ie until age 26, after which they cannot be vaccinated, so they are no longer under consideration), male partners are also in the pool for 10 years, the vaccine may not confer 100% protection and the probability of transmission from men to women is higher than the probability of transmission from women to men. Although women can still transmit the virus after age 26, we make the simplifying assumption that they will only infect male partners who do not continue to find new, younger female partners. Thus, while 50 year old men may have sexual relationships with 20 year old women (and may have multiple partners among this cohort), we assume that, after age 60, they do not begin relationships with a new cohort of 20 year old women. Parameters were also chosen so that the prevalence of HPV in the sexually active cohort was 24%, in line with Canadian estimates [4]. The maximum possible vaccination rate was set so that, if 100% of sexually active adults were vaccinated with a perfect vaccine, then the average woman would be vaccinated by 16 years 9 months. See the Table for parameters used and sample values.

A childhood-only vaccination program can eradicate targeted types of the disease, assuming 95% efficacy of the vaccine, if 81% of children are vaccinated (Figure 2A). However, if adult vaccination is included, the burden of childhood vaccination falls; if 20% of adults are vaccinated, then eradication of targeted types could be achieved with 74% childhood vaccination coverage. However, if 50% of adults are vaccinated, then eradication of targeted types could be achieved with only 55% childhood vaccination coverage. If no children were to be vaccinated, then eradication could still be achieved, assuming 81% of adults were vaccinated. As efficacy is reduced, these thresholds increase (Figure 2B). In particular, for our parameter set, the burden of childhood vaccination coverage without adult vaccination rises to 90% when the efficacy is reduced to 85%.

Figure 2
figure 2

Thresholds of eradication. Thresholds of eradication of targeted types for vaccination coverage levels in children versus adults. A. Threshold curve assuming 95% efficacy. Parameters are given in the Table. If no adult vaccination is undertaken, then a childhood vaccination program must cover 81% of the school-aged population for eradication of targeted types. If 20% of adults are vaccinated, then the burden of childhood vaccination reduces to 74%. However, if 50% of adults are vaccinated, then eradication of targeted types can be achieved with only 55% of children vaccinated. B. If the efficacy of the vaccine is reduced to 85%, the eradication threshold is increased, assuming all other parameters remain the same.

If the efficacy is reduced further, then eradication of targeted types is no longer possible (Figure 3A).

Figure 3
figure 3

Suboptimal efficacy and immunogenicity. Suboptimal values of the efficacy or the immunogenicity lead to qualitative changes in the outcome. A. There is a critical vaccine efficacy (77% in this example), such that no amount of vaccination can eradicate targeted types of the disease if the efficacy falls below this critical value. In this case, eradication of targeted types could only occur if more than 100% of children were vaccinated, which is not possible. B. If the vaccine efficacy allows for eradication of targeted types, there is a critical vaccine immunogenicity (80% in this example), such that even 100% of childhood vaccination will not eradicate targeted types of the disease. In this case, there is a minimum level of adult vaccination coverage that is required for eradication of targeted types (18% if the immunogenicity falls to 75%), even if 100% childhood coverage levels can be achieved. Adult immunogenicity was assumed to be equal to childhood immunogenicity. All other parameters as in the Table.

There is a critical vaccine efficacy ψ*, satisfying

ψ = 1 μ 4 ( μ + f ) ( α + μ C ) β M β N π M π W α , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqiYdK3aaWbaaSqabeaacqGHxiIkaaGccqGH9aqpcqaIXaqmcqGHsisljuaGdaWcaaqaaiabeY7aTnaaCaaabeqaaiabisda0aaacqGGOaakcqaH8oqBcqGHRaWkcqWGMbGzcqGGPaqkcqGGOaakcqaHXoqycqGHRaWkcqaH8oqBdaWgaaqaaiabdoeadbqabaGaeiykaKcabaGaeqOSdi2aaSbaaeaacqWGnbqtaeqaaiabek7aInaaBaaabaGaemOta4eabeaacqaHapaCdaWgaaqaaiabd2eanbqabaGaeqiWda3aaSbaaeaacqWGxbWvaeqaaiabeg7aHbaakiabcYcaSaaa@5089@

such that, if ψ <ψ*, then eradication of targeted types is not possible. See Methods for details. Thus, even if the vaccine mounts an immune response 100% of the time and we can vaccinate 100% of the population, if the efficacy is below this threshold (77% in our example), then the disease will persist. If the immunogenicity is reduced, then eradication of targeted types is not possible with a childhood-only vaccination program (Figure 3B). There is a critical immunogenicity value, ε*, satisfying

ε = 1 ψ [ 1 μ 4 ( α + μ C ) β N β M π M π W α ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqyTdu2aaWbaaSqabeaacqGHxiIkaaGccqGH9aqpjuaGdaWcaaqaaiabigdaXaqaaiabeI8a5baadaWadaqaaiabigdaXiabgkHiTmaalaaabaGaeqiVd02aaWbaaeqabaGaeGinaqdaaiabcIcaOiabeg7aHjabgUcaRiabeY7aTnaaBaaabaGaem4qameabeaacqGGPaqkaeaacqaHYoGydaWgaaqaaiabd6eaobqabaGaeqOSdi2aaSbaaeaacqWGnbqtaeqaaiabec8aWnaaBaaabaGaemyta0eabeaacqaHapaCdaWgaaqaaiabdEfaxbqabaGaeqySdegaaaGaay5waiaaw2faaOGaeiilaWcaaa@4F83@

such that, if ε <ε*, then even 100% childhood vaccination coverage will not eradicate targeted types of the disease. In this case (assuming ψ > ψ*), there is a critical proportion of adults who must be vaccinated. In our example, in the unrealistic case that 100% of children are vaccinated, we would still require 18% of adults to be vaccinated to achieve eradication of targeted types if the immunogenicity is reduced to 75% (Figure 3B).

The time course for eradication of targeted types, as with eradicating any disease, is a long one. We illustrate the temporal dynamics in Figure 4. Here, we model a community where the initial condition is the equilibrium infection among women in the sexually active cohort we are investigating. (We stress that this number will be much smaller than the total disease burden in this community.) To demonstrate the importance of adult vaccination, we model two scenarios: in the first, 65% of children are vaccinated, using a vaccine with 95% efficacy and 98% immunogenicity, but no adults are vaccinated. In this case, the number of infections is reduced by about half (Figure 4A), but is not eradicated after more than a century. However, if the same vaccination program is administered, but, additionally, 40% of adults are vaccinated, then the number of infections approaches zero (Figure 4B).

Figure 4
figure 4

Timecourse of eradication. Time course dynamics among a community where initial conditions are taken to be the equilibrium values of an existing infection, for a vaccine with 95% efficacy and 98% immunogenicity. A. If 65% of children are vaccinated, but no adult vaccination occurs, the prevalence in the sexually active cohort is reduced to approximately half of the previous disease burden after 150 years and remains so thereafter. B. If 65% of children and 40% of adults are vaccinated, the number of infections in the sexually active cohort continues to decline, eventually approaching zero. All other parameters as in Figure 4A.

The dependence of the results on parameter variation was also determined. While equation (3) is independent of the sample parameters we chose, and thus is effectively a robust sensitivity analysis unto itself, we illustrate the variation among likely parameter ranges, as given in the Table. The results are illustrated in Figure 5. We use the output variable as the proportion of adults who must be vaccinated in order to achieve eradication of targeted types, assuming that 77% of children are vaccinated [21]. If women spend less time sexually active, then the threshold may be reduced to zero (Figure 5A). Conversely, if women spend more time sexually active, then the threshold may rise considerably, although it will not exceed 100%, even if girls become sexually active as early as 13, which we argue is unlikely for the Canadian average. If the optimal age of vaccination with perfect coverage, γ/c, increases, then the threshold increases slightly (Figure 5B), but there is little variation even if a perfect vaccine can reach all "adult" women by age 13. Since parameters β N , β M , π M and π W always appear together in (3), we can consider them as a single unit. Even if the transmission probabilities or birth rates collectively double, then the threshold will only increase to 65% required adult vaccination (Figure 5C). Finally, there is little dependence on childhood mortality, as expected (Figure 5D); if life expectancy is significantly shorter than 70 years, then the threshold will be reduced, but we do not expect much variation in this parameter, given overall Canadian life expectancies.

Figure 5
figure 5

Dependence on parameter variation. Sensitivity of eradication threshold on parameter variation, assuming 77% childhood vaccination. A. If the average length of sexual activity before age 26 is 9 and a half years or less, then no adults need be vaccinated. If this time increases, then the threshold also increases, but still remains below 100%, even if children begin sexual activity at age 13. B. Variation as the optimal age of vaccination varies. This rate measures the age at which adult women are vaccinated, assuming 100% of adults are vaccinated with a perfect vaccine. In this case, there is little variation in the output. C. Variation in the transmission probabilities and birth rates. These parameters are always linked, so we examine them as a single unit. Even if these parameters collectively double, eradication could still be achieved with 65% adult vaccination. D. Variation in the childhood mortality rate. The output is relatively steady, unless mortality is very high. In this (unlikely) scenario, life expectancy of 30 years (17 years after vaccination at age 13) or less would clearly obviate the need for an adult vaccination program.


Despite human papillomavirus being one of the most prevalent sexually transmitted infectious diseases [23], eradication of targeted types is not only possible, but is feasible, under existing vaccination programs, if adults are included. While eradication of targeted types using a childhood-only vaccination program is possible for vaccines with currently estimated efficacy (90–99%) and immunogenicity (98%), it nevertheless poses an extremely high burden of coverage that must be met if only children are to be vaccinated. Crucially, sociological surveys of parents indicated that only 77% were willing to immunise their children with the vaccine [21]. Even in the best-case scenario of a vaccine with 95% efficacy and 98% immunogenicity, our results demonstrate that this proportion of children vaccinated is insufficient for eradication unless a small, but nontrivial, proportion of adults (approximately 16%) is also vaccinated (Figure 2A).

While vaccinating children is desirable, given that the majority of children can be located via the school system, our model demonstrates that any such vaccination program should be supplemented by an adult vaccination program. A widespread adult vaccination program may be unfeasible; however, our results indicate that eradication of targeted types is possible if the proportion of vaccinated children is large, while the proportion of vaccinated adults is modest. Since we do not need to locate all, or even most, sexually active young adults, the required coverage may be achieved merely by not posing obstacles to voluntary adult vaccination. Consequently, we recommend that adult vaccination in Canada be paid for by provincial healthcare, just as both childhood and adult vaccination of Hepatitis B currently is.

By determining an eradication threshold (see Methods), our model can also be used to examine scenarios outside the range of currently available data. Specifically, we investigated the effects that suboptimal efficacy or immunogenicity of the vaccine may have. Although these values have been high in the literature, such data have been taken from clinical trials, which have a relatively small sample size. With a vaccination program targeting a significant number of children, we allowed for the possibility that these effects may not be as high in the broader population. We demonstrated that there is a critical vaccine efficacy, below which eradication of targeted types is not possible, regardless of how many people are vaccinated. We also showed that there is a critical vaccine immunogenicity, below which even 100% childhood vaccination coverage will not lead to eradication of targeted types. Thus, a small but nonzero proportion of adult vaccination coverage is crucial for eradication of targeted types in this case.

A number of cohort, dynamic and population models have been developed in the literature, most of which show that vaccinating women can be cost-effective [24]. Most mathematical models have consistently predicted a useful role for vaccination, although assumptions about naturally acquired immunity and heterogeneity in risk behaviours have varied [26]. One such model, developed before the vaccine was made available, predicted a 30% reduction in epidemic prevalence if females were vaccinated [25]. However, this model assumed the vaccine would only be 75% effective, in line with our conclusions. The annual proportion of cervical cancer cases prevented was shown to be much higher when early adolescents were targeted, in a modelling study based on vaccination in Finland [27]. The authors recommended the implementation of catch-up vaccination at the start of a vaccination programme (ie vaccinating outside the target age group) in order to increase the speed with which a decrease in HPV and cervical cancer incidence is observed. Vaccinating girls and women in a catch-up program in the United States was estimated to be cost-effective, relative to other commonly accepted healthcare programs [28].

Our model has several limitations, which should be noted. First, we assume that the vaccine does not wane. Although this is consistent with current understandings of the vaccine, we only have five years of data to rely on. We assume that men who have sexual relations with women in the sexually active cohort (ie women who are eligible for adult vaccination) do not continue to find new partners in this age group as time goes on. Thus, while 50 year old men may sometimes have sexual relations with (multiple) 25 year old women, we assume that those same men stop having sexual relations with new 25 year olds after a certain period of time. In this way, the sexually active cohorts of men and women are linked only for the time (approximately 10 years) that adult women are sexually active and eligible for vaccination. We also assume that the vaccine plays no role in modifying the nature of the disease for infected vaccinated women; consequently, we assume that men infected by vaccinated women have the same probability of infection as men infected by unvaccinated women. However, we argue that if this were not the case, then the overall burden of infection would only be lowered, so our results are conservative.

It should also be noted that achieving actual eradication of targeted strains may be difficult in practice. Our model does not include stochasticity, although stochastic effects tend to hasten eradication. Furthermore, our model does not include age structure, which can influence prediction significantly, since sexual-mixing patterns tend to be highly age-dependent and heterogeneity will likely slow down transmission dynamics relative to homogenous mixing. The model also does not include type structure; this will have important consequences for replacement effects, ecological effects and other kinds of type interaction. Even eradicating types 16 and 18 may prove difficult within the context of high Canadian vaccine coverage; reintroductions from other populations with lower vaccine coverage, or vaccination heterogeneity among Canadian communities (eg rural or indigenous populations) may result in a recurring low epidemic. Such a situation currently exists with Hepatitis A, which has high vaccination coverage in Canada, but reintroductions occur from the United States [29]. Nevertheless, reducing the prevalence to the lowest levels possible, as predicted by our results, would obviously be the most desirable outcome.

A small proportion of adult vaccination can have a significant impact on eradication of targeted HPV types, when supplementing a large childhood vaccination program. This is due to the fact that childhood vaccination is a binary choice (children are either vaccinated or they are not), whereas adult vaccination is a rate, so that unvaccinated adults have a continual chance of being vaccinated, during their years in the sexually active cohort, if they are not infected first. However, there are drawbacks to adult vaccination: administering a vaccine to adults requires expensive physician office visits and compliance rates may be low. School-based vaccination is the cheapest and most efficient form of vaccination in the conventional sense of the term.

The mathematical model is a general one, with results that extend (via specific choice of parameters) to a variety of communities. We focus on Canada for specific implementation recommendations, but the model is applicable to other communities. It should also be noted that it may take decades, even centuries, to eradicate targeted types of the disease (Figure 4), but the most significant reduction should happen within 50 years. This is a remarkably fast timescale for disease eradication [30], but is possible because of the high efficacy and immunogenicity of the vaccine. However, even if eradication of targeted types is not achieved, vaccination of both children and adults is clearly desirable and efforts should still be made to vaccinate as widely as possible.

Future work will involve examining the differences between the Merck and GSK vaccines, the possibility of vaccines that wane and targeted vaccination for adults. A targeted vaccination program for adults could be implemented by offering free vaccination at STD clinics. This allows targeting of the subpopulation responsible for transmitting most sexually transmitted infections, although it should be noted that such a subpopulation is likely to already have high HPV incidence before attending the clinic. Such a model of targeted vaccination would require separate consideration for high-risk individuals. We will also generalise the model to include factors such as age-dependency in sexual partner rates, sexual activity outside of the vaccine age range among females and variable ages of sexual debut.


The model

The mathematical model that we have developed accounts for the vaccination of female children (C) and adult women (A). In particular, we assume that the pool of sexually active adults only includes women in an age range where they can still be vaccinated, as well as their male partners. Since men are not being vaccinated, we include them in our model only as age-equivalent cohort sexual partners of adult women. Thus, we account for the situation where older infected males may infect younger females who have or have not been vaccinated and who in turn infect other males. However, we assume that any men who interact with women in this cohort only do so for a finite amount of time (eg 10 years), although they may have multiple partners among this age group during that time. That is, while these men may have multiple partners among women aged 16–26, after a while they exhaust their contacts in this age group and are thus removed from the model.

We assume all vaccinated individuals are vaccinated before infection. Adult females vaccinated after infection are classified as "infected", since we assume that vaccination has no effect on this group. We consider vaccinated female children (C V ) and female adults (A V ) as those who received the vaccine and for whom the vaccine conferred immunogenicity. All other female individuals are functionally unvaccinated (C U , A U ).

We treat men as either uninfected (M) or infected (N) with infection probability β M after interaction with an infected woman. Infection of an unvaccinated woman (A U ) via contact with an infected man (N) occurs with transmission probability β N and results in an infected unvaccinated woman (I U ). We assume that the vaccine confers no disease-modifying effects once infection has occurred; thus, we assume that the probability of infection is the same for men, regardless of whether their (infected) partners are vaccinated or unvaccinated. We also assume that the vaccine does not wane over time.

If the vaccine efficacy is not perfect, then vaccinated individuals have the potential to become infected. Infection of a vaccinated woman (A V ) via contact with an infected man occurs with transmission probability β N , mediated by the efficacy (ψ) of the vaccine (β N (1 - ψ)), resulting in an infected vaccinated woman (I V ). Thus, if the vaccine is 100% effective, transmission will be zero. Conversely, if the vaccine is completely ineffective, transmission is identical to the probability of infection for unvaccinated individuals, β N .

The rate of appearance of new females is π W while that of males is π M . The mortality rate of children is μ C and the leaving rate of adults (including mortality) is μ. In this model, we do not consider disease-induced mortality, since it does not play a role in removing individuals from the sexually active pool that we are considering [8, 31]. Female children progress to adulthood at a rate α. Unvaccinated adult females who were not vaccinated as children may be vaccinated as adults via a function f denoting the rate at which unvaccinated adult females are vaccinated. This function depends on the proportion of adults vaccinated p ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiCaaNbaebaaaa@2D56@ and the immunogenicity of the vaccine in adults ε ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafqyTduMbaebaaaa@2D94@ :

f ( ε ¯ p ¯ ) = c ε ¯ p ¯ 1 ε ¯ p ¯ + γ , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOzayMaeiikaGIafqyTduMbaebacuWGWbaCgaqeaiabcMcaPiabg2da9KqbaoaalaaabaGaem4yamMafqyTduMbaebacuWGWbaCgaqeaaqaaiabigdaXiabgkHiTiqbew7aLzaaraGafmiCaaNbaebacqGHRaWkcqaHZoWzaaGccqGGSaalaaa@412D@

with c an attenuation constant and -1 the maximum possible rate of vaccination, assuming perfect coverage and immunogenicity. The proportion of adult vaccinations is thus ε ¯ p ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafqyTduMbaebacuWGWbaCgaqeaaaa@2F15@ . Since f corresponds to the rate of vaccination of adults (ie the rate at which unvaccinated adults leave the unvaccinated pool to go to the vaccinated pool), we want the rate to be zero if no adults are being vaccinated and we want the rate to be high if all adults are vaccinated (ie adults are vaccinated immediately upon entering). The model is illustrated in Figure 1.

Mathematical analysis

Mathematically, our model is represented as follows:

d C U d t = π W ( 1 ε p ) α C U μ C C U d C V d t = π W ε p α C V μ C C V d A U d t = α C U f ( ε ¯ p ¯ ) A U μ A U β N A U N d A V d t = α C V f ( ε ¯ p ¯ ) A U μ A V ( 1 ψ ) β N A V N d I U d t = β N A U N μ I U d I V d t = ( 1 ψ ) β N A V N μ I V d M d t = π M β M I U M μ M β M I V M d N d t = β M I U M μ N + β M I V M , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabGWaaaaaaKqbagaadaWcaaqaaiabdsgaKjabdoeadnaaBaaabaGaemyvaufabeaaaeaacqWGKbazcqWG0baDaaaakeaacqGH9aqpaeaacqaHapaCdaWgaaWcbaGaem4vaCfabeaakiabcIcaOiabigdaXiabgkHiTiabew7aLjabdchaWjabcMcaPiabgkHiTiabeg7aHjabdoeadnaaBaaaleaacqWGvbqvaeqaaOGaeyOeI0IaeqiVd02aaSbaaSqaaiabdoeadbqabaGccqWGdbWqdaWgaaWcbaGaemyvaufabeaaaKqbagaadaWcaaqaaiabdsgaKjabdoeadnaaBaaabaGaemOvayfabeaaaeaacqWGKbazcqWG0baDaaaakeaacqGH9aqpaeaacqaHapaCdaWgaaWcbaGaem4vaCfabeaakiabew7aLjabdchaWjabgkHiTiabeg7aHjabdoeadnaaBaaaleaacqWGwbGvaeqaaOGaeyOeI0IaeqiVd02aaSbaaSqaaiabdoeadbqabaGccqWGdbWqdaWgaaWcbaGaemOvayfabeaaaKqbagaadaWcaaqaaiabdsgaKjabdgeabnaaBaaabaGaemyvaufabeaaaeaacqWGKbazcqWG0baDaaaakeaacqGH9aqpaeaacqaHXoqycqWGdbWqdaWgaaWcbaGaemyvaufabeaakiabgkHiTiabdAgaMjabcIcaOiqbew7aLzaaraGafmiCaaNbaebacqGGPaqkcqWGbbqqdaWgaaWcbaGaemyvaufabeaakiabgkHiTiabeY7aTjabdgeabnaaBaaaleaacqWGvbqvaeqaaOGaeyOeI0IaeqOSdi2aaSbaaSqaaiabd6eaobqabaGccqWGbbqqdaWgaaWcbaGaemyvaufabeaakiabd6eaobqcfayaamaalaaabaGaemizaqMaemyqae0aaSbaaeaacqWGwbGvaeqaaaqaaiabdsgaKjabdsha0baaaOqaaiabg2da9aqaaiabeg7aHjabdoeadnaaBaaaleaacqWGwbGvaeqaaOGaeyOeI0IaemOzayMaeiikaGIafqyTduMbaebacuWGWbaCgaqeaiabcMcaPiabdgeabnaaBaaaleaacqWGvbqvaeqaaOGaeyOeI0IaeqiVd0Maemyqae0aaSbaaSqaaiabdAfawbqabaGccqGHsislcqGGOaakcqaIXaqmcqGHsislcqaHipqEcqGGPaqkcqaHYoGydaWgaaWcbaGaemOta4eabeaakiabdgeabnaaBaaaleaacqWGwbGvaeqaaOGaemOta4eajuaGbaWaaSaaaeaacqWGKbazcqWGjbqsdaWgaaqaaiabdwfavbqabaaabaGaemizaqMaemiDaqhaaaGcbaGaeyypa0dabaGaeqOSdi2aaSbaaSqaaiabd6eaobqabaGccqWGbbqqdaWgaaWcbaGaemyvaufabeaakiabd6eaojabgkHiTiabeY7aTjabdMeajnaaBaaaleaacqWGvbqvaeqaaaqcfayaamaalaaabaGaemizaqMaemysaK0aaSbaaeaacqWGwbGvaeqaaaqaaiabdsgaKjabdsha0baaaOqaaiabg2da9aqaaiabcIcaOiabigdaXiabgkHiTiabeI8a5jabcMcaPiabek7aInaaBaaaleaacqWGobGtaeqaaOGaemyqae0aaSbaaSqaaiabdAfawbqabaGccqWGobGtcqGHsislcqaH8oqBcqWGjbqsdaWgaaWcbaGaemOvayfabeaaaKqbagaadaWcaaqaaiabdsgaKjabd2eanbqaaiabdsgaKjabdsha0baaaOqaaiabg2da9aqaaiabec8aWnaaBaaaleaacqWGnbqtaeqaaOGaeyOeI0IaeqOSdi2aaSbaaSqaaiabd2eanbqabaGccqWGjbqsdaWgaaWcbaGaemyvaufabeaakiabd2eanjabgkHiTiabeY7aTjabd2eanjabgkHiTiabek7aInaaBaaaleaacqWGnbqtaeqaaOGaemysaK0aaSbaaSqaaiabdAfawbqabaGccqWGnbqtaKqbagaadaWcaaqaaiabdsgaKjabd6eaobqaaiabdsgaKjabdsha0baaaOqaaiabg2da9aqaaiabek7aInaaBaaaleaacqWGnbqtaeqaaOGaemysaK0aaSbaaSqaaiabdwfavbqabaGccqWGnbqtcqGHsislcqaH8oqBcqWGobGtcqGHRaWkcqaHYoGydaWgaaWcbaGaemyta0eabeaakiabdMeajnaaBaaaleaacqWGwbGvaeqaaOGaemyta0KaeiilaWcaaaaa@1160@

with the function f given by (2).

The disease-free equilibrium satisfies

C ¯ U = π W ( 1 ε p ) α + μ C C ¯ V = π W ε p α + μ C A ¯ U = α C U f + μ A ¯ V = α C V + f A U μ M ¯ = π M μ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabuWaaaaabaGafm4qamKbaebadaWgaaWcbaGaemyvaufabeaaaOqaaiabg2da9aqcfayaamaalaaabaGaeqiWda3aaSbaaeaacqWGxbWvaeqaaiabcIcaOiabigdaXiabgkHiTiabew7aLjabdchaWjabcMcaPaqaaiabeg7aHjabgUcaRiabeY7aTnaaBaaabaGaem4qameabeaaaaaakeaacuWGdbWqgaqeamaaBaaaleaacqWGwbGvaeqaaaGcbaGaeyypa0dajuaGbaWaaSaaaeaacqaHapaCdaWgaaqaaiabdEfaxbqabaGaeqyTduMaemiCaahabaGaeqySdeMaey4kaSIaeqiVd02aaSbaaeaacqWGdbWqaeqaaaaaaOqaaiqbdgeabzaaraWaaSbaaSqaaiabdwfavbqabaaakeaacqGH9aqpaKqbagaadaWcaaqaaiabeg7aHjabdoeadnaaBaaabaGaemyvaufabeaaaeaacqWGMbGzcqGHRaWkcqaH8oqBaaaakeaacuWGbbqqgaqeamaaBaaaleaacqWGwbGvaeqaaaGcbaGaeyypa0dajuaGbaWaaSaaaeaacqaHXoqycqWGdbWqdaWgaaqaaiabdAfawbqabaGaey4kaSIaemOzayMaemyqae0aaSbaaeaacqWGvbqvaeqaaaqaaiabeY7aTbaaaOqaaiqbd2eanzaaraaabaGaeyypa0dajuaGbaWaaSaaaeaacqaHapaCdaWgaaqaaiabd2eanbqabaaabaGaeqiVd0gaaaaaaaa@71B1@

and I ¯ U = I ¯ V = N ¯ = 0 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmysaKKbaebadaWgaaWcbaGaemyvaufabeaakiabg2da9iqbdMeajzaaraWaaSbaaSqaaiabdAfawbqabaGccqGH9aqpcuWGobGtgaqeaiabg2da9iabicdaWaaa@364C@ . We assume 0 ε , ε ¯ , p , p ¯ 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeGimaaJaeyizImQaeqyTduMaeiilaWIafqyTduMbaebacqGGSaalcqWGWbaCcqGGSaalcuWGWbaCgaqeaiabgsMiJkabigdaXaaa@3A0D@ . All other parameters are assumed to be positive. At the disease-free equilibrium, the Jacobian matrix is

J = [ μ C α 0 0 0 0 0 0 0 0 μ C α 0 0 0 0 0 0 α 0 f μ 0 0 0 0 β N A U 0 α f μ 0 0 0 ( 1 ψ ) β N A V 0 0 0 0 μ 0 0 β N A U 0 0 0 0 0 μ 0 ( 1 ψ ) β N A V 0 0 0 0 β M M β M M μ 0 0 0 0 0 β M M β M M 0 μ ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOsaOKaeyypa0ZaamWaaeaafaqabeaciaaaaaaaaeaacqGHsislcqaH8oqBdaWgaaWcbaGaem4qameabeaakiabgkHiTiabeg7aHbqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabgkHiTiabeY7aTnaaBaaaleaacqWGdbWqaeqaaOGaeyOeI0IaeqySdegabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeqySdegabaGaeGimaadabaGaeyOeI0IaemOzayMaeyOeI0IaeqiVd0gabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeyOeI0IaeqOSdi2aaSbaaSqaaiabd6eaobqabaGccqWGbbqqdaWgaaWcbaGaemyvaufabeaaaOqaaiabicdaWaqaaiabeg7aHbqaaiabdAgaMbqaaiabgkHiTiabeY7aTbqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabgkHiTiabcIcaOiabigdaXiabgkHiTiabeI8a5jabcMcaPiabek7aInaaBaaaleaacqWGobGtaeqaaOGaemyqae0aaSbaaSqaaiabdAfawbqabaaakeaacqaIWaamaeaacqaIWaamaeaacqaIWaamaeaacqaIWaamaeaacqGHsislcqaH8oqBaeaacqaIWaamaeaacqaIWaamaeaacqaHYoGydaWgaaWcbaGaemOta4eabeaakiabdgeabnaaBaaaleaacqWGvbqvaeqaaaGcbaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaeyOeI0IaeqiVd0gabaGaeGimaadabaGaeiikaGIaeGymaeJaeyOeI0IaeqiYdKNaeiykaKIaeqOSdi2aaSbaaSqaaiabd6eaobqabaGccqWGbbqqdaWgaaWcbaGaemOvayfabeaaaOqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabicdaWaqaaiabgkHiTiabek7aInaaBaaaleaacqWGnbqtaeqaaOGaemyta0eabaGaeyOeI0IaeqOSdi2aaSbaaSqaaiabd2eanbqabaGccqWGnbqtaeaacqGHsislcqaH8oqBaeaacqaIWaamaeaacqaIWaamaeaacqaIWaamaeaacqaIWaamaeaacqaIWaamaeaacqaHYoGydaWgaaWcbaGaemyta0eabeaakiabd2eanbqaaiabek7aInaaBaaaleaacqWGnbqtaeqaaOGaemyta0eabaGaeGimaadabaGaeyOeI0IaeqiVd0gaaaGaay5waiaaw2faaiabc6caUaaa@B37F@

Then det (J - λI) = (- μ C - α - λ)2 (- μ - λ)2(- f - μ - λ) det M, where

M = [ μ λ 0 β N A U 0 μ λ ( 1 ψ ) β N A V β M M β M M μ λ ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyta0Kaeyypa0ZaamWaaeaafaqabeWadaaabaGaeyOeI0IaeqiVd0MaeyOeI0Iaeq4UdWgabaGaeGimaadabaGaeqOSdi2aaSbaaSqaaiabd6eaobqabaGccqWGbbqqdaWgaaWcbaGaemyvaufabeaaaOqaaiabicdaWaqaaiabgkHiTiabeY7aTjabgkHiTiabeU7aSbqaaiabcIcaOiabigdaXiabgkHiTiabeI8a5jabcMcaPiabek7aInaaBaaaleaacqWGobGtaeqaaOGaemyqae0aaSbaaSqaaiabdAfawbqabaaakeaacqaHYoGydaWgaaWcbaGaemyta0eabeaakiabd2eanbqaaiabek7aInaaBaaaleaacqWGnbqtaeqaaOGaemyta0eabaGaeyOeI0IaeqiVd0MaeyOeI0Iaeq4UdWgaaaGaay5waiaaw2faaiabc6caUaaa@5B5A@

Thus, the largest eigenvalue for J will be the largest eigenvalue for M and so we can reduce the problem to solving

λ 3 + αλ 2 + βλ + γ = 0


α = 3μ

β = 3μ 2 - (1 - ψ)β M β N MA V + β N A U β M M

γ = μ 3 - (1 - ψ) μβ M β N MA V + μβ N A U β M M.

By the Routh-Hurwitz condition, all roots will have negative real parts if α > 0, γ > 0, and αβ - γ > 0. Clearly α > 0. We can write the third condition as

αβ - γ = 6μ 3 + 2γ,

which will be positive if γ > 0. Thus, all roots will have negative real part if

μ 2 - β N β M M[A U + (1 - ψ)A V ] > 0.

Solving for εp in terms of ε ¯ p ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafqyTduMbaebacuWGWbaCgaqeaaaa@2F15@ and substituting equilibrium values, our eradication threshold is thus

ε p = 1 ψ μ [ μ + f ( ε ¯ p ¯ ) ( 1 ψ ) μ 4 ( μ + f ( ε ¯ p ¯ ) ) ( α + μ C ) β M β N π M π W α ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqyTduMaemiCaaNaeyypa0tcfa4aaSaaaeaacqaIXaqmaeaacqaHipqEcqaH8oqBaaGcdaWadaqaaiabeY7aTjabgUcaRiabdAgaMjabcIcaOiqbew7aLzaaraGafmiCaaNbaebacqGGPaqkcqGGOaakcqaIXaqmcqGHsislcqaHipqEcqGGPaqkcqGHsisljuaGdaWcaaqaaiabeY7aTnaaCaaabeqaaiabisda0aaacqGGOaakcqaH8oqBcqGHRaWkcqWGMbGzcqGGOaakcuaH1oqzgaqeaiqbdchaWzaaraGaeiykaKIaeiykaKIaeiikaGIaeqySdeMaey4kaSIaeqiVd02aaSbaaeaacqWGdbWqaeqaaiabcMcaPaqaaiabek7aInaaBaaabaGaemyta0eabeaacqaHYoGydaWgaaqaaiabd6eaobqabaGaeqiWda3aaSbaaeaacqWGnbqtaeqaaiabec8aWnaaBaaabaGaem4vaCfabeaacqaHXoqyaaaakiaawUfacaGLDbaacqGGSaalaaa@69F1@


f ( ε ¯ p ¯ ) = c ε ¯ p ¯ 1 ε ¯ p ¯ + γ . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOzayMaeiikaGIafqyTduMbaebacuWGWbaCgaqeaiabcMcaPiabg2da9KqbaoaalaaabaGaem4yamMafqyTduMbaebacuWGWbaCgaqeaaqaaiabigdaXiabgkHiTiqbew7aLzaaraGafmiCaaNbaebacqGHRaWkcqaHZoWzaaGaeiOla4caaa@4127@

Differentiating (3), we have

p p ¯ = c ε ¯ ( 1 + γ ) ε ψ μ ( 1 ε ¯ p ¯ + γ ) 2 [ ( 1 μ ) μ 4 ( α + μ C ) β N β M π M π W α ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqGHciITcqWGWbaCaeaacqGHciITcuWGWbaCgaqeaaaacqGH9aqpdaWcaaqaaiabdogaJjqbew7aLzaaraGaeiikaGIaeGymaeJaey4kaSIaeq4SdCMaeiykaKcabaGaeqyTduMaeqiYdKNaeqiVd0MaeiikaGIaeGymaeJaeyOeI0IafqyTduMbaebacuWGWbaCgaqeaiabgUcaRiabeo7aNjabcMcaPmaaCaaabeqaaiabikdaYaaaaaWaamWaaeaacqGGOaakcqaIXaqmcqGHsislcqaH8oqBcqGGPaqkdaWcaaqaaiabeY7aTnaaCaaabeqaaiabisda0aaacqGGOaakcqaHXoqycqGHRaWkcqaH8oqBdaWgaaqaaiabdoeadbqabaGaeiykaKcabaGaeqOSdi2aaSbaaeaacqWGobGtaeqaaiabek7aInaaBaaabaGaemyta0eabeaacqaHapaCdaWgaaqaaiabd2eanbqabaGaeqiWda3aaSbaaeaacqWGxbWvaeqaaiabeg7aHbaaaiaawUfacaGLDbaacqGGUaGlaaa@6AF0@

It follows that there is a critical vaccine efficacy ψ*, satisfying

ψ = 1 μ 4 ( α + μ C ) β M β N π M π W α MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqiYdK3aaWbaaSqabeaacqGHxiIkaaGccqGH9aqpcqaIXaqmcqGHsisljuaGdaWcaaqaaiabeY7aTnaaCaaabeqaaiabisda0aaacqGGOaakcqaHXoqycqGHRaWkcqaH8oqBdaWgaaqaaiabdoeadbqabaGaeiykaKcabaGaeqOSdi2aaSbaaeaacqWGnbqtaeqaaiabek7aInaaBaaabaGaemOta4eabeaacqaHapaCdaWgaaqaaiabd2eanbqabaGaeqiWda3aaSbaaeaacqWGxbWvaeqaaiabeg7aHbaaaaa@4A00@

such that if ψ <ψ*, then eradication of targeted types is not possible. Thus, even if the vaccine mounts an immune response 100% of the time and we can vaccinate 100% of the population, if the efficacy is below this threshold, then the disease will persist. See Figure 3A.

By setting p = 1 and p ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiCaaNbaebaaaa@2D56@ = 0, it also follows from (3) that there is a critical immunogenicity value ε*, satisfying

ε = 1 ψ [ 1 μ 4 ( α + μ C ) β N β M π M π W α ] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqyTdu2aaWbaaSqabeaacqGHxiIkaaGccqGH9aqpjuaGdaWcaaqaaiabigdaXaqaaiabeI8a5baadaWadaqaaiabigdaXiabgkHiTmaalaaabaGaeqiVd02aaWbaaeqabaGaeGinaqdaaiabcIcaOiabeg7aHjabgUcaRiabeY7aTnaaBaaabaGaem4qameabeaacqGGPaqkaeaacqaHYoGydaWgaaqaaiabd6eaobqabaGaeqOSdi2aaSbaaeaacqWGnbqtaeqaaiabec8aWnaaBaaabaGaemyta0eabeaacqaHapaCdaWgaaqaaiabdEfaxbqabaGaeqySdegaaaGaay5waiaaw2faaaaa@4E99@

such that if ε <ε*, then even 100% childhood vaccination coverage will not eradicate targeted types of the disease.

From (3), when p = 1, we have

ε = ε + f ( ε ¯ p ¯ ) ψ μ [ ( 1 ψ ) μ 4 ( α + μ C ) β N β M π M π W α ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqyTduMaeyypa0JaeqyTdu2aaWbaaSqabeaacqGHxiIkaaGccqGHRaWkjuaGdaWcaaqaaiabdAgaMjabcIcaOiqbew7aLzaaraGafmiCaaNbaebacqGGPaqkaeaacqaHipqEcqaH8oqBaaGcdaWadaqaaiabcIcaOiabigdaXiabgkHiTiabeI8a5jabcMcaPiabgkHiTKqbaoaalaaabaGaeqiVd02aaWbaaeqabaGaeGinaqdaaiabcIcaOiabeg7aHjabgUcaRiabeY7aTnaaBaaabaGaem4qameabeaacqGGPaqkaeaacqaHYoGydaWgaaqaaiabd6eaobqabaGaeqOSdi2aaSbaaeaacqWGnbqtaeqaaiabec8aWnaaBaaabaGaemyta0eabeaacqaHapaCdaWgaaqaaiabdEfaxbqabaGaeqySdegaaaGccaGLBbGaayzxaaGaeiOla4caaa@5E22@


θ ψ μ ( ε ε ) μ 4 ( α + μ C ) β N β M π M π W α ( 1 ψ ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqiUdeNaeyyyIOBcfa4aaSaaaeaacqaHipqEcqaH8oqBcqGGOaakcqaH1oqzdaahaaqabeaacqGHxiIkaaGaeyOeI0IaeqyTduMaeiykaKcabaWaaSaaaeaacqaH8oqBdaahaaqabeaacqaI0aanaaGaeiikaGIaeqySdeMaey4kaSIaeqiVd02aaSbaaeaacqWGdbWqaeqaaiabcMcaPaqaaiabek7aInaaBaaabaGaemOta4eabeaacqaHYoGydaWgaaqaaiabd2eanbqabaGaeqiWda3aaSbaaeaacqWGnbqtaeqaaiabec8aWnaaBaaabaGaem4vaCfabeaacqaHXoqyaaGaeyOeI0IaeiikaGIaeGymaeJaeyOeI0IaeqiYdKNaeiykaKcaaOGaeiOla4caaa@5972@

For ψ > ψ*, and ε <ε*, θ > 0. It follows that the minimum level of adult vaccination required for eradication of targeted types, p ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiCaaNbaebadaahaaWcbeqaaiabgEHiQaaaaaa@2E72@ satisfies

p ¯ = θ ( 1 + γ ) ε ¯ ( c + θ ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiCaaNbaebadaahaaWcbeqaaiabgEHiQaaakiabg2da9KqbaoaalaaabaGaeqiUdeNaeiikaGIaeGymaeJaey4kaSIaeq4SdCMaeiykaKcabaGafqyTduMbaebacqGGOaakcqWGJbWycqGHRaWkcqaH4oqCcqGGPaqkaaGaeiOla4caaa@3F8B@

Since θ > 0, it follows that p ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiCaaNbaebadaahaaWcbeqaaiabgEHiQaaaaaa@2E72@ > 0. If the immunogenicities ε and ε ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafqyTduMbaebaaaa@2D94@ are not too small, then p ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiCaaNbaebadaahaaWcbeqaaiabgEHiQaaaaaa@2E72@ < 1 (since γ is small). Note that we are not assuming that childhood immunogenicity is necessarily the same as adult immunogenicity. Thus, if childhood immunogenicity is below ε* (but not so small that the vaccine is nonfunctional), then there is a minimum level of adult vaccination coverage that must be achieved for eradication of targeted types. See Figure 3B.

To examine sensitivity of results on parameter variation, we used the output parameter as the proportion of adults who should be vaccinated in order to eradicate targeted types. Thus, rearranging equation (3), we have

p ¯ = δ ( 1 + γ ) ε ¯ ( c + δ ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiCaaNbaebacqGH9aqpjuaGdaWcaaqaaiabes7aKjabcIcaOiabigdaXiabgUcaRiabeo7aNjabcMcaPaqaaiqbew7aLzaaraGaeiikaGIaem4yamMaey4kaSIaeqiTdqMaeiykaKcaaOGaeiilaWcaaa@3E49@


δ = β M β N π M π W α ( 1 ψ ) β M β N π M π W α μ 4 ( α + μ C ) ( ε μ ψ p + μ 5 ( α + μ C ) β M β N π M π W α μ ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqiTdqMaeyypa0tcfa4aaSaaaeaacqaHYoGydaWgaaqaaiabd2eanbqabaGaeqOSdi2aaSbaaeaacqWGobGtaeqaaiabec8aWnaaBaaabaGaemyta0eabeaacqaHapaCdaWgaaqaaiabdEfaxbqabaGaeqySdegabaGaeiikaGIaeGymaeJaeyOeI0IaeqiYdKNaeiykaKIaeqOSdi2aaSbaaeaacqWGnbqtaeqaaiabek7aInaaBaaabaGaemOta4eabeaacqaHapaCdaWgaaqaaiabd2eanbqabaGaeqiWda3aaSbaaeaacqWGxbWvaeqaaiabeg7aHjabgkHiTiabeY7aTnaaCaaabeqaaiabisda0aaacqGGOaakcqaHXoqycqGHRaWkcqaH8oqBdaWgaaqaaiabdoeadbqabaGaeiykaKcaaOWaaeWaaeaacqaH1oqzcqaH8oqBcqaHipqEcqWGWbaCcqGHRaWkjuaGdaWcaaqaaiabeY7aTnaaCaaabeqaaiabiwda1aaacqGGOaakcqaHXoqycqGHRaWkcqaH8oqBdaWgaaqaaiabdoeadbqabaGaeiykaKcabaGaeqOSdi2aaSbaaeaacqWGnbqtaeqaaiabek7aInaaBaaabaGaemOta4eabeaacqaHapaCdaWgaaqaaiabd2eanbqabaGaeqiWda3aaSbaaeaacqWGxbWvaeqaaiabeg7aHbaakiabgkHiTiabeY7aTbGaayjkaiaawMcaaiabc6caUaaa@7F41@


  1. eMedicine, Human papillomavirus. []

  2. Koutsky LA, Galloway DA, Holmes KK: Epidemiology of genital human papillomavirus infection. Epidemiol Rev. 1988, 10: 122-163.

    CAS  PubMed  Google Scholar 

  3. McKaig RG, Baric RS, Olshan AF: Human papillomavirus and head and neck cancer: epidemiology and molecular biology. Head Neck. 1998, 20: 250-265. 10.1002/(SICI)1097-0347(199805)20:3<250::AID-HED11>3.0.CO;2-O.

    Article  CAS  PubMed  Google Scholar 

  4. Richardson H, Kelsall G, Tellier P, Voyer H, Abrahamowicz M, Ferenczy A, Coutlée F, Franco EL: The Natural History of Type-specific Human Papillomavirus Infections in Female University Students. Cancer Epidemiology, Biomarkers & Prevention. 2003, 12 (6): 485-490.

    Google Scholar 

  5. Ontario Ministry of Health and Long-Term Care, Ontario's Grade 8 Vaccination program. []

  6. Winer RL, Hughes JP, Feng Q, O'Reilly S, Kiviat NB, Holmes KK, Koutsky LA: Condom use and the risk of genital human papillomavirus infection in young women. N Engl J Med. 2007, 354 (25): 2645-2654. 10.1056/NEJMoa053284.

    Article  Google Scholar 

  7. Passmore JA, Morroni C, Shapiro S, Williamson AL, Hoffman M: Papanicolaou smears and cervical inflmmatory cytokine responses. J Inamm. 2007, 4: 8-

    Google Scholar 

  8. Lowy DR, Schiller JT: Prophylactic human papillomavirus vaccines. J C Invest. 2006, 116 (5): 1167-1173. 10.1172/JCI28607.

    Article  CAS  Google Scholar 

  9. Walboomers LM, Jacobs MV, Manos MM, Bosch FX, Kummer JA, Shah KV, Snijders PJ, Peto J, Meijer CJ, Muñoz N: Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. J Pathol. 1999, 189: 12-19. 10.1002/(SICI)1096-9896(199909)189:1<12::AID-PATH431>3.0.CO;2-F.

    Article  CAS  PubMed  Google Scholar 

  10. Insinga RP, Dasbach EJ, Elbasha EH, Liaw K, Barr : Progression and regression of incident cervical HPV 6, 11, 16 and 18 infections in young women. Infectious Agents and Cancer. 2007, 2: 15-10.1186/1750-9378-2-15.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Public Health Agency of Canada, Human Papillomavirus (HPV) Prevention and HPV Vaccine: Questions and Answers. []

  12. Harper DM, Franco EL, Wheeler CM, Moscicki AB, Romanowski B, Roteli-Martins CM, Jenkins D, Schuind A, Costa Clemens SA, Dubin G, HPV Vaccination Study Group: Sustained efficacy up to 4.5 years of a bivalent L1 virus-like particle vaccine against human papillomavirus types 16 and 18: follow-up from a randomized control trial. Lancet. 2006, 367: 1247-1255. 10.1016/S0140-6736(06)68439-0.

    Article  CAS  PubMed  Google Scholar 

  13. Boot HJ, Wallenburg I, de Melker HE, Mangen MJ, Gerritsen AA, van der Maas NA, Berkhof J, Meijer CJ, Kimman TG: Assessing the introduction of universal human papillomavirus vaccination for preadolescent girls in The Netherlands. Vaccine. 2007, 25: 6245-6256. 10.1016/j.vaccine.2007.05.061.

    Article  PubMed  Google Scholar 

  14. Brisson M, van de Velde N, de Wals P, Boily MC: The potential cost-effectiveness of prophylactic human papillomavirus vaccines in Canada. Vaccine. 2007, 25: 5399-5408. 10.1016/j.vaccine.2007.04.086.

    Article  PubMed  Google Scholar 

  15. Hrgović Z, Izetbegović S: Primary prevention of cervical carcinoma. Med Arh. 2007, 61: 49-51.

    PubMed  Google Scholar 

  16. Pedersen C, Petaja T, Strauss G, Rume JC, Poder A, Richardus JH, Spiessens B, Descamps D, Hardt K, Lehtinen M, Dubin G, HPV Vaccine Adolescent Study Investigators Network: Immunization of early adolescent females with human papillomavirus type 16 and 18 L1 virus-like particle containing AS04 adjuvant. Journal of Adolescent Health. 2007, 40: 564-571. 10.1016/j.jadohealth.2007.02.015.

    Article  PubMed  Google Scholar 

  17. Signorelli C, Colzani E: Age at first intercourse and HPV immunization. J Prev Med Hyg. 2007, 48: 37-38.

    CAS  PubMed  Google Scholar 

  18. Schneider A, Schwartz TF, Hammerschmidt T, Siebert U: Impact of vaccination against oncogenic human papillomavirus on the incidence and mortality of cervical cancer in Germany. Med Klin (Munich). 2007, 102: 515-523. 10.1007/s00063-007-1065-9.

    Article  Google Scholar 

  19. Szarewski A: Prophylactic HPV vaccines. Eur J Gynaecol Oncol. 2007, 28: 165-169.

    CAS  PubMed  Google Scholar 

  20. Ontario Ministry of Health and Long-Term Care, Human Papilloma Virus Immunization Program Backgrounder. []

  21. Marshall H, Ryan P, Robertson D, Baghurst P: A cross-sectional survey to assess community attitudes to introduction of Human papillomavirus vaccine. Aust NZ J Public Health. 2007, 31: 235-242. 10.1111/j.1467-842X.2007.00054.x.

    Article  Google Scholar 

  22. Pallecaros A, Vonau B: Human papilloma virus vaccine – more than a vaccine. Current Opinion in Obstetrics & Gynecology. 2007, 19 (6): 541-546.

    Article  Google Scholar 

  23. Ho GY, Bierman R, Beardsley L, Chang CJ, Burk RD: Natural History of Cervicovaginal Papillomavirus Infection in Young Women. N Engl J Med. 1998, 338: 423-428. 10.1056/NEJM199802123380703.

    Article  CAS  PubMed  Google Scholar 

  24. Dasbach EJ, Elbasha EH, Insinga RP: Mathematical models for predicting the epidemiologic and economic impact of vaccination against human papillomavirus infection and disease. Epidemiol Rev. 2006, 28: 88-100. 10.1093/epirev/mxj006.

    Article  PubMed  Google Scholar 

  25. Hughes JP, Garnett GP, Koutsky L: The theoretical population level impact of a prophylactic human papilloma virus vaccine. Epidemiology. 2002, 13: 631-9. 10.1097/00001648-200211000-00006.

    Article  PubMed  Google Scholar 

  26. Garnett GP, Kim JJ, French K, Goldie SJ: Modelling the impact of HPV vaccines on cervical cancer and screening programmes. Vaccine. 2006, 24 (Suppl 3): 178-186. 10.1016/j.vaccine.2006.05.116.

    Article  Google Scholar 

  27. French KM, Barnabas RV, Lehtinen M, Kontula O, Pukkala E, Dillner J, Garnett JP: Strategies for the introduction of human papillomavirus vaccination: modelling the optimum age- and sex-specific pattern of vaccination in Finland. British Journal of Cancer. 2007, 96: 514-518. 10.1038/sj.bjc.6603575.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Elbasha EH, Dasbach EJ, Insinga RP: Model for Assessing Human Papillomavirus Vaccination Strategies. Emerging Infectious Diseases. 2007, 13: 28-41.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Bauch CT, Rao AS, Pham BZ, Krahn M, Gilca V, Duval B, Chen MH, Tricco AC: A dynamic model for assessing universal Hepatitis A vaccination in Canada. Vaccine. 2007, 25: 1719-1726. 10.1016/j.vaccine.2006.11.020.

    Article  CAS  PubMed  Google Scholar 

  30. Velasco-Hernandez JX, Gershengorn HB, Blower SM: Could widespread use of combination antiretroviral therapy eradicate HIV epidemics?. Lancet Infect Dis. 2002, 2: 487-493. 10.1016/S1473-3099(02)00346-8.

    Article  CAS  PubMed  Google Scholar 

  31. Schiffman M, Castle PE, Jeronimo J, Rodriguez AC, Wacholder S: Human papillomavirus and cervical cancer. Lancet. 2007, 370: 890-907. 10.1016/S0140-6736(07)61416-0.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references


ML is supported by an Ontario Graduate Scholarship. RS is supported by an NSERC Discovery grant and funding from MITACS. The authors thank Gaya Jayaraman (Public Health Agency of Canada), Tara Mac Eachern and Rad Mourtada for technical discussions and acknowledge two reviewers whose comments greatly improved the manuscript.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Robert J Smith.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors contributions

ML performed the mathematical analysis and wrote the analysis section. RS conceived of the study, performed the numerical simulations and wrote the main part of the paper.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Llamazares, M., Smith, R.J. Evaluating human papillomavirus vaccination programs in Canada: should provincial healthcare pay for voluntary adult vaccination?. BMC Public Health 8, 114 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: