Application of machine learning model to predict osteoporosis based on abdominal computed tomography images of the psoas muscle: a retrospective study
Authors: Huang Cheng-bin , Hu Jia-sen , Tan Kai , Zhang Wei , Xu Tian-hao , Yang Lei 
Published: October 2022

Background

With rapid economic development, the world's average life expectancy is increasing, leading to the increasing prevalence of osteoporosis worldwide. However, due to the complexity and high cost of dual-energy x-ray absorptiometry (DXA) examination, DXA has not been widely used to diagnose osteoporosis. In addition, studies have shown that the psoas index measured at the third lumbar spine (L3) level is closely related to bone mineral density (BMD) and has an excellent predictive effect on osteoporosis. Therefore, this study developed a variety of machine learning (ML) models based on psoas muscle tissue at the L3 level of unenhanced abdominal computed tomography (CT) to predict osteoporosis.

Methods

Medical professionals collected the CT images and the clinical characteristics data of patients over 40 years old who underwent DXA and abdominal CT examination in the Second Affiliated Hospital of Wenzhou Medical University database from January 2017 to January 2021. Using 3D Slicer software based on horizontal CT images of the L3, the specialist delineated three layers of the region of interest (ROI) along the bilateral psoas muscle edges. The PyRadiomics package in Python was used to extract the features of ROI. Then Mann–Whitney U test and the least absolute shrinkage and selection operator (LASSO) algorithm were used to reduce the dimension of the extracted features. Finally, six machine learning models, Gaussian naïve Bayes (GNB), random forest (RF), logistic regression (LR), support vector machines (SVM), Gradient boosting machine (GBM), and Extreme gradient boosting (XGBoost), were applied to train and validate these features to predict osteoporosis.

Results

A total of 172 participants met the inclusion and exclusion criteria for the study. 82 participants were enrolled in the osteoporosis group, and 90 were in the non-osteoporosis group. Moreover, the two groups had no significant differences in age, BMI, sex, smoking, drinking, hypertension, and diabetes. Besides, 826 radiomic features were obtained from unenhanced abdominal CT images of osteoporotic and non-osteoporotic patients. Five hundred fifty radiomic features were screened out of 826 by the Mann–Whitney U test. Finally, 16 significant radiomic features were obtained by the LASSO algorithm. These 16 radiomic features were incorporated into six traditional machine learning models (GBM, GNB, LR, RF, SVM, and XGB). All six machine learning models could predict osteoporosis well in the validation set, with the area under the receiver operating characteristic (AUROC) values greater than or equal to 0.8. GBM is more effective in predicting osteoporosis, whose AUROC was 0.86, sensitivity 0.70, specificity 0.92, and accuracy 0.81 in validation sets.

Conclusion

We developed six machine learning models to predict osteoporosis based on psoas muscle images of abdominal CT, and the GBM model had the best predictive performance. GBM model can better help clinicians to diagnose osteoporosis and provide timely anti-osteoporosis treatment for patients. In the future, the research team will strive to include participants from multiple institutions to conduct external validation of the ML model of this study.