Social clustering of unvaccinated children: measles-mumps-rubella vaccination coverage in schools in the Netherlands ==================================================================================================================== * Don Klinkenberg * Albert Jan van Hoek * Irene Veldhuijzen * Susan Hahné * Jacco Wallinga ## Abstract **Background** For the measles-mumps-rubella (MMR) vaccine, the WHO-recommended coverage for herd protection is 95% for measles and 80% for rubella and mumps. However, a national vaccine coverage does not reflect social clustering of unvaccinated children, e.g. in schools of Orthodox Protestant or Anthroposophic identity in The Netherlands. To fully characterise this clustering, we estimated one-dose MMR vaccination coverages at all schools in the Netherlands. **Methods** By combining postcode catchment areas of schools and school feeder data, each child in the Netherlands was characterised by residential postcode, primary and secondary school (referred to as school career). Postcode-level vaccination data were used to estimate vaccination coverages per school career. These were translated to coverages per school, stratified by school identity. **Results** Most schools had vaccine coverages over 99%, but major exceptions were Orthodox Protestant schools (63% in primary and 58% in secondary schools) and Anthroposophic schools (67% and 78%). **Conclusions** School-level vaccine coverage estimates reveal strong clustering of unvaccinated children. The school feeder data reveal strongly connected Orthodox Protestant and Anthroposophic communities, but separated from one another. This suggests that even at a national one-dose MMR coverage of 97.5%, thousands of children per cohort are not protected by herd immunity. Keywords * vaccine * vaccine-preventable disease * immunisation * herd immunity * childhood disease * infectious disease * clustering * social network ## Introduction The impact of childhood vaccination programmes on protecting children from infection can be measured in various ways. Common indicators are the remaining disease incidence and disease burden, and the national vaccination coverage [1]. For measles, WHO recommends a threshold value for the vaccination coverage of 95% for two doses [2], and for rubella and mumps of 80% [3]. Many EU countries have an explicit target coverage for the measles-mumps-rubella (MMR) vaccine of 95% [1]. However, outbreaks occur among unvaccinated children when they are clustered, even at high values for an average national vaccination coverage [4]. The Netherlands, in particular, have an average national vaccination coverage for MMR that is among the highest in Europe ([http://data.euro.who.int/cisid/](http://data.euro.who.int/cisid/)), but the substantial socio-geographic clustering of unvaccinated children results in repeated outbreaks of measles, mumps and rubella [5-7]. Clustering of unvaccinated children is often described geographically [8-13], but most contacts between children are made at school [14]. Indeed, daycare centres or schools play a major role during outbreaks in well-immunised countries [15-18]. In the Netherlands, the clustering of unvaccinated children at school is enhanced by a school system where parents are free to start and choose schools based on their own philosophy or conviction, resulting in a range of registered school identities. There are two social groups with their own school identities that are known for a lower vaccination coverage: Orthodox Protestants, who live spatially clustered in the so-called Bible Belt and have a vaccination coverage of about 60% [19, 20]; and Anthroposophics, who live spatially more dispersed and have school vaccination coverages ranging from 60% to 90% [21]. There is also a more diffuse group of people who fear side-effects and perceive the diseases as mild [22], but there is not one school identity linked to that group. Vaccination coverages by school are routinely collected in various countries, including Germany and the United States [17, 18, 23-25], but reporting is frequently done by averaging school coverages locally, by county or district [24, 26]. Whereas this is very helpful for signaling which regions are at risk of infection, an average or even a local average does not fully capture the variability between schools, which is essential to fully understand the risk of outbreaks [25]. For that, we require the complete frequency distribution of school vaccination coverages. The aim of our study was to estimate the MMR vaccination coverages for all primary and secondary schools in the Netherlands, and assess the variation in vaccination coverage in primary and secondary schools. We focussed on a single-dose rather than the WHO-recommended two-dose vaccination, because the second dose in the Netherlands is given only at nine years of age so that the one-dose coverage is a better indication of protection in primary schools. Also, a single-dose coverage better reflects willingness versus refusal to vaccinate, and shows where the completely unvaccinated children are. We discuss the implications of the observed variation for the risk of outbreaks in a population with high mean vaccination coverage. ## Methods Briefly, we assumed a steady state population of children characterised by residential postcode, primary school, and secondary school (school career). This population is constructed by use of school-level catchment data of the postcode areas in which school children live, and school feeder data of children moving from primary to secondary school. The steady-state assumption allowed us to use postcode-level coverage data for a single cohort to estimate coverages in both primary and secondary schools. We used a statistical hierarchical model, in which the vaccination status of a child is explained by the school career, with a fixed effect of the school identities (identity career) and random effects of the two schools. The school feeder data create dependencies of coverage estimates between primary and secondary schools sharing many children, which is especially relevant if those schools have different identities. ### Datasets We used five datasets to estimate school vaccination coverages (dark green in Figure 1), and two datasets for validation of estimated coverages (Utrecht and Anthroposophic school data). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/14/2022.05.12.22273875/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/05/14/2022.05.12.22273875/F1) Figure 1. Conceptual outline of the method to estimate school vaccination coverages (full details in Supplement S1). Numbers are fictional but mutually consistent. Green tables are data: dark green are original datasets, light green are linked data. Blue tables are results: dark blue are original estimates (logit-transformed in the actual analysis), light blue are derived quantities. Top (data): three population datasets are linked to build the population structure of numbers of children by school career and residential postcode area. Middle (model): regression model with vaccination data per postcode as dependent variable; the independent variables are the school careers of the children, mapped to postcode areas through function f(.), in which the population structure is used to map school career coverages to postcode area coverages, and in which the logit(coverage) of each school career is modelled as the sum of a fixed identity career effect and two random school effects. Bottom (results): the estimated coverages and deviations are used to calculate the coverage per school career, and combined with the population structure to calculate the coverages per school and per postcode. For estimation, we combined postcode-level MMR vaccination coverage data with postcode-level data of the schools that children attend (primary school and secondary school catchment data), school-level data on numbers of children moving from primary to secondary schools (school feeder data), and school identity data. The coverage data came from the Dutch vaccination registry Praeventis [27], of which we used the single MMR vaccination status in 2013 of 197,382 children of birth cohort 2003 at age 10. In the Netherlands, MMR is given at 14 months and 9 years of age, so focusing on a single-dose vaccination better reflects vaccine coverage across all ages in schools. The school catchment, feeder, and identity data came from the Education Executive Agency of the Ministry of Education [28], of which we used primary school catchment data from October 2012, secondary school catchment data from October 2013, and feeder data from September 2013. Schools were categorised into 15 groups based on their identity [28]. Most identities are related to religion, but exceptions are Municipal (“Openbaar”) which are run by the local government, and General (“Algemeen bijzonder”) which is a heterogeneous group of schools with various non-religious backgrounds (e.g. Montessori, Dalton, Jenaplan). Special types of identity were “Collaborations”, which consisted of schools categorised as collaboration (“Samenwerking”) of two or more identities among Protestant, Roman Catholic, General, and/or Municipal; and “Other”, which consisted of all schools categorized as Rest (“Overige”) in the dataset, or left blank. For validation we used data collected by the municipal health service Utrecht (Utrecht data), which covers part of the Bible Belt and consists of actual coverage data for 488 schools; and data from a study in 2014 on vaccination hesitancy on Anthroposophic schools (Anthroposophic school data), consisting of coverage data for 11 schools [21]. ### Statistical analysis Figure 1 shows how the data are used in a model to obtain the MMR coverage per school (more detail in Supplement S1). First, the school catchment data and feeder data were combined to build the population structure, in which each child is characterized by a postcode and a school career, i.e. a primary and secondary school, one of which the child actually attends and the other which it has attended or will attend. Then, this population structure is used with the vaccination data and school identities in a Bayesian regression analysis in stan ([mc-stan.org](http://mc-stan.org)), called from R statistical software [29] with the Rstan package (mc-stan.org/rstan). Conceptually, through the regression equation (Figure 1) the observed coverages in postcode areas are explained by the school careers of the children living in that postcode area, with a fixed effect due to the identities of the schools, and a random effect for each school. The regression analysis resulted in estimated vaccination coverages per school career, which were used to derive our results of interest: vaccination coverages per school. The analysis also provided estimated coverages per postcode area. Because there were 15 different school identities in the data, and therefore many more possible identity careers, we had to aggregate identities and identity careers to limit the number of variables in the regression. This was done in two preliminary analyses. The first was to define a “Rest” identity by grouping all identities without evidence of lower coverage than the major identities (Municipal, Protestant, and Roman Catholic). That resulted in five primary school identities and three secondary school identities. The second was to aggregate identity careers where possible, resulting in a final set of eight identity careers in the final analysis (see Supplement S1 for more detail). ## Results ### Aggregating school identities and identity careers We first ran two preliminary regression analyses (detail in Supplement S1) in order to reduce the number of school identities and identity careers in our final analysis. Most school identities were aggregated into one “Rest” group, but identities that remained separate were Anthroposophic and Orthodox Protestant in both primary and secondary schools, and General and Hindu in primary schools (Table 1). From these remaining identities, 15 identity careers were formed, of which eight were kept separate after a second aggregation step to obtain our final results (Table 2). View this table: [Table 1.](http://medrxiv.org/content/early/2022/05/14/2022.05.12.22273875/T1) Table 1. Identities of primary and secondary schools: numbers of schools per original identity, and identity after aggregation View this table: [Table 2.](http://medrxiv.org/content/early/2022/05/14/2022.05.12.22273875/T2) Table 2. Identity careers: numbers of children in feeder data, and aggregate careers after the second aggregation step (Ant = Anthroposophic, Gen = General, Hin = Hindu, Ort = Orthodox Protestant, Rest = Rest) ### Estimating school vaccination coverages The overall vaccination coverage in our dataset (one dose of MMR vaccine in birth cohort 2003 at the age of 10) was 97.5%. Figure 2 shows the comparison of estimated school vaccination coverages with the validation datasets. In all Utrecht schools as well as the Anthroposophic schools there is a clear correlation between true and estimated coverage, and almost all credible intervals include the true coverage. Two Rest schools with coverages below 90% have too high estimates; these are both Protestant schools. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/14/2022.05.12.22273875/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/05/14/2022.05.12.22273875/F2) Figure 2: validation plots of estimated vaccination coverages with 95% credible interval compared to actual vaccine coverage in Utrecht, 2013 (left) and on Anthroposophical schools in the Netherlands, 2012 (right). There is clear spatial clustering of unvaccinated children: most postcode areas had vaccination coverages around 99%, but about 4% of postcode areas had coverages below 95%, mainly in the Bible belt (Figures 3, 4). At school level, nonetheless, clustering was much more pronounced. The large majority of schools have have higher coverages than postcode area, above 99%, whereas specific schools have much lower coverages. In primary schools, vaccination coverage on Orthodox Protestant schools (62.8%) is slightly lower than on Anthroposophic schools (66.6%), whereas in secondary schools this difference is much larger (58.4% and 78.4%, respectively, Table 3). General primary schools have a slightly lower vaccination coverage than Rest primary schools (Table 3), but less than 1% of all these schools have a coverage below 95% (Figure 3). Coverages of Hindu schools are estimated to be lower than in the Rest or General schools, but this result is highly uncertain because of the small number of schools (Table 3). Of all identities, Orthodox Protestant schools have the largest variation in vaccine coverage (Figure 3). View this table: [Table 3.](http://medrxiv.org/content/early/2022/05/14/2022.05.12.22273875/T3) Table 3. Vaccination status by school identity: mean vaccination coverage at school level (95% credible interval), and number of unvaccinated children per cohort of 197,382 (95% credible interval across all schools). ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/14/2022.05.12.22273875/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/05/14/2022.05.12.22273875/F3) Figure 3: posterior median vaccination coverages of all schools, grouped by school identity. Displayed are medians (of the median coverages) with interquartile range of (IQR, boxes), and whiskers extending to 1.5 times the IQR or most extreme value, whichever is closest to the median. Dots are outliers. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/14/2022.05.12.22273875/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/05/14/2022.05.12.22273875/F4) Figure 4: geographic distributions of all schools by identity, with colour coded estimated vaccination coverages. The bottom right map shows the estimated MMR coverages by four-digit postcode area. Figure 4 shows maps of the school locations, partitioned by school identity, colour coded to indicate estimated vaccination coverage, and as a reference a postcode map of estimated vaccination coverages. The Orthodox Protestant schools are clearly concentrated in the Bible Belt, with most of the very-low coverage schools east from the centre. The lower-coverage Rest schools are also predominantly located in the Bible Belt area; most of these have the Protestant identity (not shown). Anthroposophic schools are more widely distributed across the country. ### Clustering of unvaccinated children In our age cohort dataset of 197,382 10-year old children in the Netherlands in 2013, 4998 (2.5%) were not vaccinated against MMR. In primary schools, 51% of unvaccinated children attend General or Rest schools most of which have a very high vaccination coverage, but the rest is clustered in schools with much lower coverages: 38% in Orthodox Protestant schools, 10% in Anthroposophic schools, and 1% in Hindu schools (Figure 3, Table 3). When these children go to secondary school, this hardly changes: 56% go to Rest, 39% to Orthodox Protestant, and 5% to Anthroposophic schools. Many children move to secondary schools of the same identity as their primary school, more so in the Orthodox Protestant community (85%) than in the Anthroposophic community (44%), showing that these children go to low-coverage schools during their whole school career. There were, however, no children moving between schools with these two identities. ## Discussion Our aim was to estimate single-dose MMR vaccination coverages on schools in the Netherlands, and assess clustering of unvaccinated children in schools. It turns out that about half of the unvaccinated children are clustered into a small minority of schools with very low vaccination coverage. These are almost all schools with an Orthodox Protestant or Anthroposophic identity. That implies that thousands of children per cohort are not protected by vaccination, even when the mean vaccination coverage of 97.5% is well above WHO recommendations for both infections. Spatial clustering of unvaccinated children in the Bible belt is clearly observed at the postcode area level (Figure 4), and has been described before. Our results, however, show that clustering is much stronger at the level of schools, which better reflects the social environment of children. Most schools have vaccination coverages higher than in most postcode areas, whereas some schools of mainly Orthodox Protestant and Anthroposophic identity have much lower coverages (Figure 3). The analysis also revealed social clustering outside the Bible belt area (Figure 4). Apart from the good agreement with the validation data, the Orthodox Protestant coverages also confirm estimates of 60% in this community from online surveys by Ruijs et al [19]. Whereas Orthodox Protestant and Anthroposophical schools are only minority identities, about half of all unvaccinated children attend these schools with low vaccination coverage (Table 3). Children not only spend much of their time at school, their social network outside school is probably very similar. That creates an infection risk for many unvaccinated children, who seem to be protected when considering their residential four-digit postcode area (Figure 5), let alone the national vaccination coverage. However, for virus introductions to cause large epidemics, there should also be sufficient contacts between the schools. Most of these contacts are probably through households, mainly between primary and secondary schools. Aside from the geographic proximity (Figure 2), the school identity careers (Table 2) show that the Orthodox Protestant community is closely knit: 85% of children attending Orthodox Protestant primary schools, also attend secondary school with that identity, and vice versa. In addition, families in the Orthodox Protestant community tend to be large. This explains the occasional large epidemics mostly within this community. Because secondary schools are essential links in the school network, large epidemics may only occur when the interepidemic period is long enough so that these schools are populated with susceptible children. This is consistent with the pattern of measles epidemics in the Bible Belt in the Netherlands, which took place in 1988, 1999, and 2013 [7, 30]. Anthroposophic schools are less well connected, firstly because they are more evenly distributed geographically (Figure 2), and secondly because more children going to either an Anthroposophic primary or an Anthroposophic secondary school but not both (Table 2). Indeed, outbreaks in Anthroposophic schools tend to stay confined to a few schools [31]. Also, the school feeder data suggest that the two low-coverage communities are socially separated, as there were no children in the dataset moving between Orthodox Protestant and Anthroposophic schools (Table 2). As a consequence, during the epidemics in the Orthodox Protestant community, Anthroposophic schools were hardly affected [7]. The estimated school vaccination coverages captured the trend of the validation datasets well, but were too unreliable to report for individual schools (Figure 2). First, they were unprecise, which is reflected by the wide credible intervals, and caused by the large number of schools compared to postcode area coverage data. Second, in a few cases estimates were inaccurate, notably for two schools of the Rest identity (Figure 2), which were Protestant by their registered identity (before aggregation of identities). Misclassification of some Protestant school is likely: the websites of these schools reveal an Orthodox identity, and most low-coverage Rest schools in the Bible Belt (Figure 4) are also Protestant. Third, confounding may have resulted in the estimated low coverages of the Hindu schools, as a lower coverage was not associated with the Hindu religious group in a large cross-sectional survey [32], in which 32 out of 35 Hindu respondents of age 1.5-20 were vaccinated (Mollema, personal communication). Fourth, we used data of a single birth cohort, which may not be representative if the vaccination coverage has increased or decreased, and which may overestimate coverage in younger children if they receive their first MMR vaccination at the age of nine, when the routine second vaccination is offered. If reliable coverages could be obtained, e.g. by coupling of school and vaccination registries, they could be used to inform parents, to prioritize supplemental vaccinations during outbreaks, or for more detailed studies into changes in vaccine hesitancy in different social groups, representated by school identity. Pockets of low vaccination coverage and resulting risk of outbreaks are better identified in social than geographic environments [33]. We assessed heterogeneity in vaccine coverage by estimating coverages at schools. The reason is that schools are good reflections of the social environment of children, especially if parents can choose a school with a clear identity like in the Netherlands. For the Netherlands, we were able to estimate school-level vaccination coverages by relating available postcode-level coverages to schools through the school catchment data, and we could use school identity to explain much of the variation between schools. In combination with the school feeder data, we could distinguish two well-connected by mutually separated clusters of schools with low vaccination coverage. For other countries, the number of postcode areas may be larger than the number of schools, which would simplify analyses, or alternative variables such as public and private schools may be more useful. Ideally, school vaccination coverage data are directly available: for instance, in the United States these are available at school level for some states [23], but often grouped by county for analysis [26]; also in Germany they are collected at school entry and then grouped by district [24]. For describing connections between schools, feeder data may be supplemented (or replaced) by individual level household data to build larger networks of schools and identify communities at risk. ## Supporting information Supplemental Methods and Figures [[supplements/273875_file03.pdf]](pending:yes) ## Data Availability Some data are privacy sensitive because of small numbers and may not automatically be shared. Please contact the authors to explore what is possible. ## Footnotes * **Conflict of interest:** The authors declare no conflicts of interest * **Funding:** No external funding has been used for this study * Received May 12, 2022. * Revision received May 12, 2022. * Accepted May 14, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Sheikh S, Biundo E, Courcier S, Damm O, Launay O, Maes E, et al. A report on the status of vaccination in Europe. Vaccine. 2018;36(33):4979–92. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) 2. 2.World Health Organization. Measles vaccines: WHO position paper - April 2017. Weekly epidemiological record. 2017;92(17):205–27. 3. 3.World Health Organization. Rubella vaccines: WHO position paper. Weekly epidemiological record. 2011;86(29):301–16. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21766537&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) 4. 4.Truelove SA, Graham M, Moss WJ, Metcalf CJE, Ferrari MJ, Lessler J. Characterizing the impact of spatial clustering of susceptibility for measles elimination. Vaccine. 2019;37(5):732–41. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) 5. 5.Hahne S, Macey J, van Binnendijk R, Kohl R, Dolman S, van der Veen Y, et al. Rubella outbreak in the Netherlands, 2004-2005: high burden of congenital infection and spread to Canada. Pediatr Infect Dis J. 2009;28(9):795–800. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/INF.0b013e3181a3e2d5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19710586&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269407500008&link_type=ISI) 6. 6.Wielders CC, van Binnendijk RS, Snijders BE, Tipples GA, Cremer J, Fanoy E, et al. Mumps epidemic in orthodox religious low-vaccination communities in the Netherlands and Canada, 2007 to 2009. Euro Surveill. 2011;16(41). 7. 7.Woudenberg T, van Binnendijk RS, Sanders EA, Wallinga J, de Melker HE, Ruijs WL, et al. Large measles epidemic in the Netherlands, May 2013 to March 2014: changing epidemiology. Euro Surveill. 2017;22(3). 8. 8.Brownwright TK, Dodson ZM, van Panhuis WG. Spatial clustering of measles vaccination coverage among children in sub-Saharan Africa. BMC Public Health. 2017;17(1):957. 9. 9.Hagemann C, Streng A, Kraemer A, Liese JG. Heterogeneity in coverage for measles and varicella vaccination in toddlers - analysis of factors influencing parental acceptance. BMC Public Health. 2017;17(1):724. 10. 10.Lieu TA, Ray GT, Klein NP, Chung C, Kulldorff M. Geographic clusters in underimmunization and vaccine refusal. Pediatrics. 2015;135(2):280–9. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBlZGlhdHJpY3MiO3M6NToicmVzaWQiO3M6OToiMTM1LzIvMjgwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDUvMTQvMjAyMi4wNS4xMi4yMjI3Mzg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 11. 11.Schlenker TL, Bain C, Baughman AL, Hadler SC. Measles herd immunity. The association of attack rates with immunization rates in preschool children. JAMA. 1992;267(6):823–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.1992.03480060069032&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1732654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992HC48400027&link_type=ISI) 12. 12.Thompson KM, Logan GE, Florida SRT. Characterization of Heterogeneity in Childhood Immunization Coverage in Central Florida Using Immunization Registry Data. Risk Anal. 2016;36(7):1418–26. 13. 13.Spencer CN, Delamater PL. Examining vaccination coverage in Germany: spatiotemporal clustering of MMR coverage, 2008-14. Eur J Public Health. 2020;30(5):993–5. 14. 14.Xiao X, van Hoek AJ, Kenward MG, Melegaro A, Jit M. Clustering of contacts relevant to the spread of infectious disease. Epidemics. 2016;17:1–9. 15. 15.Barrabeig I, Rovira A, Munoz P, Batalla J, Rius C, Sanchez JA, et al. MMR vaccine effectiveness in an outbreak that involved day-care and primary schools. Vaccine. 2011;29(45):8024–31. 16. 16.Sugerman DE, Barskey AE, Delea MG, Ortega-Sanchez IR, Bi D, Ralston KJ, et al. Measles outbreak in a highly vaccinated population, San Diego, 2008: role of the intentionally undervaccinated. Pediatrics. 2010;125(4):747–55. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBlZGlhdHJpY3MiO3M6NToicmVzaWQiO3M6OToiMTI1LzQvNzQ3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDUvMTQvMjAyMi4wNS4xMi4yMjI3Mzg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 17. 17.Sinclair DR, Grefenstette JJ, Krauland MG, Galloway DD, Frankeny RJ, Travis C, et al. Forecasted Size of Measles Outbreaks Associated With Vaccination Exemptions for Schoolchildren. JAMA Netw Open. 2019;2(8):e199768. 18. 18.Gromis A, Liu KY. Spatial Clustering of Vaccine Exemptions on the Risk of a Measles Outbreak. Pediatrics. 2022;149(1). 19. 19.Ruijs WL, Hautvast JL, van Ansem WJ, Akkermans RP, van’t Spijker K, Hulscher ME, et al. Measuring vaccination coverage in a hard to reach minority. Eur J Public Health. 2012;22(3):359–64. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eurpub/ckr081&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21715468&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000304529400014&link_type=ISI) 20. 20.Spaan DH, Ruijs WLM, Hautvast JLA, Tostmann A. Increase in vaccination coverage between subsequent generations of orthodox Protestants in The Netherlands. Eur J Public Health. 2017;27(3):524–30. 21. 21.Klomp JH, van Lier A, Ruijs WL. Vaccination coverage for measles, mumps and rubella in anthroposophical schools in Gelderland, The Netherlands. Eur J Public Health. 2015;25(3):501–5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eurpub/cku178&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25406448&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) 22. 22.Harmsen IA, Ruiter RAC, Paulussen TGW, Mollema L, Kok G, De Melker HE. Factors that influence vaccination decision-making by parents who visit an enthroposophical child welfare center: a focus group study. Advances in Preventive Medicine. 2012;2012:175694. 23. 23.Leslie TF, Street EJ, Delamater PL, Yang YT, Jacobsen KH. Variation in Vaccination Data Available at School Entry Across the United States. Am J Public Health. 2016;106(12):2180–2. 24. 24.Weigel M, Bruns R, Weitmann K, Hoffmann W. Immunization rates at the school entry in 2012. Dtsch Arztebl Int. 2014;111(46):788–94. 25. 25.Morrison M, Castro LA, Ancel Meyers L. Conscientious vaccination exemptions in kindergarten to eighth-grade children across Texas schools from 2012 to 2018: A regression analysis. PLoS Med. 2020;17(3):e1003049. 26. 26.Kluberg SA, McGinnis DP, Hswen Y, Majumder MS, Santillana M, Brownstein JS. County-level assessment of United States kindergarten vaccination rates for measles mumps rubella (MMR) for the 2014-2015 school year. Vaccine. 2017;35(47):6444–50. 27. 27.van Lier A, Oomen P, de Hoogh P, Drijfhout I, Elsinghorst B, Kemmeren J, et al. Praeventis, the immunisation register of the Netherlands: a tool to evaluate the National Immunisation Programme. Euro Surveill. 2012;17(17). 28. 28.Dienst Uitvoering Onderwijs [Education Executive Agency]. Open databestanden [Open datafiles] 2013 [Available from: [https://duo.nl/open\_onderwijsdata/databestanden/](https://duo.nl/open_onderwijsdata/databestanden/). 29. 29.Team RC. R: A language and environment for statistical computing. Vienna, Austria.: R foundation for statistical computing; 2016. 30. 30.van den Hof S, Conyn-van Spaendonck MA, van Steenbergen JE. Measles epidemic in the Netherlands, 1999-2000. J Infect Dis. 2002;186(10):1483–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/344894&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12404165&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000178891000014&link_type=ISI) 31. 31.van Velzen E, de Coster E, van Binnendijk R, Hahne S. Measles outbreak in an anthroposophic community in The Hague, The Netherlands, June-July 2008. Euro Surveill. 2008;13(31). 32. 32.Mollema L, Wijers N, Hahne SJ, van der Klis FR, Boshuizen HC, de Melker HE. Participation in and attitude towards the national immunization program in the Netherlands: data from population-based questionnaires. BMC Public Health. 2012;12:57. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2458-12-57&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22264347&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F14%2F2022.05.12.22273875.atom) 33. 33.Kauffmann F, Heffernan C, Meurice F, Ota MOC, Vetter V, Casabona G. Measles, mumps, rubella prevention: how can we do better? Expert Rev Vaccines. 2021;20(7):811–26.