收稿日期: 2026-03-06
修回日期: 2026-05-26
录用日期: 2026-05-27
网络出版日期: 2026-08-31
基金资助
教育部中国高校产学研创新基金(2024DR031)
Research on construction of osteoporotic fracture prediction models based on machine learning algorithms
Received date: 2026-03-06
Revised date: 2026-05-26
Accepted date: 2026-05-27
Online published: 2026-08-31
目的: 基于骨质疏松患者的临床资料,采用5种机器学习算法分别构建骨质疏松性骨折风险预测模型,并对模型的效能进行评价及比较。方法: 根据纳排标准,选取2020年1月至2025年12月于上海交通大学医学院附属松江医院就诊的333例骨质疏松症连续病例(男性66例,女性267例;中位年龄69岁),根据有无脆性骨折史分为骨折组(n=146)和未骨折组(n=187)。收集患者基线资料及实验室指标,缺失数据采用链式方程多重插补处理。采用R语言开展相关研究,通过逻辑回归、拉索回归筛选特征变量,随机森林进行重要性排序。按7∶3比例将数据分为训练集和测试集,采用逻辑回归、决策树(decision tree,DT)、随机森林、支持向量机和极限梯度提升(extreme gradient boosting,XGBoost)算法分别构建模型,采用受试者操作特征曲线的曲线下面积等指标评估模型预测骨质疏松性骨折风险的效能,并应用SHAP(SHapley Additive exPlanations)值对最优模型进行解释。结果: 单因素逻辑回归分析显示,性别、甲状旁腺激素(parathyroid hormone,PTH)、N端中段骨钙素(N-middle osteocalcin,N-MID-OT)、碱性磷酸酶、白蛋白、钙、平均红细胞血红蛋白含量(mean corpuscular hemoglobin,MCH)、平均红细胞血红蛋白浓度(mean corpuscular hemoglobin concentration,MCHC)、血红蛋白、磷与骨质疏松性骨折显著相关。拉索回归筛选出性别及PTH、N-MID-OT、白蛋白、钙、MCHC、磷水平共7个关键特征,进行随机森林重要性排序后,纳入性别、PTH、N-MID-OT、钙、磷、MCHC这6项指标进行机器学习建模。利用测试集数据作为输入,比较每种机器学习模型预测骨质疏松骨折中的效能,结果显示,其中XGBoost模型表现最佳,预测患者发生骨质疏松性骨折风险的曲线下面积为0.974(95%CI为0.950~0.998),其他4种机器学习模型也取得了良好的预测效能,曲线下面积介于0.891~0.964之间。SHAP结果显示,N-MID-OT、钙、PTH在模型中的权重较大,与临床认知一致。结论: 基于常规临床指标构建的XGBoost模型在骨质疏松性骨折风险预测中表现出良好的区分度和稳定性,具有潜在的临床应用价值,适用于无法完整开展骨代谢标志物检测的基层医院。
马善源 , 彭亮 , 侯彦强 , 刘维薇 . 基于机器学习算法的骨质疏松性骨折预测模型构建[J]. 诊断学理论与实践, 2026 , 25(04) : 467 -475 . DOI: 10.16150/j.1671-2870.2026.04.010
Objective Based on clinical data from patients with osteoporosis, five machine learning algorithms are used to construct osteoporotic fracture risk prediction models, and their performance is evaluated and compared. Methods A total of 333 consecutive patients with osteoporosis (66 males and 267 females; median age, 69 years) who attended Song-jiang Hospital Affiliated to Shanghai Jiao Tong University School of Medicine from January 2020 to December 2025 were selected according to inclusion and exclusion criteria. Patients were divided into a fracture group (n=146) and a non-fracture group (n=187) based on the presence or absence of a history of fragility fractures. Baseline patient data and laboratory indicators were collected, and missing data were handled using multiple imputation by chained equations. Relevant studies were conducted using the R programming language. Feature variables were screened using logistic regression and lasso regression, and variable importance was ranked using random forest. Data were split into a training set and a test set at a 7∶3 ratio. Models were constructed using logistic regression, decision tree (DT), random forest, support vector machine, and extreme gradient boosting (XGBoost) algorithms. The performance of the models in predicting osteoporotic fracture risk was evaluated using indicators such as the area under the receiver operating characteristic curve (AUC), accuracy, sensitivity, and specificity. The optimal model was interpreted using SHAP (SHapley Additive exPlanations) values. Results Univaria-te logistic regression analysis showed that sex, parathyroid hormone (PTH), N-terminal mid-fragment osteocalcin (N-MID-OT), alkaline phosphatase, albumin, calcium, mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), hemoglobin, and phosphorus were significantly associated with osteoporotic fractures. Lasso regression identified seven key features: sex, PTH, N-MID-OT, albumin, calcium, MCHC, and phosphorus. After random forest importance ranking, six variables-sex, PTH, N-MID-OT, calcium, phosphorus, and MCHC-were included in machine learning modeling. Using the test set data as input, the performance of each machine learning model in predicting osteoporotic fractures was compared. The results showed that the XGBoost model performed best with an AUC of 0.974 (95% CI: 0.950-0.998). The other four machine learning models also achieved good prediction performance, with AUCs ranging from 0.891 to 0.964. SHAP results indicated that N-MID-OT, calcium, and PTH had relatively large weights in the model, which was consistent with clinical cognition. Conclusions The XGBoost model constructed based on routine clinical indicators demonstrates good discriminative ability and stability in predicting the risk of osteoporotic fractures, has potential clinical application value, and is suitable for primary hospitals that are unable to conduct complete bone metabolism marker testing.
Key words: Machine learning; Osteoporotic fracture; Prediction model
| [1] | 中华人民共和国国家卫生健康委员会. 一图读懂:中国骨质疏松症流行病学调查结果[EB/OL]. (2018-10-20)[2026-06-30]. https://www.nhc.gov.cn/wjw/zcjd/201810/1aee4c12bec7415cbebfe7eacbc09552.shtml. |
| [2] | Wang L, Yu W, Yin X, et al. Prevalence of osteoporosis and fracture in China: The China osteoporosis prevalence study[J]. JAMA Netw Open, 2021, 4(8):e2121106. |
| [3] | Kanis J A, Johnell O, Oden A, et al. FRAX? and the assessment of fracture probability in men and women from the UK[J]. Osteoporos Int, 2008, 19(4):385-397. |
| [4] | Kanis J A, Harvey N C, Cooper C, et al. A systematic review of intervention thresholds based on FRAX: A report prepared for the National Osteoporosis Guideline Group and the International Osteoporosis Foundation[J]. Arch Osteoporos, 2016, 11:25. |
| [5] | Kanis J A, Cooper C, Rizzoli R, et al. European guidance for the diagnosis and management of osteoporosis in postmenopausal women[J]. Osteoporos Int, 2019, 30(1):3-44. |
| [6] | Schini M, Vilaca T, Gossiel F, et al. Bone turnover mar-kers: Basic biology to clinical applications[J]. Endocr Rev, 2023, 44(3):417-473. |
| [7] | Tian A, Ma J, Feng K, et al. Reference markers of bone turnover for prediction of fracture: A meta-analysis[J]. J Orthop Surg Res, 2019, 14(1):68. |
| [8] | Greener J G, Kandathil S M, Moffat L, et al. A guide to machine learning for biologists[J]. Nat Rev Mol Cell Biol, 2022, 23(1):40-55. |
| [9] | 倪康康, 董单, 李国庆, 等. 儿童及青少年骨质疏松症研究进展[J]. 中华内分泌代谢杂志, 2025, 41(5):430-434. |
| Ni Kangkang, Dong Dan, Li Guoqing, et al. Recent advances in osteoporosis in children and adolescents[J]. Chin J Endocrinol Metab, 2025, 41(5):430-434. | |
| [10] | Siris E S, Adler R, Bilezikian J, et al. The clinical diagnosis of osteoporosis: A position statement from the National Bone Health Alliance Working Group[J]. Osteoporos Int, 2014, 25(5):1439-1443. |
| [11] | Wu Y, Chao J, Bao M, et al. Predictive value of machine learning on fracture risk in osteoporosis: A systematic review and meta-analysis[J]. BMJ Open, 2023, 13(12):e071430. |
| [12] | Zhang D, Yang X, Wang F, et al. Comparative analysis of feature extraction methods and machine learning models for predicting osteoporosis prevalence[J]. J Med Syst, 2025, 49(1):72. |
| [13] | 陈虹谷, 赵国阳, 马安培, 等. 健康中老年人群血清骨转换标志物与FRAX骨折风险相关性分析[J]. 中国骨质疏松杂志, 2024, 30(2):204-209. |
| Chen Honggu, Zhao Guoyang, Ma Anpei, et al. Analysis of correlation between serum bone turnover markers and FRAX in healthy middle-aged and elderly population[J]. Chin J Osteoporos, 2024, 30(2):204-209. | |
| [14] | Komaba H, Imaizumi T, Hamano T, et al. Lower parathyroid hormone levels are associated with reduced fracture risk in Japanese patients on hemodialysis[J]. Kidney Int Rep, 2024, 9(10):2956-2969. |
| [15] | Ying P, Gu M, Jiang X, et al. Serum calcium-phosphorus product for predicting the risk of osteoporotic vertebral compression fractures in elderly patients: A retrospective observational study[J]. J Orthop Surg Res, 2022, 17(1):57. |
| [16] | Deng A F, Wang F X, Wang S C, et al. Bone-organ axes: Bidirectional crosstalk[J]. Mil Med Res, 2024, 11(1):37. |
| [17] | 刘玉洁, 孙世威, 孟霞. 中老年骨质疏松性骨折危险因素及预测模型[J]. 中国骨质疏松杂志, 2026, 32(2):171-175. |
| Liu Yujie, Sun Shiwei, Meng Xia. Risk factors and a predictive model for osteoporotic fracture in the middle age and elderly patients[J]. Chin J Osteoporos, 2026, 32(2):171-175. | |
| [18] | Qi X, Wang S, Fang C, et al. Machine learning and SHAP value interpretation for predicting comorbidity of cardiovascular disease and cancer with dietary antioxidants[J]. Redox Biol, 2025, 79:103470. |
| [19] | Liang D, Wang L, Zhong P, et al. Perspective: Global burden of iodine deficiency: Insights and projections to 2050 using XGBoost and SHAP[J]. Adv Nutr, 2025, 16(3):100384. |
| [20] | 中华医学会骨质疏松和骨矿盐疾病分会. 骨转换生化标志物临床应用指南[J]. 中华骨质疏松和骨矿盐疾病杂志, 2021, 14(4):321-336. |
| China Society of Osteoporosis and Bone Mineral Research, China Medical Association. Guidelines for clinical application of bone turnover biomarkers[J]. Chin J Osteoporos Bone Miner Res, 2021, 14(4):321-336. | |
| [21] | 肖文金, 刘彩蕊, 朱柯雨, 等. 肥胖、2型糖尿病与骨折风险的关联及健康管理:大数据研究提示新观点[J]. 中华内分泌代谢杂志, 2025, 41(10):818-822. |
| Xiao Wenjin, Liu Cairui, Zhu Keyu, et al. Obesity, type 2 diabetes and fracture risk: New insights from large-scale data studies[J]. Chin J Endocrinol Metab, 2025, 41(10):818-822. | |
| [22] | 任春蕊, 刘建凤, 安向莲, 等. 体重指数与骨质疏松症之间的因果关系:一项孟德尔随机化研究[J]. 中华内分泌代谢杂志, 2024, 40(2):108-114. |
| Ren Chunrui, Liu Jianfeng, An Xianglian, et al. Causal relationship between body mass index and osteoporosis: A Mendelian randomization study[J]. Chin J Endocrinol Metab, 2024, 40(2):108-114. |
/
| 〈 |
|
〉 |