S. pneumoniaetransmission according to inclusion in conjugate vaccines: Bayesian analysis of a longitudinal followup in schools
 Simon Cauchemez^{1, 2}Email author,
 Laura Temime^{1, 3},
 AlainJacques Valleron^{1, 2, 4},
 Emmanuelle Varon^{5},
 Guy Thomas^{1, 2, 4},
 Didier Guillemot^{6} and
 PierreYves Boëlle^{1, 2, 4}
DOI: 10.1186/14712334614
© Cauchemez et al; licensee BioMed Central Ltd. 2006
Received: 28 June 2005
Accepted: 30 January 2006
Published: 30 January 2006
Abstract
Background
Recent trends of pneumococcal colonization in the United States, following the introduction of conjugate vaccination, indicate that nonvaccine serotypes tend to replace vaccine serotypes. The eventual extent of this replacement is however unknown and depends on serotypespecific carriage and transmission characteristics.
Methods
Here, some of these characteristics were estimated for vaccine and nonvaccine serotypes from the followup of 4,488 schoolchildren in France in 2000. A Bayesian approach using Markov chain Monte Carlo data augmentation techniques was used for estimation.
Results
Vaccine and nonvaccine serotypes were found to have similar characteristics: the mean duration of carriage was 23 days (95% credible interval (CI): 21, 25 days) for vaccine serotypes and 22 days (95% CI: 20, 24 days) for nonvaccine serotypes; within a school of size 100, the Secondary Attack Rate was 1.1% (95% CI: 1.0%, 1.2%) for both vaccine and nonvaccine serotypes.
Conclusion
This study supports that, in 3–6 years old children, no competitive advantage exists for vaccine serotypes compared to nonvaccine serotypes. This is an argument in favour of important serotype replacement. It would be important to validate the result for infants, who are known to be the main reservoir in maintaining transmission. Overall reduction in pathogenicity should also be taken into account in forecasting the future burden of pneumococcal colonization in vaccinated populations.
Background
Less than 10 Streptococcus pneumoniae (S. pneumoniae) serotypes, among more than 90, have been included in the S. pneumoniae conjugate vaccine formulation. Since these serotypes account for a large part of carriage (almost 80% of all carriage in the United States [1]) and the vaccine protects against colonization [2], a reduction in overall carriage after vaccination is expected, and has indeed been observed in vaccinated populations [3]. However, recent trends of pneumococcal colonization in the United States, following the introduction of conjugate vaccination, indicate serotype replacement [4, 5], whereby a decrease in vaccine serotypes carriage is followed by an increase in carriage of nonvaccine serotypes. The eventual extent of this replacement is as yet unknown.
Mathematical models have predicted that replacement could occur in case of direct competition between serotypes in the absence of crossimmunity [6–9]. In these models, the extent of replacement is expected to be important if the degree of competition is high; and in this latter case, all the more that vaccine and nonvaccine serotypes have the same duration of carriage and transmission rate. However, this last assumption (same duration and transmission) could be challenged, considering the wide heterogeneity observed in serotypespecific prevalence, in favor of reduced capacity to colonize in the less frequent serotypes.
Ekdahl et al. [10] reported that there was no difference in the duration of carriage among 5 serogroups, 3 of which are included in the vaccine. Smith et al. [11] found that average duration of carriage ranged 6.7 to 62.5 days among 25 serotypes, 6 of which are included in the vaccine. However, these variations did not clearly identify vaccine serotypes as having longer carriage duration. Reported transmission rates of S. pneumoniae estimated from household data did not investigate serotype specific characteristics, nor inclusion in the vaccine [12, 13]. While it has been reported that nonvaccine serotypes may be less invasive [14], there is currently little information on carriage and transmission characteristics according to serotype inclusion in the vaccine. This lack of knowledge may stem from the difficulty to analyze field data on colonization, which generally consist of samples irregularly spaced in time, and lack information about the times of colonization and decolonization.
Here, using recent computationalbased statistical techniques, we analyze a prospective followup study of schoolchildren in France, with a view to estimating carriage and transmission characteristics of S. pneumoniae serotypes according to inclusion in the vaccine.
Methods
Data
A five month longitudinal study of 3–6yearold children in 81 schools, presented in detail elsewhere [15], was conducted from January to May 2000 in France. Oropharyngeal pneumococcal colonization was monitored. Swabs were collected in the schools approximately every month, for a five month period. The mean time lag between consecutive swabs was 37 days (standard deviation: 15 days). During the observation period, 9,857 swabs were collected for serotyping. Among the 4,488 3–6yearold children attending the schools (88% of the 3–6yearold children in the area under study), 2,445 (55%) gave at least one swab. Among children giving at least one swab, the mean number of swabs was four (minmax: one – five). All children attending the schools were included in the analysis, even those without a single observation during the followup.
Results of serotyping during followup of 4,488 children in 81 schools, JanuaryMay 2000, France.
Frequency (%)  

Serotypes with at least 30 swabs  
Vaccine serotypes†  19F  403 (4.1) 
6B  313 (3.2)  
23F  274 (2.8)  
14  191 (1.9)  
18C  79 (0.8)  
9V  79 (0.8)  
Nonvaccine serotypes  6A  281 (2.9) 
3  238 (2.4)  
19A  212 (2.2)  
11A  80 (0.8)  
15A  72 (0.7)  
15B  54 (0.5)  
23A  46 (0.5)  
17F  44 (0.4)  
10A  42 (0.4)  
9L  32 (0.3)  
Other serotypes  247 (2.5)  
Noncolonized  7,170 (72.7)  
Total  9,857 (100.0) 
The only currently available pneumococcal conjugate vaccine is the sevenvalent vaccine, which includes serotypes 4, 6B, 9V, 14, 18C, 19F and 23F. Hence, we divided the serotypes selected for analysis into two groups: vaccine serotypes, and nonvaccine serotypes. The two groups are detailed in table 1. Vaccine serotype 4 was not selected for analysis because it was isolated in 10 swabs only. At the time of the study, conjugate vaccine had not been introduced in France, so that children participating in the study were all unvaccinated.
Transmission model
Here, we present a dynamic model for S. pneumoniae transmission in schools. The model detailed the individual rate of colonization and decolonization by S. pneumoniae for all children attending the schools (even those that gave no swab at all). In the subsequent section "Estimation of transmission parameters", we show how model parameters may be estimated from the data.
In the model, we assumed that dual colonization was impossible (a child colonized by one serotype may not acquire another unless first clearing). Noncolonized children could be colonized within their school (see § Withinschool acquisition), or in the community (i.e. out of the school, see § Other assumptions).
Duration of carriage
We assumed that the duration of carriage of a given serotype had an Exponential distribution with mean μ _{ V }for vaccine serotypes and μ _{ U }for nonvaccine serotypes. This particular distribution is in good agreement with observed data [10].
Withinschool acquisition
We considered the school as a dynamic environment, i.e the number C _{ s }(t) of children colonized by serotype s at time t was a function of time. We assumed that, for each child who was not colonized at time t, the individual rate to acquire serotype s in the school at time t was i) proportional to C _{ s }(t), since this increased opportunities for transmission; and ii) inversely proportional to the size n of the school to allow for reduced frequency of contacts in each pair of children in larger schools [16]. With these assumptions, the individual rate to acquire serotype s at time t was β C _{ s } (t)/n, where n was the size of the school (including children that gave no swab at all) and β corresponded to pairwise childtochild transmission rate, irrespective of the size n of the school. Note that this formulation leads to an aggregate rate of colonization due to intraschool transmission of β C _{ s } (t) S(t)/n, where S(t) is the number of non colonized children, in agreement with the standard SusceptibleInfectiousSusceptible model [17]. We denoted the child to child transmission rate β _{ V }for vaccine serotypes, and β _{ U }for nonvaccine serotypes.
Other assumptions
We assumed that, during the 5months followup, a noncolonized child was exposed to a rate α _{ s }to acquire serotype s in the community, with α _{ s }constant over time. It was also necessary to define the probability that a child carried serotype s at the beginning of the followup by π _{ s }. The values of α _{ s }and π _{ s }were serotype dependent. Parameters α _{ s }and π _{ s }will be considered according to two clusters to be learnt from the data through the model.
Characterizing transmission from the model
The model was used to calculate the monthly probability to acquire a serotype in the community, which is 1exp(α _{ s } 30) for serotype s. It was also possible to calculate the Secondary Attack Rate (SAR), defined here as the probability that a colonized child transmits the bacteria to a noncolonized child of his/her school for a colonization event. The Secondary Attack Rate combines information on the mean duration of carriage μ and the child to child transmission rate β, and therefore allowed investigation of whether vaccine serotypes had 'globally' a better fitness for transmission than nonvaccine serotypes. Assuming that there is no transmission beyond the secondary case and that there is at maximum one acquisition of carriage per susceptible during the one month period, the probability of transmission between two children of a school of size n is:
$SAR={\displaystyle {\int}_{0}^{\infty}\left(1\mathrm{exp}(\beta L/n)\right)f(L)dL}$
where L is the duration of carriage with density f(L)=exp(L/μ)/μ and 1  exp( β L/n) is the probability that a child colonized during time period L transmits to a noncolonized child of the school. Eventually, the Secondary Attack Rate was equal to (1+Nμ ^{1} β ^{1})^{1} within a school of size n.
Estimation of transmission parameters
Data augmentation
Conditional on these augmented times and model parameters, the likelihood of the data is available, but since there is no unique way to choose the augmented times given the observation, a systematic exploration of the augmented times is necessary for inference. This is performed by Markov chain Monte Carlo (MCMC) sampling. In the algorithm, augmented periods of carriage may be added/suppressed, split/combined; and the times of colonization/decolonization may change [13]. Figure 2 presents two augmented trajectories that are compatible with the observation and could be explored by our algorithm.
Heterogeneity in the community acquisition rates of serotypes and in the probabilities of carriage at the beginning of followup
Heterogeneity was allowed in the community acquisition rates of serotypes and in the probabilities of carriage at the beginning of followup. More precisely, we allowed serotypes to cluster in two subgroups, not fixed in advance, for each of these parameters. In the MCMC sampling scheme, changes in serotypes allocation were proposed independently for community acquisition rate and probability of carriage at the beginning of the followup.
Bayesian hierarchical framework
The statistical framework has a Bayesian hierarchical structure [18], with 3 levels:
A) The observation level ensures that the augmented data are consistent with the observation;
B) The transmission level describes the latent transmission process (transmission model);
c) The prior level specifies the prior distributions of the parameters: For the child to child transmission rate β, the community acquisition rate α and the mean duration of carriage μ, we specified vague flat priors consisting of Exponential distributions with means 10^{5} day^{1}, 10^{5} day^{1} and 10^{5} days, respectively. The prior distribution for the probability π of carriage at the beginning of the followup was uniform from 0 to 1.
MCMC implementation
The MCMC algorithm was developed in C; the output was analyzed with R software. The seed used in the simulations was given by the computer clock. We performed 1,000,000 iterations for each run of the MCMC algorithm. The first 500,000 were discarded as the burnin period. The output was then recorded once every 5 iterations to constitute a sample from the posterior distribution. The convergence was visually assessed, and tested with the Geweke criteria [19]. We also checked that estimates were robust to a change in the initial values.
Eventually, the joint posterior distribution of augmented data, clusters of serotypes and parameters was explored by MCMC sampling, and characterized by means and equaltailed 95% credible intervals (CI). For the community acquisition rates and the probabilities of carriage at the beginning of the followup, reported results correspond to the partition of serotypes with the largest posterior probability.
Results
Posterior mean and 95% credible interval of transmission parameters of S. pneumoniae according to serotype inclusion in the vaccine. Transmission parameters consist of the mean duration of carriage μ, the child to child transmission rate β and the Secondary Attack Rate (SAR).
Vaccine serotype*  Nonvaccine serotype†  Ratio  

Mean  95% CI‡  Mean  95% CI‡  Mean  95% CI‡  
μ (days)  23  21, 25  22  20, 24  1.06  0.94, 1.18 
β (% day^{1})  4.6  4.2, 5.0  5.1  4.5, 5.6  0.91  0.80, 1.05 
SAR (%)  
n = 30 §  3.4  3.2, 3.7  3.6  3.3, 3.8  0.97  0.88, 1.06 
n = 50 §  2.1  1.9, 2.2  2.2  2.0, 2.3  0.97  0.88, 1.06 
n = 100 §  1.1  1.0, 1.2  1.1  1.0, 1.2  0.97  0.88, 1.06 
Durations of carriage
The mean duration of carriage was μ _{ V }= 23 (95% CI: 21, 25) days for vaccine serotypes and μ _{ U }= 22 (95% CI: 20, 24) days for nonvaccine serotypes, leading to a relative mean duration of carriage (μ _{ V } /μ _{ U }) equal to 1.06 (95% CI: 0.94, 1.18).
Colonization rate within a school
The child to child transmission rate was β _{ V }= 0.046 (95% CI: 0.042, 0.050) day^{1} for vaccine serotypes and β _{ U }= 0.051 (95% CI: 0.045, 0.056) day^{1} for nonvaccine serotypes, leading to a relative child to child transmission rate (β _{ V } /β _{ U }) equal to 0.91 (95% CI: 0.80, 1.05). In schools of size 30, 50, and 100, the Secondary Attack Rate was respectively 3.4% (95% CI: 3.2%, 3.7%), 2.1% (95% CI: 1.9%, 2.2%) and 1.1% (95% CI: 1.0%, 1.2%) for vaccine serotypes, and 3.6% (95% CI: 3.3%, 3.8%), 2.2% (95% CI: 2.0%, 2.3%) and 1.1% (95% CI: 1.0%, 1.2%) for nonvaccine serotypes. The relative Secondary Attack Rate (SAR _{ V } /SAR _{ U }) was 0.97 (95% CI: 0.88, 1.06) irrespective of the size of the school.
Other parameters
The posterior probability of the partition with largest support was 84% for the community acquisition rate and 92% for the probability of carriage at the beginning of the followup. Serotypes clustered identically regarding community acquisition rate and probability of carriage at the beginning of the followup. The first cluster of serotypes contained 9 serotypes (23A, 17F, 10A, 15B, 9V, 11A, 18C, 15A, 9L). It was characterized by a monthly probability of colonization in the community estimated at 0.40% (95% CI: 0.32%, 0.47%) per serotype, and a probability of carriage at the beginning of the followup equal to 0.48% (95% CI: 0.40%, 0.57%) per serotype. The second cluster consisted of 7 serotypes (19A, 14, 3, 6A, 23F, 6B, 19F). The monthly probability to acquire one of these serotypes out of the school was 1.40% (95% CI: 1.20%, 1.65%) per serotype, with a probability of carriage at the beginning of the followup equal to 2.3% (95% CI: 2.1%, 2.5%) per serotype. Serotypes that were found to have high community acquisition rate/probability of carriage at the beginning of the followup were those that were the most prevalent in the schools (191–403 swabs, as opposed to 32–80 swabs for other serotypes).
Discussion
We found that, in 3–6 years old children, the mean duration of carriage and the child to child transmission rate of S. pneumoniae was the same for vaccine and nonvaccine serotypes, from the analysis of a large cohort of schoolchildren.
Carriage was determined by oropharyngeal swabbing at each visit. This procedure is known to be less sensitive than nasopharyngeal swabbing [20]. Therefore, some samples could be falsely negative, and lead to underestimate the duration of carriage. However, our estimates of the mean durations of carriage (about 20 days for 3–6 years old children) are consistent with those reported by Ekdahl et al. [10]. Melegaro et al. [12] also found a carriage duration of 20 days for children ≥ 5 years old, but of 50 days for children <5 years old. This difference could be explained by the presence of children under 2 years old, who are known to have the largest durations of carriage [10] and were not covered by our study, and also by a coarser characterization of serotypes in Melegaro et al. [12]. A larger estimate of the duration of carriage (45 days) was found by Auranen et al. [13] for subjects older than 2 years old. Here again, a possible difference in the age distributions of the population under study might explain the differences in the estimates.
Here, we detected no difference between vaccine and nonvaccine serotypes for durations of carriage, child to child transmission rates and Secondary Attack Rates among school age children. These results are consistent with those of Ekdahl et al. [10] who found no significant difference between durations of carriage of serogroups when adjusted for age; but contrast with those of Smith et al. [11] who concluded that there was a factor of 10 between the mean durations of carriage of the most/least persistent serotypes. Swabs were collected weekly in [10]; and monthly in [11] and in our study (80% of intervals between paired nasal swab specimens were larger than ≥ 28 days in the study of Smith et al.; this proportion is even larger in our study). Our approach could be used to reanalyze the data of Smith et al [11] and investigate whether their results are robust to: i) the possible occurrence of decolonization/recolonization between consecutive positive swabs; ii) the inclusion of microepidemics in the stochastic transmission model. Since sensitivity of oropharyngeal swabbing is not likely to be serotypedependent, differences between vaccine and nonvaccine durations of carriage should not be affected if durations of carriage were underestimated. Simulation studies showed that we could have detected differences between durations of carriage of the order of one week, but also that the sampling interval (approximately 1 month in our study) was sufficient to estimate durations of carriage in the considered range (≥18 days). It would be important to investigate if it is still the case with shorter durations of carriage.
Infants are known to be the main reservoir in maintaining transmission. Consequently, the differences we found in the community acquisition rates could proceed from heterogeneity in serotypespecific prevalence in infants. Our data does not allow to judge whether serotypespecific prevalence in infants could be compatible with the absence of differences in the duration of carriage and in the child to child transmission rate (like it is for 3–6 years old children).
While the absence of difference between vaccine and nonvaccine serotypes seems at odds with the large differences observed in serotypes prevalence, mathematical models have shown that, in the context of direct competition with short term immunity [6–9] or in the context of indirect competition with long term immunity [6], the difference in serotypes prevalence could be much larger than what would be expected in the absence of competition. For example, in [6], competition between 2 serotypes with close transmission characteristics (basic reproduction number 2.2 and 1.8 respectively) led to extinction of the serotype with reduced transmissibility. Temime et al. [9] have shown that competition between serotypes was sufficient to generate large differences in prevalence according to serotypes, even in the absence of difference in transmission characteristics. While different modes of competition may influence serotypes equilibrium prevalence, replacement may occur only in case of direct competition in the absence of longterm immunity [6]. In this context, the extent of replacement depends on the degree of competition, but also on serotype specific transmission characteristics, which determine the ability of non vaccine serotypes to recolonize the ecological niche released by vaccine serotypes.
In order to explain differences in the prevalence of carriage, heterogeneity is required in at least one parameter among the duration of carriage, the child to child transmission rate and the community acquisition rate. In this respect, if duration of carriage and child to child transmission rate had been fixed to the same value for all serotypes, differences in the community acquisition rate would have naturally shown to match the differences in serotype specific prevalence in the data. However, our purpose was here to determine whether the data could indicate which differences were the most likely. To that end, we imposed that the mean duration and the child to child transmission rate should be the same among vaccine serotypes and among non vaccine serotypes, but allowed for differences between these two classes. In this framework, it is possible to obtain differences in the mean duration of carriage, or child to child transmission rate, between vaccine and nonvaccine serotypes, irrespective of a difference in the community acquisition rate. We have found, using simulations, that the data and procedure were informative and selective enough to do this: changes in the mean duration of carriage, in the presence of heterogeneity in the community acquisition rate, were identifiable in the framework described.
In this paper, we focused on pneumococcal carriage and not on pneumococcal acute otitis media or invasive disease. While our results favour the hypothesis that replacement should follow vaccination, the overall impact on pneumococcal disease is less clear. In a vaccinated population the othopathogenic capacities of the replacing nonconjugate vaccine type pneumococci have been established [21, 22]. For invasive diseases, some serotypes not included in the vaccine appear to be less invasive [14] and a reduction in incidence could follow the introduction of the vaccine despite replacement in carriage.
We considered a model in which competition was direct and prevented dual colonization. This last assumption is strong since dual colonization is possible [23, 24]. However, given the nature of our data in which at most 1 serotype may be detected per swab, this was the only sensible approach. If information on dual colonization were available in the data, the approach could easily be extended to estimate competition parameters as listed by Lipsitch [7]. Our approach, similar to that of Auranen et al. [13], amounted to consider that in case of dual colonization in one individual there was a "major" serotype that would be transmitted. The results should be robust to this assumption unless there was a significant transmission of "minor" serotypes in case of dual transmission.
In our analysis, we investigated transmission of the 16 serotypes that were the most prevalent in the dataset. We assumed that individuals carrying less prevalent serotypes were susceptible to colonization by serotypes under study (absence of direct competition). This may have an effect on the duration of carriage estimates as out competition may have occurred rather than clearance of the serotype that was originally carried. However, the effect should be limited since serotypes included in the analysis represent 93% of overall carriage.
Heterogeneity in community acquisition rates and probabilities of carriage at the beginning of the followup was explored by allowing free cluster formation according to these parameters. The number of clusters was fixed (= 2), but the composition of the clusters was determined by the MCMC algorithm. With a number of clusters fixed to 3 or 4, the estimates of the mean durations of carriage and the child to child transmission rates remained unchanged.
Conclusion
In conclusion, this study supports that, in 3–6 years old children, no competitive advantage exists for vaccine serotypes compared to nonvaccine serotypes. This is an argument in favour of important serotype replacement. It would be important to validate the result for infants, who are known to be the main reservoir in maintaining transmission. Overall reduction in pathogenicity should also be taken into account in forecasting the future burden of pneumococcal colonization in vaccinated populations.
Appendix: Simulations
Seven hundred epidemics were simulated from the model, with transmission parameters for the serotypes (child to child transmission rate β, community acquisition rate α and mean duration of carriage μ for each serotype) drawn from the posterior distribution. The probability of carriage at time 0 was not sampled from the posterior distribution. Rather, we let the simulation start in the remote past (t = 1000) with no carriage in the school. The distribution of carriage at time 0 is therefore only dependent on the transmission parameters β, α and μ. The simulations were designed so that field and simulated data had the same structure: same number of schools and same number of children per school; same number of swabs for children j of schools i of both data, collected at the same times.
Abbreviations
 CI:

Credible Interval
 SAR :

Secondary Attack Rate
 S. pneumoniae :

Streptococcus pneumoniae
 MCMC:

Markov chain Monte Carlo.
Declarations
Acknowledgements
The project was supported, in part, by grants from Aventis Company, the "Fédération Nationale de la Mutualité Française" (FNMF) and the European Community (EURIS project).
Authors’ Affiliations
References
 Butler JC, Breiman RF, Lipman HB, Hofmann J, Facklam RR: Serotype distribution of Streptococcus pneumoniae infections among preschool children in the United States, 19781994: implications for development of a conjugate vaccine. J Infect Dis. 1995, 171 (4): 885889.View ArticlePubMedGoogle Scholar
 Dagan R, GivonLavi N, Zamir O, SikulerCohen M, Guy L, Janco J, Yagupsky P, Fraser D: Reduction of nasopharyngeal carriage of Streptococcus pneumoniae after administration of a 9valent pneumococcal conjugate vaccine to toddlers attending day care centers. Journal of Infectious Diseases. 2002, 185 (7): 927936. 10.1086/339525.View ArticlePubMedGoogle Scholar
 Whitney CG, Farley MM, Hadler J, Harrison LH, Bennett NM, Lynfield R, Reingold A, Cieslak PR, Pilishvili T, Jackson D, Facklam RR, Jorgensen JH, Schuchat A, Active Bacterial Core Surveillance of the Emerging Infections Program N: Decline in invasive pneumococcal disease after the introduction of proteinpolysaccharide conjugate vaccine.[see comment]. New England Journal of Medicine. 2003, 348 (18): 17371746. 10.1056/NEJMoa022823.View ArticlePubMedGoogle Scholar
 Schutze GE, Tucker NC, Mason EOJ: Impact of the Conjugate Pneumococcal Vaccine in Arkansas. Pediatr Infect Dis J. 2004, 23 (12): 11251129.PubMedGoogle Scholar
 Pelton SI, Loughlin AM, Marchant CD: Seven valent pneumococcal conjugate vaccine immunization in two Boston communities: changes in serotypes and antimicrobial susceptibility among Streptococcus pneumoniae isolates. Pediatr Infect Dis J. 2004, 23 (11): 10151022.View ArticlePubMedGoogle Scholar
 Zhang Y, Auranen K, Eichner M: The influence of competition and vaccination on the coexistence of two pneumococcal serotypes. Epidemiol Infect. 2004, 132 (6): 10731081. 10.1017/S0950268804002894.View ArticlePubMedPubMed CentralGoogle Scholar
 Lipsitch M: Vaccination against colonizing bacteria with multiple serotypes. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94 (12): 65716576. 10.1073/pnas.94.12.6571.View ArticlePubMedPubMed CentralGoogle Scholar
 Temime L, Boelle PY, Courvalin P, Guillemot D: Bacterial resistance to penicillin G by decreased affinity of penicillinbinding proteins: a mathematical model. Emerging Infectious Diseases. 2003, 9 (4): 411417.View ArticlePubMedPubMed CentralGoogle Scholar
 Temime L, Guillemot D, Boelle PY: Short and longterm effects of pneumococcal conjugate vaccination of children on penicillin resistance. Antimicrob Agents Chemother. 2004, 48 (6): 22062213. 10.1128/AAC.48.6.22062213.2004.View ArticlePubMedPubMed CentralGoogle Scholar
 Ekdahl K, Ahlinder I, Hansson HB, Melander E, Molstad S, Soderstrom M, Persson K: Duration of nasopharyngeal carriage of penicillinresistant Streptococcus pneumoniae: experiences from the South Swedish Pneumococcal Intervention Project. Clinical Infectious Diseases. 1997, 25 (5): 11131117.View ArticlePubMedGoogle Scholar
 Smith T, Lehmann D, Montgomery J, Gratten M, Riley ID, Alpers MP: Acquisition and invasiveness of different serotypes of Streptococcus pneumoniae in young children. Epidemiol Infect. 1993, 111 (1): 2739.View ArticlePubMedPubMed CentralGoogle Scholar
 Melegaro A, Gay NJ, Medley GF: Estimating the transmission parameters of pneumococcal carriage in households. Epidemiology & Infection. 2004, 132 (3): 433441. 10.1017/S0950268804001980.View ArticleGoogle Scholar
 Auranen K, Arjas E, Leino T, Takala AK: Transmission of pneumococcal carriage in families: a latent Markov process model for binary longitudinal data. J Amer Statist Assoc. 2000, 95 (452): 10441053.View ArticleGoogle Scholar
 Brueggemann AB, Griffiths DT, Meats E, Peto T, Crook DW, Spratt BG: Clonal relationships between invasive and carriage Streptococcus pneumoniae and serotype and clonespecific differences in invasive disease potential. J Infect Dis. 2003, 187 (9): 14241432. 10.1086/374624.View ArticlePubMedGoogle Scholar
 Guillemot D, Varon E, Bernède C, Weber P, Henriet L, Simon S, Laurent C, Lecoeur H, Carbon C: Reduction of antibiotic use in the community reduces the rate of Penicillin GNonSusceptible Streptococcus pneumoniae. Clinical Infectious Diseases.
 De Jong MCM, Diekmann O, Heesterbeek H: How does transmission of infection depend on population size?. Epidemic Models: Their structure and relation to data models. Edited by: Mollison D. 1995, Cambridge , Cambridge University Press, 8494.Google Scholar
 Anderson RM, May RM: Infectious Diseases of Humans: Dynamics and Control. 1991, Oxford , Oxford University PressGoogle Scholar
 Robert C: The Bayesian choice: from decisiontheoretic motivations to computational implementation. 2001, NewYork , SpringerVerlag, 2ndGoogle Scholar
 Geweke J: Evaluating the accuracy of samplingbased approaches to calculating posterior moments. Bayesian Statistics 4. Edited by: Bernardo JM, Smith AFM, Dawid AP, Berger JO. 1992, New York , Oxford University Press, 169193.Google Scholar
 Rapola S, Salo E, Kiiski P, Leinonen M, Takala AK: Comparison of four different sampling methods for detecting pharyngeal carriage of Streptococcus pneumoniae and Haemophilus influenzae in children. Journal of Clinical Microbiology. 1997, 35 (5): 10771079.PubMedPubMed CentralGoogle Scholar
 Eskola J, Kilpi T, Palmu A, Jokinen J, Haapakoski J, Herva E, Takala A, Kayhty H, Karma P, Kohberger R, Siber G, Makela PH: Efficacy of a pneumococcal conjugate vaccine against acute otitis media. N Engl J Med. 2001, 344 (6): 403409. 10.1056/NEJM200102083440602.View ArticlePubMedGoogle Scholar
 Kilpi T, Ahman H, Jokinen J, Lankinen KS, Palmu A, Savolainen H, Gronholm M, Leinonen M, Hovi T, Eskola J, Kayhty H, Bohidar N, Sadoff JC, Makela PH: Protective efficacy of a second pneumococcal conjugate vaccine against pneumococcal acute otitis media in infants and children: randomized, controlled trial of a 7valent pneumococcal polysaccharidemeningococcal outer membrane protein complex conjugate vaccine in 1666 children. Clin Infect Dis. 2003, 37 (9): 11551164. 10.1086/378744.View ArticlePubMedGoogle Scholar
 Obaro SK, Adegbola RA, Banya WA, Greenwood BM: Carriage of pneumococci after pneumococcal vaccination. Lancet. 1996, 348 (9022): 271272. 10.1016/S01406736(05)655857.View ArticlePubMedGoogle Scholar
 Gratten M, Montgomery J, Gerega G, Gratten H, Siwi H, Poli A, Koki G: Multiple colonization of the upper respiratory tract of Papua New Guinea children with Haemophilus influenzae and Streptococcus pneumoniae. Southeast Asian Journal of Tropical Medicine & Public Health. 1989, 20 (4): 501509.Google Scholar
 The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712334/6/14/prepub
Prepublication history
Copyright
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.