PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of bmcphBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Public Health
 
BMC Public Health. 2017; 17: 957.
Published online 2017 December 15. doi:  10.1186/s12889-017-4961-9
PMCID: PMC5732449

Spatial clustering of measles vaccination coverage among children in sub-Saharan Africa

Abstract

Background

During the past two decades, vaccination programs have greatly reduced global morbidity and mortality due to measles, but recently this progress has stalled. Even in countries that report high vaccination coverage rates, transmission has continued, particularly in spatially clustered subpopulations with low vaccination coverage.

Methods

We examined the spatial heterogeneity of measles vaccination coverage among children aged 12–23 months in ten Sub-Saharan African countries. We used the Anselin Local Moran’s I to estimate clustering of vaccination coverage based on data from Demographic and Health Surveys conducted between 2008 and 2013. We also examined the role of sociodemographic factors to explain clustering of low vaccination.

Results

We detected 477 spatial clusters with low vaccination coverage, many of which were located in countries with relatively high nationwide vaccination coverage rates such as Zambia and Malawi. We also found clusters in border areas with transient populations. Clustering of low vaccination coverage was related to low health education and limited access to healthcare.

Conclusions

Systematically monitoring clustered populations with low vaccination coverage can inform supplemental immunization activities and strengthen elimination programs. Metrics of spatial heterogeneity should be used routinely to determine the success of immunization programs and the risk of disease persistence.

Keywords: Africa, MCV, Measles, Spatial analysis, Spatial heterogeneity, Spatial regression

Background

Measles is a highly contagious viral disease and is one of the leading causes of death among children in low-income countries, accounting for 114,900 deaths globally in 2014 of which 73,914 (63%) occurred in Africa [1, 2]. Measles also continues to cause epidemics in high-income countries, despite the availability of a safe and highly efficacious vaccine [3, 4].

The Measles-Rubella Initiative, spearheaded by the American Red Cross, the US Centers for Disease Control and Prevention (CDC), the World Health Organization (WHO), and others, has targeted the measles virus for global elimination. This initiative aims to reduce annual measles incidence rates (IRs) to less than five cases per million, requiring >90% coverage of at least one dose of measles containing vaccine (MCV), recommended at age 12 months, by the end of 2015 and >95% coverage by 2020 in all countries [5]. Improvement in vaccination coverage has decreased measles deaths from over half a million globally in 2000 to 114,900 in 2014 [1]. Since 2010, however, progress has stalled [1]: The 2015 vaccination goal was not met and IRs remained relatively unchanged between 2013 and 2014 [1].

Measles elimination is complicated by the high transmission rate of the measles virus. This transmission rate can be expressed as the basic reproductive rate (R0), defined as the number of infections caused, on average, by one infectious person in a fully susceptible population [6]. The R0 for measles ranges from 15 to 20 infections, which is one of the highest among all infectious diseases (e.g., influenza has an R0 around 1.5–2.0) [7]. This high R0 leads to the very high critical vaccination fraction for measles of 95%, i.e., the vaccination coverage needed for herd immunity [8]. This critical vaccination fraction assumes that vaccination coverage and population mixing are distributed homogeneously throughout a country [9]. This assumption of homogeneity is not always realistic, as recently found in Mozambique [10] and Malawi [11]. Spatial heterogeneity of vaccination coverage can increase the critical vaccination fraction required for herd immunity to a level exceeding the 95% coverage goal set by the Measles-Rubella Initiative [12, 13].

Spatial heterogeneity of vaccination coverage can delay disease elimination, as illustrated by continued measles outbreaks, even in countries with high average nationwide vaccination coverage rates [14, 15]. Substantial heterogeneity in measles vaccination coverage has been demonstrated previously in Sub-Saharan Africa [16], but drivers for this heterogeneity are poorly understood. Timely detection and targeting of low-coverage population clusters by supplemental immunization activities (SIAs) can lead to protective herd immunity and accelerate disease elimination, as demonstrated by the successful strategy used for measles elimination in the Americas [17]. We used publicly available microdata from the Demographic and Health Surveys (DHS) in Sub-Saharan Africa to determine if subpopulations with low vaccination coverage are clustered and explore possible determinants.

Methods

Clustering algorithm

We collected measles vaccination coverage data from the most recent DHS conducted in ten countries: Burundi in 2010; the Democratic Republic of the Congo (DRC) in 2013; Kenya in 2008; Madagascar in 2008; Malawi in 2010; Mozambique in 2011; Rwanda in 2010; Tanzania in 2010; Zambia in 2013–14; and Zimbabwe in 2010–11 [18]. These countries were selected based on their contiguity and data availability. We obtained approval from DHS to download and use these data for this study.

DHS are nationwide surveys that are representative of the population and detailed survey methodology has been published elsewhere [19]. In short, DHS are performed using a two-stage cluster sampling design: In the first stage, the DHS selects a random sample of clusters (groups of possible sample households in close proximity to each other) from an already existing sample frame (e.g., a population census); in the second stage, a random sample of households is selected within each cluster. The DHS also determines sample weights that should be applied to survey data to ensure that all subpopulations are equally represented [20]. We extracted, from DHS, the vaccination status (first dose MCV) of children aged 12–23 months measured from sampled households in each cluster, not differentiating between vaccine doses received from routine immunization or SIAs. In the DHS, vaccination status is obtained from vaccination cards where possible and otherwise from mothers’ reports [20]. We calculated the cluster-level vaccination coverage rate as the median of the weighted household-level vaccination coverage rates.

We estimated the spatial association of MCV coverage rates among DHS clusters using the Global Moran’s I and Anselin Local Moran’s I statistics. The Global Moran’s I ranges from −1 to 1 and is a single estimate of spatial association among all DHS clusters (spatial autocorrelation). Values close to zero indicate the absence of a spatial association (i.e., a random distribution), values close to negative one indicate strong spatial dispersion, and values close to positive one indicate strong clustering (autocorrelation). The Anselin Local Moran’s I estimates the association of vaccination coverage rates between a DHS cluster and its neighboring clusters within a specified geographical area (inter-cluster variation). The Anselin Local Moran’s I has been used previously for similar analyses to locate pockets of childhood stunting in Nigeria [21]. We used the Anselin Local Moran’s I to estimate spatial clustering of low (< 75%), high (≥ 75%), or mixed (low near high or vice versa) weighted vaccination coverage. We considered Moran’s I statistics with p-values < 0.05 to be statistically significant. DHS cluster data have been used previously for similar clustering analyses, such as examining malnutrition in Ethiopia [22], of HIV prevalence in Burundi [23], and of childhood stunting in Nigeria [24].

Determinants of low-vaccination clusters

We explored possible determinants for clustering of low-vaccination using additional information from country DHS: (1) child in possession of a health card or not (H c); (2) mother had heard of oral rehydration salts (ORS) or not (O); (3) mother is literate or not (T); (4) mother visited a health facility in the last 12 months or not (H f); (5) mother mentioned that money had been a barrier to seeking healthcare in the past or not (M). We calculated the cluster-level percent children with a health card as the median of the weighted household-level percentages of children with a health card. From the other household-level variables, we computed the cluster-level equivalents as the percent of mothers (households) that answered each question affirmatively.

We used a logistic regression model to estimate the association between the odds for a DHS cluster being part of a low-vaccination spatial cluster and the aforementioned explanatory factors. We adjusted for spatial autocorrelation (inter-cluster variation) of vaccination status among clusters with a queen contiguity weights matrix based on spatial lags. Queen contiguity calculates spatial autocorrelation of the outcome variable among all contiguous neighbors, after creating Thiessen polygons around cluster coordinates [2527]. We selected a contiguity weights matrix instead of a distance matrix due to the large variability in distances between clusters in our ten-country study area.

Our model took the following form:

lnLx1Lx=ρWlnLx1Lx+β0+β1Hc+β2O+β3T+β4Hf+β5M+ε

where Lx1Lx represented the odds of being in a low-vaccination spatial cluster, ρ represented the spatial autoregressive coefficient for the log odds of being in a low-vaccination spatial cluster, W represented the queen contiguity matrix, β 1–5 represented the regression coefficients for aforementioned covariates, and ε represented the error term.

We used SAS version 9.4 and ArcGIS version 10.4 for this analysis.

Results

Country-level vaccination coverage

We included a total of 5458 DHS clusters containing 70,092 households across all ten countries (Fig. 1, Table 1). This sample is representative of a total population of 214,339,000 people. Nationwide MCV coverage among children aged 12–23 months was below the measles critical vaccination fraction of 95% for nine out of the ten countries and ranged from a low of 69.6% for Madagascar to a high of 95% for Rwanda. The average MCV coverage for all 10 countries, weighted by population size, was 83.6%.

Fig. 1
Vaccination coverage and DHS clusters in the study area. The location of each DHS cluster is depicted as a grey circle. We computed the average vaccination coverage rate for each country from DHS cluster-level data. Both maps were created by study investigators ...
Table 1
DHS survey populations included by country

Clustering of low vaccination coverage rates

We found strong spatial heterogeneity in measles vaccination coverage across the entire ten-country region (Global Moran’s I of 0.388, p < 0.001). Based on the Anselin Local Moran’s I, we identified statistically significant spatial correlation of low vaccination coverage (< 75%) between 477 DHS clusters, of mixed coverage between 148 clusters, and of high coverage (≥ 75%) between 645 clusters (Fig. 2a). The DRC had the second-lowest nationwide vaccination coverage rate in our study region and had clustering of low coverage throughout the country. We found clustering of high coverage almost uniformly throughout Rwanda and Burundi, two countries with the highest nationwide average vaccination coverage in our sample. In other countries, clustering of low-coverage was concentrated in specific geographic areas: e.g., East Kenya, North Malawi, North Zambia, South Zimbabwe, and South Mozambique. Madagascar had the lowest average nationwide MCV coverage in our sample and had clustering of low-coverage throughout the country except in the capital region. We also found clustering of low coverage across the Kenya-Tanzania and the Malawi-Zambia borders.

Fig. 2
Spatial clustering of vaccination coverage in DHS clusters. Using the Anselin Local Moran’s I, we classified each DHS cluster as being part of a spatial cluster with low-vaccination, high-vaccination, or mixed vaccination coverage (low-vaccination ...

We found three areas with clustering of vaccination coverage that contrasted nationwide average rates: (1) the Zanzibar/Pemba island population in Kenya-Tanzania had clustering of low coverage while nationwide rates were relatively high (Fig. (Fig.2b);2b); (2) the Madagascar capital region had clustering of high coverage while the nationwide rates were low (Fig. (Fig.2c);2c); and (3) the Northern Malawi region had clustering of low vaccine coverage compared to high nationwide coverage (Fig. (Fig.2d).2d). Each of these areas has distinctive geological features that separated them from the surrounding area: Zanzibar and Pemba are Tanzanian island-populations with semi-autonomous governments; the Madagascar capital region of Antananarivo is located in the mountainous Hauts Plateaux region, separated from the rest of the country; and Northern Malawi includes most of Lake Malawi and areas of higher elevation compared to the South of the country that includes the river Shire.

Determinants of low-vaccination

We explored possible determinants for clustering of low vaccination coverage using a spatial regression model. Clustering of low vaccination coverage was associated with children not having a health card and mothers not having knowledge of ORS. Clustering of low vaccination was 4.6% less likely for each percentage point increase in children with a health card (95% CI: -0.066, −0.026, p < 0.01) and 1.7% less likely for each percentage point increase in mothers with knowledge of ORS (95% CI: -0.033, −0.001, p < 0.05). In addition, clustering of low coverage was inversely related to having financial restrictions to healthcare, where we found 1.6% less likely per percentage point increase in mothers listing financial barriers to seeking healthcare (95% CI: -0.029, −0.003, p < 0.05). Maternal literacy rates and a maternal history of visiting a health clinic were not statistically significantly associated with clustering of low vaccination coverage.

Discussion

Using publicly available DHS data from 10 countries, we found 477 geographical clusters of low measles vaccination coverage spread across Sub-Saharan Africa, many of which contrasted relatively high nationwide average vaccination coverage rates. These clusters can weaken herd immunity, cause inequity in disease risk, and delay elimination programs. Indeed, recent measles outbreaks have occurred in subpopulations with low immunization rates: Zambia had an average MCV coverage of 84.9%, and Malawi of 93%, but both countries experienced a large measles outbreak in 2010–2011 [28]; this outbreak spread from high-risk subpopulations in South Africa to Zambia, Malawi, and to high-risk subpopulations in Tanzania consistent with the clusters that we identified [29, 30]. The persistence of virus transmission due to highly connected, clustered, unimmunized subpopulations has also been demonstrated by mathematical metapopulation models [3133]. These models can be used to compute vaccination coverage goals that take into account spatial clustering of low vaccination.

We found that clustering of low vaccination coverage was more likely in populations with low health education and with limited access to healthcare. Previous studies have found similar risk factors for low immunization rates [11, 34, 35]. We also found that financial barriers to healthcare were associated with better vaccination rates, which seems counterintuitive. One possibility for this observed relationship may be that vaccination is often free of charge and may not be affected by financial barriers. In Malawi, for example, high vaccine uptake was observed despite significant cost and travel time, possibly related to high levels of trust in the effectiveness of the vaccine to prevent serious disease [36]. In high-income countries, indirect measures of wealth have been found to correlate with decreased vaccination coverage due to vaccination hesitancy [37], but income-related vaccination hesitancy has not been found (yet) in Sub-Saharan Africa.

Subpopulations with low vaccination coverage across country borders are a particular concern for measles elimination because these transient populations are often not covered by national immunization programs [10, 16, 17]. We found such subpopulations at the Kenya-Tanzania border and the Malawi-Zambia border. The Kenya-Tanzania border area includes the famous Serengeti and Kilimanjaro national parks and is inhabited by the nomadic Maasai people, who have among the lowest vaccination coverage in Tanzania due to lower use of healthcare services [38, 39]. The Malawi-Zambia border is crossed frequently by the Chewa people that reside in both countries, though this group has been found to be no less vaccination-compliant than other groups of similar socioeconomic status in the region [40, 41]. Trans-border populations with low vaccination coverage can be especially vulnerable to disease importations from one country into another. Such importations occurred during the 2010–2011 measles outbreak that spread from Malawi into Zambia [28]. Coordination of immunization activities between countries will be essential to increase coverage and eliminate measles in these cross-border populations [17].

Most countries in our sample had vaccination coverage rates well below the measles critical vaccination fraction and have already been identified by the Measles-Rubella Initiative as high priority areas for continued activities to increase immunization rates [42]. We found strong spatial heterogeneity of vaccination coverage in some of these countries, indicating that the nationwide vaccination coverage target of 95% set by the Measles-Rubella Initiative may not lead to herd immunity, but that targeted SIAs will be necessary to reach particularly vulnerable populations. For this reason, the Measles-Rubella Initiative also monitors vaccination coverage at the district level.

Conclusions

Even in countries where national average vaccination coverage rates approach the critical vaccination fraction for measles, subnational vaccination coverage can fall short, leaving subpopulations, including those that cross country borders, vulnerable to outbreaks. Systematically identifying and monitoring these low vaccination sub-populations can inform SIAs and strengthen elimination programs. In addition to monitoring average vaccination coverage statistics, metrics of spatial heterogeneity should be used to determine the success of immunization programs and the risk of disease persistence.

Acknowledgements

Not applicable.

Funding

This project was supported by the Bill and Melinda Gates Foundation (Grant 49,276, Evaluation of Candidate Vaccine Technologies Using Computational Models) and the US National Institute of General Medical Sciences (Grant 5U54GM088491, Computational Models of Infectious Disease Threats).

Availability of data and materials

The datasets generated and/or analyzed during the current study are available upon request from the DHS website, http://dhsprogram.com/data/Access-Instructions.cfm.

Approval must be obtained from DHS before use of these data.

Abbreviations

CDC
Centers for Disease Control and Prevention
DHS
Demographic and Health Surveys
DRC
Democratic Republic of the Congo
IR
Incidence rate
MCV
Measles containing vaccine
ORS
Oral rehydration salts
SIA
Supplemental immunization activity
WHO
World Health Organization

Authors’ contributions

Authors’ contributions

TKB preformed the statistical analysis, figures, and wrote the first version of the manuscript. ZMD contributed to the analysis and edited the manuscript. WGVP contributed to study design, edited the manuscript and prepared the final manuscript version. All authors read and approved the final manuscript.

Notes

Ethics approval and consent to participate

We have received written approval to use Demographic and Health Survey (DHS) data for this study from the DHS administration, via their online application system.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Tenley K. Brownwright, ude.ttip@44bet.

Zan M. Dodson, ude.ttip@nosdod.naz.

Willem G. van Panhuis, Phone: 412-624-7693, ude.ttip@01vaw.

References

1. Perry RT, Murray JS, Gacic-Dobo M, Dabbagh A, Mulders MN, Strebel PM, et al. Progress toward regional measles elimination - worldwide, 2000-2014. MMWR. Morb. Mortal. Wkly. Rep. United States; 2015;64:1246–51. [PubMed]
2. World Health Organization. Measles [Internet]. 2016 [cited 2016 Mar 1]. Available from: www.who.int/mediacentre/factsheets/fs286/en/.
3. Sudfeld CR, Navar AM, Halsey NA. Effectiveness of measles vaccination and vitamin A treatment. Int. J. Epidemiol. England; 2010;39 Suppl 1:i48–i55. [PMC free article] [PubMed]
4. Demicheli V, Rivetti A, Debalini M, Di Pietrantonj C. Vaccines for measles, mumps and rubella in children (Review). cochrane Libr. [Internet]. 2012 [cited 2016 Aug 17];2:CD004407. Available from: http://www.ncbi.nlm.nih.gov/pubmed/22336803. [PubMed]
5. World Health Organization. Global measles and rubella strategic plan 2012-2020. Geneva: World Health Organization; 2012.
6. Blumberg S, Lloyd-Smith JO. Comparing methods for estimating R0 from the size distribution of subcritical transmission chains. Epidemics Netherlands; 2013;5:131–145. [PMC free article] [PubMed]
7. Griffin DE, Pan C-H, Moss WJ. Measles vaccines. Front Biosci United States; 2008;13:1352–1370. [PubMed]
8. Scherer A, McLean A. Mathematical models of vaccination. Br Med Bull England; 2002;62:187–199. [PubMed]
9. Feng Z, Hill AN, Smith PJ, Glasser JW. An elaboration of theory about preventing outbreaks in homogeneous populations to include heterogeneity or preferential mixing. J Theor Biol. 2015;386:177–187. doi: 10.1016/j.jtbi.2015.09.006. [PMC free article] [PubMed] [Cross Ref]
10. Muloliwa AM, Camacho LAB, Verani JFS, Simões TC, Dgedge M do C. Impact of vaccination on the incidence of measles in Mozambique in the period 2000 to 2011. Cad. saúde pública.2013;29:257–69. [PubMed]
11. Abebe DS, Nielsen VO, Finnvold JE, You D, Wardlaw T, Salama P, et al. Regional inequality and vaccine uptake: a multilevel analysis of the 2007 Welfare Monitoring Survey in Malawi. BMC Public Health [Internet]. BioMed Central; 2012 [cited 2016 Sep 16];12:1075. Available from: http://bmcpublichealth.biomedcentral.com/articles/10.1186/1471-2458-12-1075. [PMC free article] [PubMed]
12. Liu F, Enanoria WTA, Zipprich J, Blumberg S, Harriman K, Ackley SF, et al. The role of vaccination coverage, individual behaviors, and the public health response in the control of measles epidemics: an agent-based simulation for California. BMC Public Health. England; 2015;15:447. [PMC free article] [PubMed]
13. Wallinga J, Heijne JCM, Kretzschmar M. A measles epidemic threshold in a highly vaccinated population. PLoS Med. United States; 2005;2:e316. [PMC free article] [PubMed]
14. Suijkerbuijk AWM, Woudenberg T, Hahne SJM, Nic Lochlainn L, de Melker HE, Ruijs WLM, et al. Economic costs of measles outbreak in the Netherlands, 2013-2014. Emerg Infect Dis United States; 2015;21:2067–2069. [PMC free article] [PubMed]
15. Zipprich J, Winter K, Hacker J, Xia D, Watt J, Harriman K. Measles outbreak--California, December 2014-February 2015. MMWR. Morb. Mortal. Wkly. Rep. United States; 2015;64:153–154. [PMC free article] [PubMed]
16. Takahashi S, Metcalf CJE, Ferrari MJ, Tatem AJ, Lessler J. The geography of measles vaccination in the African Great Lakes region. Nat. Commun England; 2017;8:15585. [PMC free article] [PubMed]
17. Castillo-Solorzano C, Marsigli C, Danovaro-Holliday MC, Ruiz-Matus C, Tambini G, Andrus JK. Measles and rubella elimination initiatives in the Americas: lessons learned and best practices. J Infect Dis Oxford University Press; 2011;204:S279–83. [PubMed]
18. ICF International. Children’s recode [multiple] [Internet]. Rockville, MD: Measure DHS; 2008-2014. Available from: dhsprogram.com
19. DHS. Methodology [Internet]. Meas. DHS. 2016 [cited 2016 May 12]. Available from: http://dhsprogram.com/What-We-Do/Survey-Types/DHS-Methodology.cfm.
20. Rutstein SO, Rojas G. Guide to DHS statistics: demographic and health surveys methodology [internet]. Calverton MD, USA; 2006. Available from: http://dhsprogram.com/pubs/pdf/DHSG1/Guide_to_DHS_Statistics_29Oct2012_DHSG1.pdf.
21. Adekanmbi VT, Uthman OA, Mudasiru OM, Adekanmbi V, Kayode G, Uthman O, et al. Exploring variations in childhood stunting in Nigeria using league table, control chart and spatial analysis. BMC Public Health [Internet]. BioMed Central; 2013 [cited 2016 Aug 17];13:361. Available from: http://bmcpublichealth.biomedcentral.com/articles/10.1186/1471-2458-13-361. [PMC free article] [PubMed]
22. Haile D, Azage M, Mola T, Rainey R. Exploring spatial variations and factors associated with childhood stunting in Ethiopia: spatial and multilevel analysis. BMC Pediatr. England; 2016;16:49. [PMC free article] [PubMed]
23. Barankanira E, Molinari N, Niyongabo T, Laurent C. Spatial analysis of HIV infection and associated individual characteristics in Burundi: indications for effective prevention. BMC Public Health. England; 2016;16:118. [PMC free article] [PubMed]
24. Adekanmbi VT, Uthman OA, Mudasiru OM, Adekanmbi V, Kayode G, Uthman O, et al. Exploring variations in childhood stunting in Nigeria using league table, control chart and spatial analysis. BMC Public Health. BioMed Central; 2013;13:361. [PMC free article] [PubMed]
25. Getis A, Aldstadt J. Constructing the Spatial Weights Matrix Using a Local Statistic. Geogr. Anal. [Internet]. Blackwell Publishing Ltd; 2004 [cited 2016 Aug 17];36:90–104. Available from: http://doi.wiley.com/10.1111/j.1538-4632.2004.tb01127.x.
26. Anselin L. Spatial econometrics: a companion to theoretical econometrics. Hoboken NJ: Blackwell Publishing Ltd; 2001.
27. Anselin L, Syabri I, Kho Y. GeoDa: An Introduction to Spatial Data Analysis. Geogr Anal. [Internet]. 2006;38:5–22. [cited 2017 Mar 15] Available from: http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.74.9014.
28. Pinchoff J, Chipeta J, Banda GC, Miti S, Shields T, Curriero F, et al. Spatial clustering of measles cases during endemic (1998–2002) and epidemic (2010) periods in Lusaka, Zambia. BMC Infect. Dis. England; 2015;15:121. [PMC free article] [PubMed]
29. World Health Organization. Global health observatory data respository [Internet]. Geneva; 2016. Available from: http://apps.who.int/gho/data/?theme=main.
30. Sartorius B, Cohen C, Chirwa T, Ntshoe G, Puren A, Hofman K. Identifying high-risk areas for sporadic measles outbreaks: lessons from South Africa. Bull World Health Organ Switzerland; 2013;91:174–183. [PubMed]
31. Metcalf CJE, Hampson K, Tatem AJ, Grenfell BT, Bjornstad ON. Persistence in epidemic metapopulations: quantifying the rescue effects for measles, mumps, rubella and whooping cough. PLoS One. United States; 2013;8:e74696. [PMC free article] [PubMed]
32. Gunning CE, Wearing HJ. Probabilistic measures of persistence and extinction in measles (meta)populations. Ecol Lett England; 2013;16:985–994. [PMC free article] [PubMed]
33. Bharti N, Djibo A, Ferrari MJ, Grais RF, Tatem AJ, McCabe CA, et al. Measles hotspots and epidemiological connectivity. Epidemiol Infect England; 2010;138:1308–1316. [PMC free article] [PubMed]
34. Russo G, Miglietta A, Pezzotti P, Biguioh RM, Bouting Mayaka G, Sobze MS, et al. Vaccine coverage and determinants of incomplete vaccination in children aged 12–23 months in Dschang, West Region, Cameroon: a cross-sectional survey during a polio outbreak. BMC Public Health. England; 2015;15:630. [PMC free article] [PubMed]
35. Okwaraji YB, Mulholland K, Schellenberg JRMA, Andarge G, Admassu M, Edmond KM. The association between travel time to health facilities and childhood vaccine coverage in rural Ethiopia. A community based cross sectional study. BMC Public Health. England; 2012;12:476. [PMC free article] [PubMed]
36. Holte JH, Maestad O, Jani J V. The decision to vaccinate a child: an economic perspective from southern Malawi. Soc. Sci. Med. England; 2012;75:384–391. [PubMed]
37. McNutt L-A, Desemone C, DeNicola E, El Chebib H, Nadeau JA, Bednarczyk RA, et al. Affluence as a predictor of vaccine refusal and underimmunization in California private kindergartens. Vaccine Netherlands; 2016;34:1733–1738. [PubMed]
38. Maasai Association. The Maasai People [Internet]. [cited 2016 Sep 15]. Available from: http://www.maasai-association.org/maasai.html.
39. Lawson DW, Borgerhoff Mulder M, Ghiselli ME, Ngadaya E, Ngowi B, Mfinanga SGM, et al. Ethnicity and child health in northern Tanzania: Maasai pastoralists are disadvantaged compared to neighbouring ethnic groups. PLoS One. United States; 2014;9:e110447. [PMC free article] [PubMed]
40. The Peoples of the World Foundation. Education for and about Indigenous Peoples [Internet]. 2016 [cited 2016 Sep 15]. Available from: http://www.peoplesoftheworld.org/hosted/chewa/.
41. Bowie C, Mathanga D, Misiri H. Poverty, access and immunization in Malawi - a descriptive study. Malawi Med J. 2006;18:19–27.
42. American Red Cross. Measles and rubella initiative [Internet]. 2016 [cited 2016 Feb 12]. Available from: measlesandrubellainitiative.org

Articles from BMC Public Health are provided here courtesy of BioMed Central