Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 10 March 2022
Sec. Stroke
This article is part of the Research Topic Leveraging Machine and Deep Learning Technologies for Clinical Applications in Stroke Imaging View all 7 articles

A Comparison of LASSO Regression and Tree-Based Models for Delayed Cerebral Ischemia in Elderly Patients With Subarachnoid Hemorrhage

  • 1Department of Neurosurgery, Renmin Hospital of Wuhan University, Wuhan, China
  • 2Department of Neurosurgery, Affiliated Hospital of Panzhihua University, Panzhihua, China
  • 3Department of Neurosurgery, Huzhou Central Hospital, Huzhou, China
  • 4Department of Neurosurgery, First Hospital of Shanxi Medical University, Taiyuan, China
  • 5Department of Neurosurgery, General Hospital of Northern Theater Command, Shenyang, China
  • 6Department of Neurosurgery, The Second Affiliated Hospital of Nanchang University, Nanchang, China

Backgrounds: As a most widely used machine learning method, tree-based algorithms have not been applied to predict delayed cerebral ischemia (DCI) in elderly patients with aneurysmal subarachnoid hemorrhage (aSAH). Hence, this study aims to develop the conventional regression and tree-based models and determine which model has better prediction performance for DCI development in hospitalized elderly patients after aSAH.

Methods: This was a multicenter, retrospective, observational cohort study analyzing elderly patients with aSAH aged 60 years and older. We randomly divided the multicentral data into model training and validation cohort in a ratio of 70–30%. One conventional regression and tree-based model, such as least absolute shrinkage and selection operator (LASSO), decision tree (DT), random forest (RF), and eXtreme Gradient Boosting (XGBoost), was developed. Accuracy, sensitivity, specificity, area under the precision-recall curve (AUC-PR), and area under the receiver operating characteristic curve (AUC-ROC) with 95% CI were employed to evaluate the model prediction performance. A DeLong test was conducted to calculate the statistical differences among models. Finally, we figured the importance weight of each feature to visualize the contribution on DCI.

Results: There were 111 and 42 patients in the model training and validation cohorts, and 53 cases developed DCI. According to AUC-ROC value in the model internal validation, DT of 0.836 (95% CI: 0.747–0.926, p = 0.15), RF of 1 (95% CI: 1–1, p < 0.05), and XGBoost of 0.931 (95% CI: 0.885–0.978, p = 0.01) outperformed LASSO of 0.793 (95% CI: 0.692–0.893). However, the LASSO scored a highest AUC-ROC value of 0.894 (95% CI: 0.8–0.989) than DT of 0.764 (95% CI: 0.6–0.928, p = 0.05), RF of 0.821 (95% CI: 0.683–0.959, p = 0.27), and XGBoost of 0.865 (95% CI: 0.751–0.979, p = 0.69) in independent external validation. Moreover, the LASSO had a highest AUC-PR value of 0.681 than DT of 0.615, RF of 0.667, and XGBoost of 0.622 in external validation. In addition, we found that CT values of subarachnoid clots, aneurysm therapy, and white blood cell counts were the most important features for DCI in elderly patients with aSAH.

Conclusions: The LASSO had a superior prediction power than tree-based models in external validation. As a result, we recommend the conventional LASSO regression model to predict DCI in elderly patients with aSAH.

Introduction

Subarachnoid hemorrhage (SAH) secondary to the ruptured aneurysm is a potentially fatal cerebrovascular disease that mainly occurs in middle-aged patients <60 years (1, 2). However, the number of elderly patients with SAH has been increasing due to improved life expectancy (3). It was reported that the annual incidence of SAH in persons over 70 years of age was estimated to exceed 25/100,000 (4). Delayed cerebral ischemia (DCI) is the most frequent complication after SAH, affecting approximately 30% of patients and often leading to poor neurology outcomes or deterioration of patients' conditions (5, 6). Nevertheless, the prognostic effect maybe worse when DCI occurs in elderly patients during hospitalization (7). The timely accurate prediction of DCI development is critical for the clinical management and prognosis of elderly patients with SAH; hence, a reliable, precise prediction model for early identifying DCI is urgently needed.

The conventional logistic regression (LR) is still the primary method for developing prediction models for clinical disease. Such as the previous studies revealed that independent risk factors were identified via LR to further construct models for predicting DCI in patients with SAH (812). Yet, the conventional LR method could not fully utilize the clinical data during the developing model process, may contributing to a low prediction power. Machine learning (ML) as a domain of artificial intelligence can solve this limitation, and recent research showed that ML algorithms outperformed traditional statistic modeling approaches (1315). Meanwhile, tree-based methods have been considered one of the best and most extensively used statistical ML methods for analyzing the complex clinical data. Tree-based models produce high accuracy and ease of interpretation of results (16). For instance, predicting long-term prognostic outcomes after SAH (1719), mortality analysis after SAH (20), and utility analysis of management strategies after SAH (21). However, after carefully reviewing the literature, we did not find any research using tree-based methods to predict DCI development in the elderly patient population after aneurysmal SAH (aSAH).

Therefore, the purpose of this study is to develop conventional regression and tree-based models and compare which model had better prediction performance for the DCI development in hospitalized elderly patients after aSAH.

Materials and Methods

Study Design and Patient Enrollment

This was a multicenter, retrospective, observational cohort study that utilized admission clinical information from the electronic health record system. This study participant consisted of all elderly patients with aSAH within 24 h of onset who were admitted in the department of neurosurgery of several medical centers from April 2019 to June 2021, such as Renmin Hospital of Wuhan University, Huizhou Third People's Hospital, Affiliated Hospital of Panzhihua University, First Hospital of Shanxi Medical University, and General Hospital of Northern Theater Command. The elderly patients were defined as those aged 60 years and older. Out of all consecutive 215 patients, 153 eligible elderly patients with aSAH were eventually enrolled in our study. Figure 1 was a flowchart that showed exclusion details and the procedure of this study. Head CT, head and neck CT angiography, or intracerebral digital subtraction angiography was used for the diagnosis of aSAH.

FIGURE 1
www.frontiersin.org

Figure 1. The flowchart that showed exclusion details and the procedure of this study.

The inclusion criteria were as follows: elderly patients aged over 60 years, admission within 24 h after onset, spontaneous SAH caused by aneurysm, head CT scan and blood laboratory tests within 24 h after admission, surgical treatment within 3 days after onset, and DCI after SAH occurred within 4–30 days.

The exclusion criteria included aSAH patients complicated with vascular malformation or intracerebral hemorrhage, postoperative state on admission, cases complicated by acute infection, permanent brain injuries or bilateral mydriasis, nonsurgical treatment, and larger missing data.

Clinical Information Collection

The patient clinical information that included sex, age, past medical history (hypertension, diabetes mellitus, coronary heart disease, smoking, and alcohol consumption), and admission state (World Federation of Neurosurgical Societies [WFNS], Hunt and Hess grade [HH], and modified Fisher scale [mFS]) was collected. In addition, aneurysmal details that included aneurysm location, number, length size, neck size, and surgical treatment method were recorded. Admission blood laboratory tests (D-dimer, glucose, white blood cell [WBC], neutrophil, lymphocyte, and monocyte counts) and CT value of subarachnoid clots and cerebral edema were also utilized in this study. The CT value evaluation method is provided in the Methods in the Data Supplement.

All hospitalized patients received standardized postoperative treatment based on the SAH guidelines (22), such as nimodipine to prevent vasospasms, anti-inflammatory drugs, hemostasis, and analgesic. A postoperative head CT scan was performed to determine the presence of intracranial rebleeding or cerebral infarction.

Delayed Cerebral Ischemia Definitions

The definition of DCI is consistent with Vergouwen et al. (23). (1) No other cause could have led to the occurrence of a permanent or temporary focal neurological impairment (such as aphasia, apraxia, hemianopia, or neglect) between 4 and 14 days after aSAH; (2) the Glasgow Coma Scale score decreased by at least two points (either on one of its components [eye opening, verbal response, and motor response] or on total score); and (3) head CT scans revealed a low-density area that was not noticeable on admission or immediately after the operation, and there were no other causes except vasospasms between 4 and 30 days after aSAH.

Sample Size Evaluation

Events per variable with a value of 10 was used to determine the effective sample size in our study (24). A total of six variables were entered into a multivariable regression model in our preliminary analysis. Hence, there should be at least 60 patients with DCI. In addition, according to the incidence rate of 30% of DCI occurrence after SAH worldwide, at least 200 patients should be enrolled in the model training cohort. Based on the limited effective sample size, the least absolute shrinkage and selection operator (LASSO) regression analysis was used to develop a conventional regression model.

Missing Data Processing

A total of five elderly patients had missing data, which accounted for <5% of the patient population. Therefore, a direct deletion method was applied to process the data (25).

Prediction Model Development

In this study, each patients with aSAH in the dataset was regarded as a single data point, clinical information measured at admission (demographic data, past medical history, WFNS grade, HH grade, mFS, aneurysm information, treatment methods, serum laboratory test, and image CT value) was used as feature input, and DCI occurred was used as the label of the algorithm. We randomly divided the multicentral data into the model development cohort and model validation cohort in a ratio of 70–30%. The training cohort of 111 patients was used to construct the conventional LASSO regression model and tree models, such as decision tree (DT), random forest (RF), eXtreme Gradient Boosting (XGBoost). We used grid search to find optimal parameters. Since the computational resource is limited, only some critical parameters are taken into account for each model. The searching range and steps of the chosen parameters for all investigated modes are listed below.

The LASSO Model

The LASSO regression, suitable for small sample size and high-dimensional data, was used to select the most informative prediction variables to construct the model. We used the “glmnet, corrplot, caret” packages of R and five-fold cross-validation to obtain the optimal λ and the variables selecting results.

The DT Model

Decision tree algorithms partition the sample data by splitting prediction features at discrete cut-points and are usually presented in the form of a tree. The DT algorithm uses the Gini index to determine each split's optimal variable and location in this study. The cost complexity parameter that penalizes more complex trees is used to control the size of the final tree. Several important parameters, such as max_depth, min_samples_split, min_samples_leaf, and max_leaf_nodes were adjusted by the 10-fold cross-validation and “rpart, partykit, caret” packages of R.

The RF Model

Random forest generates multiple DTs by sampling objects and variables, and then classifies the objects in turn to build a predictive model. Finally, to summarize the classification results of each DT, the mode category in all prediction categories is the category of the RF model prediction object. The important parameters, such as n-estimators, min_sample split, and min_sample_leaf, were determined using the 10-fold cross-validation and “randomForest” package.

An XGBoostmodel

The XGBoost is an optimized distributed gradient enhancement library whose design is efficient, flexible, and portable. An ML algorithm is implemented under the framework of gradient enhancement. XGBoost provides the promotion of parallel trees, such as gradient boosting decision trees, which can solve many data science problems quickly and accurately. The important parameters, such as gamma, subsample, nrounds, max_depth, eta, colsample_bytree, and min_child_weight were evaluated by “xgboost” package and a 10-fold cross-validation.

Evaluation of Prediction Model Performances

The area under the receiver operating characteristic (ROC) curve (AUC) with 95% CIs, a precision-recall curve, accuracy, sensitivity, and specificity were used to evaluate the model performance. Additionally, we calculated precision and recall indicators using a validation cohort. We used the optimalCutoff function to obtain the optimal threshold of the model outputs to evaluate the model performance. To better demonstrate the generalization of the above-mentioned models, we calculated those indices on both model training and validation cohorts. Furthermore, we compared the errors of the two cohorts to assess the considered models. Finally, to visualize the contribution of each clinical feature, feature importance calculated via the XGBoost method was generated to rank their relative influence on the risk of DCI development.

Statistical Analysis

The Kolmogorov–Smirnov test was used to determine the distribution type of the data before formally analyzing the dataset. Continuous variables analyzed using the Mann–Whitney U-test, or independent t-test, is presented as a median with interquartile range (IQR) or mean ± SD. Categorical variables analyzed using the chi-square test, or Fisher's exact test are expressed as numbers (percentages). All statistical analyses were two-tailed, and the values of p lower than 0.05 were considered statistically significant. All statistical analyses in this study were completed using IBM SPSS Statistics for Windows, version 26.0 (IBM Corp., Armonk NY, USA) and R software (https://www.r-project.org/).

Results

Baseline Characteristics

The mean age of elderly patients in the model training and validation cohorts was 67 years (IQR: 63, 71) and 66 years (IQR: 63, 69), respectively. We observed that elderly patients with aSAH were more likely to be women, and there were no significant distinctions in past medical history among the two cohorts. In addition, the admission state, aneurysmal information, admission laboratory results, CT value in subarachnoid clots, and cerebral edema had no significant differences between the two groups. The number of patients with DCI in the two groups was 31 (28%) and 11 (26%). Table 1 shows the baseline characteristics in training and validation cohorts. Moreover, we analyzed the baseline information of elderly patients with or without DCI in the training cohorts. Details are placed in Table 2.

TABLE 1
www.frontiersin.org

Table 1. Baseline characteristics of the elderly patients in model training and validation cohorts.

TABLE 2
www.frontiersin.org

Table 2. Baseline characteristics of model training cohort based on delayed cerebral ischemia.

LASSO and Tree Models Development

The training process and optimal parameters of the LASSO and tree-based models are demonstrated in Figure 2. In the regression model, we used the LASSO method to select the optimal predictors. An optimal λ of 0.1356784 and log (λ) of −1.997 were adopted in LASSO, and Figure 2A demonstrates that 23 features finally decreased to two features when using the above parameters. The independent predictors were CT value in subarachnoid clots (adjusted odds ration [aOR]: 1.115, 95% CI: 1.028–1.220, p = 0.011) and aneurysm treatment method (aOR: 0.196, 95% CI: 0.067–0.522, p = 0.001) after the multivariable regression analysis. Min_samples_split, min_samples_leaf, and max_leaf_nodes were set to 2, 2, 0, respectively, and Figure 2B shows that the optimal decision nodes were CT values of subarachnoid clot and WBC count during the training process of DT. When n_estimators, min_sample_leaf, and min_sample_split indiceswere set to 63, 4, 2, respectively, and Figure 2C demonstrates that the minimum error was 0.05 corresponding to the optimal tree number of 63 during the training process of RF. Figure 2D displays the training process of XGBoost, and we can obtain the best prediction power when gamma of 0.25, subsample of 0.5, nrounds of 100, max_depth of 2, eta of 0.01, colsample_bytree of 1, and min_child_weight of 1. In addition, the optimal thresholds of LASSO, DT, RF, and XGBoost were 0.3, 0.13, 0.48, and 0.43, respectively.

FIGURE 2
www.frontiersin.org

Figure 2. The training process and optimal parameters of the LASSO and tree ML models were demonstrated. (A) Demonstrates that 23 features finally decreased to two features when using an optimal λ of 0.1356784 and log(λ) of −1.997 parameters. (B) Shows that the optimal decision nodes were CT value of subarachnoid clot and WBC count during the training process of DT. The minimum error was obtained when the optima tree number was 63, the training process of RF is shown in (C). (D) Displays the training process of XGBoost, and we can obtain the best prediction power when gamma of 0.25, maximum depth of 2, and nrounds value of 100.

LASSO and Tree Models Performance and Validation

When using the training cohort to evaluate the model performance, the LASSO model had a lowest AUC-ROC value of 0.793 (0.692, 0.893) than the single DT of 0.836 (95% CI: 0.747–0.926, p = 0.15), RF of 1 (95% CI: 1–1, p < 0.05), and XGBoost of 0.931 (95% CI: 0.885–0.978, p = 0.01). Moreover, the accuracy of 80.9% of the LASSO was lower than the RF of 85.7% and the XGBoost of 83.3%. However, the LASSO scored a highest AUC value of 0.894 (95% CI: 0.8–0.989) in external verification than DT of 0.764 (95% CI: 0.6–0.928, p = 0.05), RF of 0.821 (95% CI: 0.683–0.959, p = 0.27), and XGBoost of 0.865 (95% CI: 0.751–0.979, p = 0.69). Figure 3 shows the performance and evaluation of the LASSO regression model and tree ML models.

FIGURE 3
www.frontiersin.org

Figure 3. The performance and evaluation of the LASSO regression and tree-based models. (A) ROC and AUC value of LASSO and tree-based modes in training cohort; (B) the ROC and AUC value of LASSO and tree-based modes in validation cohort. ROC, receiver operating characteristic curve; AUC, area under the curve; LASSO, least absolute shrinkage and selection operator; DT, decision tree; RF, random forest; XGBoost, extreme gradient boosting.

Table 3 illustrates the accuracy, specificity, sensitivity, precision, and recall indicators of the above models. As we can see, the RF model with an accuracy of 100% (100%) is higher than other models using the training cohort. However, its accuracy value decreased to 85.7% (83.7, 100%) in external validation cohort. On the contrary, the XGBoost model with an accuracy value of 87.4 and 83.3% had a stable performance in the two cohorts, and the DT model with an accuracy value of 81.1 and 78.5% performed the worst among all tree models. In the regression model, the LASSO model's accuracy improved by 3.5% from training to the external validation cohort. Moreover, the LASSO model had a higher precision and recall value of 62 and 90% than tree-based models in external validation. When evaluating the model performance using area under the precision-recall curve (AUC-PR), LASSO model scored a highest AUC-PR value of 0.681 than DT of 0.615, RF of 0.667, and XGB of 0.622 in external validation. Figure 4 shows the P–R curve and AUC value of all models.

TABLE 3
www.frontiersin.org

Table 3. LASSO and tree-based model performance and validation.

FIGURE 4
www.frontiersin.org

Figure 4. Area under the precision-recall curve (AUC-PR) of all prediction models. (A) the AUC-PR of LASSO; (B) the AUC-PR of DT; (C) the AUC-PR of RF; (D) the AUC-PR of XGBoost. LASSO, least absolute shrinkage and selection operator; DT, decision tree; RF, random forest; XGBoost, extreme gradient boosting.

As shown in the Supplementary Table 1, the errors of LASSO, DT, RF, and XGBoost in model training were 18.1, 20.5, 21.6, and 19.8%, respectively, while in the model validation cohort, the error percentage were 23.7, 26.9, 26.1, and 21.2%, respectively.

Feature Importance

The feature importance was scaled so that the sum added up to 1, with a higher importance score indicating a stronger impact on the occurrence of DCI. The most important three features for DCI prediction in elderly patients were CT value pf subarachnoid clots (0.239), aneurysm therapy methods (0.184), and admission WBC counts (0.132). Figure 5 shows all admission clinical feature importance.

FIGURE 5
www.frontiersin.org

Figure 5. All admission clinical feature importance. ClotCT, CT value of subarachnoid blood clot; WBC, white blood cell; LC, lymphocyte; NC, neutrophil; edemaCT, CT value of cerebral edema; HH, Hunt and Hess grade; MC, monocyte.

Discussion

In this study, we enrolled the eligible elderly patients with aSAH from five medical centers and randomly divided them into model training and external validation cohorts, and discussed whether tree-based models can improve the DCI prediction power compared with the regression model during hospitalization. Due to our limited effective sample size, the LASSO method was applied to construct one conventional regression model, and compared with three tree models. To our knowledge, this study is the first to develop tree-based models using complete admission clinical information and to systematically compare the performances of the LASSO regression model for DCI prediction in elderly patients.

The LASSO regression as a special method performs well when reducing data dimensions and multicollinearity among features. For instance, our previous study (26) used the LASSO method to select three optimal variables for establishing a dynamic nomogram for predicting an unfavorable prognosis after aSAH in the case of limited effective sample size. In this study, the optimal independent risk factor for DCI prediction in elderly patients were CT value of subarachnoid clots and aneurysmal treatment method. Previous studies have shown that a CT value of SAH more than 49.95 HU is correlated with DCI. The CT values in SAH are generally considered to represent the density of subarachnoid clots (27). It can reflect the neural inflammatory response after SAH, while the neurovascular inflammation would be a potential mechanism of early brain injury and delay cerebral vasospasm (28, 29). In our study, a CT value > 61.24 HU would be an independent predictor for DCI in elderly patients with aSAH. In addition, our study suggests that an aneurysm endovascular therapy is a vital factor for DCI prevention in elderly patients compared with the neurosurgical treatment. Montalverne et al. (30) reported that endovascular treatment should be considered as a first option for the ruptured aneurysm in elderly patients since an overall favorable prognosis can be achieved in most persons. Yue et al. (31) considered that an interventional treatment presented a better outcome than the surgical treatment for elderly patients. By fitting the above two variables, the AUC and accuracy index of the LASSO regression model demonstrated a good predictive performance, which is generally better than the tree models.

Tree learning algorithms are the most widely used supervised learning methods in clinical making-decision at present (16). In this research, we, respectively, constructed single tree, RF, and XGBoost models for DCI prediction in elderly patients, and the results indicated that classification power of the DT model was worse than other tree models. Although there were no previous research utilizing DT for predicting DCI in elderly patient population; however, Churpek et al. (32) argued that the DT model was still less capable of predicting ward deterioration than the RF and XGBoost models. We know that single tree models, while easy to construct and interpret, do not have much prediction power necessary for our attempts to solve this particular kind of outcome classification problem (33). As we can see, the RF and XGBoost models both improved the prediction ability for DCI based on a single tree model. On the whole, the XGBoost model was even better and stable when applying to training and validation cohorts. The possible explanation is that the expected result of the XGB technology is to build a series of trees, the latter trees can improve the shortcomings of the previous trees, and ultimately reduce the deviation and classification to achieve the best prediction performance (34).

In the field of predicting DCI in hospitalized elderly patients, no research has been performed comparing conventional regression and tree-based methods. Most previous studies aimed to predict the occurrence of unfavorable prognoses among the population of SAH based on the DT method. For example, a prospective cohort study of negative outcomes after aSAH by Liu et al. (11) compared the conventional regression model to the DT model. They found that the DT model had a similar predictive performance to the regression model since both achieved a high accuracy of 0.895 in the validation dataset. In addition, similar literature has been reported the field (17, 19, 35). In our study, although the prediction power of tree models is generally superior to the LASSO in the training cohort; however, the AUC value of LASSO regression was higher than tree models in external validation. Since the prediction ability of the LASSO model varies greatly between two datasets, the small sample size of the validation set may explain this phenomenon.

To visualize the contribution of each feature to the occurrence of DCI, we also calculated the importance of each feature. We can see the CT value of subarachnoid clots, aneurysm therapy methods, and WBC counts that were the three most critical features for DCI prediction in hospitalized elderly patients with aSAH. Why the first two features are most important in predicting the occurrence of DCI in elderly patients has been explained above. The number of WBCs in the peripheral blood is well-known to directly reflect the level of inflammation in the body. A large prospective observational cohort study by Al-Mufti et al. (36) considered that WBCs of more than 12.1 × 109/L were the most important predictor for DCI prediction in patients with good-grade after aSAH. In our study, the WBCs >12.8 × 109/L in the peripheral blood, similar to the previous study's result, was deemed as an important factor for DCI prediction in hospitalized elderly patients. Ruptured aneurysms events in elderly patients often result in poor-grade at admission. Clinical signs of the early pro-inflammatory cytokine cascade caused by aSAH are almost double in poor-grade patients (36, 37). This may interpret the difference in WBC counts between our results and the previous study.

Our research has several points of clinical value. For elderly patients with aSAH as a particular cohort, there is currently no literature on the early prediction of DCI in elderly patients during hospitalization. Based on this, we created and compared the regression and tree models to predict the DCI performance of elderly patients with aSAH. Although DT, RF, and XGBoost had a better prediction performance than the LASSO regression in the training cohort. However, the LASSO model demonstrated a superior generalization ability than all tree-based models in external validation cohort. These results need to be further validated in the future. Second, we found the three most important clinical predictive features: CT value of subarachnoid clots, WBCs in the peripheral blood, and aneurysmal therapy method.

However, the limitation of this study deserves attention. The object of this research was a special elderly population with aSAH, and the time of primary SAH must be guaranteed within 24 h. This has led to the fact that although we have carried out a multi-center study, the sample size was relatively limited. In the future, large-sample prospective clinical studies are still needed to verify our results. Second, the CT value of the subarachnoid blood clot is measured by manually drawing on the ROI, and we should pay attention to the measurement error. Therefore, after agreeing with the radiologist, the CT value was obtained by two clinicians without knowing the patient's clinical information to reduce errors.

Conclusions

For the early prediction of DCI in elderly patients with aSAH, the LASSO model had a superior prediction power than tree-based models in external validation. As a result, we recommend the conventional LASSO regression model to predict DCI in elderly patients with aSAH. However, these results need to be further validated in the future.

Data Availability Statement

The raw data the supporting the findings of this study are available from the corresponding author upon reasonable request.

Ethics Statement

The studies involving human participants were reviewed and approved by the Medical Ethics Committees of five Medical Centers, including Renmin Hospital of Wuhan University (Principal Affiliation Center; WDRM2021-K022), First Hospital of Shanxi Medical University (2021-Y6), Affiliated Hospital of Panzhihua University (202105002), Huzhou Central Hospital (202108005-01), and the General Hospital of Northern Theater Command (Y2021060). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

PH and QC: study design. PH, YLi, TY, and QS: literature search. PH, HZ, YLiu, GG, ZS, and XG: data acquisition. PH, YQ, LY, YX, YLiu, and JC: data analysis and statistical analysis. PH, YLi, TY, GD, YLiu, and QC: manuscript preparation, editing, and review. All authors read and approved the final manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 82001311).

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2022.791547/full#supplementary-material

References

1. van Gijn J, Rinkel GJ. Subarachnoid haemorrhage: diagnosis, causes and management. Brain J Neurol. (2001) 124:249–78. doi: 10.1093/brain/124.2.249

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Ziu E, Mesfin FB. Subarachnoid Hemorrhage. Treasure Island, FL: StatPearls, StatPearls Publishing (2021).

Google Scholar

3. Kaminogo M, Yonekura M. Trends in subarachnoid haemorrhage in elderly persons from Nagasaki, Japan: analysis of the Nagasaki SAH Data Bank for cerebral aneurysm, 1989-1998. Acta Neurochirurgica. (2002) 144:1133–8; discussion 1138-9. doi: 10.1007/s00701-002-1026-2

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Rooij NKde, Linn FH, van der Plas JA, Algra A, Rinkel GJ. Incidence of subarachnoid haemorrhage: a systematic review with emphasis on region, age, gender and time trends. J Neurol Neurosurg Psychiatry. (2007) 78:1365–72. doi: 10.1136/jnnp.2007.117655

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Macdonald RL. Delayed neurological deterioration after subarachnoid haemorrhage. Nature reviews. Neurology. (2014) 10:44–58. doi: 10.1038/nrneurol.2013.246

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Francoeur CL, Mayer SA. Management of delayed cerebral ischemia after subarachnoid hemorrhage. Critical Care (London, England). (2016) 20:277. doi: 10.1186/s13054-016-1447-6

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Darkwah Oppong M, Iannaccone A, Gembruch O, Pierscianek D, Chihi M, Dammann P., et al. Vasospasm-related complications after subarachnoid hemorrhage: the role of patients' age and sex. Acta Neurochirurgica. (2018) 160:1393–400. doi: 10.1007/s00701-018-3549-1

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Duan W, Pan Y, Wang C, Wang Y, Zhao X, Wang Y, Liu L. Risk factors and clinical impact of delayed cerebral ischemia after aneurysmal subarachnoid hemorrhage: analysis from the china national stroke registry. Neuroepidemiology. (2018) 50:128–36. doi: 10.1159/000487325

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Al-Mufti F, Amuluru K, Smith B, Damodara N, El-Ghanem M, Singh IP., et al. Emerging markers of early brain injury and delayed cerebral ischemia in aneurysmal subarachnoid hemorrhage. World Neurosurg. (2017) 107:148–59. doi: 10.1016/j.wneu.2017.07.114

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Fang YJ, Mei SH, Lu JN, Chen YK, Chai ZH, Dong X., et al. New risk score of the early period after spontaneous subarachnoid hemorrhage: For the prediction of delayed cerebral ischemia. CNS Neurosci Ther. (2019) 25:1173–81. doi: 10.1111/cns.13202

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Liu J, Xiong Y, Zhong M, Yang Y, Guo X, Tan X., et al. Predicting long-term outcomes after poor-grade aneurysmal subarachnoid hemorrhage using decision tree modeling. Neurosurgery. (2020) 87:523–9. doi: 10.1093/neuros/nyaa052

PubMed Abstract | CrossRef Full Text | Google Scholar

12. de Rooij NK, Greving JP, Rinkel GJ, Frijns CJ. Early prediction of delayed cerebral ischemia after subarachnoid hemorrhage: development and validation of a practical risk chart. Stroke. (2013) 44:1288–94. doi: 10.1161/STROKEAHA.113.001125

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Beam AL, Kohane IS. Big data and machine learning in health care. JAMA. (2018) 319:1317–8. doi: 10.1001/jama.2017.18391

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chen JH, Asch SM. Machine learning and prediction in medicine—beyond the peak of inflated expectations. N Engl J Med. (2017) 376:2507–9. doi: 10.1056/NEJMp1702071

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Goldstein BA, Navar AM, Carter RE. Moving beyond regression techniques in cardiovascular risk prediction: applying machine learning to address analytic challenges. Eur Heart J. (2017) 38:1805–14.

PubMed Abstract | Google Scholar

16. Banerjee M, Reynolds E, Andersson HB, Nallamothu BK. Tree-based analysis. Circ Cardiovasc Qual Outcomes. (2019) 12:e004879. doi: 10.1161/CIRCOUTCOMES.118.004879

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Toledo Pde, Rios PM, Ledezma A, Sanchis A, Alen JF. A Lagares Predicting the outcome of patients with subarachnoid hemorrhage using machine learning techniques. IEEE Trans Inf Technol Biomed. (2009) 13:794–801. doi: 10.1109/TITB.2009.2020434

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Chen G, Lu M, Shi Z, Xia S, Ren Y, Liu Z., et al. Development and validation of machine learning prediction model based on computed tomography angiography-derived hemodynamics for rupture status of intracranial aneurysms: a Chinese multicenter study. Eur Radiol. (2020) 30:5170–82. doi: 10.1007/s00330-020-06886-7

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hostettler IC, Muroi C, Richter JK, Schmid J, Neidert MC, Seule M., et al. Decision tree analysis in subarachnoid hemorrhage: prediction of outcome parameters during the course of aneurysmal subarachnoid hemorrhage using decision tree analysis. J Neurosurg. (2018) 129:1499–510. doi: 10.3171/2017.7.JNS17677

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Rau CS, Wu SC, Chien PC, Kuo PJ, Chen YC, Hsieh HY, et al. Prediction of mortality in patients with isolated traumatic subarachnoid hemorrhage using a decision tree classifier: a retrospective analysis based on a trauma registry system. Int J Environ Res Public Health. (2017) 14:1420. doi: 10.3390/ijerph14111420

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Wu X, Kalra VB, Durand D, Malhotra A. Utility analysis of management strategies for suspected subarachnoid haemorrhage in patients with thunderclap headache with negative CT result. Emerg Med J. (2016) 33:30–6. doi: 10.1136/emermed-2015-204634

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Connolly ES, Rabinstein AA, Carhuapoma JR, Derdeyn CP, Dion J, Higashida RT., et al. Guidelines for the management of aneurysmal subarachnoid hemorrhage: a guideline for healthcare professionals from the American Heart Association/American Stroke Association. Stroke. (2012) 43:1711–37. doi: 10.1161/STR.0b013e3182587839

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Vergouwen MD, Vermeulen M, van Gijn J, Rinkel GJ, Wijdicks EF, Muizelaar JP., et al. Definition of delayed cerebral ischemia after aneurysmal subarachnoid hemorrhage as an outcome event in clinical trials and observational studies: proposal of a multidisciplinary research group. Stroke. (2010) 41:2391–5. doi: 10.1161/STROKEAHA.110.589275

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Peduzzi P, Concato J, Kemper E, Holford TR, Feinstein AR. A simulation study of the number of events per variable in logistic regression analysis. J Clin Epidemiol. (1996) 49:1373–9. doi: 10.1016/S0895-4356(96)00236-3

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Eekhout I, de Boer RM, Twisk JW, de Vet HC, Heymans MW. Missing data: a systematic review of how they are reported and handled. Epidemiology (Cambridge, Mass). (2012) 23:729–32. doi: 10.1097/EDE.0b013e3182576cdb

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hu P, Xu Y, Liu Y, Li Y, Ye L, Zhang S., et al. An Externally Validated Dynamic Nomogram for Predicting Unfavorable Prognosis in Patients With Aneurysmal Subarachnoid Hemorrhage. Front Neurol. (2021) 12:683051. doi: 10.3389/fneur.2021.683051

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kanazawa T, Takahashi S, Minami Y, Jinzaki M, Toda M, Yoshida K. Early prediction of clinical outcomes in patients with aneurysmal subarachnoid hemorrhage using computed tomography texture analysis. J Clin Neurosci. (2020) 71:144–9. doi: 10.1016/j.jocn.2019.08.098

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Kubo Y, Ogasawara K, Kakino S, Kashimura H, Tomitsuka N, Sugawara A, et al. Serum inflammatory adhesion molecules and high-sensitivity C-reactive protein correlates with delayed ischemic neurologic deficits after subarachnoid hemorrhage. Surg Neurol. (2008) 69:592–6; discussion 596. doi: 10.1016/j.surneu.2008.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Helbok R, Schiefecker AJ, Beer R, Dietmann A, Antunes AP, Sohm F., et al. Early brain injury after aneurysmal subarachnoid hemorrhage: a multimodal neuromonitoring study. Critical Care (London, England). (2015) 19:75. doi: 10.1186/s13054-015-0809-9

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mont'alverne F, Musacchio M, Tolentino V, Riquelme C, Tournade A. Endovascular management for intracranial ruptured aneurysms in elderly patients: outcome and technical aspects. Neuroradiology. (2005) 47:446–57. doi: 10.1007/s00234-005-1345-0

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Yue Q, Liu Y, Leng B, Xu B, Gu Y, Chen L., et al. A prognostic model for early post-treatment outcome of elderly patients with aneurysmal subarachnoid hemorrhage. World Neurosurg. (2016) 95:253–61. doi: 10.1016/j.wneu.2016.08.020

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Churpek MM, Yuen TC, Winslow C, Meltzer DO, Kattan MW, Edelson DP. Multicenter Comparison of Machine Learning Methods and Conventional Regression for Predicting Clinical Deterioration on the Wards. Crit Care Med. (2016) 44:368–74. doi: 10.1097/CCM.0000000000001571

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Zhang Z. Decision tree modeling using R. Ann Transl Med. (2016) 4:275. doi: 10.21037/atm.2016.05.14

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Ogunleye A, Wang QG. XGBoost Model for Chronic Kidney Disease Diagnosis. IEEE/ACM Trans Comput Biol Bioinform. (2020) 17:2131–40. doi: 10.1109/TCBB.2019.2911071

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Otite F, Mink S, Tan CO, Puri A, Zamani AA, Mehregan A., et al. Impaired cerebral autoregulation is associated with vasospasm and delayed cerebral ischemia in subarachnoid hemorrhage. Stroke. (2014) 45:677–82. doi: 10.1161/STROKEAHA.113.002630

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Al-Mufti F, Misiolek KA, Roh D, Alawi A, Bauerschmidt A, Park S., et al. White blood cell count improves prediction of delayed cerebral ischemia following aneurysmal subarachnoid hemorrhage. Neurosurgery. (2019) 84:397–403. doi: 10.1093/neuros/nyy045

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Provencio JJ, Vora N. Subarachnoid hemorrhage and inflammation: bench to bedside and back. Semin Neurol. (2005) 25:435–44. doi: 10.1055/s-2005-923537

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: delayed cerebral ischemia, subarachnoid hemorrhage, aneurysm, LASSO, tree model

Citation: Hu P, Liu Y, Li Y, Guo G, Su Z, Gao X, Chen J, Qi Y, Xu Y, Yan T, Ye L, Sun Q, Deng G, Zhang H and Chen Q (2022) A Comparison of LASSO Regression and Tree-Based Models for Delayed Cerebral Ischemia in Elderly Patients With Subarachnoid Hemorrhage. Front. Neurol. 13:791547. doi: 10.3389/fneur.2022.791547

Received: 08 October 2021; Accepted: 31 January 2022;
Published: 10 March 2022.

Edited by:

Adriano Pinto, University of Minho, Portugal

Reviewed by:

Zhengbing Yan, Wenzhou University, China
Alex Jung, Aalto University, Finland
Lorenzo Camponovo, University of Applied Sciences and Arts of Southern Switzerland (SUPSI), Switzerland

Copyright © 2022 Hu, Liu, Li, Guo, Su, Gao, Chen, Qi, Xu, Yan, Ye, Sun, Deng, Zhang and Chen. 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: Hongbo Zhang, aG9uZ2Jvemhhbmc5OSYjeDAwMDQwOzE2My5jb20=; Qianxue Chen, Y2hlbnF4NjY2JiN4MDAwNDA7d2h1LmVkdS5jbg==

These authors have contributed equally to this work

Disclaimer: 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.