Study design and selection of study patients
This study was a retrospective analysis of a prospective cohort from 23 EDs, namely the Injury Surveillance Cohort, which was generated by the Korea Center for Disease Control and Prevention (KCDC) from 2011 to 2018. This registry comprised prospectively collected data on the epidemiology and outcome variables of injury patients who presented at an ED [8]. The registry included poisoning cases as a type of injury. We selected poisoning patients from this cohort. This selected registry included the baseline characteristics of poisoning patients: age; sex; time-related factors, such as ED presenting time and poison exposure time; poisoning-related variables, such as the intent of poisoning, route of exposure, type of substance (7 categories and 44 types of substance); and initial vital signs at ED presentation, such as systolic blood pressure (SBP), heart rate (HR), respiration rate (RR), body temperature (BT), and AVPU scale of mental status. The registry also contained outcome-related variables, such as mortality at the ED or after hospital admission.
Patients who were transferred from the initial ED to another hospital; those who had incomplete data of poisoning-related variables, initial physiological condition-related variables, or outcome-related variables; and those who died on arrival at the ED were excluded from this study (Figure 1).
The selected study population was divided into two groups, namely the derivation group for the prediction of in-hospital mortality and the validation group for the external validation of the developed prediction model (Figure 1).
The Institutional Review Board of the Korea University hospital approved this study (IRB No. 2020AN0195).
Data analysis
The primary outcome was in-hospital mortality. We compared the characteristics of the poisoning patients between the derivation and validation groups (Table 1). Age, sex, time from exposure to ED presentation, classes of substance, intent of poisoning, route of exposure, vital signs of the patient at ED presentation, and in-hospital mortality were analyzed (Table 1). For analysis, the variables related to poisoning characteristics were categorized as follows: intent of poisoning: 1) unintentional, 2) intentional, and 3) unknown; route of exposure: 1) dermal, ocular, or injection; 2) oral; and 3) inhalation; and toxic substance included 44 kinds of substances that were classified into eight categories from A to H. For categorization of substances, we considered the classification in the types of substances. And then we categorized the substances in the same classification according to the mortality index (MI) of each substance: A) pharmaceutical agents with MI of less than 0.5%, B) pharmaceutics with MI 0.5 – 5%, C) artificial toxic substances with MI of less than 1.0%, D) artificial toxic substances or pesticides with MI of 1.0 – 10.0%, respectively, E) artificial toxic substances or pesticides with MI of 11.0 – 20.0%, respectively, F) paraquat with MI of 52.5%, G) gases with MI of less than 1.0%, H) natural toxic substances with MI of less than 1.0% (Table 2)(An additional file 1 shows this in more detail [see Additional file 1]). The patient’s physiological variables included age, SBP, HR, RR, BT, and mental status (AVPU scale), in accordance with the predictors in SAPS-II [9]. However, because SAPS-II does not include RR score, we categorized RR according to the normal range (12–24 breaths/min) [10, 11].
Development of the new poisoning mortality scoring system
We developed the new-Poisoning Mortality Scoring system (new-PMS) to generate a prediction model for the derivation group (2011–2017 data of the KCDC cohort) (Figure 1). First, we compared poisoning- and physiological condition-related variables between the patients who survived and were discharged (survivor subgroup) and those who died at the hospital (in-hospital death subgroup) among the derivation group (Table 3). We selected variables that had statistical and clinical significance in acute poisoning as predictors for developing the new-PMS [12]. Points for each category of each predictor were computed using multivariable logistic regression, in which the regression coefficient for each category of each predictor was converted into points by dividing the smallest regression coefficient in the model (Table 4) [13]. The sum of these points of each category in each predictor was the new-PMS.
Performance evaluation of the new-PMS
We analyzed the performance of the new-PMS in predicting mortality using receiver operating characteristic (ROC) curves in both the derivation and validation groups. The validation group (2018 data of the KCDC cohort) was subjected to external temporal validation.
For simple interpretation in a clinical setting, we created the following risk groups: very low, low, intermediate, and high risk according to the quartile range of the new-PMS. Real mortalities were investigated in the derivation and validation groups, respectively [12]. Moreover, we generated a mortality curve according to the new-PMS in the derivation group.
Statistical analysis
Continuous variables were reported as the median with interquartile ranges (IQR). Differences in the medians were compared using the Mann-Whiney U-test. Categorical variables were compared using the chi-square test. All statistical analyses were performed using SPSS version 20.0 (IBMSPSS Inc., Chicago IL, USA). Two-tailed p values < 0.05 were considered statistically significant.