Selecting Optimal Long Short-Term Memory (LSTM) Architectures for Online Estimation of Mean Arterial Pressure (MAP) in Patients Undergoing General Anesthesia
Abstract
:1. Introduction
2. Methodology
2.1. Clinical Trial Protocol
2.2. Model Design and Architecture
2.3. Performance Metrics
- (i)
- Mean Squared Error (MSE): MSE measures the average squared difference between the actual and predicted values and is calculated as
- (ii)
- Mean Absolute Error (MAE): MAE computes the average absolute difference between the predicted and actual values, providing a more interpretable measure:
- (iii)
- Coefficient of Determination (): quantifies the proportion of the variance in the dependent variable that is predictable from the independent variable(s), and is given by
2.4. Population-Specific Method
2.5. Patient-Specific Method
Algorithm 1 Best LSTM architecture search |
|
2.6. Category-Specific Method
3. Results
3.1. Population-Specific Method Results
3.2. Patient-Specific Method Results
3.3. Category-Specific Method Results
3.4. Comparative Study
4. Discussion
- The current study method provides estimation, not forecasting, which means it estimates the MAP value at a particular moment rather than forecasting future values. Although the proposed method offers the advantage of having MAP data that can be sampled as desired, the method could be improved in future research to forecast MAP within a precise prediction horizon. Forecasting MAP values can significantly enhance monitoring capabilities and early warning systems for adverse events. By predicting MAP values, it is possible to enable early detection of potential issues such as hypotension (low blood pressure) [36] and hypertension (high blood pressure) [37], allowing timely interventions to prevent complications. This predictive approach improves patient outcomes by providing healthcare professionals with advanced notice of potential adverse events, facilitating prompt and appropriate responses.
- This approach needs an offline process at the start of a TIVA trial to categorize patients based on their correlation and to identify the best model architecture for each category. Investigating the factors underlying these correlations could improve this approach. Factors such as the medical history of the patient (e.g., heart disease), biometric data (e.g., weight, age, sex), and the type of surgery may all influence the correlation between anesthetic drug concentrations (Propofol/Remifentanil) and MAP. Understanding these factors could facilitate the development of a more efficient categorization process. In our database, patients are categorized into four groups to examine the influence of BMI on the average correlations between the inputs (CeRemi and CeProp) and the output MAP.In Table 5, the number of patients found in each BMI category is shown. The figure below illustrates the distribution of average correlations between MAP and the inputs (CeRemi and CeProp) across different BMI categories.The boxplot in Figure 8 shows a significant difference (p-value = 0.0233) in the average correlation based on BMI categories. This indicates that BMI significantly affects the correlation between drug concentrations and MAP. This finding suggests that BMI can be a factor used in the development of a more efficient categorization process.
- The dataset used in this paper consists of a relatively small sample size of 70 patients from a single population, recognizing the non-generalizability of our approach, which relies on a single database. Future research should expand upon these findings, exploring alternative deep learning models, e.g., feed-forward neural networks (FFNs) and recurrent neural networks (RNNs).
5. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
Appendix A. LSTM Model for MAP Prediction
- Input gate (): The input gate controls the flow of new information into the cell state. It is computed using the sigmoid activation function and is defined as follows:
- Forget gate (): The forget gate controls the retention of past information in the cell state and acts as a weighting factor of past-to-new data. It is computed similarly to the input gate and is defined as follows:
- Candidate cell state (): The candidate cell state represents the new information that could be stored in the cell state. It is computed using the hyperbolic tangent activation function, defined as follows:
- Cell state update (): The cell state is updated by combining the previous cell state with the new information from the input gate and candidate cell state:
- Output gate (): The output gate controls what information from the cell state should be used to compute the output. It is defined as follows:
Appendix B. Patient Database
Index | Gender | Age (Years) | Height (cm) | Weight (kg) | BMI (km/m2) |
---|---|---|---|---|---|
1 | F | 31 | 163 | 49 | 18.4 |
2 | F | 44 | 168 | 57 | 20.2 |
3 | F | 42 | 168 | 70 | 24.8 |
4 | F | 71 | 158 | 69 | 27.6 |
5 | F | 38 | 182 | 87.6 | 26.4 |
6 | F | 41 | 175 | 135 | 44.1 |
7 | F | 68 | 172 | 67 | 22.6 |
8 | M | 59 | 183 | 76 | 22.7 |
9 | M | 50 | 186 | 96 | 27.7 |
10 | M | 73 | 181 | 83 | 25.3 |
11 | X | 22 | 165 | 87 | 32 |
12 | M | 29 | 183 | 92 | 27.5 |
13 | F | 54 | 160 | 48 | 18.8 |
14 | X | 19 | 168.7 | 59 | 20.7 |
15 | X | 26 | 155 | 61 | 25.4 |
16 | F | 54 | 163 | 58 | 21.8 |
17 | F | 50 | 174 | 64 | 21.1 |
18 | F | 30 | 176 | 74 | 23.9 |
19 | F | 57 | 164 | 70 | 26 |
20 | X | 33 | 162 | 80 | 30.5 |
21 | F | 62 | 168 | 88 | 31.2 |
22 | F | 48 | 155 | 56 | 23.3 |
23 | M | 62 | 183 | 85 | 25.4 |
24 | F | 36 | 168 | 63 | 22.3 |
25 | M | 58 | 179 | 94 | 29.3 |
26 | F | 65 | 162 | 87 | 33.2 |
27 | F | 49 | 167 | 86 | 30.8 |
28 | M | 46 | 187 | 97 | 27.7 |
29 | M | 68 | 176 | 85 | 27.4 |
30 | F | 50 | 167 | 70 | 25.1 |
31 | M | 54 | 175 | 90 | 29.4 |
32 | F | 64 | 164 | 96 | 35.7 |
33 | F | 46 | 170 | 82 | 28.4 |
34 | F | 58 | 157 | 53.5 | 21.7 |
35 | F | 41 | 167 | 65 | 23.3 |
Index | Gender | Age (Years) | Height (cm) | Weight (kg) | BMI (kg/m2) |
---|---|---|---|---|---|
36 | F | 40 | 159 | 63 | 24.9 |
37 | F | 30 | 160 | 53 | 20.7 |
38 | F | 18 | 163 | 64 | 24.1 |
39 | F | 64 | 158 | 58 | 23.2 |
40 | M | 63 | 171 | 104 | 35.6 |
41 | F | 30 | 156 | 44 | 18.1 |
42 | F | 19 | 168 | 54 | 19.1 |
43 | F | 51 | 168 | 58 | 20.5 |
44 | F | 59 | 171 | 75 | 25.6 |
45 | M | 73 | 174 | 67 | 22.1 |
46 | M | 23 | 173 | 66 | 22.1 |
47 | M | 25 | 181 | 74.5 | 22.7 |
48 | M | 31 | 170 | 85 | 29.4 |
49 | M | 32 | 186 | 100 | 28.9 |
50 | F | 40 | 157 | 85 | 34.5 |
51 | F | 27 | 157 | 58 | 23.5 |
52 | F | 71 | 174 | 71 | 23.5 |
53 | M | 32 | 180 | 78 | 24.1 |
54 | F | 21 | 166 | 63 | 22.9 |
55 | X | 48 | 170 | 84 | 29.1 |
56 | M | 63 | 168 | 82 | 29.1 |
57 | M | 71 | 175 | 66 | 21.6 |
58 | F | 48 | 165 | 76 | 27.9 |
59 | F | 71 | 167 | 83 | 29.8 |
60 | X | 21 | 160 | 60 | 23.4 |
61 | X | 21 | 171 | 73 | 25 |
62 | X | 34 | 172 | 90 | 30.4 |
63 | X | 26 | 170 | 50 | 17.3 |
64 | F | 33 | 172 | 94 | 31.8 |
65 | X | 22 | 172 | 69 | 23.3 |
66 | F | 70 | 165 | 62 | 22.8 |
67 | X | 21 | 168 | 56 | 19.8 |
68 | F | 35 | 160 | 82 | 32 |
69 | X | 21 | 165 | 70 | 25.7 |
70 | X | 20 | 155 | 49 | 20.4 |
Appendix C. MAP Variability among Different Patients and the Safety Intervals
References
- Hino, H.; Matsuura, T.; Kihara, Y.; Tsujikawa, S.; Mori, T.; Nishikawa, K. Comparison between hemodynamic effects of propofol and thiopental during general anesthesia induction with remifentanil infusion: A double-blind, age-stratified, randomized study. J. Anesth. 2019, 33, 509–515. [Google Scholar] [CrossRef]
- Menacho, S.T.; Floyd, C. Current practices and goals for mean arterial pressure and spinal cord perfusion pressure in acute traumatic spinal cord injury: Defining the gaps in knowledge. J. Spinal Cord Med. 2021, 44, 350–356. [Google Scholar] [CrossRef] [PubMed]
- Cava, C.; Salvatore, C.; Castiglioni, I. Pan-cancer classification of gene expression data based on artificial neural network model. Appl. Sci. 2023, 13, 7355. [Google Scholar] [CrossRef]
- Zhang, C.; Zhang, L.; Tian, Y.; Bao, B.; Li, D. A machine-learning-algorithm-assisted intelligent system for real-time wireless respiratory monitoring. Appl. Sci. 2023, 13, 3885. [Google Scholar] [CrossRef]
- De Grauwe, T.; Ghita, M.; Copot, D.; Ionescu, C.M.; Neckebroek, M. Artificial intelligence for pain classification with the non-invasive pain monitor Anspec-Pro. Acta Anaesthesiol. Belg. 2022, 73, 45–52. [Google Scholar] [CrossRef]
- Shamshirband, S.; Fathi, M.; Dehzangi, A.; Chronopoulos, A.T.; Alinejad-Rokny, H. A review on deep learning approaches in healthcare systems: Taxonomies, challenges, and open issues. J. Biomed. Inform. 2021, 113, 103627. [Google Scholar] [CrossRef] [PubMed]
- Xu, Y.; Lu, X.; Cetiner, B.; Taciroglu, E. Real-time regional seismic damage assessment framework based on long short-term memory neural network. Comput.-Aided Civ. Infrastruct. Eng. 2021, 36, 504–521. [Google Scholar] [CrossRef]
- Shah, S.Y.; Yuan, Z.; Lu, S.; Zerfos, P. Dependency analysis of cloud applications for performance monitoring using recurrent neural networks. In Proceedings of the 2017 IEEE International Conference on Big Data (Big Data), Boston, MA, USA, 11–14 December 2017; IEEE: Piscataway, NJ, USA, 2017; pp. 1534–1543. [Google Scholar]
- Hewamalage, H.; Bergmeir, C.; Bandara, K. Recurrent neural networks for time series forecasting: Current status and future directions. Int. J. Forecast. 2021, 37, 388–427. [Google Scholar] [CrossRef]
- Ionescu, C.M.; Keyser, R.D.; Copot, D.; Yumuk, E.; Ynineb, A.; Othman, G.B.; Neckebroek, M. Model Extraction From Clinical Data Subject to Large Uncertainties and Poor Identifiability. IEEE Control Syst. Lett. 2024, 1. [Google Scholar] [CrossRef]
- Pham, T.D. Time–frequency time–space LSTM for robust classification of physiological signals. Sci. Rep. 2021, 11, 6936. [Google Scholar] [CrossRef]
- Heynen, J.; Copot, D.; Ghita, M.; Ionescu, C.M. Using convolutional neural network online estimators for predicting pain-level variability enables predictive control of anesthesia. In Proceedings of the 2021 25th International Conference on System Theory, Control and Computing (ICSTCC), Iași, Romania, 20–23 October 2021; IEEE: Piscataway, NJ, USA, 2021; pp. 194–199. [Google Scholar]
- Gundu, V.; Simon, S.P. PSO–LSTM for short term forecast of heterogeneous time series electricity price signals. J. Ambient Intell. Humaniz. Comput. 2021, 12, 2375–2385. [Google Scholar] [CrossRef]
- Khorram, A.; Khalooei, M.; Rezghi, M. End-to-end CNN+ LSTM deep learning approach for bearing fault diagnosis. Appl. Intell. 2021, 51, 736–751. [Google Scholar] [CrossRef]
- Schnider, T.W.; Minto, C.F.; Gambus, P.L.; Andresen, C.; Goodale, D.B.; Shafer, S.L.; Youngs, E.J. The influence of method of administration and covariates on the pharmacokinetics of propofol in adult volunteers. J. Am. Soc. Anesthesiol. 1998, 88, 1170–1182. [Google Scholar] [CrossRef]
- Minto, C.F.; Schnider, T.W.; Egan, T.D.; Youngs, E.; Lemmens, H.J.; Gambus, P.L.; Billard, V.; Hoke, J.F.; Moore, K.H.; Hermann, D.J.; et al. Influence of age and gender on the pharmacokinetics and pharmacodynamics of remifentanil: I. Model development. J. Am. Soc. Anesthesiol. 1997, 86, 10–23. [Google Scholar] [CrossRef]
- Karnika, S. ISO 14155: Clinical Investigation of Medical Devices for Human Subjects. In Medical Device Guidelines and Regulations Handbook; Springer: Berlin/Heidelberg, Germany, 2022; pp. 1–18. [Google Scholar]
- Global, S.P.; Windisch, F.; Zimmermann, N.; Knoll, V.; Christodoulou, M.; Habimana, K.; Piccoli, G.; Del Prete, M. Study on the Implementation of Article 17 of Regulation (EU) 2017/745 on Medical Devices on the EU Market; European Union: Brussels, Belgium, 2024. [Google Scholar]
- Ionescu, C.M.; Copot, D.; Yumuk, E.; De Keyser, R.; Muresan, C.; Birs, I.R.; Ben Othman, G.; Farbakhsh, H.; Ynineb, A.R.; Neckebroek, M. Development, Validation, and Comparison of a Novel Nociception/Anti-Nociception Monitor against Two Commercial Monitors in General Anesthesia. Sensors 2024, 24, 2031. [Google Scholar] [CrossRef] [PubMed]
- Ahmed, I.; Ahmad, M.; Chehri, A.; Jeon, G. A heterogeneous network embedded medicine recommendation system based on LSTM. Future Gener. Comput. Syst. 2023, 149, 1–11. [Google Scholar] [CrossRef]
- Kaushik, S.; Choudhury, A.; Sheron, P.K.; Dasgupta, N.; Natarajan, S.; Pickett, L.A.; Dutt, V. AI in healthcare: Time-series forecasting using statistical, neural, and ensemble architectures. Front. Big Data 2020, 3, 4. [Google Scholar] [CrossRef] [PubMed]
- Goodfellow, I.; Bengio, Y.; Courville, A. Deep Learning; MIT Press: Cambridge, MA, USA, 2016. [Google Scholar]
- Smagulova, K.; James, A.P. A survey on LSTM memristive neural network architectures and applications. Eur. Phys. J. Spec. Top. 2019, 228, 2313–2324. [Google Scholar] [CrossRef]
- Mohamed, M.A.; Hassan, H.A.; Essai, M.H.; Esmaiel, H.; Mubarak, A.S.; Omer, O.A. Modified gate activation functions of Bi-LSTM-based SC-FDMA channel equalization. J. Electr. Eng. 2023, 74, 256–266. [Google Scholar] [CrossRef]
- Nagabushanam, P.; Thomas George, S.; Radha, S. EEG signal classification using LSTM and improved neural network algorithms. Soft Comput. 2020, 24, 9981–10003. [Google Scholar] [CrossRef]
- Yu, Y.; Si, X.; Hu, C.; Zhang, J. A review of recurrent neural networks: LSTM cells and network architectures. Neural Comput. 2019, 31, 1235–1270. [Google Scholar] [CrossRef] [PubMed]
- Beale, M.H.; Hagan, M.T.; Demuth, H.B. Deep Learning Toolbox; User’s Guide; The MathWorks, Inc.: Natick, MA, USA, 2018. [Google Scholar]
- Jin, G.; Yi, X.; Zhang, L.; Zhang, L.; Schewe, S.; Huang, X. How does weight correlation affect generalisation ability of deep neural networks? Adv. Neural Inf. Process. Syst. 2020, 33, 21346–21356. [Google Scholar]
- Chen, J.; Liu, Z.; Yin, Z.; Liu, X.; Li, X.; Yin, L.; Zheng, W. Predict the effect of meteorological factors on haze using BP neural network. Urban Clim. 2023, 51, 101630. [Google Scholar] [CrossRef]
- Sharma, J.; Soni, S.; Paliwal, P.; Saboor, S.; Chaurasiya, P.K.; Sharifpur, M.; Khalilpoor, N.; Afzal, A. A novel long-term solar photovoltaic power forecasting approach using LSTM with Nadam optimizer: A case study of India. Energy Sci. Eng. 2022, 10, 2909–2929. [Google Scholar] [CrossRef]
- Chang, Z.; Zhang, Y.; Chen, W. Effective adam-optimized LSTM neural network for electricity price forecasting. In Proceedings of the 2018 IEEE 9th International Conference on Software Engineering and Service Science (ICSESS), Beijing, China, 23–25 November 2018; IEEE: Piscataway, NJ, USA, 2018; pp. 245–248. [Google Scholar]
- Vijayaprabakaran, K.; Sathiyamurthy, K. Towards activation function search for long short-term model network: A differential evolution based approach. J. King Saud Univ.-Comput. Inf. Sci. 2022, 34, 2637–2650. [Google Scholar]
- Salman, A.G.; Heryadi, Y.; Abdurahman, E.; Suparta, W. Single layer & multi-layer long short-term memory (LSTM) model with intermediate variables for weather forecasting. Procedia Comput. Sci. 2018, 135, 89–98. [Google Scholar]
- Abbasimehr, H.; Shabani, M.; Yousefi, M. An optimized model using LSTM network for demand forecasting. Comput. Ind. Eng. 2020, 143, 106435. [Google Scholar] [CrossRef]
- Ionescu, C.M.; Neckebroek, M.; Ghita, M.; Copot, D. An Open source patient simulator for design and evaluation of computer based multiple drug dosing control for anesthetic and hemodynamic variables. IEEE Access 2021, 9, 8680–8694. [Google Scholar] [CrossRef]
- Frassanito, L.; Sonnino, C.; Piersanti, A.; Zanfini, B.A.; Catarci, S.; Giuri, P.P.; Scorzoni, M.; Gonnella, G.L.; Antonelli, M.; Draisci, G. Performance of the hypotension prediction index with noninvasive arterial pressure waveforms in awake cesarean delivery patients under spinal anesthesia. Anesth. Analg. 2022, 134, 633–643. [Google Scholar] [CrossRef]
- Kandil, H.; Soliman, A.; Alghamdi, N.S.; Jennings, J.R.; El-Baz, A. Using mean arterial pressure in hypertension diagnosis versus using either systolic or diastolic blood pressure measurements. Biomedicines 2023, 11, 849. [Google Scholar] [CrossRef]
- Khanna, A.K.; Kinoshita, T.; Natarajan, A.; Schwager, E.; Linn, D.D.; Dong, J.; Ghosh, E.; Vicario, F.; Maheshwari, K. Association of systolic, diastolic, mean, and pulse pressure with morbidity and mortality in septic ICU patients: A nationwide observational study. Ann. Intensive Care 2023, 13, 9. [Google Scholar] [CrossRef] [PubMed]
Parameter | Role | MATLAB Code Example |
---|---|---|
LSTM layers | Determines the depth of the neural network. | layers = [lstmLayer(16), lstmLayer(32), …] |
Neurons per layer | Defines the number of basic processing units in each LSTM layer. | neurons = [16, 64]; |
State activation function | Activation function for updating the cell and hidden state in the LSTM layer. | stateActivation = {tanh, softsign}; |
Gate activation function | Activation function for the gates in the LSTM layer. | gateActivation = {sigmoid, hard-sigmoid}; |
Number of epochs | Controls the number of training iterations through the entire dataset. | epochs = [50, 100, 200, …]; |
Optimizer | Adjusts model parameters during training to minimize the error. | optimizer = {adam, rmsprop, sgd, …}; |
Mini-batch size | Determines the number of training examples used in each update. | miniBatchSize = [16, 32, 64]; |
Categories | Correlation Range | Patients |
---|---|---|
1 Weak or no correlation | 0 to 0.10 | Patient 15 |
2 Weak correlation | 0.10 to 0.20 | Patient 59 |
3 Mild correlation | 0.20 to 0.40 | Patients 10, 28, 30, 52, 60, 62, 65, 68 |
4 Moderate correlation | 0.40 to 0.70 | Patients 2, 4, 6, 12, 13, 19, 20, 21, 23, 26, 27, 29, 31, 32, 33, 34, 37, 40, 41, 42, 44, 45, 46, 47, 49, 50, 51, 53, 55, 57, 58, 61, 67, 69 |
5 Strong correlation | 0.70 to 1.00 | Patient 1, 3, 5, 7, 8, 9, 11, 14, 16, 17, 18, 22, 24, 25, 35, 36, 38, 39, 43, 48, 54, 56, 63, 64, 66, 70 |
Population-Specific | Category-Specific | Patient-Specific | |
---|---|---|---|
MAE | 0.4139 ± 0.05 | 0.38115 ± 0.04 | 0.36325 ± 0.03 |
MSE | 0.332 ± 0.04 | 0.280 ± 0.03 | 0.222 ± 0.02 |
0.594 ± 0.07 | 0.606 ± 0.06 | 0.722 ± 0.05 |
Method | Time Consumed to Find the Architecture |
---|---|
Population-specific | Less than 1 h |
Patient-specific | 130 h |
Category-specific | 53 h |
BMI Category | BMI Range | Number of Patients |
---|---|---|
Underweight | <18.5 | 3 |
Normal weight | 18.5–24.9 | 31 |
Overweight | 25–29.9 | 24 |
Obesity | ≥30 | 12 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Ben Othman, G.; Copot, D.; Yumuk, E.; Neckebroek, M.; Ionescu, C.M. Selecting Optimal Long Short-Term Memory (LSTM) Architectures for Online Estimation of Mean Arterial Pressure (MAP) in Patients Undergoing General Anesthesia. Appl. Sci. 2024, 14, 5556. https://doi.org/10.3390/app14135556
Ben Othman G, Copot D, Yumuk E, Neckebroek M, Ionescu CM. Selecting Optimal Long Short-Term Memory (LSTM) Architectures for Online Estimation of Mean Arterial Pressure (MAP) in Patients Undergoing General Anesthesia. Applied Sciences. 2024; 14(13):5556. https://doi.org/10.3390/app14135556
Chicago/Turabian StyleBen Othman, Ghada, Dana Copot, Erhan Yumuk, Martine Neckebroek, and Clara M. Ionescu. 2024. "Selecting Optimal Long Short-Term Memory (LSTM) Architectures for Online Estimation of Mean Arterial Pressure (MAP) in Patients Undergoing General Anesthesia" Applied Sciences 14, no. 13: 5556. https://doi.org/10.3390/app14135556
APA StyleBen Othman, G., Copot, D., Yumuk, E., Neckebroek, M., & Ionescu, C. M. (2024). Selecting Optimal Long Short-Term Memory (LSTM) Architectures for Online Estimation of Mean Arterial Pressure (MAP) in Patients Undergoing General Anesthesia. Applied Sciences, 14(13), 5556. https://doi.org/10.3390/app14135556