Background: Immunization coverage in many parts of Nigeria is far from optimal, and far from equitable. Methods: Nigeria accounts for half of the deaths from Measles in Africa, the highest prevalence of circulating wild poliovirus in the world, and the country is among the ten countries in the world with vaccine coverage below 50 percent. Studies focusing on community-level determinants therefore have serious policy implications. Results: Multilevel multivariable regression analysis was used on a nationally-representative sample of women aged 15-49 years from the 2003 Nigeria Demographic and Health Survey. Multilevel regression analysis was performed with children (level 1) nested within mothers (level 2), who were in turn nested within communities (level 3). Conclusion: Results show that the pattern of full immunization clusters within families and communities, and that socio-economic characteristics are important in explaining the differentials in full immunization among the children in the study. At the individual level, ethnicity, mothers’ occupation, and mothers’ household wealth were characteristics of the mothers associated with full immunization of the children. At the community level, the proportion of mothers that had hospital delivery was a determinant of full immunization status. Significant community-level variation remaining after having controlled for child- and mother-level characteristics is indicative of a need for further research on community-levels factors, which would enable extensive tailoring of community-level interventions aimed at improving full immunization and other child health outcomes. © 2009 Antai; licensee BioMed Central Ltd.
Data on the health and mortality of children in Nigeria were collected as part of the Nigeria Demographic and Health Survey (DHS). This study uses data from the 2003 edition of this survey, which is a nationally-representative probability sample, collected using a stratified two-stage cluster sampling procedure. Sampling of women was performed according to the list of enumeration areas developed from the 1991 Population Census sampling frame. The initial sampling stage involved selecting 365 clusters, also known as primary sampling units (PSUs) with a probability proportional to the size. The size, in this case, is the number of households in the cluster. Subsequent sampling involved systematically selecting households from the already selected clusters. This resulted in a probability sample of 7864 households, from which data was collected by face-to-face interviews from 3725 women aged 15 to 49 years. These women contributed a total of 6029 live born children born to the survey. Information collected included birth histories, in-depth demographic and socio-economic information on illnesses, medical care, immunizations, and anthropometric details of children [20]. Immunization status of a child was determined from vaccination cards shown to the DHS interviewer. In the absence of vaccination cards, mothers were asked to recall whether the child had received BCG, Polio, DPT (including the number of doses for each) and Measles vaccinations. The outcome variable is the likelihood of a child 12 months of age and older having received all of the eight required vaccinations (full immunization). Eight additional child- and mother-level variables of interest were examined: i) sex of the child, assessed as: male and female; ii) birth order and interval between births, created by merging “birth order” and “preceding birth interval” classified as: first births, birth order 2-4 with short birth interval (<24 months), birth order 2-4 with medium birth interval (24-47 months), birth order 2-4 with long birth interval (48+ months), birth order 5+ with short birth interval (<24 months), birth order 5+ with medium birth interval (24-47 months), and birth order 5+ with long birth interval (48 months); iii) mothers' age, grouped as: 15-18, 19-23, 24-28, 29-33, and 34 years and older; iv) marital status, grouped as: single, married, and divorced; iv) ethnicity, categorized as: a) Hausa/Fulani/Kanuri (grouped on the basis that these ethnic groups either speak a common language or dialect, share a common sense of identity, cohesion and history; or have a single set of customs and behavioural rules as in marriage, clothing, diet, taboos); b) Igbo; c) Yoruba; and d) Others (a merger of various other minority ethnic groups from the more than 374 identifiable ethnic groups in Nigeria); v) vi) mothers' education, categorized as: no education, primary, and secondary or higher education; vii) mothers' occupation, categorized as: professional/technical/managerial, clerical/sales/services/skilled manual, agricultural self-employed/agricultural employee/household & domestic/unskilled manual occupations, and not working; and viii) mothers' household wealth index, categorized into five quintiles as: poorest, poorer, middle, richer and richest. Primary sampling units or clusters are administratively-defined areas used as proxies for "neighbourhoods" or "communities" [21,22], and are relevant when the hypothesis involves policies. Primary sampling units are small and designed to be fairly homogenous units with respect to population socio-demographic characteristics, economic status and living conditions, and consist of one or more enumeration areas (EAs), which are the smallest geographic units for which census data are available in Nigeria. Each cluster was made up of a minimum of 50 households; in the case of less than 50 households, a contiguous enumeration area was added [20]. Four community-level variables were assessed. Community prenatal care by doctor was assessed because prenatal care directly increases the chances that mothers would access subsequent health care services for their child, such as institutional delivery and immunization [23,24]. Community hospital delivery was included because the proportion of mothers that delivered in a hospital setting is a predictor of child immunization uptake. Hospital delivery is one of the most important preventive measures against maternal and child health outcomes, and an important determinant of full immunization [25,26]. Community mother's education was assessed because higher levels of maternal education are associated with better child health outcomes, such as child immunization rates [23,24]. These community-level variables were: i) community mother's education, defined as the percentage of mothers with secondary or higher education in the primary sampling unit, and categorized as: low, middle, and high (cut-off at median value in all primary sampling units combined; "middle" referring to the proportion at the median value, "low" referring to the proportion below the median value, and "high" referring to the proportion above the median value); ii) community hospital delivery, defined as the percentage of mothers who delivered their child in the hospital, and categorized as: low, middle, and high (cut-off at median value in all primary sampling units combined); iii) Community prenatal care by doctor, defined as the percentage of mothers who received prenatal care by a doctor and categorized as: low, and high (cut-off at 13% in all primary sampling units combined); and iv) mother's region of residence, categorized according to the six geo-political zones in Nigeria, as: North Central, North East, North West, South East, South South, and South West. Community-level variables were estimated at the level of the primary sampling unit (n = 365). The distribution of the children and mothers in the sample by full immunization status was assessed. Normalized sample weights provided in the DHS data were used for all analyses using Stata 10 software package [27], so as to adjust for non-response and enable generalization of findings to the general population. A three-level multilevel logistic regression model was applied in order to account for the hierarchical structure of the DHS data [28]. Children (level 1), were nested within mothers (level 2), who were in turn nested within communities (level 3). Four models containing variables of interest were fitted. Model 0 (empty model) contained no exposure variable and only focused on decomposing the total variance into its mother and community components. Model 1 contained child-level variables (sex of the child, birth order/birth interval of the children) and Model 2 included mother-level variables (mothers' age, marital status, ethnicity, mothers' education, mothers' occupation, and mothers' household wealth index). Model 3 contained community-level variables (community mother's education, community hospital delivery, community prenatal care by doctor, and mothers' region of residence). The three-level multilevel model is written as follows: where πijk is the probability of dying for the ith child of the jth mother in the kth community, eijk is a child-level error term distributed as Bernoulli constant, Xijk is a vector of covariates corresponding to the ith child of the jth mother in the kth community including mother's ethnicity, and educational background, β0 is a vector of unknown parameters, u0jk is the random effect at the mother level, and v0k is the random effect at the community level. The intercept or average probability of being fully immunized is assumed to vary randomly across mothers and communities. The fixed effects (measures of association) are expressed as odds ratio (OR) and 95% confidence intervals (95% CI). The random effects (measures of variation) are expressed as Variance Partition Coefficient (VPC) and proportional change in variance (PCV). We appraised the precision by the standard error (SE) of the explanatory variables, and tested parameters using the Wald statistic i.e. the ratio of the estimated variance to its standard error [29], and we calculated p-values. MLwiN software package 2.0.2 [30] was used for the multilevel analyses, with Binomial, Penalized Quasi-Likelihood (PQL) procedures [31]. Missing data were excluded from the analysis. This study is based on analysis of secondary data with all participant identifiers removed. The survey was approved by the National Ethics Committee in the Federal Ministry of Health, Nigeria and the Ethics Committee of the Opinion Research Corporation Macro International, Incorporated (ORC Macro Inc.), Calverton, USA. Informed consent was obtained from the participants prior to participation in the survey, and data collection was done confidentially. Permission to use the DHS data in this study was obtained from ORC Macro Inc.
N/A