This retrospective observational case-control study was conducted at an 858-bed tertiary care teaching hospital in southern Thailand from January 2015 to December 2020. All procedures adhered to relevant laws and institutional guidelines. The study protocol was approved by the Human Research Ethics Unit, Faculty of Medicine, Prince of Songkla University (approval number 63-560-8-1), which granted a waiver for informed consent due to its retrospective design. The privacy rights of all participants were rigorously observed. As this study was retrospective and observational in nature and did not involve any prospective intervention, clinical trial registration was not required.
Patient selection
The inclusion comprised patients undergoing general anaesthesia with endotracheal intubation for elective or emergency surgery who encountered difficult airway management, defined as requiring three or more intubation attempts. Patients who were already on preoperative mechanical ventilation or had anticipated difficult airways were excluded.
Matching procedure
A rigorous matching algorithm was implemented to minimize selection bias and control for potential confounding variables. Cases of difficult airway management were matched with controls who underwent routine airway management based on age (± 5 years), surgery type, and year, in a 1:3 ratio. For each difficult airway case, three controls were selected: one with a preoperatively anticipated difficult airway and two with normal airway assessment.
Definition of unanticipated difficult intubation
Unanticipated difficult intubation was defined as three or more attempts at tracheal tube intubation using direct or video laryngoscopy in a patient not previously identified as having a difficult airway. Intubations were performed by anaesthesia staff with at least one year of experience, including anaesthesiology residents, nurse anaesthetists, and attending anaesthesiologists.
Potential risk factors and confounding variables
Potential predictors included patient-, surgery-, and anaesthesia-related factors such as sex, age, weight, height, body mass index (BMI), American Society of Anesthesiologists (ASA) physical status, preoperative airway assessment, history of difficult intubation and ventilation, comorbidities (e.g., medical conditions, congenital heart disease, obstructive sleep apnoea, head and neck radiation, anatomical abnormalities due to infection, trauma, tumours, or burns, and abnormal laboratory values), type of surgery, laryngoscopic view grades, intubation attempts, intubation devices, incidence of emergency tracheostomy, first intubation attempt by attending staff, and intubation experience.
Sample size determination
A sample size calculation estimated that 70 unanticipated difficult airway cases and 280 controls were required to detect an odds ratio of 2.5, assuming a 15% prevalence of exposure among controls, with 80% power and a 0.05 significance level. Based on the institutional annual incidence of 10–15 cases, a 6-year study period was required.
Statistical analysis
Statistical analyses were performed using R version 4.3.1 (R Foundation, Vienna, Austria). Descriptive statistics are reported as medians with interquartile ranges for continuous variables and as frequencies with percentages for categorical variables. Associations between categorical variables and difficult intubation were assessed using the chi-squared or Fisher’s exact test, while continuous variables were evaluated with Student’s t-tests or Mann–Whitney U tests, depending on data distribution. Collinearity diagnostics and bivariate correlation matrices were evaluated for all variables. In cases of multicollinearity, only one variable was retained for multivariate analysis. Variables with a p-value < 0.25 in the univariate analysis were considered for inclusion in the initial multivariate logistic regression model. The final regression model was derived using backward selection, retaining all significant variables. The optimal cut off point was identified using Youden’s index. Statistical significance was defined as a p-value < 0.05.
Score derivation and validation
The risk prediction score system was developed using predictors derived from the multivariate logistic regression model. Risk scores were calculated by assigning weights to regression coefficients and scaling the total to 5.
In the final model, the total predictor score was used to estimate the likelihood of unanticipated difficult airways. Youden’s index determined the cutoff value that maximized specificity and sensitivity. The performance of the final model was reported as the area under the receiver operating characteristic curve (AUC), with an AUC greater than 80 indicating excellent predictive accuracy.