Figures
Abstract
Background
Despite recent success towards controlling poliovirus transmission, Nigeria has struggled to achieve uniformly high routine vaccination coverage. A lack of reliable vaccination coverage data at the operational level makes it challenging to target program improvement. To reliably estimate vaccination coverage, we conducted district-level vaccine coverage surveys using a pre-existing infrastructure of polio technical staff in northern Nigeria.
Methods
Household-level cluster surveys were conducted in 40 polio high risk districts of Nigeria during 2014–2015. Global positioning system technology and intensive supervision by a pool of qualified technical staff were used to ensure high survey quality. Vaccination status of children aged 12–23 months was documented based on vaccination card or caretaker’s recall. District-level coverage estimates were calculated using survey methods.
Results
Data from 7,815 children across 40 districts were analyzed. District-level coverage with the third dose of diphtheria-pertussis-tetanus vaccine (DPT3) ranged widely from 1–63%, with all districts having DPT3 coverage below the target of 80%. Median coverage across all districts for each of eight vaccine doses (1 Bacille Calmette-Guérin dose, 3 DPT doses, 3 oral poliovirus vaccine doses, and 1 measles vaccine dose) was <50%. DPT3 coverage by survey was substantially lower (range: 28%–139%) than the 2013 administrative coverage reported among children aged <12 months. Common reported reasons for non-vaccination included lack of knowledge about vaccines and vaccination services (50%) and factors related to access to routine immunization services (15%).
Conclusions
Survey results highlighted vaccine coverage gaps that were systematically underestimated by administrative reporting across 40 polio high risk districts in northern Nigeria. Given the limitations of administrative coverage data, our approach to conducting quality district-level coverage surveys and providing data to assess and remediate issues contributing to poor vaccination coverage could serve as an example in countries with sub-optimal vaccination coverage, similar to Nigeria.
Citation: Gunnala R, Ogbuanu IU, Adegoke OJ, Scobie HM, Uba BV, Wannemuehler KA, et al. (2016) Routine Vaccination Coverage in Northern Nigeria: Results from 40 District-Level Cluster Surveys, 2014-2015. PLoS ONE 11(12): e0167835. https://doi.org/10.1371/journal.pone.0167835
Editor: Ray Borrow, Public Health England, UNITED KINGDOM
Received: April 22, 2016; Accepted: November 21, 2016; Published: December 9, 2016
This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by the US Centers for Disease Control and Prevention. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Vaccine preventable diseases cause significant mortality among children aged under 5 years [1]. In 2012, the World Health Assembly endorsed the Global Vaccine Action Plan, which includes an aim of ensuring delivery of universal access to immunization with associated targets of reaching 90% national vaccination coverage and at least 80% vaccination coverage in every district [2]. To support this goal and to respond to the 2012 World Health Assembly’s declaration of polio eradication as a public health emergency, Nigeria is working to achieve polio-free status [3], and continues to strengthen routine immunization (RI) service delivery to achieve high and equitable coverage [4].
In Nigeria, the Expanded Program on Immunization (EPI) was launched in 1979 and provides routine vaccinations to all children aged <12 months and pregnant women [5]. Core vaccines provided to infants as part of EPI include Bacille Calmette-Guérin vaccine (BCG), diphtheria-pertussis-tetanus vaccine (DPT), oral poliovirus vaccine (OPV), and measles vaccine. Though centrally managed by the National Primary Health Care Development Agency (NPHCDA) of the Federal Ministry of Health (FMOH), RI services are organized and implemented at subnational levels, including states and Local Government Areas (LGAs), which are equivalent to districts. Given the variation in sociodemographic characteristics of the population over a large number of states (36 plus the Federal Capital Territory, or FCT) and LGAs (774), RI service delivery must be tailored to suit local needs and assure high, uniform vaccination coverage.
During 2010–2014, national coverage with the third dose of DPT vaccine (DPT3), a key indicator of RI program performance, i.e., the percentage of children aged <12 months reported to have received the vaccine dose, ranged from 57%–74% by administrative reporting [6]. Multiple coverage surveys suggest even lower vaccination coverage than reported by administrative data and demonstrate that national estimates can often mask subnational immunity gaps at state and LGA levels [7]. Data from the 2013 National Demographic and Health Survey (DHS) show that while the national estimate for DPT3 coverage was 38%, state-level coverage ranged from 3%–76% across 20 states in the northern part of the country [8].
A lack of reliable vaccination coverage data at the operational-level (i.e., LGA), makes it challenging to monitor and remediate RI service delivery to achieve uniformly high vaccination coverage. To fill this information gap, we used a pre-existing infrastructure of polio technical field staff to conduct RI coverage surveys in 40 polio high-risk LGAs across eight states in northern Nigeria during 2014–2015. The overall survey objectives of providing LGA-level vaccine coverage estimates and RI service delivery information to LGA immunization staff for use in planning and program improvement were successfully achieved.
Materials and Methods
Survey Setting
The National Stop Transmission of Polio (NSTOP) program was collaboratively established in 2012 by the NPHCDA, the Nigeria Field Epidemiology and Laboratory Training Program (NFELTP) and U.S. Centers for Disease Control and Prevention (CDC), to place staff at national, state and LGA levels to strengthen RI service delivery with the goal of aiding polio eradication efforts. Senior NSTOP staff placed at the central and state level are graduates of the 2-year NFELTP program, including classroom training on polio and RI services and supervisory experience during OPV supplemental immunization activities (SIAs) [9]. At the LGA level, 100 NSTOP staff were hired in 100 LGAs deemed to be at high risk for polio based on the CDC-Global Good risk analysis algorithm [10]. The LGA-level coverage surveys provided an opportunity to strengthen capacity of staff and obtain data about RI service delivery performance. As neither NSTOP nor NFELTP staff are involved with service delivery at the health facility or community levels, they provided a pool of independent, qualified field staff and a supervisory structure that was well-suited for implementing high quality coverage surveys within a short timeframe.
Survey Population
The target population was children aged 12–23 months living in 100 polio high-risk LGAs of northern Nigeria where NSTOP staff had been assigned. Forty high-risk LGAs were selected for the survey using the following algorithm. We excluded the insecure states of Borno and Yobe, and also Kano state, where ongoing LGA-level coverage surveys were being conducted (Dale Rhoda, KANRICS Task Force, and Biostat Global Consulting, personal communication). In the remaining states, all LGAs were included except in states with more than six total polio high-risk LGAs, where three LGAs with the highest administrative DPT3 coverage from 2013 and three LGAs with the lowest DPT3 coverage were chosen. The exception was Kaduna state, where two additional LGAs were included for the survey pilot phase. Thus, included LGAs were from Bauchi (which had four polio high risk LGAs out of the 20 total LGAs in the state), Jigawa (six out of 27), Kaduna (eight out of 23), Katsina (six out of 34), Kebbi (three out of 21), Sokoto (six out of 23), and Zamfara (six out of 14) states, and the FCT (one).
Survey Design
We conducted a two-stage household-level cluster survey in each LGA using the 2005 World Health Organization (WHO) survey methodology [11]. The sampling frame of primary sampling units (clusters) was a list of enumeration areas (EA) obtained from the National Population Commission (NPC). In the first stage, 30 clusters in each LGA were systematically sampled probability proportional to estimated size, using 2006 census data. PDF maps for selected clusters were obtained from NPC. In the second stage, a standardized protocol was used to select households, defined as groups of persons living and eating together under the same roof. To reduce potential selection bias in the field, pre-assignment of random start points for each cluster was done at the central level. From the start point, teams followed a standardized path through the cluster (moving clockwise and turning right at every opportunity). All housing structures along the path were visited to determine if they contained eligible children 12–23 months of age. Only one eligible child per household was included in the survey; if multiple eligible children lived in the household, one child was selected using a random number table. Teams continued visiting households until they enrolled a total of seven households with eligible children within each cluster. Households that could not be contacted or who declined to participate were replaced. The expected sample size for each LGA survey was 210 children (8,400 in 40 LGAs), based on an expected coverage of 50%, desired precision of +/- 10% with 95% confidence intervals, and a design effect of 2.0 [11]. This survey was approved by the National Health Research Ethics Committee of Nigeria and was determined to be a non-research, public health program evaluation by the Human Research Protection Office of the U.S. CDC, according to its human subjects’ procedures.
Survey Implementation
Data collection tools were tested in March 2014. Implementation of data collection followed in three phases: a pilot phase was conducted in two LGAs in Kaduna (May 2014), followed by Phase 1 (October 2014) and Phase 2 (January 2015) in 19 LGAs each that were delayed due to human resource limitations secondary to ongoing Ebola outbreak response in Nigeria [12]. A central command center for the survey was established in Abuja which included a NSTOP central coordinator, NSTOP data managers, and a GPS specialist. Senior supervisors were central and state-level NSTOP staff that supervised multiple LGAs. Three survey teams were allocated per LGA, with each team consisting of one team supervisor (mostly NFELTP residents and several staff from NPC, FMOH, and NPHCDA) and two interviewers (who were recruited locally and required to have a university degree, be fluent in English and Hausa, and preferably have experience with SIAs). An excess number of interviewers was recruited within the respective LGA and trained on the survey, with final selection for field work dependent on their performance during the training and post-training evaluation.
During the pilot phase, senior and team supervisors were trained centrally by CDC staff and central NSTOP staff over 3 days in Abuja, followed by training of interviewers over 2 days in Kaduna. Lessons learned from the pilot phase were used to adapt survey implementation for subsequent phases: the variable skill-level of LGA level staff required a robust supervisory structure for high-quality survey implementation; long distances and difficult terrain between clusters led to reduction in the number of clusters assigned to teams per day; and GPS navigation was so useful in helping survey teams locate assigned clusters that a full time GPS specialist was hired to be part of the control center team in Abuja for phases 1 and 2. In October 2014, prior to phases 1 and 2, a 3-day refresher training of trainers was conducted in Abuja. This was followed by a 2-day cascade training for interviewers in each LGA conducted by master trainers who had participated in the previous phase.
Each survey team completed 10 clusters within 5 days (an average of two clusters per team per day), for a total of 30 clusters completed by three teams per LGA. Clusters found to lack settlements, contain fewer than seven eligible children, or be inaccessible due to security risk were not replaced or combined with other clusters. Households were visited up to two times in attempts to enroll eligible children, and reasons for non-response were documented. After obtaining verbal consent, the pre-tested, standardized paper questionnaire (S1 Appendix) was administered to the caretaker (or a knowledgeable adult, if the mother was unavailable) of the child in the selected household. Questionnaires were written in both English and Hausa, and interviews were conducted in the appropriate language for each household. Questions included socio-demographic information, RI vaccination history, awareness of RI service opportunities, use of RI services (any vaccination at health clinic or health outreach program), and reasons for non-vaccination. Vaccination status was first assessed based on caretaker’s recall, and second, based on vaccination card.
A high supervisory ratio (i.e., one team supervisor per one interview team, and one senior supervisor per three interview teams) was employed to ensure survey implementation was of the highest quality. Standardized monitoring forms were used by all supervisors in the field to assess adherence to survey protocol and data quality. To reduce selection bias of households in the field, each team was equipped with an Android™ device loaded with base maps and global positioning system (GPS) coordinates of their pre-assigned cluster start points for offline navigation. GPS coordinates were also taken at the beginning and end of fieldwork in each cluster to verify that the cluster was visited.
Data Analysis
Data from completed paper forms were singly entered into an electronic database in Abuja using Census and Survey Processing System (CSPro) software (Version 5.0, U.S. Census Bureau). Both data cleaning and all data management were conducted using CSPro and SAS (Version 9.3, SAS Institute, Cary, North Carolina). Data cleaning included reviewing of the completeness and validity of data entered, making necessary corrections based on review of paper questionnaires, and excluding ineligible children based on date of birth (birthdates within 4 weeks of the eligible dates were allowed to account for the possibility of inaccurate documentation). For each phase, data entry and cleaning began during the survey period, and was completed immediately after field work ended; data analysis was completed at the end of the pilot phase, and then after Phase 2.
Within each LGA, data were analyzed using SAS to account for the cluster sampling, but the sample was assumed to be self-weighting. Vaccination coverage estimates and 95% confidence intervals (CIs) were calculated; modified Wilson (score) CIs are reported for estimates ≥80% and ≤20% and Wald CIs are reported for all others. Graphics were produced in R v3.2 [13]. Since LGAs were selected purposively within states, LGA data were not pooled to obtain state-level estimates. To provide a succinct summary of results, the LGA-level coverage median and range across all LGAs are presented. All coverage estimates represent a combination of data from vaccine cards and caretakers’ recall. OPV and measles coverage did not include vaccine doses given during SIAs. “Complete vaccination coverage” was defined in accordance with NPHCDA policy as receipt of the following eight vaccine doses: BCG, OPV 1/2/3, DPT 1/2/3, and measles. Reasons for not receiving vaccinations through RI services were grouped into categories of similar response choices. As part of our analysis, we compared DPT3 coverage estimates from our surveys of children aged 12–23 months during October 2014 and January 2015 to administrative DPT3 vaccination coverage for 2013 among children aged <12 months reported through the District Vaccine Data Management Tool (DVD-MT) at WHO.
Results
The pilot phase was conducted in two LGAs in Kaduna state in May 2014. Phase 1 was conducted in October 2014 and included survey implementation in 19 LGAs in four states (Katsina, Jigawa, Sokoto and the FCT); and Phase 2 was conducted in January 2015 in 19 LGAs in four additional states (Bauchi, Kaduna, Kebbi, and Zamfara). Of the 1,200 clusters visited during the surveys in 40 LGAs, 26 clusters (2%) were excluded due to local security issues, and 5 clusters (<1%) were excluded because they lacked household structures (settlement had been relocated post-cartography). In addition, 43 clusters (4%) were included which had fewer than the target of seven households with eligible children (number of eligible children ranged 1–6; median of 5 per cluster). The average number of respondents per cluster in each LGA is described in S2 Appendix. Clusters reported to have fewer than seven eligible children were accepted only after verification of number of household structures by central staff using GPS information and satellite imagery. After excluding children based on age ineligibility (190 children), 7,815 children from 1,169 clusters in 40 LGA surveys were included in the analysis. A total of 573 households had more than one eligible child, which required survey teams to use a random number table to select one child for inclusion in the survey. Revisits were conducted to 74 households found to have no caretaker at home during the initial visit.
As shown in Table 1, the median percentage of male children was 49% (range: 43%–59%). For the highest level of education among mothers/caretakers and heads of households, Quranic schooling had the highest median percentage (57% and 48%, respectively). Residence in urban or rural areas varied widely within each LGA, but the median percentage of households residing in rural areas was 82% compared to 22% for urban areas across all LGAs.
Coverage estimates of routine vaccinations ranged widely across all LGAs (Fig 1 and S3 Appendix). Median coverage for all vaccines, as well as the proportion of children with complete vaccination was below 50%. Median OPV1 coverage was highest at 44% (range: 10%–98%), and median DPT3 coverage was lowest at 14% (range: 1%-–63%). Coverage estimates also ranged widely between LGAs within the same state, as shown for DPT3 (Fig 2 and S3 Appendix).
Estimates of vaccine coverage do not include OPV or measles doses given during supplemental immunization activities. Complete vaccination coverage is defined as receiving eight antigens [Bacille Calmette-Guérin (BCG), Oral Poliovirus (OPV) 1/2/3, Diphtheria-pertussis-tetanus (DPT) 1/2/3, and measles). For boxplot interpretation: the dark line is the median; “hinges” are the top and bottom of the box. The upper and lower "hinges" correspond to the first and third quartiles (25th and 75th percentiles). The upper and lower whiskers represent +/- 1.5 * IQR, where IQR is the inter-quartile range, or distance between the first and third quartiles. Data beyond the end of the whiskers are outliers and plotted as points (as specified by Tukey).
LGAs are grouped by state to illustrate variability in coverage across LGAs within the same state. This data is not representative of state-level coverage since LGAs were purposefully selected. Black vertical lines depict 95% confidence intervals. For comparison, administrative coverage for each LGA (source: World Health Organization, Nigeria; DVD-MT data, 2013) is represented with a black dot.
The reported use of RI services also varied widely; by LGA, median use was 47%, with a range of 13% to 98% (Table 2). Among those who ever used RI services, the median proportion of mothers/caretakers who reported receiving a vaccine card was 87%; among those who reported receiving a card, the median proportion with available cards was 45%. For children not receiving all the recommended vaccinations through RI services, the highest median responses were reasons attributed to lack of knowledge or education (50%), followed by mother forgetting or being too busy (16%), and suboptimal access to RI services (15%) (Table 2).
The highest median response by LGA for main information source about RI services was loud speaker/town announcer (43%), followed by health worker (21%) and radio (15%) (Table 3). We also asked about the main type of information needed to help mothers/caretakers decide to vaccinate children through RI services (Table 3). Knowing the safety of the vaccine and the importance of vaccination had the highest median responses (37% and 23%, respectively), but also the widest ranges. Knowing which diseases the vaccine protected against had a 10% median response.
When comparing administrative DPT3 coverage data reported for 2013 to the DPT3 coverage estimates obtained during our survey by respective LGA, we found that the administrative estimates were higher than the survey estimates for all 40 LGAs (Fig 2 and Table 4). The difference in survey vs. administrative coverage ranged from 9% in Abuja Metropolitan Area Council (AMAC) LGA, Abuja FCT to 139% in Birnin Kudu LGA, Jigawa state. Reported administrative coverage fell within the 95% CI of our survey estimate in only 2 LGAs; AMAC LGA, Abuja FCT and Itas/Gadau, Bauchi state.
Discussion
Our survey findings highlight gaps in vaccination coverage at the LGA level in areas of northern Nigeria that are high risk for polio. In the LGAs surveyed, median DPT3 coverage was 14%, with all 40 LGAs having coverage below the district-level target of 80% and all at or below the national reported administrative coverage of 63% [2, 5]. While some LGAs were estimated to have coverage >80% for BCG (4 LGAs), DPT1 (2), OPV0 (2), OPV1 (3), and measles (2), the median coverage for all individual antigens across all LGAs was still <50%. Substantial differences in DPT3 coverage by survey compared with administrative reporting were observed in 38 (95%) LGAs, consistent with the findings from state-level coverage surveys [7, 8]. Our results reinforce the limitation of using administrative vaccination coverage to assess RI service delivery performance in these settings [7].
We observed a wide range in routine vaccination coverage for all antigens across the LGAs surveyed (e.g., DPT3 coverage ranging 1%—63%) and even among LGAs within the same state, similar to the results of a 2014 coverage survey completed in Kano state where DPT3 coverage by LGA ranged from 9%–80% (Dale Rhoda, KANRICS Task Force and Biostat Global Consulting, personal communication). These findings emphasize the need for LGA-level program evaluation to identify specific programmatic issues and appropriately target improvement in the respective LGAs.
Despite varying results among LGAs, the overall survey findings point to both factors of supply and demand in contributing to low vaccination coverage in Nigeria, as seen in other weak health systems [14]. Lack of knowledge or education about vaccines and vaccination services was the most commonly reported reason for non-vaccination, and knowledge about the safety and importance of vaccines were identified as the most important types of information needed to help mothers/caretakers decide to vaccinate their children. These findings highlight the need for targeted social mobilization strategies to increase the demand for RI services. The low reported use of RI services and the common identification of issues relating to access to RI services as a reason for non-vaccination (i.e., inconvenient location and time of vaccination, vaccine stock outs, inconsistent outreach schedules, and inability to pay for vaccination services) suggest more can be done to improve vaccine delivery mechanisms.
Significant strengths of this survey are that LGA-level data was gathered, analyzed, and shared for program action in LGAs designated as high-risk for polio. Using an efficient and rapidly deployable human resource infrastructure that exists within NSTOP and NFELTP, we were able to assure high quality data collection while conducting multiple concurrent LGA surveys. Local NSTOP and government immunization staff had the first-hand opportunity to assess issues contributing to poor vaccination coverage within their LGAs and inform prioritization and intensification of resources for RI service delivery. Detailed reports of survey results for individual LGAs were disseminated within 6 months to each LGA immunization program team for use in work plans for targeted improvement (S4 Appendix).
The survey analyses had limitations. Since we purposively sampled high-risk LGAs within eight states, we expected low vaccination coverage and cannot generalize results to larger geographic areas. Additionally, the sampling of children within clusters was not random, therefore enrolled households may not be representative of the cluster. The selection process, and the assumption of a self-weighting sample likely yields biased point estimates and an underestimate of the standard errors. However, it is unlikely that bias would explain the large discrepancy between administrative coverage and the survey estimates. Challenges during field operations that led to an overall decrease in sample size included: clusters with fewer than seven eligible households, clusters with insecurity issues prohibiting field work (including early withdrawal of survey teams before completing some clusters), and inaccurate estimation of children’s age during enrollment that led to exclusion based on ineligible age during analysis. The surveys were not powered for sub-population analysis, so inferential associations cannot be made about factors associated with low coverage. Overall card availability was low (median value of 45%), and reliance on maternal recall in the absence of a vaccination card can be unreliable and lead to over- or under-estimation of coverage estimates [15]. Future surveys could incorporate verification of vaccination status from health facility records, if available, to improve the documentation of vaccination status.
Conclusions
Vaccine preventable diseases continue to contribute significantly to the morbidity and mortality burden in Nigeria. Given the limitations of using administrative vaccination coverage data and the heavy resource demand required to conduct regular household coverage surveys, practical means of gathering high quality RI data are needed to inform program planning at the LGA level. Our approach could serve as an example of efficient human resource utilization for conducting quality LGA-level coverage surveys, particularly in countries with sub-optimal vaccine coverage and heterogeneous populations similar to Nigeria. We recommend using coverage surveys as capacity building exercises for immunization staff where possible, while incorporating the improved sampling methodology outlined in the recently updated WHO cluster survey guidelines [16], to achieve high quality results.
Supporting Information
S1 Appendix. Nigeria Routine Immunization Coverage Survey Questionnaire, Phase 1 For reference, this is an example of a blank questionnaire used during Phase 1.
https://doi.org/10.1371/journal.pone.0167835.s001
(PDF)
S2 Appendix. Data tables of antigen coverage estimates, by LGA
https://doi.org/10.1371/journal.pone.0167835.s002
(PDF)
S3 Appendix. Graphs of routine immunization coverage estimates for each LGA, grouped by antigen All coverage estimates combine maternal recall + vaccine card data; complete coverage = 8 antigens (BCG, OPV 1, DPT 1, OPV 2, DPT 2, OPV 3, DPT 3, Measles); does not include OPV or measles doses from SIAs.
The order of the LGAs within a state is based on DPT3 coverage. LGAs are grouped by state to illustrate variability in coverage across LGAs within the same state. This data is not representative of state-level coverage since LGAs were purposefully selected. Graphs of routine immunization coverage estimates, grouped by state. All coverage estimates combine maternal recall + vaccine card data; complete coverage = 8 antigens (BCG, OPV 1, DPT 1, OPV 2, DPT 2, OPV 3, DPT 3, Measles); does not include OPV or measles doses from SIAs.
https://doi.org/10.1371/journal.pone.0167835.s003
(PDF)
S4 Appendix. Survey results shared with Maru LGA immunization staff.
This is an example of the slide sets that were distributed to each LGA after data analysis was completed.
https://doi.org/10.1371/journal.pone.0167835.s004
(PDF)
Acknowledgments
We would like to thank the NFELTP residents and graduates who participated during field work implementation; Dr Namadi M Lawal from NPHCDA for his support during the planning phases; Dr Frank Mahoney and Dr Steve McCracken from U.S. CDC for their assistance during the planning and field work implementation phases.
Author Contributions
- Conceptualization: IUO HMS HE CJO MM PN JFV.
- Data curation: RG IUO OJA KAW AR.
- Formal analysis: RG IUO OJA KAW AR.
- Funding acquisition: IUO CJO PN JFV.
- Investigation: RG IUO OJA BVU NEW.
- Methodology: IUO HMS KAW HE CJO.
- Software: OJA AR.
- Supervision: RG IUO NEW JFV.
- Validation: AR.
- Writing – original draft: RG.
- Writing – review & editing: RG IUO OJA HMS KAW.
References
- 1.
World Health Organization. Immunization surveillance, assessment and monitoring: Vaccine-preventable diseases. Available: http://apps.who.int/immunization_monitoring/diseases/en/. Accessed July 2015.
- 2.
World Health Organization. Immunization, Vaccines and Biologicals. Available: http://www.who.int/immunization/global_vaccine_action_plan/en/. Accessed July 2015.
- 3. Etsano A, Gunnala R, Shuaib F, Damisa E, Mkanda P, Ticha J, et al. Progress toward poliomyelitis eradication–Nigeria, January 2014-July 2015. MMWR Morb Mortal Wkly Rep. 2015; 64:878–882. pmid:26292207
- 4.
National Primary Healthcare Development Agency. Nigeria polio eradication emergency plan 2015. Abuja, Nigeria (2015). Available: http://www.polioeradication.org/Portals/0/Document/Aboutus/Governance/IMB/12IMBMeeting/4.2_12IMB.pdf
- 5.
Federal Ministry of Health, Nigeria. Comprehensive EPI Multi-Year Plan, 2011–2015. Abuja, Nigeria (2011). Available: http://www.nationalplanningcycles.org/sites/default/files/planning_cycle_repository/nigeria/cmyp_2011-2015_country_final.pdf
- 6.
World Health Organization. Immunization, Vaccines and Biologicals; Data, Statistics and Graphics. Available: http://www.who.int/immunization/monitoring_surveillance/data/nga.pdf. Accessed February 2016.
- 7. Dunkle SE, Wallace AS, MacNeil A, Mustafa M, Gasasira A, Ali D, et al. Limitations of using administratively reported immunization data for monitoring routine immunization system performance in Nigeria. J Infect Dis. 2014 Nov 1; 210 Suppl 1:S523–30.
- 8.
National Population Commission (NPC) [Nigeria] and ICF International. Nigeria Demographic and Health Survey 2013. Abuja, Nigeria, and Rockville, Maryland, USA (2014): NPC and ICF International.
- 9. Waziri NE, Ohuabunwo CJ, Nguku PM, Ogbuanu IU, Gidado S, Biya O, et al. Polio eradication in Nigeria and the role of the National Stop Transmission of Polio program, 2012–2013. J Infect Dis. 2014 Nov 1; 210 Suppl 1:S111–S117.
- 10. Upfill-Brown AM, Lyons HM, Pate MA, Shuaib F, Baig S, Hu H, et al. Predictive spatial risk model of poliovirus to aid prioritization and hasten eradication in Nigeria. BMC Med. 2014; 12:92 pmid:24894345
- 11.
World Health Organization (WHO) Department of Immunization, Vaccines and Biologicals. Immunization coverage cluster survey—Reference manual. Geneva, Switzerland (2005). Available: http://apps.who.int/iris/bitstream/10665/69087/1/WHO_IVB_04.23.pdf
- 12. Shuaib F, Gunnala R, Musa EO, Mahoney FJ, Oguntimehin O, Nguku PM, et al. Ebola virus disease outbreak—Nigeria, July–September 2014. MMWR Morb Mortal Wkly Rep. 2014; 63:867–872. pmid:25275332
- 13.
R Core Team (2014). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available at http://www.R-project.org/
- 14. Johri M, Pérez MC, Arsenault C, Sharma JK, Pant Pai N, Pahwa S, et al. Strategies to increase the demand for childhood vaccination in low- and middle-income countries: a systematic review and meta-analysis. Bulletin of the World Health Organization. 2015; 93:339–346C. pmid:26229205
- 15. Miles M, Ryman TK, Dietz V, Zell E, Lulman ET. Validity of vaccination cards and parental recall to estimate vaccination coverage: a systematic review of the literature. Vaccine. 2013; 31: 1560‒1568. pmid:23196207
- 16.
World Health Organization (WHO). Vaccine coverage cluster surveys: Reference manual. Geneva, Switzerland (2015). Available: http://www.who.int/immunization/monitoring_surveillance/Vaccination_coverage_cluster_survey.pdf