Study design
This study used data from a health database in Korea, the National Emergency Department Information System (NEDIS), between 2017 and 2019. The NEDIS is a nationwide ED-based database for evaluating the emergency care system in Korea, established in accordance with Article 15 of the Emergency Medical Service Act. To achieve this goal, the NEDIS collects ED visit-level data, including demographic, clinical, and administrative information. Each visit-level datum also has the corresponding hospital identifier and hospital characteristics, such as total staffed beds, level of ED, and region. All patient-related information was anonymised and electronically submitted to the central processing facility, which was examined both manually and using computerised algorithms to detect data inconsistencies. Between 2017 and 2019, the participation rate of nationwide EDs in the NEDIS was 99.3% (413/416) in 2017, 99.5% (399/401) in 2018, and 99.8% (401/402) in 2019. The design and variables of the NEDIS database have been described elsewhere [33-35].
From the NEDIS database, we identified all patients admitted to the ICU directly from the ED between 1 January 2017 and 31 December 2019 based on the date of presentation to the ED. This operational definition of critically ill patients was adopted from previous studies [4, 6, 29]. Patients with missing age or sex information, those <18 years old, and those with missing days and times of ED presentation and departures were excluded.
Outcomes and variables
The primary outcome of this study was prolonged EDLOS, which was defined as an EDLOS of 6 h or more. This decision was based on existing evidence suggesting that an EDLOS of 6 h or more is associated with increased mortality risk and influences the quality of care in critically ill patients in the ED [23, 36, 37]. The secondary outcome was in-hospital mortality.
We identified patient and hospital variables a priori as potential predictors of prolonged EDLOS and in-hospital mortality risk in critically ill patients. Potential predictors were selected based on a review of the academic literature and data available in the NEDIS database [23-28, 38-41].
Patient variables included age, sex, insurance type, injury code, emergency ambulance attendance, transferred-in, date and time of ED presentation, initial triage score, mechanical ventilation in the ED, diagnosis codes during hospitalisation, Charlson comorbidity index (CCI), and discharge status. The initial triage was scored according to the Korean Triage and Acuity Scale (KTAS), which prioritizes patients according to the five ordinal scales reflecting clinical severity and acuity as follows: resuscitation, 1; emergent, 2; urgent, 3; less urgent, 4; non-urgent, 5 [42]. The date and time of ED presentation were categorised according to the year, season (spring, March–May; summer, June–August; fall, September–November; winter, December–February) and ED shift time (day, 07:00–14:59; evening, 15:00–22:59; night, 23:00–06:59). Diagnostic codes used during hospitalisation were identified based on codes defined in the International Classification of Diseases, Tenth Revision (ICD-10). The CCI score was calculated based on diagnostic codes used during hospitalisation by applying the methods proposed in previous studies, which showed good-to-excellent discriminant power in predicting in-hospital mortality risk [43, 44].
Hospital variables were hospital staffed beds (1,000 or more, 800–999, 600–799, 300–599, and < 300), type of ED (levels 1, 2, and 3), and location (metropolitan city versus provincial area).
Statistical analysis
We calculated the proportion of the study population with overall ED presentations and overall adult ED presentations, as well as the annual incidence/100,000 adult ED presentations.
Descriptive analyses were performed to compare the patient and hospital characteristics between critically ill patients with an EDLOS of 6 h or more and critically ill patients with an EDLOS of <6 h. Categorical variables were reported as frequencies and proportions and were compared between patient groups using Pearson’s Chi-squared test. Continuous variables were described as the median and interquartile range (IQR) and were tested using the Wilcoxon rank-sum test. The median EDLOS with IQR and percentages of in-hospital mortality for each of the most common primary diagnoses were calculated.
We constructed multivariate logistic regressions to model the effects of patient and hospital variables as predictors of prolonged EDLOS (both 6 h an 12 h) and in-hospital mortality risk. To evaluate the potential differential associations of hospital characteristics with prolonged EDLOS vs. in-hospital mortality, we performed a stratified analysis with the highest hospital staffed bed category (1,000 or more) or type of ED (level 1) as the reference in the same logistic regression model.
All analyses were performed using SAS version 9.4 (SAS Inc., Cary, NC, USA) and R version 4.1.3 (R Development Core Team, https://cran.r-project.org/). All tests were two-tailed, and a p-value <0.05 was considered statistically significant.