Abstract
IMPORTANCE: Decreasing mortality in the field of pediatric critical care medicine has shifted practicing clinicians' attention to preserving patients' neurodevelopmental potential as a main objective. Earlier identification of critically ill children at risk for incurring neurologic morbidity would facilitate heightened surveillance that could lead to timelier clinical detection, earlier interventions, and preserved neurodevelopmental trajectory. OBJECTIVES: To develop machine-learning models for identifying acquired neurologic morbidity in hospitalized pediatric patients with critical illness and assess correlation with contemporary serum-based, brain injury-derived biomarkers. DESIGN, SETTING, AND PARTICIPANTS: This prognostic study used data from all children admitted to a quaternary pediatric intensive care unit in a large, freestanding children's hospital in Western Pennsylvania between January 1, 2010, and December 31, 2022. External model validation used data from children admitted between January 1, 2018, and December 31, 2023, to a quaternary pediatric intensive care unit in a large, freestanding children's hospital that serves as a referral center for the 5-state region of Washington, Wyoming, Alaska, Montana, and Idaho. EXPOSURES: Critical illness. MAIN OUTCOMES AND MEASURES: The outcome was neurologic morbidity, defined according to a computable, composite definition at the development site or an order for neurocritical care consultation at the validation site. Models were developed using varying time windows for temporal feature engineering and varying censored time horizons between the last feature and the identified neurologic morbidity. A generalizable model created at the development site was optimized and assessed at an external validation site. Correlation was assessed between development site model predictions and measurements of brain biomarkers from a convenience cohort. RESULTS: After exclusions, there were 18 568 encounters from 2010 to 2022 in the development site generalizable model cohort (median age, 70 [IQR, 18-161] months; 8325 [45%] female). There were 6825 encounters from 2018 to 2021 at the external validation site (median age, 96 [IQR 18-171] months; 3159 [46%] female). A generalizable extreme gradient boosted model with a 24-hour time horizon and 48-hour feature engineering window demonstrated an F1 score of 0.37 (95% CI, 0.33-0.40), area under the receiver operating characteristics curve of 0.81 (95% CI, 0.78-0.83), and number needed to alert of 4 at the validation site. After recalibration at the validation site, the Brier score was 0.04. Serum levels of the brain injury biomarker glial fibrillary acidic protein significantly correlated with model output (rs = 0.34; P = .007). CONCLUSIONS AND RELEVANCE: This prognostic study of prediction models for detecting neurologic morbidity in critically ill children demonstrated a well-performing ensemble of models with biomolecular corroboration. Prospective assessment and refinement of biomarker-coupled risk models in pediatric critical illness are warranted.