Establishing a nomogram to predict refracture after percutaneous kyphoplasty by logistic regression
Aiqi Zhang
Hongye Fu2 - 1The Second Clinical Medical College of Zhejiang Chinese Medical University, Hangzhou, Zhejiang, China
- 2The Fourth Clinical Medical College of Zhejiang Chinese Medical University, Hangzhou, Zhejiang, China
- 3Department of Orthopedics, The Second Affiliated Hospital of Zhejiang Chinese Medical University, Hangzhou, Zhejiang, China
Introduction: Several studies have examined the risk factors for post-percutaneous kyphoplasty (PKP) refractures and developed many clinical prognostic models. However, no prior research exists using the Random Forest (RF) model, a favored tool for model development, to predict the occurrence of new vertebral compression fractures (NVCFs). Therefore, this study aimed to investigate the risk factors for the occurrence of post-PKP fractures, compare the predictive performance of logistic regression and RF models in forecasting post-PKP fractures, and visualize the logistic regression model.
Methods: We collected clinical data from 349 patients who underwent PKP treatment at our institution from January 2018 to December 2021. Lasso regression was employed to select risk factors associated with the occurrence of NVCFs. Subsequently, logistic regression and RF models were established, and their predictive capabilities were compared. Finally, a nomogram was created.
Results: The variables selected using Lasso regression, including bone density, cement distribution, vertebral fracture location, preoperative vertebral height, and vertebral height restoration rate, were included in both the logistic regression and RF models. The area under the curves of the logistic regression and RF models were 0.868 and 0.786, respectively, in the training set and 0.786 and 0.599, respectively, in the validation set. Furthermore, the calibration curve of the logistic regression model also outperformed that of the RF model.
Conclusion: The logistic regression model provided better predictive capabilities for identifying patients at risk for post-PKP vertebral fractures than the RF model.
1 Introduction
Osteoporotic vertebral compression fractures are one of the most prominent etiological factors that contribute to back pain, reduced daily activity levels, and increased bed rest durations in patients. This condition diminishes the quality of life and is associated with elevated mortality rates in more severe cases (Lyritis et al., 1989; Goldstein et al., 2015). Therapeutic interventions for osteoporotic vertebral compression fractures involve conservative and surgical approaches. Although conservative management provides symptom alleviation, it is burdened by the risk of malunion and the progression of kyphotic deformities, culminating in symptom recurrence (Soultanis et al., 2021). Percutaneous kyphoplasty (PKP) has garnered widespread recognition for its efficacy in addressing osteoporotic vertebral compression fractures. It has also demonstrated analgesic effects exceeding 90%, facilitating the early mobilization of patients and mitigating the incidence of complications (Li et al., 2015). Scholarly attention has been drawn to post-PKP cases with the evolution of PKP techniques, presenting with new vertebral compression fractures (NVCFs). Reports have cited NVCF occurrence rates ranging from 2 to 23% (Chen et al., 2021). Furthermore, several variables have been implicated in NVCF etiology, including age, bone mineral density (BMD), body mass index, the volume of administered bone cement, and instances of cement leakage (Bian et al., 2022; Zhang et al., 2023).
Numerous investigations have scrutinized the risk factors for post-PKP refractures and developed several clinical prognostic models (Li et al., 2021; Bian et al., 2022; Gao et al., 2023). However, most of these studies predominantly employed univariate or multivariate analyses, and the importance of the variables related to risk factors was usually not elucidated. Therefore, this ambiguity may impede clinical decision-making regarding the focal points of preventive measures post-PKP surgery. The Random Forest (RF) algorithm is a favored tool for model development and is distinguished as one of the most proficient techniques available (Iverson et al., 2008). It is also invaluable for ranking the relative importance of risk factors beyond model development, thereby assisting clinicians in pinpointing pivotal determinants for NVCF prevention. However, no prior research utilizing the RF model to predict NVCF occurrence exists. Therefore, this study aimed to compare the predictive performance of logistic regression and RF prediction models for NVCF. Specifically, this study’s objective was to select the most suitable model for predicting NVCF, providing a reference for tailoring individualized treatment plans for the early prevention of post-PKP fractures.
2 Materials and methods
2.1 Study population
We conducted a comprehensive retrospective analysis of clinical data obtained from patients who underwent PKP for single-level osteoporotic vertebral compression fractures in our healthcare institution between January 2018 and December 2021. Our inclusion criteria comprised the following: (1) the occurrence of a single-level vertebral body fracture resulting from low-energy trauma, such as a fall, which subsequently led to lower back pain or restricted mobility and subsequent PKP intervention; (2) radiographic evaluations, including X-rays and computed tomography scans, indicating the presence of vertebral compression fractures, corroborated by magnetic resonance imaging confirming recent vertebral fractures; and (3) the availability of complete clinical data. The exclusion criteria were as follows: (1) vertebral fractures caused by tumors, infections, or tuberculosis; (2) vertebral fractures resulting from high-energy trauma, such as a fall from a height or a motor vehicle accident; (3) spinal cord compression with neurological symptoms; (4) poor cardiopulmonary function rendering the patient unfit for surgery; and (5) secondary osteoporosis.
A total of 349 patients who met the inclusion criteria were included in the study. The fractured segments were distributed as follows: T5–T9, T10–L2, and L3–L5 in 31, 248, and 70 cases, respectively. The enrolled patients were randomly categorized into the training and validation groups in a 7:3 ratio. The training and validation groups were used for variable selection and model development and to assess the model’s performance, respectively.
This study was reviewed and approved by the Ethics Committee of Zhejiang Provincial People’s Hospital (QT2023051). The need for obtaining informed consent was exempted by the ethics committee because of the retrospective nature of this study.
2.2 Observation indicators
In a retrospective analysis, we collected data on age, sex, diabetes, hypertension history, body mass index, BMD, the vertebral segment of the fracture, the average volume of bone cement used, bone cement leakage, bone cement dispersion, bone cement distribution (Figure 1), contact of bone cement with the endplate, osteoporosis treatment, spinal curvature, preoperative and postoperative Cobb angle, and preoperative and postoperative anterior vertebral height (AVH) (Table 1). Adequate dispersion of bone cement on spinal X-rays was defined as the cement crossing the midline of the vertebral body; otherwise, it was considered poor dispersion. The AVH recovery ratio was defined as the difference between postoperative and preoperative AVH divided by the average of the adjacent vertebrae to the injured vertebra. Furthermore, the Cobb angle was defined as the angle formed by the upper and lower endplates of the fractured vertebra, while the Cobb angle recovery rate was the preoperative Cobb angle divided by the postoperative Cobb angle in percentiles.
Figure 1. (A) Bone cement in contact with the upper edge of the vertebral endplate. (B) Bone cement in contact with the lower edge of the vertebral endplate. (C) Bone cement in contact with both the upper and lower edges of the vertebral endplate. (D) Bone cement does not contact either the upper or lower edges of the vertebral endplate.
Table 1. Basic characteristics of clinical data.
2.3 Surgical methods
Patients were placed in the prone position, and a C-arm X-ray machine facilitated anteroposterior and lateral fluoroscopic localization of the bilateral pedicles of the affected vertebrae. Following standard disinfection and draping procedures, local infiltration anesthesia was administered at the puncture sites of the bilateral pedicles. The procedure was initiated on the left side with an approximately 1 cm longitudinal incision at the left pedicle’s projection. Next, a puncture was carefully performed at the 10 o’clock position after dissecting the skin and subcutaneous tissues, with attention to sagittal and angulation angles. Under fluoroscopic guidance using the C-arm X-ray machine in the lateral position, confirmation of the puncture needle’s central alignment within the pedicle was ensured. Subsequently, the puncture needle was withdrawn, and a guiding needle was introduced. A working cannula was placed approximately 0.5 cm anterior to the vertebral body’s posterior margin, replacing the guiding needle. A tamp was inserted, followed by an inflated balloon with a contrast medium to the appropriate pressure. Next, the balloon was extracted after the contrast medium was removed. The procedure was repeated on the right side, with a puncture made at the 2 o’clock position. Slow injection of bone cement was conducted under fluoroscopic guidance, establishing an effective mechanical column. Finally, the cannula was removed following cement solidification, and cement injection was promptly halted in cases of cement leakage during the procedure.
2.4 Postoperative management and refracture assessment
All patients received calcium supplements and active vitamin D combined with zoledronic acid or denosumab as standard treatment postoperatively. Thoracolumbar braces were worn for 3 months, and postoperative spinal X-rays, including anteroposterior and lateral views, were obtained. Spinal magnetic resonance imaging was employed for patients suspected of having vertebral refractures. The diagnostic criteria for post-PKP refractures encompassed the following: (1) reappearance of lower-back pain or restricted lumbar mobility post-PKP; (2) spinal magnetic resonance imaging indicating a new fracture characterized by low and high signals on T1- and T2-weighted imaging, respectively; and (3) the distinct identification of the fractured vertebra from the preoperative one.
2.5 Statistical analysis
Data analysis and graphical representation utilized IBM SPSS Statistics for Windows, version 26 (IBM Corp., Armonk, N.Y., United States) and R 4.22 software. Continuous variables were expressed as means ± standard deviations and compared using independent samples t-tests. Categorical variables were presented as percentages and compared using chi-square tests. Lasso regression aided in the selection of risk factors. The parameters selected by Lasso regression informed the construction of logistic regression and RF predictive models. Lasso regression analysis was conducted using the “glmnet” package. Given its non-intuitive output, variable importance was employed to elucidate the significance of each predictive factor in the RF model. The “randomForest” package was employed for the RF model and variable importance analysis. Additionally, the receiver operating characteristic curve was constructed, and the area under the curve value was calculated to assess sensitivity and specificity. The receiver operating characteristic curve was generated using the “pROC” package. A calibration plot was generated to examine the performance characteristics of the predictive model. Furthermore, the calibration curve was plotted using the “ggplot2” package, and statistical significance was considered at p < 0.05.
3 Results
A cohort of 349 patients was examined, with 48 (13.7%) experiencing refractures. The sex distribution in the control group was 56 (18.4%) males and 245 (81.6%) females, with a mean age and BMD of 73.85 ± 9.68 years and − 2.46 ± 1.27, respectively. In contrast, the refracture group included 7 (14.6%) males and 41 (85.4%) females, with a mean age and BMD of 78.04 ± 9.11 years and − 3.23 ± 1.31, respectively. Moreover, osteoporosis treatment was more prevalent in the control group than in the refracture group. The vertebral height recovery rate was significantly lower in the control group than in the refracture group. Table 1 summarizes the additional baseline characteristics of the enrolled patients.
3.1 Key variables
We observed coefficient variations for each variable using Lasso regression for variable selection (Figure 2A). Additionally, we identified an optimal model with superior performance and minimal variables at λ = 0.051 (ln λ = −2.98) employing tenfold cross-validation (Figure 2B). The selected variables comprised BMD, fracture location, bone cement shape, preoperative vertebral height, and vertebral height recovery rate.
Figure 2. Variable selection using Lasso regression. (A) Coefficient fluctuations for variables. (B) Lasso regression λ determined through tenfold cross-validation.
3.2 Establishment of logistic regression and RF models
We further established a logistic regression model based on the parameters identified via Lasso regression (Table 2), and this model demonstrated training set metrics of accuracy, specificity, sensitivity, and F1 score at 0.780, 0.800, 0.777, and 0.859, respectively. The training set area under curve, accuracy, specificity, sensitivity, and F1 score were 0.825, 0.846, 0.822, and 0.892, respectively (Table 3).
Table 2. Parameters for the logistic regression model.
Table 3. Metrics for the training and validation sets.
We constructed an RF model in the training set using the same set of variables. With an ntree value of 28, the out-of-bag error rate was low at 11.38% (Figure 3A). The training set exhibited accuracy, specificity, sensitivity, and F1 score metrics of 0.943, 0.971, 0.938, and 0.954, respectively, whereas the validation set metrics showed accuracy, specificity, sensitivity, and F1 score of 0.757, 0.615, 0.778, and 0.687, respectively (Table 3) (Figure 3B).
Figure 3. (A) Random Forest error rate. (B) Out-of-bag variable importance ranking. AVH, anterior vertebral height; AVHRR, anterior vertebral height recovery ratio; BMD, bone mineral density.
3.3 Performance of logistic regression and RF models
The performance of the logistic regression and RF models was evaluated using receiver operating characteristic and calibration curves. The area under the curve for the logistic regression and RF models within the training set were 0.868 and 0.786, respectively, but were 0.786 and 0.599 in the validation set, respectively (Figure 4). The logistic regression model consistently outperformed the RF model in the training and validation sets. Calibration curves were subsequently constructed, and the logistic regression model outperformed the RF model in both the training and validation sets (Figure 5).
Figure 4. Area under the curves (AUCs) for the logistic regression model in the (A) training and (C) validation sets. AUCs for the Random Forest model in the (B) training and (D) validation sets.
Figure 5. Calibration curves for the Logistic regression model in the (A) training and (C) validation sets. Calibration curves for the Random Forest model in the (B) training and (D) validation sets.
3.4 Nomogram development
We selected the logistic regression model to predict the recurrence risk of NCVF post-PKP. Next, we transformed the complex mathematical model into a graphical format for improved clinical applicability. For instance, if a patient is male, has a T9 vertebral fracture, a BMD of −3.0, vertebral height recovery >14%, a bone cement distribution shape of D, and preoperative vertebral height of 20 mm, their total score would be approximately 169, indicating a risk of approximately 64% for developing NCVF post-PKP (Figure 6).
Figure 6. Nomogram for predicting post-PKP refractures. AVH, anterior vertebral height; AVHRR, anterior vertebral height recovery ratio; BMD, bone mineral density; PKP, percutaneous kyphoplasty.
4 Discussion
PKP provides rapid relief from severe pain that is difficult to manage conservatively and effectively reduces cement leakage due to the presence of an expandable balloon (Chen et al., 2018; Hu et al., 2018). It allows for the rapid restoration of spinal lordosis, thereby improving biomechanical stability (Lee et al., 2014; Cao et al., 2020). Furthermore, PKP enables early mobilization, reducing bed rest time and the risk of complications. Common post-PKP complications include NCVF, pulmonary embolism, superior vena cava injury, and postoperative back pain, with NCVF being particularly common and severely affecting patients’ quality of life (Li et al., 2021; Yu et al., 2021; Zou et al., 2021). However, some patients experienced multiple NVCFs, which significantly undermined the confidence in these patients’ recovery and increased the economic and time costs associated with their treatment. Therefore, analyzing the risk factors for NVCF post-PKP and establishing a predictive model to assess this risk is essential. This can enhance patient satisfaction during the treatment process and contribute to the effective use of healthcare resources.
This study involving 349 patients resulted in the development of a nomogram based on several variables selected using Lasso regression. The predictive model incorporates five clinical assessment indicators as follows: BMD, bone cement distribution (shape D), T5–T9 vertebral fractures, preoperative vertebral height, and vertebral height recovery rate. Preoperative vertebral height and BMD evidently play a crucial role in the occurrence of post-PKP fractures (Figure 3B). Therefore, this may become a key focus for clinicians in the management of high-risk patients post-PKP in the future. These indicators allow for the evaluation of individual NCVF risk, assigning scores based on each variable’s contribution to the outcome and summing these scores. Drawing a vertical line at the corresponding total score position intersects with the line representing the risk of recurrence.
Lasso regression can address the issue of multicollinearity among variables, which is why this study selected it for variable selection. The logistic regression and RF models showed similar performance on the training dataset; however, the RF model exhibited poorer performance on the validation dataset. Therefore, a nomogram was developed in this study based on the logistic regression model. Previous research has also compared the predictive capabilities of both models for binary classification issues, and the results consistently indicate that the RF model performs poorly (Wang et al., 2022). This aligns with our study findings. Therefore, this study does not recommend using the RF model to predict the risk of post-PKP refracture.
Osteoporosis is a skeletal disorder characterized by reduced bone mass and microstructural deterioration, increasing bone fragility and fracture susceptibility (Sambrook and Cooper, 2006). BMD reflects the severity of osteoporosis, exhibiting a strong correlation with fracture occurrence (Black et al., 2020). Reduced BMD is typically associated with decreased bone formation capacity, vertebral height loss, and necrotic bone formation (Qi et al., 2021). Furthermore, changes in BMD effectively predict vertebral fractures, with increased BMD almost certainly lowering the risk of vertebral fractures (Bouxsein et al., 2019). Bisphosphonates and denosumab are first-line treatments for osteoporosis, effectively increasing BMD (Reid and Billington, 2022). Additionally, anti-osteoporosis treatment significantly reduces fracture risk, and the distribution of bone cement within the vertebra can impact NCVF occurrence. However, poor cement distribution may provide inadequate support, leading to NCVF. Ideally, bone cement should be in contact with the upper and lower vertebral endplates, which better restores vertebral strength and reduces the risk of NCVF (Tan et al., 2020). Our study found that patients with T5–T9 vertebral fractures were more prone to NCVF than those with T10–L2 fractures. T10–L2 is susceptible to compression fractures because it represents the junction between thoracic kyphosis and lumbar lordosis, concentrating trunk mobility stress. Clinical data review revealed that patients with T5–T9 vertebral fractures had lower BMD levels (−3.04 ± 1.24) than those with T10–L2 and L3–L5 fractures (−2.54 ± 1.26 and − 2.42 ± 1.45, respectively) (Figure 7). Furthermore, patients with T5–T9 vertebral fractures had lower BMD, suggesting a potential association between lower BMD and NCVF occurrence.
Figure 7. Distribution of vertebral fracture sites. BMD, bone mineral density.
Previous studies have emphasized the significance of the AVH recovery ratio in post-PKP refracture (Niu et al., 2015; Li et al., 2018). Specifically, Lin et al. (2008) reported a significant association between AVH recovery ratio and post-PKP refracture. Excessive restoration of vertebral height may lead to increased tension in the adjacent soft tissues, resulting in elevated mechanical loads on the vertebrae or increased instability in the fractured segment (Kim et al., 2004, 2010; Niu et al., 2015) and can elevate intra-vertebral pressure, promoting bone necrosis (Heo et al., 2009). Furthermore, our model incorporated the preoperative AVH as a risk factor. The risk of post-PKP refracture was found to increase with decreasing AVH. Low AVH may also indicate significant trabecular compression within the fractured vertebra, leading to increased internal pressure, which could hinder the adequate penetration of cement (Yu et al., 2018; Zhang et al., 2023). Notably, the primary goal of PKP is pain relief in vertebral fractures; therefore, the excessive pursuit of vertebral height restoration may not always be necessary.
Previous literature has reviewed the risk factors for refracture and established predictive nomograms with excellent performance, providing valuable guidance for clinical practice (Li et al., 2021; Ju and Liu, 2023; Zhang et al., 2023). Indeed, it should be emphasized that the use of Lasso regression and multivariate logistic regression for variable selection in these studies is one of the contributing factors to the disparity in selected variables. The dataset size used to build the nomogram and the variables of interest for researchers also play a role. Consequently, it is anticipated that different datasets and risk factors may yield diverse sets of refracture risk factors in future studies. We recommend that researchers make every effort to employ consistent methods for variable selection and use the same set of variables to mitigate the impact of the factors mentioned above on study outcomes. Notably, this approach can enhance the credibility and comparability of research findings.
This study had some limitations. First, it was a retrospective study, which may have introduced selection bias and data loss. Second, the nomogram has not been externally validated in other regions or multi-center settings. Therefore, we aim to validate the nomogram’s performance with larger samples and multi-center data in future research.
In conclusion, we established a nomogram for predicting post-PKP refracture based on a logistic regression model. This nomogram can assist in screening high-risk patients for post-PKP refracture. Therefore, clinicians can use the predicted risk of recurrence to develop personalized follow-up and treatment plans.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the Ethics Committee of Zhejiang Provincial People's Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
AZ: Data curation, Methodology, Writing – original draft, Writing – review & editing. HF: Data curation, Investigation, Methodology, Writing – review & editing. JW: Methodology, Software, Supervision, Writing – original draft. ZC: Writing – review & editing. JF: Funding acquisition, Investigation, Methodology, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by a grant from the National Nature and Science Foundation of China (81873128).
Acknowledgments
The authors express their sincere gratitude to all the patients who participated in this study and to all the staff members in our department for their dedication and contribution.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
AVH, anterior vertebral height; PKP, percutaneous kyphoplasty; BMD, bone mineral density; \NVCFs, new vertebral compression fractures; RF, Random Forest
References
Bian, F., Bian, G., Zhao, L., Huang, S., Fang, J., and An, Y. (2022). Risk factors for recollapse of new vertebral compression fractures after percutaneous kyphoplasty in geriatric patients: establishment of a nomogram. BMC Musculoskelet. Disord. 23:458. doi: 10.1186/s12891-022-05409-3
Black, D. M., Bauer, D. C., Vittinghoff, E., Lui, L. Y., Grauer, A., Marin, F., et al. (2020). Treatment-related changes in bone mineral density as a surrogate biomarker for fracture risk reduction: meta-regression analyses of individual patient data from multiple randomised controlled trials. Lancet Diabetes Endocrinol. 8, 672–682. doi: 10.1016/S2213-8587(20)30159-5
Bouxsein, M. L., Eastell, R., Lui, L. Y., Wu, L. A., de Papp, A. E., Grauer, A., et al. (2019). Change in bone density and reduction in fracture risk: a meta-regression of published trials. J. Bone Miner. Res. 34, 632–642. doi: 10.1002/jbmr.3641
Cao, Z., Wang, G., Hui, W., Liu, B., Liu, Z., and Sun, J. (2020). Percutaneous kyphoplasty for osteoporotic vertebral compression fractures improves spino-pelvic alignment and global sagittal balance maximally in the thoracolumbar region. PLoS One 15:e0228341. doi: 10.1371/journal.pone.0228341
Chen, X., Guo, W., Li, Q., Ou, Z., Lao, Z., Liu, Y., et al. (2018). Is unilateral percutaneous kyphoplasty superior to bilateral percutaneous kyphoplasty for osteoporotic vertebral compression fractures? Evidence from a systematic review of discordant meta-analyses. Pain Physician 21, 327–336. doi: 10.36076/ppj.2018.4.327
Chen, Z., Song, C., Lin, H., Sun, J., and Liu, W. (2021). Does prophylactic vertebral augmentation reduce the refracture rate in osteoporotic vertebral fracture patients: a meta-analysis. Eur. Spine J. 30, 2691–2697. doi: 10.1007/s00586-021-06899-w
Gao, W., Chen, Y., Wang, X., Liu, G., Cui, K., Guo, J., et al. (2023). Establishment and verification of a predictive nomogram for new vertebral compression fracture occurring after bone cement injection in middle-aged and elderly patients with vertebral compression fracture. Orthop. Surg. 15, 961–972. doi: 10.1111/os.13655
Goldstein, C. L., Chutkan, N. B., Choma, T. J., and Orr, R. D. (2015). Management of the elderly with vertebral compression fractures. Neurosurgery 77, S33–S45. doi: 10.1227/NEU.0000000000000947
Heo, D. H., Chin, D. K., Yoon, Y. S., and Kuh, S. U. (2009). Recollapse of previous vertebral compression fracture after percutaneous vertebroplasty. Osteoporos. Int. 20, 473–480. doi: 10.1007/s00198-008-0682-3
Hu, K. Z., Chen, S. C., and Xu, L. (2018). Comparison of percutaneous balloon dilation kyphoplasty and percutaneous vertebroplasty in treatment for thoracolumbar vertebral compression fractures. Eur. Rev. Med. Pharmacol. Sci. 22, 96–102. doi: 10.26355/eurrev_201807_15370
Iverson, L. R., Prasad, A. M., Matthews, S. N., and Peters, M. (2008). Estimating potential habitat for 134 eastern US tree species under six climate scenarios. For. Ecol. Manag. 254, 390–406. doi: 10.1016/j.foreco.2007.07.023
Ju, G., and Liu, X. (2023). A nomogram prediction model for refracture in elderly patients with osteoporotic vertebral compression fractures after percutaneous vertebroplasty. Eur. Spine J. 32, 3919–3926. doi: 10.1007/s00586-023-07843-w
Kim, S. H., Kang, H. S., Choi, J. A., and Ahn, J. M. (2004). Risk factors of new compression fractures in adjacent vertebrae after percutaneous vertebroplasty. Acta Radiol. 45, 440–445. doi: 10.1080/02841850410005615
Kim, Y. J., Lee, J. W., Kim, K. J., Chung, S. K., Kim, H. J., Park, J. M., et al. (2010). Percutaneous vertebroplasty for intravertebral cleft: analysis of therapeutic effects and outcome predictors. Skelet. Radiol. 39, 757–766. doi: 10.1007/s00256-009-0866-8
Lee, J. H., Lee, D. O., Lee, J. H., and Lee, H. S. (2014). Comparison of radiological and clinical results of balloon kyphoplasty according to anterior height loss in the osteoporotic vertebral fracture. Spine J. 14, 2281–2289. doi: 10.1016/j.spinee.2014.01.028
Li, Y. X., Guo, D. Q., Zhang, S. C., Liang, D., Yuan, K., Mo, G. Y., et al. (2018). Risk factor analysis for re-collapse of cemented vertebrae after percutaneous vertebroplasty (PVP) or percutaneous kyphoplasty (PKP). Int. Orthop. 42, 2131–2139. doi: 10.1007/s00264-018-3838-6
Li, Q., Long, X., Wang, Y., Fang, X., Guo, D., Lv, J., et al. (2021). Development and validation of a nomogram for predicting the probability of new vertebral compression fractures after vertebral augmentation of osteoporotic vertebral compression fractures. BMC Musculoskelet. Disord. 22:957. doi: 10.1186/s12891-021-04845-x
Li, L., Ren, J., Liu, J., Wang, H., Wang, X., Liu, Z., et al. (2015). Results of vertebral augmentation treatment for patients of painful osteoporotic vertebral compression fractures: a meta-analysis of eight randomized controlled trials. PLoS One 10:e0138126. doi: 10.1371/journal.pone.0138126
Li, Q. D., Yang, J. S., Gong, H. L., Wei, J. M., Wang, X. F., Zou, P., et al. (2021). Can additional facet joint block improve the clinical outcome of kyphoplasty for acute osteoporotic vertebral compression fractures. Pain Physician 24, 283–291. doi: 10.36076/ppj.2021/24/283
Lin, W. C., Lee, Y. C., Lee, C. H., Kuo, Y. L., Cheng, Y. F., Lui, C. C., et al. (2008). Refractures in cemented vertebrae after percutaneous vertebroplasty: a retrospective analysis. Eur. Spine J. 17, 592–599. doi: 10.1007/s00586-007-0564-y
Lyritis, G. P., Mayasis, B., Tsakalakos, N., Lambropoulos, A., Gazi, S., Karachalios, T., et al. (1989). The natural history of the osteoporotic vertebral fracture. Clin. Rheumatol. 8, 66–69. doi: 10.1007/BF02207237
Niu, J., Zhou, H., Meng, Q., Shi, J., Meng, B., and Yang, H. (2015). Factors affecting recompression of augmented vertebrae after successful percutaneous balloon kyphoplasty: a retrospective analysis. Acta Radiol. 56, 1380–1387. doi: 10.1177/0284185114556016
Qi, H., Qi, J., Gao, J., Sun, J., and Wang, G. (2021). The impact of bone mineral density on bone metabolism and the fracture healing process in elderly chinese patients with osteoporotic vertebral compression fractures. J. Clin. Densitom. 24, 135–145. doi: 10.1016/j.jocd.2020.11.003
Reid, I. R., and Billington, E. O. (2022). Drug therapy for osteoporosis in older adults. Lancet 399, 1080–1092. doi: 10.1016/S0140-6736(21)02646-5
Sambrook, P., and Cooper, C. (2006). Osteoporosis. Lancet 367, 2010–2018. doi: 10.1016/S0140-6736(06)68891-0
Soultanis, K., Thano, A., and Soucacos, P. N. (2021). Outcome of thoracolumbar compression fractures following non-operative treatment. Injury 52, 3685–3690. doi: 10.1016/j.injury.2021.05.019
Tan, L., Wen, B., Guo, Z., and Chen, Z. (2020). The effect of bone cement distribution on the outcome of percutaneous Vertebroplasty: a case cohort study. BMC Musculoskelet. Disord. 21:541. doi: 10.1186/s12891-020-03568-9
Wang, Q., Qiao, W., Zhang, H., Liu, B., Li, J., Zang, C., et al. (2022). Nomogram established on account of lasso-cox regression for predicting recurrence in patients with early-stage hepatocellular carcinoma. Front. Immunol. 13:1019638. doi: 10.3389/fimmu.2022.1019638
Yu, D., Liu, Z., Wang, H., Yao, R., Li, F., Yang, Y., et al. (2021). Analysis on the effect of different surgical methods on the treatment of senile osteoporotic spinal compression fractures and the influencing factors of complications. Evid. Based Complement. Alternat. Med. 2021, 1–8. doi: 10.1155/2021/1599470
Yu, W., Xu, W., Jiang, X., Liang, D., and Jian, W. (2018). Risk factors for recollapse of the augmented vertebrae after percutaneous vertebral augmentation: a systematic review and meta-analysis. World Neurosurg. 111, 119–129. doi: 10.1016/j.wneu.2017.12.019
Zhang, A., Lin, Y., Kong, M., Chen, J., Gao, W., Fan, J., et al. (2023). A nomogram for predicting the risk of new vertebral compression fracture after percutaneous kyphoplasty. Eur. J. Med. Res. 28:280. doi: 10.1186/s40001-023-01235-y
Keywords: vertebral compression fractures, percutaneous kyphoplasty, logistic regression, random forest model, anterior vertebral height
Citation: Zhang A, Fu H, Wang J, Chen Z and Fan J (2023) Establishing a nomogram to predict refracture after percutaneous kyphoplasty by logistic regression. Front. Neuroinform. 17:1304248. doi: 10.3389/fninf.2023.1304248
Edited by:
Rashid Mehmood, King Abdulaziz University, Saudi ArabiaReviewed by:
Chang Qing Zhao, Shanghai Jiao Tong University, ChinaGuangtao Xu, Jiaxing University Medical College, China
Copyright © 2023 Zhang, Fu, Wang, Chen and Fan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jiajun Fan, Apoptoxin4869@outlook.com