Objective This study aims to identify the individual and contextual factors consistently associated with utilisation of essential maternal and child health services in Nigeria across time and household geolocation. Design, setting and participants Secondary data from five nationally representative household surveys conducted in Nigeria from 2003 to 2018 were used in this study. The study participants are women and children depending on essential maternal and child health (MCH) services. Outcome measures The outcome measures were indicators of whether participants used each of the following essential MCH services: antenatal care, facility-based delivery, modern contraceptive use, childhood immunisations (BCG, diphtheria, tetanus, pertussis/Pentavalent and measles) and treatments of childhood illnesses (fever, cough and diarrhoea). Methods We estimated generalised additive models with logit links and smoothing terms for households’ geolocation and survey years. Results Higher maternal education and households’ wealth were significantly associated with utilisation of all types of essential MCH services (p<0.05). On the other hand, households with more children under 5 years of age and in poor communities were significantly less likely to use essential MCH services (p<0.05). Except for childhood immunisations, greater access to transport was positively associated with utilisation (p<0.05). Households with longer travel times to the most accessible health facility were less likely to use all types of essential MCH services (p<0.05), except modern contraceptive use and treatment of childhood fever and/or cough. Conclusion This study adds to the evidence that maternal education and household wealth status are consistently associated with utilisation of essential MCH services across time and space. To increase utilisation of essential MCH services across different geolocations, interventions targeting poor communities and households with more children under 5 years of age should be appropriately designed. Moreover, additional interventions should prioritise to reduce inequities of essential MCH service utilisation between the wealth quantiles and between education status.
This study uses secondary data from national representative cross-sectional surveys of geolocated households conducted in Nigeria from 2003 to 2018. We combined publicly available data from four Nigeria Demographic Health Surveys (DHS) in 2003, 2008, 2013 and 2018 and Multiple Indicator Cluster Surveys (MICS) in 2016–2017. DHS data were extracted from IPUMS DHS.16 Detailed methodologies of the four DHS surveys are published elsewhere.17–21 All surveys employed stratified two-stage or three-stage cluster sampling techniques. The primary sampling unit (PSU) for DHS 2003 was defined as one or more enumeration areas (EAs) used for Population and Housing Census 1991, while the PSU for DHS 2008, DHS 2013, DHS 2018 and MICS 2016/17 was defined as one or more EAs used for the Population and Housing Census 2006. The counts of households interviewed in DHS 2003, DHS 2008, DHS 2013, DHS 2018 and MICS 2016/17 are shown in table 1. In DHS 2018, 11 of 27 local government areas in Borno State were excluded from the sampling frame due to insecurity in those districts. Likewise, in MICS 2016/17, a total of 101 EAs in Borno, Yobe and Adamawa states were not surveyed due to insecurity. Total number of households interviewed in DHS 2003, 2008, 2013 and 2018 and MICS 2016/17 DHS, Demographic Health Survey; MICS, Multiple Indicator Cluster Survey. To protect the confidentiality of PSU geolocations, the Global Positioning System coordinates of urban PSU locations were randomly displaced within a 2 km buffer, and rural PSUs were displaced within a 5 km buffer (and in 1% of cases, a 10 km buffer). The direction and distance of the displacement for each PSU was randomly selected using a uniform distribution.22 23 Prior research found that the effect of random displacement across a 10 km2 grid to be negligible for estimating measles vaccination coverage.24 Geolocation data for 16 of 3533 PSUs (0.5%) were missing across the four DHS. Similarly, geolocation data for 1 of 2239 PSUs (0.0004%) was missing in MICS 2016/17. After initial random displacement, 14 PSUs (1 in DHS 2008 and 13 in MICS 2016/17) were ‘located’ either in the sea or out of country’s boundaries. We resampled the random displacement of those PSUs until their displaced positions lay inside the relevant boundaries (using a 5 km buffer if possible, and a 10 km buffer if necessary). Of these 14 PSUs, 8 were successfully resampled and 6 cases that could not be appropriately displaced across 10 000 attempts were discarded. The essential maternal and child health services considered in this study consist of ANC, facility-based delivery, modern contraceptive use, childhood immunisations (BCG, first and third diphtheria, tetanus, pertussis/pentavalent, measles) and treatments for childhood illnesses. The target group for ANC and facility-based delivery was women aged 15–49 years having given a live birth in the last 23 months, while that for modern contraceptive use was women aged 15–49 years not having wanted to have more children. Children aged 12–23 months and aged 0–59 months were the target groups for immunisation and treatments of childhood illnesses, respectively. Table 2 provides further details on the definitions of and study populations for these essential services. Definitions and target populations of essential health services DPT, diphtheria, tetanus, pertussis. Independent variables across essential maternal and child health services were selected based on three earlier studies.25–27 We considered five types of explanatory variables that might influence health seeking behaviours: (i) individual characteristics; (ii) the built environment; (iii) neighbourhood demographics; (iv) the social environment and (v) the healthcare environment. Maternal and households’ characteristics include the explanatory variables of maternal age, household head, education level, marital status, possession of television and radio, possession of means of transport and household’s wealth index. Possession of television and radio was categorised into three groups: (i) households possessing both a television and a radio; (ii) households possessing either of them and (iii) households possessing neither of them. Possession of means of transport means was generated using possession of car, motorcycle and bike and categorised into three groups: (i) no means of transport; (ii) one means of transport and (iii) two to three means of transport. The household wealth index was the first principal component estimated by a principal component analysis on the household assets, sources of drinking water, sanitation facilities, type of fuel for cooking and materials of floors for housing units. Gridded estimates of population density provided by WorldPop were used as a proxy for the built environment.28 The proportion of households in a PSU living under the poverty line was used as a proxy for neighbourhood demographics and the social environment. As proxy for the healthcare environment, we measured each PSU’s travel time to the most accessible health facility using the friction surface developed by the Malaria Atlas Project.29 Geolocations of health facilities managed by government, community-based organisations and faith-based organisations in Nigeria were provided by the Nature Scientific database, which records locations of health facilities as of 2018.30 We assume these facilities were present in all years surveyed by DHS; however, because Nature Scientific does not provide the date on which each facility was established, some health facilities may not have existed at the time of some of the five surveys, a possible limitation of our analysis. Finally, the number of health facilities within a 20 km buffer around each PSU was employed as the proxy for the healthcare environment, indicating the availability of accessible health facility options. In analyses of the utilisation of childhood immunisations and treatments of childhood illnesses, we added additional explanatory variables related to child characteristics (ie, age, sex and birth month). Children’s ages were rounded to whole months. In addition to descriptive analyses, we estimated generalised additive models (GAMs) with logit links to identify factors associated with the utilisation of essential service v by the ith individual in the dth PSU in the jth state at the t year. The systematic component of the model of vidjt is: s(longidj , latidj ): smooth function of longitude and latitude using isotropic smooths on the sphere to account for spatial autocorrelation. s(montht ): smooth function of time trends. DHSdjt : binary indicator recording 1 if the data source for year t is DHS and 0 if not. Headidjt : binary indicator recording 1 if the household head in year t was a mother and 0 if not. Ageidjt : maternal age of the mother in a household in year t. U5 idjt : the number of children under 5 years of age in a household in year t. Educationidjt : the education level of the mother in year t. Wealthidjt : the wealth quantile of the household in year t. Maritalidjt : the marital status of the mother in year t. Mediaidjt : possession of TV and/or radio by the household in year t. Transportidjt : possession of means of transport by the household in year t. Povertydjt : proportion of the households living below the poverty line in the dth PSU of the jth state in year t. Accessdj : travel time in minutes from the household’s PSU to the most accessible health facility. Choicedj : the number of health facilities within 20 km from the household’s PSU. PopDensitydjt : population density in the dth PSU of the jth state in year t. We included childhood covariates ChildAgeijt and ChildSexijt (child’s age and sex, respectively) in the models of childhood immunisation and care seeking for common childhood illnesses. We log-transformed population density, the number of health facilities and the number of children under 5 to improve model fit and to account for diminishing marginal effects of these variables. Two of these variables—the count of health facilities within 20 km and the count of children under 5 years of age—could in some cases have a value of precisely zero, posing a problem for taking logs. Rather than adding an arbitrary positive quantity to these count variables, we directly estimate the effect of having zero children (or zero health facilities) by including dummy variables in the model to indicate cases where each is precisely zero. In turn, and without loss of generality, before logging the count of health facilities (or children), we replaced zeros with ones, so that cases in which there are zero health facilities within 20 km (or no children under 5 years of age) affect the outcome only through the dummy variable for that case.31 32 We listwise deleted missing data, which accounted for <2% of total cases. Because the guidelines for DHS 2018 recommend against using weights for estimating relationships, we do not use sampling weights in estimating the GAMs.33 However, sampling weights were used for estimating health service coverage reported in table 3. Essential health service coverage from 2003 to 2018 in Nigeria Finally, we estimated an additional eight models for each outcome as sensitivity analyses to check the robustness of our findings: Patients and/or the public were involved in the design, or conduct, or reporting, or dissemination plans of this research.