Jurnal #2 IKA12
Jurnal #2 IKA12
Jurnal #2 IKA12
ORIGINAL ARTICLE
to the intensive care unit
Department of Pediatrics, Guangdong Second Provincial General Hospital, No. 466 Middle Xingang Road, Zhuhai District, Guangzhou,
Guangdong 510317, P.R. China
∗ Corresponding author: Tel: +86 020–61322569; E-mail: [email protected]
Received 31 July 2020; revised 30 November 2020; editorial decision 21 February 2021; accepted 1 March 2021
Background: To explore the influencing factors for in-hospital mortality in the neonatal intensive care unit (NICU)
and to establish a predictive nomogram.
Methods: Neonatal data were extracted from the Medical Information Mart for Intensive Care III (MIMIC-III)
database. Both univariate and multivariate logit binomial general linear models were used to analyse the factors
influencing neonatal death. The area under the receiver operating characteristics (ROC) curve was used to assess
the predictive model, which was visualized by a nomogram.
Results: A total of 1258 neonates from the NICU in the MIMIC-III database were eligible for the study, including
1194 surviving patients and 64 deaths. Multivariate analysis showed that red cell distribution width (RDW) (odds
ratio [OR] 0.813, p=0.003) and total bilirubin (TBIL; OR 0.644, p<0.001) had protective effects on neonatal in-
hospital death, while lymphocytes (OR 1.205, p=0.025), arterial partial pressure of carbon dioxide (PaCO2 ; OR
1.294, p=0.016) and sequential organ failure assessment (SOFA) score (OR 1.483, p<0.001) were its independent
risk factors. Based on this, the area under the curve of this predictive model was up to 0.865 (95% confidence
interval 0.813 to 0.917), which was also confirmed by a nomogram.
Conclusions: The nomogram constructed suggests that RDW, TBIL, lymphocytes, PaCO2 and SOFA score are all
significant predictors for in-hospital mortality in the NICU.
© The Author(s) 2021. Published by Oxford University Press on behalf of Royal Society of Tropical Medicine and Hygiene. This is an Open Access
article distributed under the terms of the Creative Commons Attribution-NonCommercial License (http://creativecommons.org/licenses/by-
nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For
commercial re-use, please contact [email protected]
1
X. Huang et al.
Variables Survival group (n=1194) Death group (n=64) Total (N=1258) Statistic p-Value
Values are presented as median (IQR; quartile 1–quartile 3) unless stated otherwise.
Methods red blood cell [RBC] count, platelet [PLT] count, red cell distribu-
tion width [RDW], neutrophils, lymphocytes, sodium, total biliru-
Patients bin [TBIL], fraction of inspired oxygen [FiO2 ], arterial partial pres-
The MIMIC-III database was used to extract data. The MIMIC-III sure of carbon dioxide [PaCO2 ], arterial partial pressure of oxygen
is a large, free database that contains de-identified health data [PaO2 ]), sequential organ failure assessment (SOFA) score and the
for >40 000 ICU patients from the Beth Israel Deaconess Med- Simplified Acute Physiology Score II (SAPS II). These data could be
ical Center (Boston, MA, USA) between 2001 and 2012.13 It not downloaded from several sources, including archives from crit-
only contains demographics, vital signs, medications and labo- ical care information systems, hospital electronic health record
ratory tests, but also includes observations and notes charted databases and Social Security Administration Death Master File.14
by care providers, procedures and diagnostic codes, fluid bal- Neonatal deaths in this study were defined as death occur-
ance, imaging reports, length of stay and survival data. According ring within 30 d following admission to the NICU. The SOFA score
to the MIMIC-III, 7874 neonates were admitted to the NICU.14 was developed to assess the acute morbidity of critical illness
Neonates with a gestational age of 23–31 weeks were included at the patient level based on six criteria—cardiovascular, respira-
in our study, while those with missing data, including birthweight, tory, nervous, hepatic, renal and haematological systems—each
vital signs and laboratory measurements, were excluded. scored on a scale of 0–4. The higher the SOFA score, the worse the
Since the data used in this study were accessed from MIMIC- organ dysfunction.15 The SAPS II provides an estimate of the risk
III, an openly available dataset, the study did not get approval of death without specifying a primary diagnosis. It contains 17
from the Institutional Review Board of Guangdong Second Provin- variables: age, type of admission, 12 physiology variables and 3
cial General Hospital. underlying disease variables.16 PaO2 and FiO2 are commonly used
for diagnosis of lung injury and TBIL reflects liver metabolism.17,18
Data collection
Data from the MIMIC-III used for analysis in this study included Statistical analysis
gender, birthweight, heart rate, laboratory parameters (bicarbon- Categorical data were presented as number and percentage,
ate, haematocrit, haemoglobin, white blood cell [WBC] count, while the χ 2 test or Fisher’s exact test was used for compar-
2
International Health
Gestational age (weeks) 0.971 (0.853 to 1.104) 0.649 0.968 (0.888 to 1.056) 0.464
ison. Continuous variables were expressed as medians and 1345 patients, laboratory parameters 3992 patients and scor-
interquartile ranges (IQRs). Differences in baseline information ing systems 584 patients), 1258 neonates (733 males [58.27%]
of included neonates and clinical characteristics, as well as lab- and 525 females [41.73%]) were eligible for the study. In-hospital
oratory parameters between groups, were compared according death occurred in 64 patients. All included subjects were included
to in-hospital survival outcomes for neonates in the NICU. With in the death group (n=64) and 1194 were included in the survival
neonatal in-hospital death or survival as the dependent variable group. The baseline characteristics of the two groups are com-
and demographic characteristics and laboratory indicators as pared in Table 1. It can be observed that except for gender, heart
independent variables, a univariate logit binomial general linear rate, PaCO2 , PaO2 and bicarbonate, the differences were all pro-
model (GLM) was performed. Independent variables with sig- nounced between the death group and the survival group regard-
nificant difference in the univariate analysis were included in a ing other variables, including birthweight (p<0.001), WBC count
multivariate logit binomial GLM. In view of the multicollinearity (p=0.030), RBC count (p<0.001), PLT count (p<0.001), haema-
among independent variables, a stepwise regression analysis tocrit (p=0.001), haemoglobin (p=0.002), RDW (p=0.004), neu-
was conducted to identify the risk factors and establish a pre- trophils (p=0.029), lymphocytes (p=0.026), TBIL (p<0.001), FiO2
dictive model. The receiver operating characteristics (ROC) curve (p<0.001), sodium (p=0.008), SOFA score (p<0.001) and SAPS II
and the area under the curve (AUC) were used to assess the (p<0.018).
predictive model, which was visualized by a nomogram. The
nomogram is one method used to visualize the predictive model
with graphic scoring.19 All statistical tests were two-sided and
p-values ≤0.05 were considered to be statistically significant. All Univariate and multivariate analysis of in-hospital
statistical analyses were carried out using R statistical software death in neonates
(R Foundation for Statistical Computing, Vienna, Austria).
Univariate and multivariate analysis of in-hospital death in
neonates is presented in Table 2. As it shows, birthweight (odds
ratio [OR] 0.335, p<0.001), RBC count (OR 0.459, p<0.001), PLT
Results count (OR 0.995, p=0.001), haematocrit (OR 0.957, p=0.006),
haemoglobin (OR 0.880, p=0.007), RDW (OR 0.845, p=0.037),
Baseline characteristics of included neonates neutrophils (OR 0.982, p=0.023), TBIL (OR 0.524, p<0.001) and
Of 7874 neonates admitted to the NICU, there were 7800 with bicarbonate (OR 0.892, p=0.003) were all associated with a re-
a gestational age of 23–31 weeks. After excluding 6542 patients duced risk of in-hospital death in neonates, whereas lymphocytes
with missing data (birthweight 621 patients, survival outcomes (OR 1.180, p=0.024), FiO2 (OR 1.011, p=0.004), PaCO2 (OR 1.262,
3
X. Huang et al.
p=0.039), sodium (OR 1.089, p<0.001) and SOFA score (OR 1.533,
p<0.001) related to an increased risk.
Through multivariate analysis, RDW (OR 0.813, p=0.003) and
TBIL (OR 0.644, p<0.001) were found to have protective ef-
fects on neonatal in-hospital death, while the risk of in-hospital
mortality was increased 1.205-fold and 1.294-fold for each in-
crease in lymphocytes (OR 1.205, p=0.025) and PaCO2 (OR 1.294,
p=0.016), respectively. In addition, the SOFA score (OR 1.483,
p<0.001) was also identified as an independent risk factor for
neonatal in-hospital mortality.
4
International Health
5
X. Huang et al.
predictor in oncology children mechanically ventilated >3 d, and a systematic analysis for the Global Burden of Disease Study 2013.
it was strongly associated with mortality. Lancet. 2014;384(9947):957–79.
As a prognostic device, the nomogram has been extensively 2 Orsido TT, Asseffa NA, Berheto TM. Predictors of neonatal mortality in
used in oncology and medicine to predict the probability of clini- neonatal intensive care unit at referral hospital in southern Ethiopia:
cal events. It can help clinicians make quick and sound decisions a retrospective cohort study. BMC Pregnancy Childbirth. 2019;19(1):
due to the characteristics of better accuracy, straightforward 83.
digital interfaces and readily comprehensible prognosis.30,31 To 3 Liu L, Johnson HL, Cousens S, et al. Global, regional, and na-
Ethical approval: Not required. 18 Woreta TA, Alqahtani SA. Evaluation of abnormal liver tests. Med Clin
North Am. 2014;98(1):1–16.
Data availability: The data that support the findings of this study are 19 Tunthanathip T, Udomwitthayaphiban S. Development and validation
available from the first author and corresponding author upon reasonable of a nomogram for predicting the mortality after penetrating trau-
request. matic brain injury. Bull Emerg Trauma. 2019;7(4):347–54.
20 Martin SL, Desai S, Nanavati R, et al. Red cell distribution width and its
association with mortality in neonatal sepsis. J Matern Fetal Neonatal
Med. 2019;32(12):1925–30.
21 Ellahony DM, El-Mekkawy MS, Farag MM. A study of red cell distribution
References width in neonatal sepsis. Pediatr Emerg Care. 2020;36(8):378–83.
1 Wang H, Liddell CA, Coates MM, et al. Global, regional, and national 22 Lawn JE, Cousens S, Zupan J, et al. 4 million neonatal deaths: when?
levels of neonatal, infant, and under-5 mortality during 1990–2013: where? why? Lancet. 2005;365(9462):891–900.
6
International Health
23 Said AS, Spinella PC, Hartman ME, et al. RBC distribution width: with minimal change nephropathy. Clin Ther. 2013;35(3):286–
biomarker for red cell dysfunction and critical illness outcome? Pe- 93.
diatr Crit Care Med. 2017;18(2):134–42. 28 Aperstein Y, Cohen L, Bendavid I, et al. Improved ICU mortality pre-
24 Thompson BL, Wyckoff SL, Haverstick DM, et al. Simple, reagentless diction based on SOFA scores and gastrointestinal parameters. PLoS
quantification of total bilirubin in blood via microfluidic phototreat- One. 2019;14(9):e0222599.
ment and image analysis. Anal Chem. 2017;89(5):3228–34. 29 Ha EJ, Kim S, Jin HS, et al. Early changes in SOFA score as a prognostic
25 Dani C, Poggi C, Fancelli C, et al. Changes in bilirubin in infants with factor in pediatric oncology patients requiring mechanical ventilatory