Next Article in Journal
Addition of ROCK Inhibitors Alleviates Prostaglandin-Induced Inhibition of Adipogenesis in 3T3L-1 Spheroids
Next Article in Special Issue
BU-DLNet: Breast Ultrasonography-Based Cancer Detection Using Deep-Learning Network Selection and Feature Optimization
Previous Article in Journal
Zr-Based Metal-Organic Frameworks for Green Biodiesel Synthesis: A Minireview
Previous Article in Special Issue
Estimation of Caenorhabditis Elegans Lifespan Stages Using a Dual-Path Network Combining Biomarkers and Physiological Changes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Toward Evaluating Critical Factors of Extubation Outcome with XCSR-Generated Rules

1
Institute of Computer Science and Engineering, College of Computer Science, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan
2
Division of Respiratory Therapy, Kaohsiung Medical University Hospital, Kaohsiung 80708, Taiwan
3
Division of Pulmonary and Critical Care Medicine, Kaohsiung Medical University Hospital, Kaohsiung 80708, Taiwan
4
School of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
5
Department of Computer Science, College of Computer Science, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan
6
Institute of Biomedical Engineering, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan
7
Department of Internal Medicine, Kaohsiung Municipal Cijin Hospital, Kaohsiung 80544, Taiwan
8
Division of Respiratory Therapy, College of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 21 October 2022 / Revised: 14 November 2022 / Accepted: 15 November 2022 / Published: 17 November 2022

Abstract

:
Predicting the correct timing for extubation is pivotal for critically ill patients with mechanical ventilation support. Evidence suggests that extubation failure occurs in approximately 15–20% of patients, despite their passing of the extubation evaluation, necessitating reintubation. For critically ill patients, reintubation invariably increases mortality risk and medical costs. The numerous parameters that have been proposed for extubation decision-making, which constitute the key predictors of successful extubation, remains unclear. In this study, an extended classifier system capable of processing real-value inputs was proposed to select features of successful extubation. In total, 40 features linked to clinical information and variables acquired during spontaneous breathing trial (SBT) were used as the environmental inputs. According to the number of “don’t care” rules in a population set, Probusage, the probability of the feature not being classified as above rules, can be calculated. A total of 228 subjects’ results showed that Probusage was higher than 90% for minute ventilation at the 1st, 30th, 60th, and 90th minutes; respiratory rate at the 90th minute; and body weight, indicating that the variance in respiratory parameters during an SBT are critical predictors of successful extubation. The present XCSR model is useful to evaluate critical factors of extubation outcomes. Additionally, the current findings suggest that SBT duration should exceed 90 min, and that clinicians should consider the variance in respiratory variables during an SBT before making extubation decisions.

Graphical Abstract

1. Introduction

In the intensive care unit (ICU), intubation with mechanical ventilators necessarily saves the lives of critically ill patients immediately fighting for their lives and gaining extra time to treat and recover. However, clinical reports provide evidence that prolonged mechanical ventilation can lead to the impairment of diaphragmatic or lung function [1,2,3]. Similarly, respiratory failure requiring reintubation occurs in approximately 15–20% of patients after extubation [4]. Delayed extubation or reintubation tends to increase mortality risk in critically ill patients. Therefore, it is critical to accurately determine whether patients still require mechanical ventilation support. The timing of weaning from mechanical ventilation also constitutes a crucial concern.
Patients who are intubated and undergoing machine ventilation support require frequent physician visits and checkups. After the stabilization of original cause, there is a procedure for determining the patient’s health status called the ventilator weaning protocol. The first step is to observe the patient’s daily assessment of readiness to wean, and to meet and check the weaning profiles, including the rapid shallow breathing index (RSBI), maximal inspiratory pressure (Pimax), cuff leak test, etc. [5,6,7,8]. The second step was to conduct a spontaneous breathing trial (SBT) with 120 min of low pressure support ventilation. The third step of airway protection and cough function assessment was performed. After all of them have passed, the decision to turn off mechanical ventilation and remove the endotracheal tube is made by a trained and certified physician. The strictly refined process still produces 15~20% extubation failure rate; therefore, many studies have introduced multivariate analysis, artificial intelligence, and other methods to find the appropriate factors. During SBT, respiratory distress may occur [9] and changes in the respiratory variables may also be associated with successful extubation [10]. It is, however, still not possible to completely exclude the possibility of failed extubation. The search for other key variables, the development of novel features, and the establishment of accurate predictive models for extubation based on physician decision rules are the major research and development in this study.
A learning classifier system (LCS), which is a rule-based machine learning method, interacts with the environment to generate an optimal policy, modifying the classifier parameter values according to the environmental feedback to formulate more effective rules. A popular type of LCS is the extended classifier system (XCS), for which numerous applications have been proposed [11,12]. The input data for the original XCS must be encoded into binary format. However, most problems are more easily addressed by encoding with real values. Therefore, Wilson modified the XCS for the adoption of real-value inputs; an XCS that can process such inputs is called an XCSR [13]. One feature of XCS is the “don’t care” bits in the rules. This feature enables users to determine which variables are less important in some situations but more important in other situations. A study applied this feature to identify the critical variables for detecting Internet addiction in patients [14]. Thus, the XCSR is suitable for the present objective: determining and evaluating SBT-related variables are more important for successful extubation than others. The aim of this study was not only to construct a model for predicting extubation outcomes but also to evaluate the critical SBT effect with other variables.

2. Experimental and Computational Details

2.1. Subjects Description

This retrospective study of the electronic medical records (EMR) database has been approved by the Institutional Review Board of Kaohsiung Medical University Chung-Ho Memorial Hospital (the approval number: KUMHIRB-F(1)-20200033). A total of 262 deidentified patient data were obtained from the above hospital. The recruited patient is on a ventilator for the first time and for more than 48 h and has stable vital signs before extubation. Since the literature has indicated that SBT with limited pressure support is more effective than T-tubes [15,16], the data of a 6 cmH2O pressure-supported SBT is collected in this study.
Figure 1 presents the extubation procedure used in the ICU of the internal medicine department. The intensivist decides to stop ventilation support and execute the extubation procedure after each index meting the relevant criteria. Within 48 h of extubation, the patient who experienced respiratory failure or died is defined as the extubation failure group, and the patient who did not is defined as the extubation success group. Descriptions of the variables obtained from the EMR are presented in the following section.

2.2. Variable

For evaluating critical factors, 40 variables for each patient were selected as the input variables. Table 1 lists the selected 19 variables of clinical information, consisting of basic patient information, two severity scales, mechanical ventilator-related information, hemodynamics situation before doing SBT, maximal inspiratory pressure, and disease history (Table 1). Table 2 presents another 21 variables of acquired information when doing SBT, including the hemodynamics and ventilation level at the 1st, 30th, 60th, and 90th minute. Regarding the output variable, the patient in the extubation failure group is encoded as 0 and the one in the success group is encoded as 1.

2.3. XCS and XCSR

The learning classification system (LCS) is a rule-based machine learning method that combines exploration components (evolutionary computation, EC) and learning components (reinforcement learning, RL). The development of LCS has been decades-long, and the most famous descendant in the field of the LCS is the extended classifier system (XCS). The XCS retains the main components of the LCS, but still achieves similar performance or even better. The streamlined architecture makes the XCS a benchmark for future development in the LCS field.
The architecture of the XCS is a standard RL. The XCS can be regarded as an agent interacting with the environment to achieve the desired goal. The environment can be a database to be explored or a problem to be solved. The XCS uses a set of cooperative rules to express the solutions it has learned. The rules are easily interpretable, and the expression of the rules is “IF condition THEN action”. In addition, each rule also contains 3 parameters to evaluate the quality of the rule: (1) prediction, predict the reward after the action is executed; (2) prediction error, calculate the error between predicted and actual reward; (3) fitness, evaluate the quality and reliability of rule. The reward received from the environment is used to update the rule parameters using Q-learning and the Widrow–Hoff delta rule. Introducing EC can help the XCS reorganize and explore better rules. The XCS has also been proven to be able to find the most general and accurate rules.
The operation process of the XCS is described as follows: (1) Environmental input is converted to binary encoding through the detector. The set of rules is called population [P] in the XCS. The condition part of the rule will correspond to the environmental input and use ternary-alphabet {0, 1, #} for encoding. The # symbol is treated as “don’t care”, which means it can match regardless of whether the corresponding environmental input is 0 or 1. The action part of the rule adopts binary encoding. The XCS will search for all the rules in [P] and move the rules whose condition part matches the environmental input to the matching set [M]. If [M] is empty, the XCS will start the covering mechanism. In the covering mechanism, the XCS will generate a new rule that matches the condition part with the environmental input and randomly give a binary code to the action part. In the condition part, there is a certain probability to replace the regular bit with the # symbol. The remaining parameters are given random initial values. The generated rule will be placed in [P] and rescan [P] to generate [M]; (2) The rules in [M] may contain different actions. To calculate the expected reward for each action, the rules in [M] with the same action are used to execute the fitness weighted average as the expected reward of the action. The prediction array [PA] can be built after calculating the expected reward of each action. The XCS has two modes to select action for execution, namely exploration and exploitation. In the exploration mode, the XCS adopts the roulette-wheel method, which means that the larger the expected reward, the easier to be picked. In the exploitation mode, the XCS directly selects the action with the most expected reward to execute. The selected action will be executed through the effector, and the environment will give reward. The rules in [M] whose action is the same as the selected action will be moved to action set [A]; (3) The reward obtained from the environment will be used to update the prediction, prediction error, and fitness parameters of rules in [A] through Q-learning and the Widrow–Hoff delta rule; (4) The parameters of the rule corresponding to [A] in [P] will be updated; (5) The EC method in the XCS is to use genetic algorithm (GA). When the average time of the rule in [A] from the last GA execution is greater than the threshold, GA will execute. GA will use fitness as the probability to select two rules for crossover and mutation to generate two new rules. In the XCS, the alleles of condition in the two rules will be exchanged with probability χ, and the alleles will be randomly flipped and mutated with the probability μ; (6) The XCS also has subsumption mechanism to assist in learning the most general and most accurate rules. The XCS will determine two rules corresponding to the same environmental input, and whether the more general rule has enough experience and reliability to accommodate another specific rule; (7) The XCS will repeat the above steps until it converges or reaches the maximum number of iterations.
The variables in real-world data are often continuous values. A new encoding is proposed in the condition part of the rule to allow the XCS to handle input data other than binary. The modified XCS version is called XCS with real-valued input (XCSR) [13]. The schematic of the XCSR and the corresponding operation procedure are shown in Figure 2 and Figure 3, respectively. The XCSR mainly converts the encoding of the condition of the rules from ternary-alphabet to interval predicates. The expression of interval predicate is inti = (ci, si), where inti represents the i-th interval predicate corresponding to the i-th input variable, ci represents the center, and si represents the spread. When an i-th input variable xi can match the i-th interval predicate inti, it means (ci − si   xi < ci + si). The operating mechanism of XCSR is the same as XCS, except that the mechanism is slightly changed in covering and GA. When covering generates a new rule, it will generate inti corresponding to the xi in the condition of the rule, where the value of ci is set to xi, and the value of si is set to a random number ranging from 0 to sr, where sr is defined by the user. In the part of GA, the inti of the two rules will be exchanged with probability χ. ci and si will randomly add or subtract the value of mi with probability μ, where mi is defined by the user.

2.4. Data Analysis

This study will import clinical data of 228 patients with 40 variables as input data for the XCSR. The rule condition of the XCSR contains 40 interval predicates, and the output action is the success or failure of extubation. Table 3 shows the parameter settings of XCSR used in this study. The detailed description of XCSR can be referred to [13]. After the values of the 40 variables of the 228 patients were input into the XCSR, it generated 154 rules for extubation outcome prediction after the model reached stability. Therefore, according to the number of “don’t care” bits of each variable in the rules, the usage probability Probusage can be calculated. A larger Probusage in the ith variable indicates that more model rules use this variable to determine the extubation outcomes. The Probusage of the ith variable is calculated as follows:
Prob usage ( i ) = t o t a l   r u l e s   n u m b e r d o n t   c a r e   b i t   n u m b e r ( i ) t o t a l   r u l e s   n u m b e r × 100 %
We constructed the model 10 times to verify the stability of the Probusage of each variable. Therefore, the results are presented as means ± standard deviations of Probusage.

2.5. Statistical Analysis

Because the two groups were of unequal size, the nonparametric Mann–Whitney U test was performed to examine the distribution between extubation success and failure. The continuous variables are expressed as means ± standard deviations, and the discrete variables are expressed as numbers. A p value of <0.01 indicated a significant difference between two groups in two-tailed tests. The area under the receiver operating characteristic (ROC) curve (AUC) was used to test the discriminatory ability of each variable, with results expressed as the AUC values. The analyses were conducted and the XCSR was coded using the LabVIEW platform (National Instruments Corporation, Austin, TX, USA).

3. Results and Discussion

After exclusions were made for 34 of the 262 patients because of missing data, the data of 228 patients remained for analysis. Table 4 presents the clinical information of the final sample. Height and ideal body weight (IBW) differed significantly between the two groups. The AUC of these two variables were both 0.63.
Table 5 lists the variables measured during the SBTs of the 228 patients. A significant between-group difference was only observed for the heart rate (HR) in the 60th minute. The AUC of this variable was 0.63.
Figure 4 presents the performance of one of the XCSR models (a total of 10 models were built, one of them for an example). A problem refers to an iteration of a prediction that has been created. For the model to approach 100% accuracy, only approximately 125 problems were required. Though the possibility of over fitting of the XCSR model cannot be excluded, the analysis of the rules created by the XCSR can still be useful information for clinical doctors to make a decision.
Figure 5 presents the “don’t care” numbers of each variable in all rules after sorting. The first six variables had significantly fewer “don’t care” numbers, meaning that they had higher usage probabilities. Table 6 lists both these high-performing variables and some others with lower rankings.
Variables ranked 1–6 had a substantially higher Probusage than did other model variables. The Probusage of all variables ranked ≥7 was lower than 36%.
Although height, IBW, and HR in the 60th minute typically differ significantly between the extubation failure and success groups in relevant studies, the AUC values remain lower than 0.63. Few studies using the XCSR have applied the ROC analysis. However, in clinical research, ROC analysis is the most common method for comparing different models and variables. Therefore, in the present study, the rewards in the prediction array were applied to the XCSR model for the ROC analysis. The AUC of 0.993 was significantly higher than that of individual variables. Moreover, the accuracy of the XCSR model in predicting extubation outcomes approached 100% after training. For comparison with the XCSR model, the accuracy of height, IBW, and HR in the 60th minute was determined as 59%, 65%, and 65%, respectively, by using the cutoff point from the ROC analysis to predict extubation outcomes. The prediction accuracy for individual variables was relatively low because the ratio between the number of extubation failures and successes was approximately 0.83, meaning that the accuracy would be 83% if the model predicted that all extubation outcomes would be successful. Nevertheless, these results demonstrate how powerful the present XCSR model is. In the population rule set, 130 and 24 rules predicted extubation success and failure, respectively. Overall, the rules predicted the extubation outcomes with high accuracy.
After the model was constructed 10 times, although the Probusage rank of some variables differed, the first six positions remained filled by the same variables. Moreover, the Probusage of these six variables was higher than 95%, with small standard deviations. Overall, the results indicate that the respiratory rate (RR) and minute ventilation (VE) measured during the SBTs, as well as body weight (BW), were more important than the other variables. In particular, VE was critical at all-time points during the SBT. Studies have demonstrated that some SBT-related respiratory variables are risk factors for extubation failure [10,17,18,19]. One study noted that VE is a parameter that is easy to measure during SBTs that can guide decisions regarding the discontinuation of mechanical ventilation [17]. Another study reported that including all variables (especially the variables in 90th minutes) assessed during SBTs in a training model resulted in higher performance [20]. The current results also suggest that VE during SBT is a key predictor of successful extubation.
Patients who are able to complete an SBT are likely to have more stable breathing than those who are not. This may explain why respiratory rate in the 90th minute constitutes a major factor in extubation decision-making. Notably, BW can also be a critical factor for successful extubation. Because the rules generated by the XCSR are presented in a range, the values may not be continuous. However, the reason BW is typically used for rules in XCSR models remains unclear. We speculate that patients may feel more comfortable and breathe more stably if their BW is regulated within a certain range. Further research is warranted to confirm this supposition.
Some studies have also applied machine learning to predict extubation success or failure. Zeggwagh used multivariable logistic regression (MLR) and three variables (vital capacity, RSBI, and maximal expiratory pressure) to create a prediction model [21]. Mueller (2004) observed that a model created using an artificial neural network (ANN) outperformed one created using MLR [22] in predicting extubation outcomes in infants. Another study revealed that an ANN model more accurately predicted extubation failure [23]. Mueller (2013) tested five machine learning methods, namely MLR, ANNs, boosted decision trees, naïve Bayesian classifiers, and support vector machines, to predict extubation outcomes in a large sample of premature infants [24]. Two common statistical models, MLR and partial least squares (PLS) regression, were also used to test the same dataset of present study. The AUCs of MLR (0.660) and PLS (0.812) were significantly lower than that of the XCSR model (0.993). Table 7 presents a summary of the AUC results obtained in previous studies and the current study. Because each study applied different database and variables for model training, the summary of Table 7 is just for reference instead of comparison. Interpreting the nonlinear models and determining which variables were more important than others are challenging, especially when ANNs are used. Relevant studies have typically employed the trial-and-error approach to fine-tune models and variables. Though some method had been proposed for NN model interpreting, the XCSR model is easier to interpret. In the present study, each rule created by the XCSR was readable (each variable in a rule would be a range of value or “don’t care”). In general, interpreting an XCSR model is easier than a model formulated using other nonlinear machine learning methods. The usage probabilities calculated in the present study may be useful for identifying the importance of each variable. In addition, the distributions of each variable in different rules also can be useful information for finding the physiological meaning of patients with fail extubation and patients with successful extubation. These rules may serve as a valuable reference for further research.
Recently, a study also proposed a method to visualize rule-based machine learning [25]. In future work, it is possible to apply the method for identifying the importance of each variable to compare with present study and to explore detailed information from the XCSR. In addition, the original XCSR had been proposed for almost two decades. Many modified versions of the XCSR were proposed in these years. To study which one of the modified XCSRs is more suitable for the purpose of this study is also a future work.
The present XCSR system is already built in the system of Kaohsiung Medical University Chung-Ho Memorial Hospital for providing the extubation suggestion. In the future work, the prospective research of this XCSR system for extubation suggestion will be investigated.

4. Conclusions

The XCSR model trained in the present study had a prediction accuracy of almost 100%. We identified critical factors of extubation success and failure from among 40 variables, observing that the respiratory variables assessed during SBTs were more important. The generated rules used over 95% of the VE measured in the 1st, 30th, 60th, and 90th minutes, as well as the respiratory rate in the 90th minute and BW. These parameters, especially VE, are the key to determining extubation success or failure. These findings suggest that the duration of an SBT should be longer than 90 min, and that extubation decisions should be made in consideration of the variability in such respiratory parameters during an SBT. Overall, relative to those of other methods, the AUC of the XCSR model was higher, facilitating the evaluation of the importance of each variable. Furthermore, each generated rule may constitute a useful reference for future investigations.

Author Contributions

Conceptualization, P.-H.H., L.-Y.C., T.-C.H. and J.-R.T.; methodology, P.-H.H. and L.-Y.C.; software, P.-H.H. and L.-Y.C.; validation, P.-H.H., L.-Y.C., T.-C.H. and J.-R.T.; formal analysis, P.-H.H. and L.-Y.C.; investigation, P.-H.H., L.-Y.C. and W.-C.C.; resources, C.-C.S., T.-C.H. and J.-R.T.; data curation, W.-C.C.; writing—original draft preparation, P.-H.H. and L.-Y.C.; writing—review and editing, T.-C.H. and J.-R.T.; visualization, P.-H.H. and L.-Y.C.; supervision, C.-C.S.; project administration, T.-C.H. and J.-R.T.; funding acquisition, T.-C.H. and J.-R.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC were funded by the project “Medical artificial intelligence joint project of Kaohsiung Medical University/National Chiao Tung University” (grant number: NCTUKMU 108-AI-08). This work was also partially supported by the Higher Education Sprout Project of the National Yang Ming Chiao Tung University and Ministry of Education (MOE), Taiwan. This work was supported in part by the National Science and Technology Council (Taiwan) under grant number MOST109-2221-E-009-117-MY3 and MOST111-2221-EA49-011-MY3.

Institutional Review Board Statement

This retrospective study of electronic medical records (EMR) database has been approved by the Institutional Review Board of Kaohsiung Medical University Chung-Ho Memorial Hospital (the approval number: KUMHIRB-F(1)-20200033).

Informed Consent Statement

Patient consent was waived due to only deidentified patient data were obtained.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to participates’ privacy.

Acknowledgments

This work was supported by grant NCTUKMU108-AI-08-2 from Kaohsiung Medical University and National Chiao Tung University, the Higher Education Sprout Project of the National Yang Ming Chiao Tung University and Ministry of Education (MOE), Taiwan, and by the grant number MOST109-2221-E-009-117-MY3 and MOST111-2221-EA49-011-MY3 from National Science and Technology Council, Taiwan.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tobin, M.J. Advances in mechanical ventilation. N. Engl. J. Med. 2001, 344, 1986–1996. [Google Scholar] [CrossRef]
  2. Liao, P.; Whitehead, T.; Evans, T.; Griffiths, M. Ventilator-associated lung injury. Lancet 2003, 361, 332–340. [Google Scholar]
  3. Pham, T.; Brochard, L.J.; Slutsky, A.S. Mechanical ventilation: State of the art. Mayo Clin. Proc. 2017, 92, 1382–1400. [Google Scholar] [CrossRef] [Green Version]
  4. Thille, A.W.; Richard, J.-C.M.; Brochard, L. The decision to extubate in the intensive care unit. Am. J. Respir. Crit. Care Med. 2013, 187, 1294–1302. [Google Scholar] [CrossRef]
  5. Yang, K.L.; Tobin, M.J. A prospective study of indexes predicting the outcome of trials of weaning from mechanical ventilation. N. Engl. J. Med. 1991, 324, 1445–1450. [Google Scholar] [CrossRef]
  6. Bien, U.D.S.; Souza, G.F.; Campos, E.S.; de Carvalho, E.F.; Fernandes, M.G.; Santoro, I.; Costa, D.; Arena, R.; Sampaio, L.M.M. Maximum inspiratory pressure and rapid shallow breathing index as predictors of successful ventilator weaning. J. Phys. Ther. Sci. 2015, 12, 3723–3727. [Google Scholar] [CrossRef] [Green Version]
  7. Nemer, S.N.; Barbas, C.S.; Caldeira, J.B.; Guimarães, B.; Azeredo, L.M.; Gago, R.; Souza, P.C.P. Evaluation of maximal inspiratory pressure, tracheal airway occlusion pressure, and its ratio in the weaning outcome. J. Crit. Care 2009, 24, 441–446. [Google Scholar] [CrossRef]
  8. Conti, G.; Montini, L.; Pennisi, M.A.; Cavaliere, F.; Arcangeli, A.; Bocci, M.G.; Proietti, R.; Antonelli, M. A prospective, blinded evaluation of indexes proposed to predict weaning from mechanical ventilation. Intensive Care Med. 2004, 30, 830–836. [Google Scholar] [CrossRef] [PubMed]
  9. Esteban, A.; Alía, I.; Tobin, M.J.; GIL, A.; Gordo, F.; Vallverdú, I.; Blanch, L.; Bonet, A.; Vázquez, A.; de Pablo, R.; et al. Effect of spontaneous breathing trial duration on outcome of attempts to discontinue mechanical ventilation. Am. J. Respir. Crit. Care Med. 1999, 159, 512–518. [Google Scholar] [CrossRef] [PubMed]
  10. Segal, L.N.; Oei, E.; Oppenheimer, B.W.; Goldring, R.M.; Bustami, R.T.; Ruggiero, S.; Berger, K.I.; Fiel, S.B. Evolution of pattern of breathing during a spontaneous breathing trial predicts successful extubation. Intensive Care Med. 2010, 36, 487–495. [Google Scholar] [CrossRef] [PubMed]
  11. Akbar, M.A.; Farooq, M. Application of evolutionary algorithms in detection of sip based flooding attacks. In Proceedings of the 11th Annual Conference on Genetic and Evolutionary Computation (GECCO’09), Association for Computing Machinery, New York, NY, USA, 8–12 July 2009; pp. 1419–1426. [Google Scholar] [CrossRef]
  12. Passaro, A.; Baronti, F.; Maggini, V. Exploring relationships between genotype and oral cancer development through XCS. In Proceedings of the 7th Annual Workshop on Genetic and Evolutionary Computation (GECCO’05), Association for Computing Machinery, New York, NY, USA, 25–26 June 2005; pp. 147–151. [Google Scholar] [CrossRef]
  13. Wilson, S.W. Get real! XCS with continuous-valued inputs. In Learning Classifier Systems; Lecture Notes in Computer Science; Springer: Berlin/Heidelberg, Germany, 2000; Volume 813. [Google Scholar] [CrossRef]
  14. Ji, H.-M.; Chen, L.-Y.; Hsiao, T.-C. Real-time detection of Internet addiction using reinforcement learning system. In Proceedings of the Genetic and Evolutionary Computation Conference Companion (GECCO’19). Association for Computing Machinery, New York, NY, USA, 13–17 July 2019; pp. 1280–1288. [Google Scholar] [CrossRef]
  15. Pellegrini, J.A.S.; Moraes, R.B.; Maccari, J.G.; de Oliveira, R.P.; Savi, A.; A Ribeiro, R.; Burns, K.E.; Teixeira, C. Spontaneous breathing trials with T-Piece or pressure support ventilation. Respir. Care 2016, 61, 1693–1703. [Google Scholar] [CrossRef] [Green Version]
  16. Ladeira, M.T.; Vital, F.; Andriolo, R.B.; Andriolo, B.N.G.; Atallah, N.; Peccin, M.S. Pressure support versus T-tube for weaning from mechanical ventilation in adults (Review). Cochrane Database Syst. Rev. 2014, 5, 1–58. [Google Scholar] [CrossRef]
  17. Martinez, A.; Seymour, C.; Nam, M. Minute ventilation recovery time a predictor of extubation outcome. Chest 2003, 123, 1214–1221. [Google Scholar] [CrossRef] [PubMed]
  18. Khamiees, M.; Raju, P.; DeGirolamo, A.; Amoateng-Adjepong, Y.; Manthous, C.A. Clinical investigations in critical care predictors of extubation outcome in completed a spontaneous breathing trial. Chest 2001, 120, 1262–1270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Frutos-Vivar, F.; Ferguson, N.D.; Esteban, A.; Epstein, S.K.; Arabi, Y.; Apezteguía, C.; González, M.; Hill, N.S.; Nava, S.; D’Empaire, G.; et al. Risk factors for extubation failure in patients following a successful spontaneous breathing trial. Chest 2006, 130, 1664–1671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Chung, W.; Sheu, C.-C.; Hung, J.; Hsu, T.; Yang, S.; Tsai, J. Novel mechanical ventilator weaning predictive model. Kaohsiung J. Med. Sci. 2020, 36, 841–849. [Google Scholar] [CrossRef] [PubMed]
  21. Zeggwagh, A.A.; Abouqal, R.; Madani, N.; Zekraoui, A.; Kerkeb, O. Weaning from mechanical ventilation: A model for extubation. Intensive Care Med. 1999, 25, 1077–1083. [Google Scholar] [CrossRef] [PubMed]
  22. Mueller, M.; Wagner, C.L.; Annibale, D.J.; Hulsey, T.C.; Knapp, R.G.; Almeida, J.S. Predicting extubation outcome in preterm newborns: A comparison of neural networks with clinical expertise and statistical modeling. Pediatr. Res. 2004, 56, 11–18. [Google Scholar] [CrossRef]
  23. Hsieh, M.-H.; Chen, C.-M.; Hsieh, C.-C.; Chao, C.-M.; Lai, C.-C. An artificial neural network model for predicting successful extubation in intensive care units. J. Clin. Med. 2018, 7, 240. [Google Scholar] [CrossRef] [Green Version]
  24. Mueller, M.; Almeida, J.S.; Stanislaus, R.; Wagner, C.L. Can machine learning methods predict extubation outcome in premature infants as well as clinicians? J. Neonatal Biol. 2013, 2, 1000118. [Google Scholar] [CrossRef] [Green Version]
  25. Liu, Y.; Browne, W.N.; Xue, B. Visualizations for rule-based machine learning. Nat. Comput. 2021, 1–22. [Google Scholar] [CrossRef]
Figure 1. Ventilator weaning protocol in ICU of Kaohsiung Medical University Chung-Ho Memorial Hospital.
Figure 1. Ventilator weaning protocol in ICU of Kaohsiung Medical University Chung-Ho Memorial Hospital.
Bioengineering 09 00701 g001
Figure 2. Schematic of the XCSR.
Figure 2. Schematic of the XCSR.
Bioengineering 09 00701 g002
Figure 3. Flowchart of the XCSR in stage n.
Figure 3. Flowchart of the XCSR in stage n.
Bioengineering 09 00701 g003
Figure 4. XCSR performance.
Figure 4. XCSR performance.
Bioengineering 09 00701 g004
Figure 5. Distribution of the number of rules that use each variable.
Figure 5. Distribution of the number of rules that use each variable.
Bioengineering 09 00701 g005
Table 1. Selected clinical variables.
Table 1. Selected clinical variables.
Variable (Unit)Definition
Basic information
 Gender (a.u.)Subject gender (1: male; 2: female)
 Age (year)Subject age
 Height (cm)Height
 BW (kg)Body weight
 IBW (kg)Ideal body weight
Severity scale
 APACHEII (a.u.)Acute physiology and chronic health evaluation II
 GCS_enter_ICU (a.u.)Glasgow coma scale when entering ICU
 GCS_before_SBT (a.u.)Glasgow coma scale before doing SBT
Mechanical ventilator
 Endo_size (mm)Endotracheal tube diameter
 cuff_leak (mL)Leak of cuff of endotracheal tube
 Reason_using_MV (a.u.)Reason of using mechanical ventilator
Hemodynamicsbefore SBT
 SBP_before_SBT (mmHg)Systolic blood pressure before SBT
 DBP_before_SBT (mmHg)Diastolic blood pressure before SBT
 HR_before_SBT (bpm)Heart rate before SBT
 SpO2_before_SBT (%)Pulse oxygen saturation before SBT
Ventilation
 Pimax (cmH2O)Maximal inspiratory pressure
Disease history(0: no; 1: yes)
 HX_Res (a.u.)History of respiratory disease
 HX_CV (a.u.)History of cardiovascular disease
 HX_NM (a.u.)History of neuromuscular disease
Table 2. SBT-related cardiovascular and respiratory variables.
Table 2. SBT-related cardiovascular and respiratory variables.
VariableDefinition
 SBT_times (a.u.)How many times of doing SBT
Hemodynamicsduring SBT
 mX_HR (bpm)Xth minute heart rate
 mX_SpO2 (%)Xth minute pulse oxygen saturation
Ventilation during SBT
 mX_VE (mL)Xth minute ventilation volume
 mX_SponVT (mL)Xth minute spontaneous tidal volume
 mX_RR (bpm)Xth minute respiratory rate
Table 3. Parameter setting.
Table 3. Parameter setting.
ParameterValueDefinition
Epoch100Number of dataset iterations
N1500Maximum number of rules in [P]
β0.2Learning rate for updating the parameters
α0.1Decline rate of fitness
ν0.1Exponent of fitness
ε010Acceptable range of prediction error threshold
θGA25The threshold for GA executing
χ0.8Probability of GA performing crossover
μ0.04Probability of GA performing mutation
θdel20The threshold for deleting rule
δ0.1The fraction of the mean fitness in [P]
θsub20The threshold for subsumption
P#0.5When covering, the range of inti can cover the maximum and minimum values of xi with a certain probability
pi10Initial value of prediction
εi0Initial value of prediction error
Fi10Initial value of fitness
θmna2Maximum number of actions
Table 4. Characteristics of the final sample of 228 patients.
Table 4. Characteristics of the final sample of 228 patients.
Clinical CharacteristicsExtubation Failure (n = 40)Successful Extubation (n = 188)Total (n = 228)AUC
Basic information
 Gender (male/female)18/22118/70136/92-
 Age (year)72.62 ± 14.0567.53 ± 15.7668.42 ± 15.570.60
 Height (cm)158.5 ± 9.05162.60 ± 08.56 *161.89 ± 8.760.63
 BW (kg)60.73 ± 13.2260.94 ± 14.5460.90 ± 14.290.51
 IBW (kg)55.51 ± 6.3658.27 ± 06.03 *57.78 ± 6.160.63
Severity scale
 APACHEII24.43 ± 8.1921.29 ± 08.3321.84 ± 8.380.61
 GCS_enter_ICU6.30 ± 3.416.84 ± 03.526.74 ± 3.280.55
 GCS_before_SBT9.55 ± 2.0010.10 ± 01.6310.00 ± 1.710.60
Mechanical ventilator
 Endo_size
 7.0 # (%)9 (3.9)31 (13.6)40 (17.5)-
 7.5 # (%)31 (13.6)156 (68.4)188 (82.5)-
 cuff_leak203.47 ± 124.36245.65 ± 107.31 *283.25 ± 111.360.63
 Reason_using_MV21/1993/95114/114-
Hemodynamics before SBT
 SBP_before_SBT144.25 ± 26.34140.94 ± 26.47141.52 ± 26.420.53
 DBP_before_SBT71.40 ± 17.3970.89 ± 15.8670.98 ± 16.100.52
 HR_before_SBT92.95 ± 16.9688.15 ± 15.9188.99 ± 16.160.57
 SpO2_before_SBT97.50 ± 2.1898.22 ± 01.8798.09 ± 1.950.60
Ventilation
 Pimax33.25 ± 11.9934.82 ± 12.1934.55 ± 12.140.54
Disease history
 HX_Res (yes/no)12/2844/14458/170-
 HX_CV (yes/no)30/20104/84134/94-
 HX_NM (yes/no)4/3631/15735/193-
* p < 0.01 between extubation failure and successful extubation.
Table 5. SBT variables of the 228 patients.
Table 5. SBT variables of the 228 patients.
Variables during SBTExtubation Failure (n = 40)Successful Extubation (n = 188)Total (n = 228)AUC
 SBT_times2.58 ± 1.412.49 ± 1.49 2.51 ± 1.560.54
Hemodynamics
 m1_HR94.15 ± 16.2889.07 ± 15.7789.96 ± 15.940.59
 m1_SpO297.70 ± 2.2998.02 ± 2.1397.96 ± 2.150.54
 m30_HR93.90 ± 15.9788.69 ± 17.5589.60 ± 17.360.60
 m30_SpO297.40 ± 3.1097.71 ± 2.4697.65 ± 2.580.51
 m60_HR96.40 ± 17.0488.12 ± 20.33 *89.57 ± 20.010.63
 m60_SpO297.58 ± 2.7796.12 ± 12.4796.38 ± 11.390.52
 m90_HR92.17 ± 21.8685.94 ± 20.3689.04 ± 20.720.61
 m90_SpO287.78 ± 29.7392.12 ± 23.1191.36 ± 24.380.53
Ventilation
 m1_VE7.93 ± 2.977.94 ± 2.567.94 ± 2.630.52
 m1_SponVT420.40 ± 206.83463.26 ± 222.54455.74 ± 220.030.58
 m1_RR19.60 ± 6.8918.73 ± 6.6118.88 ± 6.660.53
 m30_VE7.53 ± 2.467.90 ± 2.787.83 ± 2.720.53
 m30_SponVT442.27 ± 244.51531.98 ± 706.51516.24 ± 650.110.60
 m30_RR18.65 ± 6.1917.37 ± 5.6817.60 ± 5.780.57
 m60_VE8.20 ± 3.467.84 ± 2.627.90 ± 2.780.52
 m60_SponVT451.43 ± 208.80472.49 ± 174.87468.79 ± 180.950.55
 m60_RR18.87 ± 6.2217.57 ± 6.1717.80 ± 6.190.57
 m90_VE7.65 ± 2.637.81 ± 2.457.82 ± 2.430.53
 m90_SponVT466.93 ± 191.02478.64 ± 187.76478.68 ± 185.680.55
 m90_RR17.65 ± 4.0117.52 ± 6.1217.62 ± 5.690.53
* p < 0.01 between extubation failure and successful extubation.
Table 6. Usage probabilities of the top 10 variables. The first 6 variables have significant higher probabilities.
Table 6. Usage probabilities of the top 10 variables. The first 6 variables have significant higher probabilities.
RankingVariableProbusage (%)
1m90_RR96.74 ± 1.27
2m60_VE96.58 ± 1.92
3m90_VE96.07 ± 1.67
4m30_VE95.90 ± 1.49
5m1_VE95.71 ± 1.54
6BW95.61 ± 2.70
7age36.10 ± 5.23
8Height35.29 ± 3.33
9SpO2_before_SBT34.93 ± 4.77
10DBP_before_SBT34.76 ± 4.43
40HX_NM26.03 ± 5.28
Table 7. AUC results of the current and previous studies (bold text: the result of our group).
Table 7. AUC results of the current and previous studies (bold text: the result of our group).
Ref.MethodTypeExtubation Failure n (%)AUC
Zeggwagh, et al. (1999) [21]MLRCalibration22 (36.67)0.913
MLRValidation22 (37.29)0.855
Mueller, et al. (2004) [22]MLRCalibration22 (16.92)0.810
MLRValidation13 (24.53)0.750
ANNCalibration22 (16.92)0.810
ANNValidation13 (24.53)0.870
Mueller, et al. (2013) [24]MLRCalibration59 (12.14)0.880
MLRValidation59 (12.14)0.762
ANNCalibration59 (12.14)0.930
ANNValidation59 (12.14)0.753
BDTCalibration59 (12.14)1.000
BDTValidation59 (12.14)0.513
NBCCalibration59 (12.14)0.610
NBCValidation59 (12.14)0.626
SVMValidation59 (12.14)0.493
Hsieh, et al. (2018) [23]ANNCalibration185 (5.14)0.850
Chung, et al. (2020) [20]MLRCalibration28 (16.57)0.889
This study (2022)MLRCalibration40 (17.54)0.660
PLSCalibration40 (17.54)0.812
XCSRCalibration40 (17.54)0.993
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, P.-H.; Chen, L.-Y.; Chung, W.-C.; Sheu, C.-C.; Hsiao, T.-C.; Tsai, J.-R. Toward Evaluating Critical Factors of Extubation Outcome with XCSR-Generated Rules. Bioengineering 2022, 9, 701. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering9110701

AMA Style

Huang P-H, Chen L-Y, Chung W-C, Sheu C-C, Hsiao T-C, Tsai J-R. Toward Evaluating Critical Factors of Extubation Outcome with XCSR-Generated Rules. Bioengineering. 2022; 9(11):701. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering9110701

Chicago/Turabian Style

Huang, Po-Hsun, Lian-Yu Chen, Wei-Chan Chung, Chau-Chyun Sheu, Tzu-Chien Hsiao, and Jong-Rung Tsai. 2022. "Toward Evaluating Critical Factors of Extubation Outcome with XCSR-Generated Rules" Bioengineering 9, no. 11: 701. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering9110701

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop