Modeling Determinants of Time-To-Death in Premature Infants Admitted to Neonatal Intensive Care Unit in Jimma University Specialized Hospital

Annals of Data Science, Mar 2017

Preterm birth is the term used to define births that occur before 37 completed weeks or 259 days of gestation. The aim of this study is to model survival probability of premature infants who were under follow-up and identify significant risk factors for mortality. Recorded hospital data were obtained for a cohort of 490 infants at Jimma University Specialized Hospital, Ethiopia. The infants have been under follow-up from January 2013 to December 2015. The non-parametric, semi-parametric and parametric survival models are used to estimate the survival time as well as examine the association between the survival time with different demographic, health and risk behavior variables. The analysis shows that most factors significantly contribute to a shorter survival time of premature infants. These factors include having prenatal Asphyxia, hyaline membrane disease, sepsis, jaundice, low gestational age, respiratory distress syndrome and initial temperature. It is therefore recommended that people ought to be cognizant on the burden of these risk factors and well informed about the prematurity.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://link.springer.com/content/pdf/10.1007%2Fs40745-017-0107-2.pdf

Modeling Determinants of Time-To-Death in Premature Infants Admitted to Neonatal Intensive Care Unit in Jimma University Specialized Hospital

Million Wesenu 0 1 2 Sudhir Kulkarni 0 1 2 Tafere Tilahun 0 1 2 0 Department of Statistics, College of Natural Sciences, Jimma University , Jimma , Ethiopia 1 Department of Statistics, College of Computing and Informatics, Haramaya University , Dire Dawa , Ethiopia 2 Modeling Determinants of Time-To-Death in Premature Infants Admitted to Neonatal Intensive Care Unit in Jimma University Specialized Hospital Preterm birth is the term used to define births that occur before 37 completed weeks or 259 days of gestation. The aim of this study is to model survival probability of premature infants who were under follow-up and identify significant risk factors for mortality. Recorded hospital data were obtained for a cohort of 490 infants at Jimma University Specialized Hospital, Ethiopia. The infants have been under followup from January 2013 to December 2015. The non-parametric, semi-parametric and parametric survival models are used to estimate the survival time as well as examine the association between the survival time with different demographic, health and risk behavior variables. The analysis shows that most factors significantly contribute to a shorter survival time of premature infants. These factors include having prenatal Asphyxia, hyaline membrane disease, sepsis, jaundice, low gestational age, respiratory distress syndrome and initial temperature. It is therefore recommended that people ought to be cognizant on the burden of these risk factors and well informed about the prematurity. Premature infant; Time to death; Cox proportional hazards model; Log-logistic regression model 1 Introduction Premature infants come early into the world and they are born fragile, small and weighing less than full term infants. Many of the babies who survive face greater risks of significant health problems and disability throughout their lives (i.e, learning disabilities, visual and hearing problems, chronic lung disease and other long-term diseases) which translate into significant increased costs to healthcare, the economy and the broader society [27]. Premature is one of the major causes of infants’ death which is not an acute disease and compared to term infants experience more difficulty with feeding, blood glucose control, jaundice, temperature instability, respiratory distress and sepsis either singly or in combination [8]. Globally, an estimated 13 million infants are born before 37 completed weeks of gestation annually. Rates are generally highest in low and middle income countries and increasing in some middle and high income countries [17]. More than 1 in 10 of the world’s babies born in 2010 were born prematurely, making an estimated 15 million preterm births, of which more than 1 million died as a result of their prematurity [13]. Preterm birth accounts for 3.1% of all Disability Adjusted Life Years (DALYs) in the Global Burden of Disease, more than for HIV and malaria [21, 27]. Also, deaths constitute 28% of the 4 million annual new born deaths with 99% of these deaths occurring in developing countries [10]. Morbidity, mortality and prolonged hospital stay of preterm babies result in significant cost to the health sector, parents and the society [20]. In Ethiopia, according to report of United Nations of children fund [26], one of the main causes of neonatal death is preterm birth accounts for 23% of all other causes of neonatal death. Also, Ethiopian Demographic and Health Survey in 2011, high rate of neonatal mortality (37 deaths per 1000 live births) is reported and preterm birth is believed to be a major and direct cause of neonatal mortality. One in every 13 babies born in Ethiopia did not survive to celebrate its first birth day and one in every eight children died before its fifth birth day by Belaynew et al. [3]. Preterm birth has multiple factors whose solutions will not come through a single discovery but rather from an array of discoveries addressing multiple biological, clinical, and social behavioral risk factors. Causal factors linked to preterm birth include medical conditions of the mother or fetus, genetic influences, environmental exposure, infertility treatments, behavioral and socio-economic factors as well as iatrogenic prematurity [23]. Approximately 45–50% of preterm births are idiopathic, 30% are related to preterm rupture of membranes and another 15–20% result from medically indicated or elective preterm deliveries [11]. It is currently known that the study of risk factors for infant mortality is very important, as, particularly in the newborn, it can be considered one of the best quality indicators for health care, as well as an indicator for population social and economic welfare [24]. Neonatal death rate was overall 27.4%, which was significantly higher in gestational age subgroup of less than 28 weeks compared with other gestational age subgroups. The most prevalent etiologies of neonatal death were respiratory distress syndrome (73.8%), congenital abnormalities (13.8%) and sepsis (5.4%), respectively. Over 80% of all cases of death in the world are result of neonatology premature birth and the complications which are result of preterm birth are significant socioeconomic problem and the hyaline membrane disease (HMD) is in the leading place, in the field of prenatal medicine and neonatology are applied preventive treatments in order to reduce the risk of death of preterm newborns [1]. Even if premature birth is not an acute disease, it is one of the major causes of infants’ death and it continues to be significant public health burden. The average cost of medical care for a premature and low birth-weight baby for the first year of life is high in developing country like Ethiopia. These high medical expenses could burden the parents and family. In other case insurance coverage may provide additional benefits in the health plan in order to cover the medical cost incurred during the first year of preemies’ life and hence it helps to reduce the burden of the family. The first 4 weeks of life carries one of the highest risks of death of any 4 weeks period in the human lifespan. Of the 130 million babies born every year, about 4 million die in the first 4 weeks of life in the wide world. Reducing Neonatal morbidity and mortality are now a major focus of child health strategies [16]. The study conducted in Fawzy Moaz Hospital in Egypt, were reported that 48% death of preterm infants admitted to NICU [9]. According to Ethiopian Demographic and Health Survey in 2011, in Ethiopia, high rate of neonatal mortality (37 deaths per 1000 live births)is reported and preterm birth is believe to be a major and direct cause of neonatal mortality [7] . A total of 225 neonates were admitted during the study period of Jan 2012 to Dec 2012 in JUSH. One of the main causes of admissions were found to be prematurity (30.7%) and the hospital neonatal mortality rate was 15.9%.The average length of hospital stay was 9.5 days with SD of 8.2 days using binary logistic regression analysis [25]. Parents’ awareness of premature birth and its consequences is still at a low level. Survival rates pattern or trend over the years is very important to give a sign to the parents whether their premature baby’s life is at a risky level or otherwise [27]. This study was intended to model the survival time of premature death of infants admitted to neonatal intensive care units (NICU) in Jimma University specialized hospital (JUSH) using survival analysis framework. The primary variable in survival analysis is survival time, time to death of infants admitted to NICU for this study. It was investigated the major risk factors of preterm death which will help to guide health professionals and health policy makers to identify indicators for monitoring preterm birth strategy and applying necessary preventive and appropriate measures to decrease preterm birth, and may ultimately it will helps to reduce infant mortality rate. The survival approach used for this, were non-parametric, semi-parametric and parametric survival with (exponential, weibull and log-logistic distribution) and the model best fit the data were selected using Akaikie information criteria (AIC). 2 Study Design and Procedure 2.1 Source of Data All medical records of preterm infants those who were admitted to Neonatal Intensive Care Unit at Jimma University Specialized Hospital from January, 2013 to December, 2015 were retrospectively reviewed by medical professionals. The data consists of infants that admitted to Jimma University Specialized Hospital (Jimma, Ethiopia) with prematurity case. The total numbers of preterm infants admitted to NICU during this study were 552. Among the total of 552 of preterm registered in the given year, only 490 premature infants whose card had full information, satisfy inclusion criteria and hence are included in this study. Thus, Preterm infants admitted to NICU with a gestational age of 26 weeks or greater and less than 37 completed weeks were included and Term infants who born with in gestational age of greater than 37 weeks are excluded. 2.2 Study Population and Period The study population includes all premature infants. All preterm infants those who were admitted to Neonatal Intensive Care Unit (NICU) at Jimma University Specialized Hospital from January, 2013 to December, 2015 G.C. were eligible to be included in the study 2.3 Data Collection Procedure This study has an ethical approval from college of Natural Sciences, Jimma University, Research and Ethical Review board with number Ref. no RPG/132/08. Ethical permission was obtained from the neonatology department of JUSH. The hospital based data collected by trained enumerator and principal investigator. So this study incorporates secondary data and Data collection process was carried out in the time interval of 25/07/2008 to 15/08/2008 E.C. 2.4 Variables of the Study The response (dependent) variable is continuous and describes the length of hospital stay time in days. The response variable for the ith individual is represented by Yi and it measures duration to event and it is defined by status variable (event or censoring variable). Survival time measures the follow-up of time from a defined starting point to the occurrence of a given event. This observation time has two components, the beginning point of the study time and the observation of time to the end. In survival analysis, the outcome of interest (death in this study) is the duration of time until death occurs measured in days. 3 Method of Data Analysis 3.1 Non-parametric Survival Model Kaplan–Meier Estimator The Kaplan–Meier estimator of the survivorship function (or survival probability) S(t) = P(T ≥ t) is defined as: S(t ) = n j − d j where d j is the number of individuals who experience the event at time t j and n j is the number of individuals who have not yet experienced the event at that time [12]. When comparing groups of subjects, it is always a good idea to begin with a graphical display of the data in each group. The figure in general shows the pattern of one survivorship if the function lying above another, it means the group defined by the upper curve lived longer, or had a more favorable survival experience than the group defined by the lower curve. The statistical question is whether the observed difference seen in the figure is significant. The general form of this test statistic is given by: Q = time t(i). In group 0, n1i is the number at risk at observed survival time t(i) in the group 1, d0i is the number of observed deaths in group 0, d1i is the number of observed deaths in group 1, ni is the total number of individuals at risk before time t(i) and di is the total number of deaths at t(i). If weights (wi ) equal to 1, i.e. wi = 1, the Cochran– Mantel–Haenszel Log Rank test (QLR) is used; otherwise if wi = ni , the Generalized Wilcoxon test (QGWt) is used. Q L R = and QGW t = 3.2 Semi-parametric Survival Model 3.2.1 The Cox Proportional Hazard Model The Cox Proportional Hazard (PH) Model is a multiple regression method and is used to evaluate the effect of multiple covariates on the survival [6]. This model gives an expression for the hazard at time t for an individual with a given specification of a set of explanatory variables denoted by X and it is generally given by: h(t, xt , β) = h0(t ) exp(β Xt ) where h0(t ) is the baseline hazard function that characterizes how the hazard function changes as a function of survival time, Xi is the vector of values of the explanatory variables for the i th individual at time t and β is the vector of unknown regression parameters that are assumed to be the same for all individuals in the study, which measure the influence of covariate on the survival experience. The cumulative hazard function is given by:H (t ) = H0(t ) exp(β X ). The survivor function is S(t, X, β) = S0(t ) exp(β X ). Where S0(t ), is the baseline survival function. 3.2.2 Fitting Proportional Hazard Model The maximum likelihood estimates of the cox model parameters are derived by maximizing a likelihood function usually denoted as L. The likelihood function is a mathematical expression which describes the joint probability of obtaining the data actually observed on the subjects in the study as a function of the unknown parameters (the β s) in the model being considered. L is sometimes denoted as L(β) where β denotes the collection of unknown parameters. The partial likelihood can be written as the product of several likelihoods, one for each of, say k failure times. Thus, at the j th failure time, l j denotes the likelihood of failing at this time, given survival up to this time. Note that the set of individuals at risk at the j th failure time is called the “risk set,” R(t( j)), and this set may changes and actually gets smaller in size as the failure time increases. j=1 i=1 In general sense, the partial likelihood is given by the expression where the summation in the denominator is over all subjects in the risk set at time ti denoted by R(ti ), when there are no tied times, and it is often modified to exclude terms when ci = 0, yielding j∈R(t( j)) where the product is over the m distinct ordered survival time and X(i) denotes the value of the covariance for the subject with ordered survival times t(i). 3.3 Parametric Regression Modeling In previous topics it was focused entirely on the use of semi-parametric model of proportional hazards Cox regression model, in the analysis and prediction of the survival time of infants with prematurity. The basis of this method was to avoid having to specify the hazard function completely. However, there may be settings in which the distribution of the survival time is in specific parametric distribution that justifies the use of a fully parametric model to better address the goal of the analysis. For the time data and skewed to the right with exponential distribution, the time of survival for a set of covariate X, which is called, accelerated failure time is expressed as: where, ε , is the error component. The exponential model (t ∼ exp(α)) is the simplest parametric model and assumes a constant risk or hazard over time, which reflects the property of the distribution appropriately called “lack of memory” because the hazard function h(t) = A does not depend on time. Hence the probability of failure in a time interval (t, t + δ t) does not depend on previous interval. The survivorship function may be obtained by expressing in terms of time as: S(t, X, β) = exp(−t e−β X ) and the hazard function of the exponential regression model is h(t, X, β) = e−(β X). For the exponential regression survival model, the hazard ratio for the dichotomous covariate is H R(x = 1, x = 0) = eβ1 . 3.3.2 Weibull Regression Model Survival time t is a positive random variable with weibull probability density function that can be expressed as: lative hazard function becomes: H (t ) = μt α depending on the value of α, the hazard function can increase or decrease with increasing survival time. Hence the weibull model can yield an accelerated failure time model. Independent observations (ti , δi ), i = 1, . . ., n with survival time ti and censoring indicator δi which has value of one if i th observation is not censored and zero when the i th observation is censored and let β be the unknown parameter. The likelihood function is: and the cumui=1 n ⎧ ⎨ i=1 ⎩ exp − Re-parameterizing the weibull distribution using λ = μ−α then h0(t ) = λαt α−1 will be the baseline hazard function. Now incorporate covariates X in the hazard function, the weibull regression models become: h(t, X, β) = λαα−1 exp(Xβ ) −σγ j , α = σ −1 and μ = exp(υ). Multiple covariate log-logistic accelerated failure time may be expressed as: where σ is the scale parameter and ε is residual (unexplained) variation in the transformed survival times [5]. The survivorship function for the model (13) is: S(t, X, β, σ ) = !1 + exp(z)"−1 where z is the standardized log-time outcome variable that is; z = (y−β0σ−β1 X) and y = ln(t). OR(X = 1, X = 0) = = exp This is independent of time. Then the hazard rate is given as follows: λpt p−1 exp (Xβ) h(t, x, β) = 1 + λpt p exp (Xβ) 4 Result and Discussion The medical cards of 490 premature infants were reviewed, 319 (65.1%) were discharged at the end of the follow up and 171 (34.9%) premature infants were died. Out of which 45.3% were female having death proportion of 15.5% lower than males (19.4%). 4.1 Comparison of Survival Experience The mean survival time of premature infants who had Hyaline membrane disease were 12.319 days which is lower than infants who hadn’t hyaline membrane disease ( p < 0.000). Both log-rank and Generalized Wilcoxon test for survival difference were highly significant. Log-rank showed that there is a significant difference of survival experience among groups of prenatal Asphyxia, hyaline membrane disease, sepsis, jaundice, Antenatal care visit, gestational age, respiratory distress syndrome and weight of infant. Similarly the results of Breslow (Generalized Wilcoxon) test also show that there were significant differences among premature infants survival experience of groups prenatal Asphyxia, hyaline membrane disease, sepsis, jaundice, gestational age, respiratory distress syndrome and weight of infant except for Antenatal care visit which means there is no differences in survival experience of premature infants whose mother had Antenatal care visit and had no Antenatal care visit have at the earlier phases where the number at risk is large (Tables 1, 2). Survival of premature infants was significantly related with having sepsis, jaundice, prenatal asphyxia, hyaline membrane disease, respiratory distress syndrome, gestational Age at birth and temperature during admission. The values of the Wald-statistic for individual β coefficients support that the estimated values βi s are significantly Table 1 Summary results of premature infants’ death by different demographic, health and risk behavior variables, at JUSH, during 2013–2015 Weight at birth Table 1 continued Category Death (%) Censored (%) Total (%) Yes 125 (25.2) 101 (20.6) Minimum Maximum Median 0.13 144 48 [36.9, 38] 36.91 0.11 [36.69, 37.13] [15,…] 21.23 0.024 [19.22, 23.24] PNA prenatal asphyxia, HMD hyaline membrane disease, RDS respiratory distress syndrome different from zero at 5% level of significance. Hazard ratio having 95% CI for premature infants admitted to NICU who had Sepsis, Jaundice, Prenatal asphyxia, Hyaline membrane disease and Respiratory distress syndrome compared to those infants who hadn’t Sepsis, Jaundice, prenatal asphyxia,Hyaline membrane disease and Respiratory distress syndrome were 1.7763 (1.2234,2.5790), 2.3950 (1.7071,3.3603), 1.6399 (1.0495,2.5625), 1.8772 (1.2927,2.7261) and 2.2237 (1.5166,3.2605), respectively. That is, the risk of death for premature infants those who had Jaundice and Respiratory distress syndrome were 2.3950 and 2.2237 times higher than those infants without jaundice and Respiratory distress syndrome respectively. Premature infants who had Sepsis, Prenatal asphyxia and Hyaline membrane disease were 77.63, 63.99 and 87.72% more likely to die than those infants without Sepsis, Prenatal asphyxia and Hyaline membrane disease respectively. Also, the hazard ratio (95% CI) for premature infants admitted to NICU those who were born in between gestational age (28–30], (30–32], (32–34] and (34–37) weeks compared to those who were born at gestational age (26–28] were 0.598 (0.3232,1.1097), 0.5186 (0.2766,0.9725), 0.3787 (0.2042,0.7021) and 0.5528 (0.2940,1.0395), respectively. That is, Premature infants who were born in between gestational age (28–30], (30–32], (32–34] and (34–37) weeks were 40.2, 48.14,62.13 and 44.72% less likely to die than to those infants who were born in the interval (26–28] weeks of gestational age respectively. Moreover, by letting other covariates constant, the hazard ratio (95% CI) of temperature at admission for premature infants were 0.8553 (0.7552,0.9687) at NICU. This means, the hazard ratio for a one unit increase in temperature is around 85.53%, so that increasing infants temperature,such that infants temperature goes up by one leads to a reduction in risk of death of 14.47% among premature infants admitted to NICU survivors. The remaining variables which were used in the single covariate analysis such as antenatal care visit of mother, weight of infants and Age at admission found to be non-significant. This implies that the covariates antenatal care visit of mother, weight of infants and Age at admission were no jointly effect on the time to death of premature infants admitted to NICU rather individual effect on the survival of preterm (Table 3). m e ] r 2 ] p p ) 6 8 f u s . 5 o y 4 . o r a I 2 3 e g d , 1 cn re in C .46 ,3 ie v ( % 0 .5 r o e 5 2 6 eexp ison ltim 9 [ [ liav tcun iavv v f r ru ro su 4 e 3 9 s v n n 1 5 s 4 1 f i a a 0 a 4 3 o rv e e .54 . se .7 .2 no su M M 2 10 id 2 1 is f e r o a n ap ty ix ra oCm liauq syhp ebm e a ] ] ] ] ] ] ] ] ) 8 1 5 3 8 4 1 2 ] sy .4 .3 .4 .9 .9 .7 .1 .0 12 a 5 7 3 9 5 8 5 2 . pu dn IC ,27 ,28 ,28 ,19 ,24 ,16 ,28 ,34 ,13 rog i( % .97 .98 .57 .44 .10 .52 .40 .57 .02 r e 5 1 1 1 1 2 1 2 2 9 e m 9 [ [ [ [ [ [ [ [ [ v it o l s a n iv tio rv c u n s 4 5 0 8 0 3 7 4 rfu ean ean .72 .14 .51 .20 .06 .65 .79 em 779 .116 o 1 3 0 7 3 5 o . iv M M 2 2 2 1 2 1 22 rd 9 1 v n r y u s s ) ) ) ) ) ) ) ) ) ) R 90 03 97 25 21 95 25 61 87 05 H 7 6 0 7 0 3 6 2 6 6 fro ,.25 ,.33 ,.11 ,.09 ,.07 ,.10 ,.25 ,.27 ,.09 ,.32 I 4 1 2 6 2 0 5 7 2 6 C 32 70 32 67 40 49 94 29 55 61 % .2 .7 .3 .2 .2 .2 .0 .2 .7 .5 59 (1 (1 (0 (0 (0 (0 (1 (1 (0 (1 4 5745. 8734. 5126. 6565. 9711. 5927. 4946. 6298. 1563. 7992. ,oIC 0 0 0 0 0 0 0 0 0 0 ti − − − − − ra The relationship between covariates and survival probability of premature infants modeled by log logistic regression model are presented in Table 4 below. As can be seen from this table, survival of the infants was significantly related with having sepsis, Jaundice, prenatal Asphyxia, hyaline membrane disease, respiratory distress syndrome, with gestational age at birth and temperature at admission. Using the regression model of Eq. (15) and with the parameters found, the survival time of premature infants admitted to NICU have log logistic distribution. The log-logistic model has two parameter λ is the scale parameter and ρ is the shape parameter, which can be expressed as time ∼Log-logistic (λ, p). The default output for parametric provide maximum likelihood estimates of intercept v, and scale parameter δ, associated with the logistic distribution. The parameters of the underlying log logistic distribution are the functions of these extreme value parameters, λ = exp(−μ/p) = 0.952245 and p = 1/ α = 1.248439 where μ = exp(v) = 0.06081. Then time ∼log-logistic (0.952245,1.248439) have hazard rate of ho(t) = λ1ρ+tλρ−tρ1 = 1+10.1.985828224t05.t214.8244389439 . The log-logistic regression model that predicts the survival of premature infants admitted to neonatal intensive care unit with identical data settings were: 1.18882t 0.248exp (Xβ) h(t, x, β) = 1 + 0.952245t 1.248exp(Xβ) 1.18882t 0.248 h0(t) = 1 + 0.952245t 1.248 4.4 Discussions Results obtained in this study were found to be analogous with literature on the topic. The first factor that affects survival time of premature infant admitted to neonatal care unit is prenatal asphyxia of neonate. As it was indicated both in cox proportional hazard models and log-logistic regression models the hazard rate of premature infants admitted to NICU who had prenatal Asphyxia is about HR 1.6399 and OR 2.479 times higher than premature infants who had no prenatal Asphyxia respectively. This result is in accordance with the studies by Lawn et al. [15]. Other studies like studies of Black et al. [4]. Early onset of disease (Sepsis) of premature infants is a prognostic factor that significantly predicts survival of time to death of premature infants. The multiple covariate analysis of cox proportional hazard showed that, hazard rate of having Sepsis is 77.63% higher than premature infants who had no Sepsis and log logistic parametric Table 4 Parameter estimates, standard errors and the hazard ratios in the final log logistic regression model at JUSH, during 2013–2015 95% CI for OR (1.239,4.9585) (1.5971,4.3517) (1.7182,4.361) (1.2415,3.4587) (0.2078,1.2543) (0.1376,0.8217) (0.0993,0.5886) (0.1794,1.0549) (0.6826,0.9642) (2.033,5.3148) (1.1324,1.5374) (2.012,2.4571) model also indicates that the proportional odds of having sepsis is much higher for infants those had Sepsis (OR 2.072). The result is comparable with earlier study [4, 18] at global level. Sime et al. [25]. The results of this study suggested that hyaline membrane disease was significant predictive factor for time-to-death of premature infants. Infants who had hyaline membrane disease have hazard ratio 1.8772 (1.2927,2.7261), that is the risk of death for premature infants who had hyaline membrane disease was 87.72% higher than infants who hadn’t hyaline membrane disease. The log-logistic parametric regression model also suggested that the proportional odds of survival were 2.636 times higher odds of dying than premature infants who had no hyaline membrane disease (or 37.94% less likely of risk of dying for premature infants who had no hyaline membrane disease). This result supported by study carried out in ST. George, Germany on prevention of hyaline membrane disease in preterm infant conclude that Hyaline membrane disease is the leading cause of death for infants [1]. Jaundice disease is an important predictor of time to death of premature infants. Both cox proportional hazard and log logistic parametric models suggest that the jaundice disease have significantly associated with premature infants’ death. Premature infants who had jaundice have hazard rate of 2.395 (CI 1.7071,3.3603), higher than Premature infants who hadn’t jaundice disease. This is similar to reports by Khan et al. [14] in Karachi, Pakistan who reported jaundice and sepsis as the commonest morbidities in their preterm patients. However, Onwuanaku et al. [22] in Jos University Teaching Hospital Nigeria, reported sepsis as the commonest morbidity, followed by jaundice. They recommended that is an urgent need for the prevention and adequate management of jaundice in this vulnerable group. In our study in both models, cox proportional and log-logistic regression models, the initial body temperature of infants was significantly associated with reducing the survival probability of premature infants. Cox proportional hazard and loglogistic parametric regression model were suggested the hazard ratio and odds ratio of 0.8553 and 0.811 respectively. This indicates that one unit increase in initial temperature of premature infants reduces the risk of death by 14.47 and 18.9% respectively. This result confirmed with study carried out by Merertu et al. [19] reported that Temperature at admission has strong association with survival in a multivariate model adjusted for gestational age and birth weight of neonatal and With the neonates whose temperature at admission was less than or equal to 33 ◦C have a 5.43 times higher risk of death compared to those with T◦ at admission between 36.50–37.5 ◦C. Respiratory distress syndrome is one of the potential risk factors that affect the survival time of premature infants. Results of cox proportional hazard and log-logistic parametric regression model showed that the respiratory distress syndrome disease have significantly associated with survival time of premature infants admitted to NICU. This study confirmed with the study carried out by Behnaz [2]; concluded that the death rate of respiratory distress syndrome (73.8%) is higher in premature infants who had the disease, which predict occurrences of death in premature infants. 5 Conclusions and Recommendations 5.1 Conclusions Cox proportional hazard analysis showed that the major factors that affect the survival time of preterm are prenatal Asphyxia, Sepsis, hyaline membrane disease, jaundice, Gestational age, temperature and respiratory distress syndrome for time to death of premature infants. Preterm infants having prenatal asphyxia, Sepsis, Jaundice, Hyaline membrane and respiratory distress syndrome have higher death rate. Similarly, preterm with poor health indicators like lower gestational age (26–28) weeks and initial temperature, were less likely to survive. To predict and model the survival time of premature infants, various baseline parametric regression models were applied. Among these, parametric model with loglogistic baseline distribution is best fitted to predict the survival time of the premature infants. 5.2 Recommendations According to the results of this study the main predictive factors for the survival time of premature infants are more of clinical variables. So, health workers should be cautious when mother’s born preterm which has prenatal asphyxia, Sepsis, Jaundice, hyaline membrane disease and respiratory distress syndrome. The log-logistic regression model provides better predictions to the survival probability of premature infants’. So, future researchers could make use of this model. Future studies also need to assess the level of awareness, treatment and control of these risk factors. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Appendix 1: Checking for the Linearity of Continuous Covariates See Fig. 1. plot of martingale residual for covariate tempereture Temperature Fig. 2 Plots of score residuals for temperature of infants at admission to detect the existence of influential observation in Cox proportional hazards model Temperature Fig. 1 Plots of martingale residuals against the covariate Temperature Appendix 2: Identification of Influential and Poorly Fit Subjects See Fig. 2. See Fig. 3. Appendix 3: Cox–Snell Residual Plots for Model Assessment Exponential Weibull Cox-Snell residual Cox-Snell residual Hazard Cox-Snell residual Cox-Snell residual Fig. 3 Plots of parametric survival models to examine models that fit the data better Appendix 4: Residual Plot for Model Assessment See Fig. 4. Test PH assumption by scheonfeld Test PH assumption by scheonfeld bandwidth = .8 Time in days bandwidth = .8 Time in days Test PH assumption by scheonfeld Test PH assumption by scheonfeld bandwidth = .8 Time in days bandwidth = .8 Time in days Test PH assumption by scheonfeld bandwidth = .8 Time in days Fig. 4 Plots of scaled Schoenfeld residuals against length of hospital stay for each covariate in Cox proportional hazards model fit of premature infants 1. Anna M , Ekaterina U , Nikoleta P , Mariana P ( 2016 ) Prevention of hyaline membrane disease (HMD) in preterm infants . St. George “University Hospital”, Plovdiv 2. Behnaz B , Farzaneh E , Maryam Sh ( 2015 ) Neonatal Mortality and its Main Determinants in Premature Infants Hospitalized in NICU in Fatemieh Hospital , Hamadan, Iran, Iranian society of pediatrics . J Compr Ped 6 ( 3 ): e26965 3. Belaynew W , Teumay A , Getachew G , Mohamed K ( 2015 ) Effects of inter pregnancy interval on preterm birth and associated factors among postpartum mothers who gave birth at Felege Hiwot referral hospital . World J Pharm Pharm Sci 4 ( 4 ): 12 - 25 4. Black RE , Cousens S , Johnson HL et al ( 2010 ) Global, regional, and national causes of child mortality in 2008: a systematic analysis . Lancet 375 : 1969 - 1987 5. Collett D ( 2003 ) Modelling survival data in medical research, 2nd edn . Chapman and Hall/CRC, London 6. Cox D ( 1972 ) Regression models and life-table . J R Stat Soc B 34 : 187 - 220 7. EDHS ( 2011 ) Central Statistical Agency , ( 2011 ). Ethiopian Demographic and Health Survey , Addis Ababa 8. Engle WA , Kominiarek MA ( 2008 ) Late preterm infants, early term infants, and timing of elective deliveries . Clin Perinatol 35 ( 2 ): 325 - 341 9. Fakher M , Shaaban W , Monein AA , Hassan Z , Fikry MM ( 2005 ) Statistical study of preterm infants admitted to NICU In Fawzy Moaz . Alex J Pediatr 19 : 155 - 158 10. Federal Ministry of Health (FMOH) (2011) Saving newborn lives in Nigeria. Newborn Health I the Context of the Integrated Maternal, Newborn and Child Health Strategy . Revised 2nd Edition 11. Goldenberg RL , Culhane JF , Iams JD , Romero R ( 2008 ) Epidemiology and causes of preterm birth . Lancet 371 : 75 - 84 12. Hosmer DJ , Lemeshow S ( 1999 ) Applied survival analysis: regression modeling of time to event data, 2nd edn . Wiley 13. Howson CP , Kinney MV , Mcdougall L ( 2013 ) Born too soon: preterm birth matters . Reprod Health J 10 : 1 - 9 14. Khan MR , Maheshwari PK , Shamim H , Ahmed S , Ali SR ( 2012 ) Morbidity pattern of sick hospitalized preterm infants in Karachi, Pakistan . J Pak Med Assoc 62 : 386 - 388 15. Lawn JE , Cousens S , Zupan J ( 2005 ) 4 million neonatal deaths: when? Where? Why? Lancet 365 : 891 - 900 16. Lawn JE , Wilczynska-Ketende K , Cousens SN ( 2006 ) Estimating the causes of 4 million neonatal deaths in the year 2000 . Int J Epidemiol 5 ( 3 ): 706 - 718 17. Lawn JE , Gravett MG , Nunes T , Rubens E , Stanton C et al ( 2010 ) Global report on preterm birth and stillbirth (1 of 7): definitions, description of the burden and opportunities to improve data . BMC Pregnancy Child 10 : 1 - 22 18. Liu L , Johnson HL , Cousens S , Perin J , Scott S et al ( 2012 ) Global, regional, and national causes of child mortality in 2000-2010: an updated systematic analysis . Lancet 379 : 2151 - 2161 19. Merertu T , Bogale W , Yonas R , Amha M ( 2011 ) Survival Of preterm infants admitted to Tikur Anbessa Hospital Nicu , Addis Ababa 20. Mohammad T , Khattak AA , Rehman Shafiq-Ur ( 2009 ) Mortality and morbidity pattern in small-forgestational age and appropriate-for-gestational age very preterm babies: a hospital based study . J Ayub Med Coll Abbottabad 21 : 16 - 21 21. Offiah I , O'Donoghue K , Kenny L ( 2012 ) Clinical risk factors for preterm birth , vol 3. INTECH Open Access Publisher , pp 74 - 91 22. Onwuanaku CA , Okolo SN , Ige KO , Okpe SE , Toma BO ( 2011 ) The effects of birth weight and gender on neonatal mortality in North Central Nigeria . BMC Res Notes 4:562 23. Pennell CE , Jacobsson B , Williams SM , Buus RM , Muglia LJ , Dolan SM ( 2007 ) Genetic epidemiologic studies of preterm birth: guidelines for research . Am J Obstet Gynecol 196 : 107 - 118 24. Risso SDP , Nascimento LFC ( 2010 ) Risk factors for neonatal death in neonatal intensive care unit according to survival analysis . Rev Bras Ter Intensiva 22 : 9 - 26 25. Sime H , Workneh N , Girma E ( 2013 ) Morbidity and mortality of neonates admitted in Jimma University Specialized Hospital Paediatrics Neonatal Ward: a one year retrospective analysis . Ethiopian J Pediatr Child Health, X(10) 26. UNICEF ( 2012 ) Maternal, newborn and child survival, country profile . Ethiopia, Statistics and Monitoring Section/Policy and Practice 27. World Health organization (WHO) (2012) Born too soon: the global action report on preterm birth . In: Howson, Kinney M , Lawn E (eds). Geneva


This is a preview of a remote PDF: http://link.springer.com/content/pdf/10.1007%2Fs40745-017-0107-2.pdf

Modeling Determinants of Time-To-Death in Premature Infants Admitted to Neonatal Intensive Care Unit in Jimma University Specialized Hospital, Annals of Data Science, 2017, DOI: 10.1007/s40745-017-0107-2