Background. Childhood vaccination continues to increase dramatically. In spite of the success of immunization programs to date, millions of children continued to die each year, and sub-Saharan Africa (SSA) accounted for the world’s highest neonatal deaths. Childhood vaccination was designed as one of the most effective ways to reduce child mortalities from fatal vaccine-preventable diseases. Therefore, this study is aimed at investigating the individual- and community-level determinants of childhood complete vaccination in Ethiopia. Methods. A secondary data analysis was done based on the 2016 Ethiopian Demographic and Health Survey (EDHS). A total weighted sample of 1,984 children aged 12-23 months was included for analysis. Considering the hierarchical nature of EDHS data, a two-level multilevel analysis for assessing individual- and community-level determinants of childhood complete vaccination was done. The intraclass correlation coefficient (ICC), Median Odds Ratio (MOR), Proportional Change in Variance (PCV), and deviance (-2LL) were used for model comparison and for checking model fitness. Variables with p value < 0.2 in the bivariable multilevel analysis were considered for the multivariable multilevel analysis. In the multivariable multilevel logistic regression analysis, the Adjusted Odds Ratio (AOR) with 95% Confidence Interval (CI) was reported to declare significant determinants of complete childhood vaccination. Results. Overall complete vaccination status among children aged 12-23 months was 39% (95% CI: 36.8, 41.2). In the multilevel analysis, secondary or above educated mothers (AOR=2.48; 95% CI: 1.41, 4.36), richest wealth status (AOR=2.24; 95% CI: 1.16, 4.32), ≥four ANC visits (AOR=2.77; 95% CI: 1.90-4.02), employed mothers (AOR=1.66; 95% CI: 1.26, 2.18), urban residence (AOR=1.84; 95% CI: 1.00, 3.51), and children in city administration (AOR=2.66; 9% CI: 1.53, 4.62) were positively associated with vaccination status. On the other hand, children with a female household head (AOR=0.68; 95% CI: 0.48, 0.96) were negatively associated. Conclusion. Overall, childhood full vaccination status was low compared with the WHO targets. Maternal education, wealth status, ANC visit, maternal occupation, residence, region, and sex of household head were significant predictors of childhood complete vaccination. As a result, it is better to design a compensation mechanism to the costs associated with childhood vaccination for the poor households and strengthen awareness creation for rural residents to improve the access, utilization, and continuum of vaccination service.
A secondary data analysis was done based on the 2016 Ethiopian Demographic and Health Survey (EDHS). The 2016 EDHS is the fourth DHS survey in Ethiopia collected in the nine administrative regions, namely, Tigray, Afar, Amhara, Benishangul-Gumuz, Gambela, Harari, Oromia, Somali, and Southern Nations, Nationalities, and People's Region (SNNP), and two city administrative regions (Addis Ababa and Dire Dawa) (Figure 1). Most of the population of Ethiopia is an agrarian society, and about 43 percent of the Gross Domestic Product (GDP) of the country has been accounted for by agriculture, and 84% of the population lives in rural areas. More than 80 percent of the country's total population lives in the regional states of Amhara, Oromia, and SNNP [24]. In addition, Ethiopia is the thirteenth in the world and the second in Africa's most populous countries with a 4.46 fertility rate. Ethiopia has followed 3 tiers of preventive healthcare system approaches. These are primary-level healthcare comprising of a primary hospital, health center, and health post; secondary-level healthcare (general hospital); and tertiary-level healthcare (specialized hospital). The number of hospitals varies from region to region in relation to the size of the population. The Oromia region has the highest number of hospitals (30), and only one hospital is found in the Gambela region [25]. Map of the nine regions and two city administrations of Ethiopia, 2016. All children aged 12-23 months in Ethiopia were the source population, while children aged 12-23 months in the selected enumeration areas before five years of the survey were the study population. In EDHS, a stratified two-stage cluster sampling technique was employed to select the study participants using the 2007 Population and Housing Census (PHC) as a sampling frame. Stratification was achieved by separating each region into urban and rural areas. A total of 21 sampling strata have been created because the Addis Ababa region is entirely urban. In the first stage, 645 enumeration areas (EAs) were selected proportionally to the size of EAs with independent selection from the sampling stratum. Consecutively, a complete list of the households was carried out in all selected EAs before the actual data collection period, and 28 households were selected using a systematic sampling technique. Of these, 18,008 households and 16,583 eligible women were included, and the detailed sampling procedure was presented in the full 2016 EDHS report [6]. The full vaccination status of children in Ethiopia was the outcome variable for this study. The EDHS data about vaccination were collected from verbal reports of the mother and data extraction from the childhood immunization card. The detailed report is found in the 2016 EDHS report. A child that received one dose of BCG, 3 doses of pentavalent, 3 doses of polio, two doses of rota, three doses of PCV, and one dose of measles was considered fully vaccinated and categorized as “yes,” and the remaining was categorized as “no.” The response variable for the ith child is represented by a random variable Yi with two possible values coded as 1 and 0. As a result, the response variable of the ith child Yi was measured as a dichotomous variable with possible values of “Yi = 1” if the ith child was fully vaccinated and “Yi = 0” if the child was not fully vaccinated. Accordingly, consistent with the objective of the study and the hierarchical structural nature of the EDHS data, the women were nested within the cluster/community. As a result, we considered two levels of the independent variables. These are the individual-level (level one) factors including individual sociodemographic and economic factors, such as age, maternal education, paternal education, media exposure, wealth index, maternal occupation, health insurance, and sex of head of the household, and maternal obstetric-related factors, such as ANC visit, parity, preceding birth interval, and birth order. On the other hand, the community-level (level two) factors include the characteristics of the community such as the region, residence, community media exposure, community women education, region, and place of residence. This might help us to see whether the cluster-level variables had an effect on full childhood vaccination status. Community-level variables had two sources. These are direct community- and aggregated community-level variables used in the analysis of community-level variables. Some of the community-level variables such as community education and media exposure were an aggregate result of the individual data and categorized as low or high using the median value since the EDH data were not normally distributed. The data were weighted using sampling weight, primary sampling unit, and strata before any statistical analyses to restore the representativeness of the survey and to get reliable statistical estimates. Cross-tabulations and summary statistics were conducted to describe the study population using STATA version 14 software. In the EDHS data, children are nested within a cluster, and children within the same cluster might be more similar to each other than children in the rest of the country. This violates the assumption of a simple binary logistic regression model such as the independence of observations and equal variance across clusters. Therefore, a multilevel logistic regression model (both fixed and random effect) was fitted to take into account the clustering effect. Model comparison was done based on deviance since the models were nested. The likelihood ratio, intraclass correlation coefficient (ICC), Median Odds Ratio (MOR), and Proportional Change in Variance (PCV) were computed to measure the variation between clusters. The intraclass correlation coefficient (ICC) quantifies the degree of heterogeneity of childhood full vaccination between the clusters, or it is the proportion of the total observed individual variation in childhood full vaccination attributed to between-cluster variations. ICC = б2/(б2 + π2/3) [26], but MOR is quantifying the variation or heterogeneity in outcomes between clusters and defined as the median value of the odds ratio between the cluster at high likelihood of full vaccination and cluster at lower likelihood when randomly picking out two clusters or EAs. MOR=exp 2∗∂2∗0.6745~MOR=exp 0.95∗∂ [27]. ∂2 indicates the cluster variance, and PCV measures the total variation attributed to individual- and community-level factors in the multilevel model as compared to the null model: A two-level, at individual and community (cluster) levels, multilevel multivariable logistic regression was used to analyze factors associated with childhood full vaccination. Furthermore, four models were constructed for the multilevel logistic regression analysis. The first model was an empty model without any explanatory variable/s to determine the extent of cluster variation on full vaccination. The second model was adjusted with individual-level variables; the third model was adjusted for community-level variables, while the fourth was fitted with both individual- and community-level variables simultaneously. As a result, a model with the lowest deviance was chosen. Multicollinearity was also checked using the Variance Inflation Factor (VIF), and VIF < 10 and tolerance greater than 0.1 were used to declare the absence of multicollinearity. Variables with p values < 0.2 in the bivariable analysis for both individual- and community-level factors were fitted in the multivariable model. Adjusted Odds Ratio (AOR) with a 95% Confidence Interval (CI) and p values < 0.05 in the multivariable model were used to declare a significant association with childhood full vaccination.
N/A