Background: Ethiopia is among the countries with the highest neonatal mortality with the rate of 37 deaths per 1000 live births. In spite of many efforts by the government and other partners, non-significant decline has been achieved in the last 15 years. Thus, identifying the determinants and causes are very crucial for policy and program improvement. However, studies are scarce in the country in general and in Jimma zone in particular.
This study was a community-based prospective follow up conducted in Jimma Zone from September 2012-December 2013. Jimma Zone is one of the 17 Zones of the Oromia Regional State of Ethiopia having a total of 17 rural districts and two town administrations. According to the 2007 national population and housing census, the Zone has a total population of 2.6 million, of which 88.7% are rural residents [6], [7]. The minimum required sample size for this study was determined by using Epi-Info V.3.5.1 by considering two sample comparisons of proportions based on the following assumptions. The outcome variable was neonatal mortality. Among all the determinants of neonatal mortality considered, educational status of mothers was found to give the largest sample size. Based on this, the prevalence of neonatal mortality among mothers having educational status of secondary or above was estimated to be 4.0% (P1 = 0.040) and among those who didn’t attend secondary education was to be 8.1% (P2 = 0.081) [8]; 95% level of confidence and 80% power were considered. A ratio of 1∶3 was used (r = 3). As multistage-clustered sampling method was used, a design effect of 2 was considered. Finally, 10% was added for non-responses and miss-to-follow up and the final sample size became 3604. Multistage-clustered sampling technique was used to identify a cohort of pregnant women to be enrolled in the follow up for the study. At first stage, the Zone was stratified as rural districts (17 in number) and town administrations (2 in number, Jimma and Agaro). Then, by considering time and logistics, 5 districts (30%) were selected by simple random sampling from the 17 districts. At second stage, all the selected 5 districts were clustered by ‘Kebeles’ (A ‘kebele’ is the smallest administrative unit having 5000 population in average) and stratified in to urban and rural ‘Kebeles’. Then, by simple random sampling method, 9 rural ‘Kebeles’ and 2 urban ‘Kebeles’ were selected from each selected district. This number of clusters (‘kebeles’) was determined based on expected number of pregnant women per ‘Kebele’. Jimma town administration and Agaro town administration have 13 and 5 ‘Kebeles’, respectively and all were included purposefully. With this, a total of 73 Clusters (‘Kebeles’) were included in the study. Then, for all selected ‘kebeles’, pregnant women were enumerated by using house-to-house visit and all obtained were enrolled in the study (Figure S1). The dependent variable for this study was neonatal mortality and the independent variables were divided into two levels. Level 2 (higher-level variables) included community or cluster level variables such as place of residence, access to health centers and access to hospitals. Level 1 (lower-level variables) included individual and household characteristics such as: socio-demography, wealth quintiles, maternal obstetric factors, maternal health care use, conditions of labor, characteristics of the neonates and neonatal care practices. The detail descriptions and measurements are given below (Table 1). The data were collected by using pre-tested interviewer administered structured questionnaires which were adapted from different literatures. The indicators for the wealth index were adapted from Ethiopian Demographic and Health Survey (EDHS) [5]. Indicators to measure birth preparedness and complication readiness (BP & CR) were adapted from the safe motherhood questionnaires developed by maternal and neonatal health program of Johns Hopkins Program for International Education in Gynecology and Obstetrics (JHPIEGO) [9]. Indicators for neonatal care practices were adapted from the World Health Organization (WHO) minimum neonatal care packages [10]. Data on causes of neonatal death were collected by using structured verbal autopsy questionnaire adapted from the standard VA questionnaire developed and validated by WHO, Johns Hopkins University (JHU) and London School of Hygiene and Tropical Medicine [11]. All the questionnaires were prepared in English, then translated to local languages ‘Afan Oromoo’ and Amharic and used to collect the data after back translating to English by different experts to check its consistency. As this was prospective follow up study, data were collected in three phases. First, home-to-home visit was made to enumerate pregnant women from the selected 73 clusters. Then, all the identified pregnant women were enrolled in the study as a cohort. At a baseline, data on basic socio-demography, economy and birth preparedness and complication readiness were collected. Then, just at the end of neonatal period, maternal service use (antenatal care (ANC), delivery place and attendant and postnatal care), conditions of labor, neonatal characteristics and neonatal care practices were collected. For died neonates, VAs were conducted within 15–30 days of death. Females, who had completed 10th grade or above were recruited, trained and collected the data. The VAs were conducted by two experienced females. The data collection process was supervised strictly by trained supervisors and principal investigators. To control the quality of data, in addition to training, pretest, supervision and use of local languages, the inter-item consistency of the indicators to measure the composite score of wealth index, BP & CR and neonatal care practices were checked by using Chronbach-alpha at 0.7 cut-off points. The collected data were coded and entered into Epidata V.3.1 to minimize logical errors and design skipping patterns. Then, the data were exported to SPSS for windows version 20.0 for cleaning, editing and analysis. Descriptive analysis was done by computing proportions and summary statistics. Socioeconomic quintiles were determined by using Principal Component Analysis (PCA). Birth preparedness and complication readiness was computed by composite indicator of five items. Similarly, neonatal care practice was determined by composite variable of 12 items by using PCA. As Jimma and Agaro town administrations were both purposefully included, the status of neonatal mortality was estimated by calculating weighted percentage based on the complex sample survey procedure to avoid underestimation. Bivariate analysis was done by using cross-tabulation to see associations between the dependent and independent variables. Then, all variables having P-value 10 considered as existence of multicollinearity) before interpreting the final output. However, only skill of delivery attendant (VIF = 10.9) had multicollinearity with place of delivery (VIF = 9.1, reduced to 1.8 when delivery attendant was dropped). As a result, they were included in the model alternatively by dropping the other. For the rest of the variables, the VIF was 0.05 for each). The VAs were interpreted by two independent pediatricians and third pediatrician interpreted in case of disagreements. Ethical approval was obtained from the Institutional Review Board (IRB) of College of Health Sciences of Addis Ababa University as well as IRB of Oromia Regional State Health Bureau. Following this, formal letters and permissions were secured from all respective local administrators. Written informed consent was obtained from each respondent before actual data collection. Issues of confidentiality were maintained by removing any identifiers from the questionnaire. To protect vulnerable group, data collectors were trained to maintain confidentiality and provide necessary health information based on the need of the participants and arrange referral to health facilities for sick neonates.