Abstract

Background. Premature birth occurs before 37 completed weeks of gestation. It has a greater risk of developmental disabilities, health, and growth problems than full birth. It is the second leading cause of morbidity and mortality among under-five children. Therefore, the aim of this study was to identify determinants of the survival time of premature neonates admitted to the neonatal intensive care unit (NICU) at Shambu General Hospital. Methods. A retrospective study design was used. Data were collected from medical records of premature neonates from January 2018 to March 2021. A total of 361 premature neonates were included in the study. Descriptive statistics, Kaplan–Meier (KM) curve and log-rank test were computed. The survival time of preterm neonates were compared for different categorical covariates. The Cox’s proportional hazard model was fitted. The fitness and statistical assumptions of the model were checked. Parametric regression models were compared. Weibull regression model was fitted for premature data to identify the predictors of death time of the premature neonates. Results. The proportion of premature neonatal death was 23.3%. Gestational age, neonatal sex, place of residence, hemoglobin (Hb) level, hypertension status, HIV status, antenatal care, mode of delivery, birth weight, multiple pregnancies, perinatal asphyxia, and parity greater than 1 were significantly associated with the death time of premature neonates. Conclusion. Percentage of premature neonatal death in this study was 23.3. Improving mothers’ Hb level through routine iron supplementation, encouraging mothers to have regular antenatal follow-up at health institution were recommended.

1. Introduction

Premature births are born fragile, small, and weigh less than full term birth. It is major determinant of neonatal morbidity and mortality. It can be divided into late preterm (34–37) weeks, moderate preterm (32–34) weeks, very preterm (28–32) weeks, and extremely preterm (less than 28) weeks [1, 2].

Preterm birth commonly occurs during the third trimester of pregnancy [3]. It is the major cause of perinatal morbidity [4] and is the second leading cause of neonatal death [5]. It accounts 75% of the deaths among premature neonates in the modern and developing countries [6] resulting in 80% of deaths of neonates without congenital abnormalities. Preterm birth accounts for 3.1% of all disability adjusted life years [5].

Globally, 15-million babies born before 37 completed weeks of gestation every year and more than 1-million die due to complications related to preterm birth [7, 8]. Rates are highest in low and middle income countries [9]. Sub-Saharan Africa and south Asia account for over 60% of preterm birth worldwide. In Ethiopia, preterm birth accounts for 28% of all other causes of neonatal death. From 320,000 babies born too soon each year 24,000 children under 5-year die due to preterm complication in Ethiopia [10].

Preterm neonates born within preterm category shares similar risk for death, weight, size, and are at higher risk for health and developmental problems compared to the neonates born full term. They experience difficulty in feeding, blood glucose control, jaundice, temperature instability, apnoea, respiratory distress and sepsis. Preterm neonates are at a higher risk of cognitive and behavioral disorders compared to the other newborns. These complications are associated with genetic influences, infertility treatments, multiple pregnancies, infections, and chronic conditions such as diabetes and high blood pressure, [11] environmental exposure, medical conditions of the mother or fetus, behavioral and socioeconomic factors [12].

2. Methodology

A retrospective study design was used. All medical records of premature neonates who were admitted to NICU at Shambu General Hospital from January 2018 to March 2021 were collected by medical professional. Data were extracted for 361 neonates using a structured checklist. The dependent variable is time from the beginning of follow-up until death occurs; the study ends or the participant is lost follow-up. The independent variables were place of residence, gestational age, neonatal sex, Hb level of mother, hypertension status of mother, antenatal care of mother during pregnancy, jaundice, HIV status of mother, mode of delivery, multiple pregnancy, perinatal asphyxia (PNA), weight of neonate during birth, parity, hypothermia, sepsis, gravidity, and breast feed initiated within 1 hr. Data were analyzed by R statistical package.

2.1. Survival Analysis

Survival model is a method of data analysis where the outcome variable is time from the beginning of follow-up until an event occurs [13]. It takes into account when some subjects are lost to follow-up or when the period of observation is finite and certain subjects may not experience the event of interest over the study period.

2.1.1. The Survivor Function S (t)

Let T be a random variable associated with the survival times of premature neonates and f (t) be the probability density function of the survival time. The cumulative distribution function F(t) represents the probability that premature neonate selected at random will have a survival time less than t is given by:

The survivor function is the probability that the survival time of a randomly selected premature neonate is greater than or equal to some specified time t.

2.1.2. The Hazard Function

The hazard function is the risk of death at time t, and is obtained from the probability that a preterm neonate dies at time t given that it has survived up to time t. The hazard function is given as follows:

The hazard function can be expressed in terms of probability density function and the survivor function as follows:

2.1.3. Kaplan–Meier (KM) Estimator

Suppose that be the survival time of n independent premature neonates and be the m distinct ordered death times. The KM estimates of survival time t is given by [14]:with the convention that for where is the number of premature neonates who are at the risk of dying at time t and is the number of premature neonates died at time t. The variance of the KM estimate is [15]:

2.2. Comparison of Survival Curves

The survival time of premature neonates for the different groups can be compared by plotting the KM estimators of the groups on the same axes. The graph shows that the survival curve lying above had more survival experience than the group defined by the lower curve. The test statistics for the comparison of survival time between groups can be defined as follows:where m is the number of rank-ordered survival times. is the observed number of death in Group 1 at failure time .is the expected number of deaths corresponding in Group 1 at time .is the variance of the number of deaths in Group 1 at time .

The test statistic Q has χ2 distribution when the total number of observed events and sum of expected number of events are large and assuming that the censoring experience is independent of group [16, 17].

2.2.1. Cox Proportional Hazard Model

Cox proportional hazard model gives a hazard at time t for a premature neonate with a given specification of a set of explanatory variables denoted by X and it is generally given by:where is a baseline hazard function which is obtained when all X’s are set to zero, is a vector of explanatory variables associated with the premature neonate and is a vector of unknown regression parameters that are assumed to be the same for premature neonates, which measures the influence of the covariate.

The Cox model can also be regarded as a linear model for the logarithm of the relative hazard that is

The logarithm of the hazard ratio for two individuals having two distinct covariate values xj and xi can be expressed as follows:

2.3. Parametric Regression Modeling

A parametric survival model assumes that the survival time follows a known distribution. Many models using different distributions have been developed.

2.4. The Weibull Regression Model

The survival time t is a positive random variable having Weibull probability density function can be expressed as follows:where and the hazard function of the distribution becomes yielding in survivor function and cumulative hazard function Now incorporating covariates X in the hazard function, the Weibull regression models become:

The model assumes that individual i and j with covariates and have proportional hazard function of the form:

The quantity exp (β) can be interpreted as hazard ratios.

2.5. The Exponential Regression Model

For the time data and skewed to the right, with distribution of the time is exponential, the survival of time for a single covariate , which is called, accelerated failure time, expressed as follows:

This model can be linearized by taking the natural log of each side of the equation as follows:where is the error component. The exponential model (tExp(α)) is the simplest parametric model and assumes a constant hazard over time, which reflects the “memory less” property of the exponential distribution. The survivorship function may be obtained by expressing in terms of time as follows:And the hazard function of the exponential regression model is:

The exponential regression model for the covariates and ith individual preterm neonate is expressed as follows:

For exponential regression survival model the hazard ratio for the dichotomous covariate is

2.6. The Log-Logistic Regression Model

Single covariate log-logistic accelerated failure time may be expressed as follows:

The survivorship function for the model is:where z is the standardized log-time outcome variable, that is;

The odds of a survival time of at least are, assumes that the covariate is dichotomous and coded 0 or 1. The odds-ratio at time from the ratio the odds of a survival time evaluated at x = 0 and x = 1 is:which is independent of time.

3. Results

Among 361 preterm neonates admitted to Shambu General Hospital, 76.7% were discharged and 23.3% premature neonates died. The mean and median length of hospital stay was 16.8 (95% CI: 16.12, 17.49) and 18 days, respectively. The minimum and maximum time of hospital stay was 1 and 28 days, respectively. The result in Table 1 revealed that about 64.8% and 62.3% of premature neonates were female and rural resident, respectively. The proportion of premature neonatal death born in gestational age less than 28, (28–32), (32–34), and (34–37) weeks were 59%, 23%, 11.8%, and 6.7%, respectively. From 50.4% of premature neonates weighting greater than or equal to 1,600 g at the time of birth 35.7% were died. More than 37% and 24% of premature neonates delivered from HIV positive and HIV negative mothers were died.

Majority of premature neonates (70.6%) were delivered from mother who had Hb level greater than or equal to 11 g/dl. About 35.8% and 23.5% of premature neonates delivered from mothers who had Hb level less than and more than 11 g/dl were died, respectively. Majority of mothers (81.4%) followed antenatal visit during pregnancy. About 26.9% premature neonates delivered from mothers who did not have antenatal care visit during pregnancy were died. Rate of neonatal death varies for the different type delivery method. Among preterm neonates, 78.9% were born via spontaneous vaginal delivery (SVD). About 27.3% and 27% of premature neonates delivered by vaginal and Caesarian section (C/S) were died, respectively. Majority of premature neonates (87%) were delivered singleton. 27.1 and 27.7% of premature neonates delivered singleton and delivered twins were died, respectively. About 29.7% of premature neonates delivered from hypertensive mothers were died.

3.1. Comparison of Survivor Function

KM curve is a decreasing step functions as time increases illustrated Figure 1. It gives the probability that the survival time of premature neonate exceeds the specified day. During the first day of hospital stay the maximum survival probability of 0.994 (95% CI: 0.987–1.000) was observed with a standard error of 0.00391, at the 26th days of hospital stay the probability of survival time of premature neonates was 0.409 (95% CI: 0.299, 0.519) with a standard error of 0.05625. Separate KM curve were shown for different categorical covariates as shown in Figures 210.

KM curves of neonatal sex, Hb level of mothers, multiple pregnancies, HIV status of mother, weight of neonate, PNA, mode of delivery, hypertension status of mother, and gestational age were shown in Figures 210. The survival time of premature neonates who are female, had weight ≥1,600 g, high-gestational age, whose mother had Hb level ≥11 g/dl, had PNA, born singleton, delivered by C/S, delivered from HIV negative, and nonhypertensive mothers were consistently greater than the male premature neonates, weight less 1,600 g and lower gestational age, Hb level less than 11 g/dl, did not have PNA, born multiple, delivered with SVD, delivered from HIV positive, and hypertensive mothers, respectively.

Log-rank test evaluates whether KM curves of categories of covariates are statistically equivalent or not. This is used to test the hypothesis that survival time of premature neonates for the different categories of predictor variables are equal. Both log-rank and Generalized Wilcoxon test given in Table 2 revealed that there is a significant difference in survival experience among categories of place of residence, gestational age, neonatal sex, Hb level of mother, hypertension status of mother, antenatal care, jaundice, HIV status of mother, mode of delivery, multiple pregnancy, PNA, weight of neonate during birth, and parity at 5% level of significance.

This mean that the hypothesis of equal survival time is rejected and we have enough evidence to say that the KM curves of different categories of covariates are significantly different at 5%.

But there is no difference in survival experience of premature neonates among groups of hypothermia, sepsis, gravidity, and breast feed initiated within 1 hr.

3.2. Cox’s Proportional Hazard Model

The Cox’s proportional hazard model was used to identify factors associated with the survival time of preterm neonates. The results of multivariate Cox Proportional hazard model in Table 3 revealed that gestational age, neonatal sex, place of residence, Hb level of mother, hypertension status of mother, HIV status mother, antenatal care, mode of delivery, birth weight, multiple pregnancy, parity (2–4), and PNA were significantly related with the survival time of premature neonates. However, jaundice, hypothermia, sepsis, parity > 5, gravidity, and breast feed initiated within 1 hr were not significant at 5%.

The hazard ratio (95% CI) for premature neonates who are female and lives in rural were 0.498 (0.272–0.912) and 1.733 (1.044–2.878) compared to the male premature neonates and premature who lives in urban, respectively. The risk of death of premature neonates who lives in rural is 1.733 times higher than premature neonates who reside in urban. Female premature neonates were 50.2% times less likely die compared to the male premature neonates.

The hazard ratio (95% CI) for premature neonates who were born in between gestational age (28–32), (32–34), and (34–37) weeks compared to those who were born at gestational age less than or equal to 28 weeks were 0.217 (0.078–0.601), 0.300 (0.159–0.569), and 0.285 (0.098–0.829), respectively. That is, premature neonates who were born in between gestational age (28–32), (32–34), and (34–37) weeks were 78.3, 70, and 71.5% less likely die compared to the premature neonates who were born at less than or equal to 28 weeks of gestation, respectively.

The hazard ratio (95% CI) for premature neonates who were born from mothers who had antenatal care, hypertension, Hb level greater than or equal to 11 g/dl and HIV Positive were 0.363 (0.206–0.637), 2.789 (1.456–5.343), 0.433 (0.273–0.688), and 2.527 (1.356–4.707), respectively, compared to the neonates who were born from mothers who did not have antenatal care, not hypertensive, Hb level less than 11 g/dl, and HIV negative. This means that premature neonates who were born from mothers who had antenatal follow-up were 63.7% less likely die compared to the premature neonates whose mother did not have antenatal care. The risk of death of premature neonates who were born from mothers who had hypertension was 2.789 times higher than premature neonates who were born from mothers who did not have hypertension. Premature neonates who were born from mothers who had Hb level greater than or equal to 11 g/dl were 56.7% less likely die than neonates whose mother had Hb less than 11 g/dl. Premature neonates who were born from mothers with HIV positive were 52.7% more likely die than neonates born from HIV negative mothers.

The hazard ratio (95% CI) for premature neonates who had weight less than 1,600 g, PNA, born multiple and parity (2–4) were 1.775 (1.088–2.898), 1.914 (1.125–3.256), 2.872 (1.488–5.546), and 2.113 (1.044–4.279) compared to the premature neonates who had weight greater than or equal to 1,600 g, who did not have PNA, born single and parity 1, respectively. That is, the risk of death of premature neonates who had weight less than 1,600 g, who had PNA, twin and parity 1 were 1.775, 1.94, 2.872, and 2.133 times higher than premature neonates who weights greater than or equal to 1,600g, who did not have PNA, singleton and parity 1, respectively.

Hazard ratio (95% CI) for premature neonates delivered by spontaneous vertex were 0.394 (0.217–0.714) compared to the premature neonates delivered by C/S. This means, premature neonates delivered by spontaneous vertex were 60.6% less likely die compared with the premature neonates delivered by C/S.

3.3. Assessment of Model Assumption

The assumptions of proportional hazard were checked by χ2 test based on the Schoenfeld residuals for each explanatory variables and globally (χ2 = 20.41, p-value = 0.495). In addition, the assumption of proportionality was assessed graphically by plotting the Schoenfeld residuals of each covariate against log time. There were no covariates which show a pattern with time indicating the hazard ratio was constant for time.

3.4. Model Comparison Criteria

Based on the Akaike information criterion (AIC), Weibull model (AIC = 781.1) was more efficient than exponential (AIC = 903.4), log logistic (AIC = 799.1), and lognormal (AIC = 828.3) models.

3.5. Multivariate Weibull Regression Model

The Weibull regression model in Table 4 showed that gestational age, neonatal sex, Hb level of mother, hypertension status of mother, HIV status of mother, antenatal care, mode of delivery, multiple pregnancy, birth weight, and PNA were significant determinants of the survival time premature neonatal death. Breast feed within 1 hr, place of residence, jaundice, sepsis, gravidity, parity, and hypothermia were not significant at 5%.

The hazard ratio (95% CI) for premature neonates who had gestational age between (28–32), (32–34), and (34–37) weeks and whose mother had Hb level greater than or equal to 11 g/dl compared to the premature neonates who were born at gestational age less than or equal to 28 weeks and whose mother had Hb level less than 11 g/dl were 0.758 (0.544–1.056), 0.944 (0.772–1.156), 0.790 (0.585–1.0660, and 0.816 (0.707–0.941), respectively. Premature neonates who had gestational age between (28–32), (32–34), and (34–37) weeks and whose mother had Hb level greater than or equal to 11 g/dl were 24.2, 5.6, 21, and 18.4% less likely die compared to the premature neonates born at gestational age less than or equal to 28 weeks and whose mother had Hb level less than 11 g/dl, respectively.

The hazard ratio (95% CI) for preterm neonates who were born from HIV positive mothers and whose mother had antenatal care compared to the premature neonates who were born from HIV negative mothers and whose mother did not have antenatal care were 3.284 (2.678–4.026) and 0.661 (0.558–0.784), respectively. This mean that premature neonate who were born from HIV positive mother were 3.284 more likely die than premature neonates who were born HIV negative mothers. Premature neonates who were born from mothers who had antenatal follow-up were 33.9% less likely to die compared to premature neonates who were born from mothers who did not have antenatal care. The hazard ratio (95% CI) for premature neonates who had PNA, whose mother have hypertension and born twin compared to premature neonates who did not have PNA, whose mother did not have hypertension and single were 2.145 (1.809–2.544), 3.678 (2.988–4.527), and 4.351 (3.5–5.408), respectively. That is, the risk of death of premature neonates who had PNA, whose mother had hypertension and born twin were 2.145 and 3.678; and 4.351 times higher than premature neonates did not have PNA, whose mother did not have hypertension and singleton, respectively.

The hazard ratio (95% CI) for female premature neonate and premature neonates who weight less than 1,600 g were 0.823 (0.685–0.990) and 2.59 (2.223–3.018), respectively, compared to the male premature neonate and premature neonates whose weights greater than or equal to 1,600 g. That is, the risks of death of female premature neonates were decreased by 17.7% compared with male premature neonates.

The risk of death of premature neonates weighting less than 1,600g was increased by 59% compared with the premature neonates whose weight greater than or equal to 1,600 g.

4. Discussion

During the study period, 361 premature neonates were admitted to the NICU at Shambu General Hospital.

The proportion of premature neonatal death was 23.3%. This finding is in line with studies conducted in northern Ethiopia 25.2% [18], in Addis Ababa 25.3% [19], in urban Pakistan 34% [20], in Tigray region 34% [21], and in Jimma 34.9% [22]. The overall mean (95% CI) and median length of hospital stay was 16.8 (95% CI: 16.12, 17.49) and 18 days, respectively.

The log-rank test and KM curves show that categories of place of residence, gestational age, neonatal sex, Hb level of mothers, hypertension status of mothers, jaundice, HIV status of mothers, mode of delivery, multiple pregnancy, PNA, antenatal care, birth weight, and parity were statistically difference in experiencing death event at 5% level of significance. But sepsis, breast feed within 1-hr hypothermia, and gravidity were not clearly experiencing significance differences in the death of premature neonates.

Cox proportional hazard and Weibull regression models were used to estimate the risk of death of premature neonates. Gestational age, neonatal sex, place of residence, Hb level of mother, hypertension status of mother, HIV status of mother, antenatal care visit, mode of delivery, birth weight, multiple pregnancy, parity 2–5, and PNA were identified as significant determinants of the survival time premature neonates in both Cox proportional hazard and Weibull regression models.

Birth weight affects survival time of premature neonates. Both Cox’s proportional hazard model and Weibull regression model show that the risk of death of premature neonates who had weight less than 1,600 g is about 1.775 and 2.59 higher than the premature neonates who had weight greater than 1,600 g, respectively. This result is in accordance with the previous studies conducted in northern west and central parts of the country [19, 23].

Gestational age is prognostic factor that significantly predicts the survival time of premature neonates. The risk of death of premature neonates with gestation (28–32), (32–34), and (34–37) weeks were 0.217, 0.300, 0.285 and 0.758, 0.944, 0.79 times lower than the hazard of death premature neonates with gestational age less than 28 weeks, respectively, in both Cox’s proportional hazard and Weibull regression models. This result is comparable with earlier studies conducted in Jimma University [22] Gondar University [18] and Moi University Hospital in Kenya [24]. This is because as fetal maturity increases risk of premature neonatal death will decrease.

The risk of death of premature neonates who had PNA is about 1.914 and 2.145 times higher than premature neonates who did not have PNA, respectively, in both Cox’s proportional hazard and Weibull regression models. This result is in accordance with the studies done by Yehuala et al. [18] and Wesenu et al. [22]. This is because PNA is causes of premature neonatal death where neonatal care not adequate.

Mother’s HIV status is predictor of death time of premature neonates. This study revealed that the risk of death of premature neonates whose mother had HIV was 2.527 and 3.284 higher than premature neonates whose mother was HIV negative, respectively, in both Cox’s and Weibull regression models. The present result concords with earlier study conducted in Ethiopian [18, 25] and in Uganda [26].

The survival time of premature neonate whose mother had received antenatal care were 63.7 and 33.9% lower than premature neonates whose mother did not have received antenatal care in both Cox’s and Weibull regression models, respectively. This finding is consistent with the report of the previous study having antenatal care visit significantly reduces premature neonatal death [18, 24].

Hypertension is related to the time to death of premature neonates. In both Cox’s and Weibull regression models hazard of death of premature neonate whose mother had hypertension were 2.789 and 2.988 higher than premature neonate whose mother did not have hypertension, respectively. This is similar to reports by Gebreslasie [25]. This might be because, the complications of pregnancy induced hypertension can cause vascular damage to placenta.

Hb level of mother significantly predicts the survival time of premature neonates. The hazard of death of premature neonate whose mother had Hb level less than 11 g/dl is about 56.7% and 18.4% higher than premature neonates whose mother had Hb greater than 11 g/dl in both Cox’s proportional hazard and Weibull regression model, respectively. This result is in accordance with the previous study [18].

Survival prognosis of premature neonate was lower for twins. The risk of death of premature twins were 2.872 and 4.351 higher compared to singleton, respectively, in both Cox’s and Weibull parametric models. This find was in line with studies conducted in northern Ethiopia [18].

This study indicates that hazard of death for premature neonates delivered by SVD were 0.394 and 0.782 compared to premature neonates delivered by C/S in both Cox’s proportional hazard and Weibull models. The finding was supported by [24, 27]. But this finding was contrary to the other studies where C/S was associated with better survival [18]. The variation in the results of our study is probably related to different environments of the countries that affected the variables.

There is a significant relationship between neonatal sex and survival time. Female preterm neonates had 0.498 and 0.823 times lower risk of death compared to the male premature neonates, respectively in both Cox’s proportional hazard and Weibull regression models. This find was in line with study conducted by Nalini et al. [28].

The hazard rate of premature neonates who live in rural was 1.733 times higher than the premature neonates who reside in urban in Cox’s proportional model. This result is in accordance with the study conducted in iran [28]. But the relationship between neonatal sex and death time of premature neonate was insignificant at 5% in Weibull regression model.

5. Conclusions

The main factors of preterm death were gestational age, being HIV positive, hypertension, Hb level, antenatal care, PNA, place of delivery, multiple pregnancy, and birth weight of neonates. Efforts have to be made to decrease the prevalence of preterm death. Health workers should work on encouraging mothers to have regular antenatal follow-up, improving mothers’ Hb level, controlling of hypertension, providing quality of healthcare may decrease the rate of preterm birth and its consequences.

Data Availability

Data were extracted from neonates’ records using a pretested structured checklist.

Additional Points

Limitation. There was a lack some important variables affecting the outcome variables and missing values.

Conflicts of Interest

The author declares that there is no conflicts of interest.

Acknowledgments

We are grateful to Wollega University for providing financial support.