Objective This study aimed to develop and evaluate machine learning (ML) models for predicting survival in patients with rare and heterogeneous sinonasal malignancies, identify significant prognostic factors, and compare their performance with traditi...
Objective This study aimed to develop and evaluate machine learning (ML) models for predicting survival in patients with rare and heterogeneous sinonasal malignancies, identify significant prognostic factors, and compare their performance with traditional statistical methods. The secondary objectives were to extensively integrate large language models (LLMs) into the ML pipeline to assess their feasibility in reliably building ML models and to explore their practical applicability in medical research. Methods We retrospectively analyzed a cohort of 621 patients with sinonasal malignancies treated at a single tertiary center between 1990 and 2024. LLMs were incorporated into the analytic workflow for data preprocessing, feature engineering, and the development of multiple ML models, including logistic regression, random forest, XGBoost, LightGBM, support vector classifier (linear SVC), and multilayer perceptron for 5-year survival classification, as well as Cox proportional hazards regression, random survival forest, XGBoost-Survival, LightGBM-Survival, and DeepSurv for survival analysis. All the steps of model training, hyperparameter tuning, and variable visualization were initially attempted within the LLM environment; when execution was not feasible, the code generated by the LLM was externally verified and run by investigators. Subgroup analyses were performed for squamous cell carcinoma and epithelial and non-epithelial tumors, and additional models were constructed for recurrence prediction. Results For 5-year survival prediction in the overall cohort, random forest achieved the highest area under the receiver operating characteristic curve (AUROC) on the V1 set (0.861), while LightGBM yielded the best F1-score (0.718). The V2 feature set notably improved linear models, with the AUROC of logistic regression increasing to 0.854 (Δ +0.025). In the subgroup analyses, the performance was particularly strong in SCC, where the V2 set increased the AUROC of logistic regression to 0.906. For the time-to-death analysis, the tuned XGBoost-Survival model demonstrated superior performance (C-index 0.777), significantly outperforming the traditional Cox proportional hazards model (C-index 0.752). Across diverse models and subgroups, residual disease after initial treatment, distant metastasis, and age at diagnosis consistently emerged as the most powerful prognostic factors. The LLM-integrated workflow proved feasible; however, certain models including LightGBM and all survival analysis algorithms required external execution due to library constraints. Conclusions ML models provided more accurate survival prediction than conventional methods and consistently identified key prognostic factors including residual disease, distant metastasis, and age. By integrating LLMs into the analytic workflow, this study presents a practical framework that links advanced computational methods with clinical applicability, enhancing prognostic accuracy and supporting the adoption of AI-based prognostic modeling for rare sinonasal cancers in otolaryngology.