Impact Factors on Posterior Modified Transfacet Debridement for Thoracic Spinal Tuberculosis Basing on Regression and Classification Analysis

Posterior transfacet approach has been proved to be a safe and effective access to treat thoracic disc herniation. However, the influencing factors of posterior modified transarticular debridement for thoracic tuberculosis have not been reported in the clinical literature. From 2009 to 2014, 37 patients with TST underwent a posterior modified transfacet debridement, interbody fusion following posterior instrumentation, under the cover of 18 months of antituberculosis chemotherapy. The patients were evaluated preoperatively and postoperatively in terms of Frankel Grade, visual analog scale (VAS) pain score, kyphotic Cobb angle, and bone fusion. Blood loss (positive correlation) and focal debridement (positive correlation) could affect operative time. Operative time (positive correlation) could affect blood loss. While, age (positive correlation), PostE (negative correlation), and T_FocalDebridement (positive correlation) could affect bone fusion. The accuracy of naive bayes classifier model is 86.11%. Our preliminary results show that blood loss and focal debridement could affect operative time; operative time could affect blood loss; age, PostE, and T_FocalDebridement could affect bone fusion; the naive Bayes classifier model can predict the KirkaldyWillis accurately.


Introduction
Tuberculosis (TB), an infectious disease caused by mycobacterium tuberculosis (MTB), is the ninth leading cause of global death, ranking above acquired immune deficiency syndrome (AIDS). e spine is the most frequently involved part of extrapulmonary tuberculosis, which is more common in thoracolumbar vertebrae, 30.3% to 55.8% of which are thoracic tuberculosis and about 67% of thoracic tuberculosis patients are accompanied by tuberculosis [1,2]. Vertebral body destruction will lead to spinal instability, spinal cord damage, and deformity. erefore, in addition to standardized antituberculosis treatment, surgery is also an important treatment method for thoracic tuberculosis treatment [3].
As an important surgical intervention in the treatment of spinal tuberculosis, surgical treatment aims to effectively remove the lesion, completely relieve the neurospinal compression, restore the height of the affected vertebra as far as possible, and rebuild the stability of the spine. After continuous exploration and development, the current surgical treatment methods are classified according to the approach, mainly including anterior approach, posterior approach, and combined approach. (1) Anterior approach: compared with other methods, anterior approach is the most convenient method for removing lesions because it can directly reach the lesion and has a good operating field of vision [4]. (2) Combined approach: both lesion clearance and deformity correction. However, anterior or anterior combined approaches have a higher rate of primary activity and mortality [5,6]. (3) Posterior approach: posterior approaches such as posterolateral, pedicle, or posterior approaches and their variants have been proven to be a safe and effective method for the treatment of thoracic disc herniation (TDH) with a relatively low incidence. In addition, posterior approach can better correct kyphosis, which is the most common complication of TST [7]. Some studies believed that, in terms of operation time, the combined approach was longer than the anterior approach, but there was no statistical difference with the posterior approach. In terms of blood loss, anterior and posterior combined approaches were more common than anterior or posterior approaches. ere was no clinical literature to report the influencing factors of posterior modified transarticular debridement for thoracic tuberculosis. is paper considers the influencing factors of operative time, blood loss, bone fusion, and the classification of the KirkaldyWillis.

Methods
Between 2009 and 2014, the authors treated 37 consecutive patients with TST via a modified transfacet approach [7].
is study was approved by the Ethics Committee of Affiliated Hospital of Fujian Medical University.

Preoperative Surgery.
Patients were treated with standard chemotherapy regimen of isoniazid (H), rifampicin (R), ethambutol (E), and pyrazinamide (Z) (HREZ) at least 4 weeks before surgery, including isoniazid (300 mg/d), rifampicin (450 mg/d), ethambutol (750 mg/d), and pyrazinamide (750 mg/d). ESR was 44.7 ± 23.3 mm/h. And the mean kyphosis angle was 29.4 ± 10.9°. Surgery is performed when ESR is significantly reduced (<40 mm/h). Preoperative antituberculosis treatment can reduce mycobacterium tuberculosis in focus and improve surgical safety. Bony spinal fusion was assessed according to the criteria defined by Lee et al. [8], which was graded as shown in Table 1. And in our study, all patients obtained definitive bone fusion.

Statistical Principle.
In order to examine the factors that might affect operative time, the multivariable linear regression model is adopted. Operative time is considered as a continuous and dependent variable, whereas blood loss, PreK, FFUK, and focal debridement are constituted the independent variables. e multivariable linear regression model is also utilized for investigating whether operative time could affect blood loss. Blood loss is considered as a constant and dependent variable, whereas operative time, bone fusion, and PreK constituted the independent variables. e naive Bayes classifier model is also utilized for investigating what affects the KirkaldyWillis grade. All the statistical analyses were performed using the R Project (version 3.5.1, 2018, e R Foundation for Statistical Computing), the R package e1071, corrplot needs to be loaded. e results were considered statistically significant at p < 0.05.
Multiple Linear Regression [9] attempts to find a linear equation with known data to describe the relationship between two or more independent variables and dependent variable and use this linear equation to predict the results.
Let the dependent variable be Y, and the k independent variables affecting the dependent variable be X 1 , X 2 , . . . , X K , respectively. Let us assume that the influence of each independent variable on the dependent variable Y is linear, that is, under the condition that other independent variables remain unchanged, the mean value of Y changes uniformly with the change of independent variable X i : is is called the population regression model, and are the regression parameters. e principle of naive Bayes classifier [10] is a statistical method. Bayes theorem is the theoretical basis of classifier modeling. It uses conditional probability principle, prior information, and sample data information to determine the probability of event occurrence. In addition, naive Bayes model is a kind of generating model, which directly models the joint probability to obtain the target probability. e algorithm is as follows.

Result
We find that blood loss (positive correlation) and focal debridement (positive correlation) could affect operative time, but not bone fusion, gender, PreA, As for the influencing factor model of operative time, F test is performed on the whole model (Table 2), its F-statistic is 57.76 on 2 and 34 DF, and p value is 1.161e−11, indicating that the model is significant on the whole, that is, the model is significant highly. e goodness of fit R 2 is 0.7726, indicates that the model has a good fitting effect and can better explain the influencing factors of operative time. In patients with the same T_FocalDebridement, operative time increased by an average of 0.11575 if blood loss increases 1 ml, its p value is 0.000168, which indicates that at the significance level of 0.001, and blood loss has a significant effect on the operative time. Also, operative time increased by an average of 23.61903 as T_FocalDebridement increased by 1 segment when the blood loss is the same, its p value is 0.035237, which indicates that at the significance level of 0.05, and T_FocalDebridement factor has a significant effect on the operative time.
As for the influencing factor model of blood loss, F test is performed on the whole model (Table 3), and its F-statistic is 99.85 on 1 and 35 DF, and p value is 8.648e−12, indicating that the model is significant highly on the whole, that is, the model is significant. e goodness of fit R 2 is 0.7404, indicates that the model has a good fitting effect and can better explain the influencing factors of blood loss. Blood loss increased by an average of 4.4878 if operative time increased 1 min, its p value is 8.65e−12, which indicates that at the significance level of 0.001, and the operative time factor has a significant effect on the blood loss.

(7)
As for the influencing factor model of bone fusion, F test is performed on the whole model (Table 4), its F-statistic is 13.12 on 3 and 31 DF, and p value is 8.342e−06, indicating that the model is significant on the whole, that is, the model is significant highly. e goodness of fit R 2 is 0.5026, indicates that the model has a good fitting effect and can better explain the influencing factors of bone fusion. In patients with the same PostE and T_FocalDebridement, bone fusion increased by an average of 0.04326 if age increases 1, its p value is 0.00286, which indicates that at the significance level of 0.001, age has a significant effect on the bone fusion. Also, bone fusion decreased by an average of 0.08810 as PostE increased by 1 when the age and T_FocalDebridement are the same, its p value is 0.03415, which indicates that at the significance level of 0.05, PostE factor had a significant effect on the bone fusion. Also, bone fusion increased by an average of 1.53847 as T_FocalDebridement increased by 1 when the age and PostE are the same, its p value is 3.99e−06 indicates that at the significance level of 0.001, and T_FocalDebridement factor had a significant effect on the bone fusion.
After deleting one sample in which KirkaldyWillis is two, only the samples of KirkaldyWillis � 3 and KirkaldyWillis � 4 are trained.
e correlation coefficient matrices of gender, bone fusion, follow-up, FFUK, FFUA, and FFUE are calculated, as shown in Table 5.
As shown in Figure 1, we can find that the interior color of each circle is lighter except the diagonal line, indicating that the correlation coefficient between variables is low and the degree of correlation is not high. e above variables can be used for naive Bayes classification modeling. e confusion matrix is obtained in Table 6. cientific Programming Because naive Bayes model is a kind of generating model, it directly models the joint probability to obtain the target probability. rough modeling, we can only get the classification probability of each sample, so as to determine the category, that is, no explicit model can be obtained. e classification accuracy is 86.11%, which indicates the percentage of all correctly predicted samples to all samples; recall is 96.6%, which indicates the correct percentage of all samples predicted with KirkaldyWillis � 4; precision is 87.5%, which indicates the proportion of all KirkaldyWillis samples predicted to be KirkaldyWillis � 4.

Discussion
Operative data of 37 patients with TST were recorded with multiple indicators of the operative level. Multiple linear regression is also suitable for analyzing the influence of multiple surgical factors on one operative variable. e number of vertebral segments debridement is the most important factor affecting the operation time. When other surgical factors remain unchanged, the greater the amount of blood loss, the longer the operation time. Bone fusion time is the most important factor affecting blood loss. When other surgical factors remain unchanged, the longer the operation time, the longer the bone fusion time or the greater the preoperative the KirkaldyWillis score, the greater the blood loss. And both models have high R 2 values. e posterior transfacet-modified approach has been reported to be a safe and effective method for the treatment of thoracic spinal tuberculosis (TST) with low morbidity in the precious report [7]. However, the impact factors of the surgical approach for TST have not been intensively analyzed. Multiple linear regressions are indeed applicable for analyzing the impact of multiple factors on one operative variable. e aim    e Naive Bayes classification model could be used because the correlation coefficient between variables is low and the degree of correlation is not high. e classification accuracy, recall, and precision for KirkaldyWillis scores are high. erefore, the Naive Bayes classification model is reasonable to classify the KirkaldyWillis scores.

Data Availability
e data used to support the findings of the study are available from the corresponding author upon request.

Ethical Approval
Between 2009 and 2014, the authors treated 37 consecutive patients with TST via a modified transfacet approach [7].
is study was approved by the Ethics Committee of Affiliated Hospital of Fujian Medical University.

Conflicts of Interest
e authors declare no conflicts of interest.