Comparative Analysis of Stress Responses in Medical Students Using Virtual Reality Versus Traditional 3D-Printed Mannequins for Pericardiocentesis Training =========================================================================================================================================================== * Alberto Rubio-López * Rodrigo García Carmona * Laura Zarandieta Román * Alejandro Rubio Navas * Ángel González Pinto * Pablo Cardinal-Fernández ## ABSTRACT **Background** As medical education evolves, innovative methods like virtual reality (VR) and 3D-printed mannequins are increasingly used to simulate high-stress medical scenarios realistically. This study investigates the effectiveness of VR and 3D-printed mannequins in replicating stress levels during pericardiocentesis training, comparing their impact on the emotional and physiological responses of learners. **Methods** We enrolled 108 final-year medical students who were randomized to train with both VR and 3D-printed mannequins. Heart rate variability (HRV) analysis was employed to assess stress responses. Additionally, a secondary analysis examined the influence of demographic factors, lifestyle, medication use, and academic stress on these responses. **Results** Both VR and traditional mannequin-based training methods proved equally effective in simulating the stress levels encountered in real medical procedures. Our findings indicate significant interactions between stress markers and demographic factors, which highlights the complex nature of stress responses in medical education and underscores the necessity for personalized training approaches. **Conclusion** The study validates the use of VR as a viable alternative to traditional mannequins, capable of simulating the technical skills and emotional pressures of medical procedures such as pericardiocentesis. Incorporating VR into medical training programs may enhance learning outcomes and accessibility, particularly in settings constrained by resources. Keywords * Virtual Reality (VR) * Medical Education * Pericardiocentesis Training * Stress Levels * Heart Rate Variability (HRV) * 3D-Printed Mannequin * Biometric Data * Real-life Scenarios * Demographic Factors * Personalized Educational Approaches ## INTRODUCTION The field of medical education has been revolutionized by the adoption of simulation-based learning, which aligns with the “learning by doing” [1] philosophy. This approach enhances patient safety by allowing students to practice clinical procedures in safe environments, offering reproducible scenarios for consistent learning outcomes while accommodating various educational levels. Despite its widespread benefits, the application of simulation-based learning is often hampered by the high costs of advanced simulators and the steep learning curves for instructors, which can exacerbate educational disparities, particularly in economically constrained settings [2]. In response to these challenges, recent advancements in 3D printing and virtual reality (VR) technologies have begun to bridge the gap [3]. 3D printing provides a cost-effective means to create detailed medical models [4], thereby broadening the scope and flexibility of simulations. Similarly, VR offers immersive and realistic environments that enhance both sensory and emotional training experiences [5], making these tools accessible across different educational settings and platforms such as Oculus and HTC. Despite these technological advances, a critical gap persists in evaluating the emotional impact of these simulations, especially their ability to evoke authentic stress responses during high-stakes procedures [6]. This study addresses this gap by focusing on pericardiocentesis—a procedure notorious for its complexity and high stress levels [7]. We aim to compare the efficacy of traditional 3D-printed mannequins and VR simulations in eliciting realistic stress responses as measured by heart rate variability (HRV), a validated stress biomarker. Through this comparison, we seek to ascertain which simulation method better prepares medical students for the real-world pressures of medical emergencies, particularly in resource-limited settings [8]. ## METHODS ### Study Design This validation study builds on our pilot research to compare two simulation methods for pericardiocentesis training: a 3D-printed mannequin and a virtual reality (VR) setup using the Unity game engine. The primary objective was to evaluate and contrast the ability of these methods to simulate realistic emotional and physiological responses encountered in actual medical procedures. ### Randomization and Participant Grouping We enrolled 132 final-year medical students from the CEU University School of Medicine in Madrid. Before participation, students completed demographic and health behavior questionnaires. Participants were randomly assigned into two groups using a computer-generated sequence, with one group starting with VR training followed by the mannequin training, and the other group starting in reverse order. This counterbalanced design aimed to mitigate any order effects on learning outcomes and stress responses. ### Development and Validation of the Training Models The 3D-printed mannequin was designed using Tinkercad and Adobe Fusion 360 to ensure anatomical precision. It was printed with MeshMixer software, Ultimaker Cura for slicing, and an Ender 3 printer. Concurrently, the VR model was developed in Unity to mirror the mannequin’s design, ensuring a consistent and immersive training experience. Both models were rigorously validated by clinical experts from HM Montepríncipe University Hospital in Madrid. ### Ethical considerations and recruitment The study received approval from the HM Montepríncipe University Hospital Research Ethics Committee (Code: 18.12.1339.GHM) and spanned from November 7, 2023, to February 12, 2024. All participants provided written informed consent, and their data were anonymized to ensure confidentiality. ### Inclusion and Exclusion Criteria All consenting final-year medical students were eligible, except those on cardiac medication, consuming excessive caffeine, or simpaticomimetics. Students with incomplete data or improperly recorded signals were excluded from the analysis. Participants received standardized training through an instructional video on pericardiocentesis before the simulations. ### Stress response measurement To assess stress responses, students were equipped with three electrodes on the torso to continuously monitor heart rate using the Biosignal Plux system. This setup was established prior to the initiation of simulation training to ensure robust data collection on physiological stress indicators. The collected data were time-stamped, allowing for the separate analysis of each part of the corresponding procedure. Data acquisition, processing, and conversion into values for statistical analysis were conducted using OpenSignals digital signal processing software Heart rate variability (HRV) analysis was employed to assess stress, as it serves as a sensitive marker for dysregulation in the Autonomic Nervous System (ANS). HRV is defined as the temporal variation in the intervals between consecutive heartbeats over a predefined period, providing a quantitative measure of autonomic nervous system activity that is critical under stressful conditions. The parameters used in our HRV analysis include: * **Frequency-domain parameters**: * Low Frequency (LF): Primarily reflects sympathetic activation but also includes parasympathetic influences. It is used to assess the balance of autonomic input. * High Frequency (HF): Acts as a specific indicator of parasympathetic nervous system (PNS) activity and is linked to stress relief mechanisms. * LF/HF Ratio: Serves as an indicator of the balance between sympathetic and parasympathetic activity, with higher values indicating increased stress levels. * **Time-domain parameters**: * Root Mean Square of Successive Differences (rMSSD): Measures the short-term variability in heart rate and is closely related to the parasympathetic nervous system’s control over heart rate. * Percentage of Successive NN Intervals that Differ by More than twenty ms (PNN20) and fifty ms (PNN50): These parameters provide additional insight into the variability of heart rate, correlating strongly with parasympathetic activity. * **Non-linear parameters**: * SD1/SD2 Ratio: Derived from the Poincaré plot, which visually represents the temporal intervals between consecutive heartbeats. This ratio correlates with the LF/HF ratio and provides insight into the complexity of heart rate variability patterns, reflecting stress levels. These measures collectively provide a comprehensive profile of the students’ stress responses during simulation training, offering valuable insights into how different simulation modalities impact the autonomic nervous system under stress-inducing conditions ### Secondary Analysis of Stress Indicators We conducted a secondary analysis to examine the effects of external factors such as medication intake, lifestyle habits, and academic pressures on stress responses. This analysis helped contextualize the physiological data, offering deeper insights into the multifactorial influences on stress levels during clinical training ### Data Collection Phases * **Baseline Data Collection**: Initial stress data were collected in a resting state to establish baseline levels for each participant. * **First Training Scenario**: Students underwent training with either the mannequin or VR, followed by a resting phase to mitigate carryover effects. * **Second Training Scenario**: Participants then trained with the alternate method, facilitating a comparative analysis of stress responses ### Statistical analysis Stress parameters between groups (VR-First and Mannequin-First) were compared using IBM SPSS software version 20. Analysis included: * **Normality Testing**: Kolmogorov–Smirnov test to assess data distribution. * **Comparative Analysis**: Friedman test for intra-group changes over conditions (baseline, simulation, rest, VR). * **Paired Comparisons**: Friedman test and Wilcoxon signed-rank test for specific condition comparisons. * **Secondary Analysis**: Repeated-measures ANOVA to explore relationships between stress responses and demographic factors, with post hoc analyses using Bonferroni and Tukey corrections. ## RESULTS ### Demographic and lifestyle characteristics Our study initially enrolled 132 medical students, which after screening for eligibility criteria, was reduced to 108 participants. These students were equally divided into two groups of fifty-four. The cohort consisted of female students (87%) with an average age of 23.6 years. A significant portion of participants (78%) reported no prior experience with VR headsets. Comprehensive assessments of lifestyle factors—such as living arrangements, relationship status, income levels, familial responsibilities, and work-study balance—were conducted to understand the external influences on stress (Table 1). The impact of medication usage (e.g., anxiolytics, antidepressants) and lifestyle choices (smoking, alcohol consumption, caffeine intake, physical activity) were also examined (Table 2). Additionally, academic factors such as study hours, perceived stress levels, and mental health concerns were evaluated (Table 3) View this table: [Table 1.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T1) Table 1. Lifestyle factors View this table: [Table 2.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T2) Table 2. Medication and Substance Use View this table: [Table 3.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T3) Table 3. Academic and mental health factors. ### Biological stress analysis Heart rate variability (HRV) analysis was employed to assess stress levels across four distinct states: basal, during VR simulation, resting, and during mannequin simulation. The Friedman test was used to analyze differences in biometric parameters across these states, revealing significant changes in HRV parameters (p < 0.01) from the baseline rest state to the active simulation states. This significant variation across time, frequency, and non-linear domains of HRV data underscores the efficacy of our biometric markers in detecting stress under different conditions, as detailed in Table 4 and Figures 1,2. View this table: [Table 4.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T4) Table 4. Biometric parameters at distinct stages of simulation (Friedman test) ![Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/22/2024.08.22.24312406/F1.medium.gif) [Fig. 1.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/F1) Fig. 1. rMSDD values boxplot in the four scenarios (basal, mannequin, rest, and VR) ![Fig. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/22/2024.08.22.24312406/F2.medium.gif) [Fig. 2.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/F2) Fig. 2. SD1/SD2 values boxplot in the four scenarios (basal, mannequin, rest, and VR) Furthermore, we conducted a Wilcoxon signed-rank test to compare the HRV parameters between the two rest phases (basal and resting) and, more critically, between the VR and mannequin simulation phases. The analysis showed no significant differences between the stress responses elicited by the VR and mannequin simulations, indicating that both modalities equally effectively simulate the stress levels experienced in clinical scenarios. This equivalence in stress induction by both VR and traditional mannequin-based training highlights VR’s potential as a viable alternative in medical education, as documented in Table 5 View this table: [Table 5.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T5) Table 5. Biometric parameters paired (Wilcoxon signed-rank test) ### Bias and Sequence Effect Analysis A detailed analysis was conducted to address potential selection biases and the effects of the sequence of simulation exposure. The data confirmed that the order of engagement—whether participants started with VR or with the mannequin—did not significantly impact stress outcomes. This robustness supports the validity of our experimental design and substantiates the equivalence of both training modalities in terms of their stress-inducing capabilities (Tables 6 and 7). View this table: [Table 6.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T6) Table 6. Biometric parameters (First VR Session, Friedman test). View this table: [Table 7.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T7) Table 7. Biometric (First Mannequin Session, Friedman Test) ### Secondary Analysis of the Stress Parameters This secondary analysis delved into the intricate relationships between stress responses and a range of demographic and lifestyle factors, employing repeated-measures ANOVA with Bonferroni and Tukey post hoc tests to rigorously examine these interactions. #### Age and Gender Influences We observed significant interactions between the HRV parameter RMSSD and the combined factors of age and gender (Table 8). This highlights the importance of demographic variations in assessing physiological stress responses in medical education, suggesting that different age groups and genders may experience stress differently during simulations. View this table: [Table 8.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T8) Table 8. Multivariate analysis of variance results for RMSSD by age and sex. #### Lifestyle dynamics Our findings indicated that lifestyle variables, particularly marital status, living conditions, and family income, significantly influenced stress responses. Notably, the interaction between family income and work-study balance was strongly correlated with the LF/HF ratio, a marker of sympathetic activity, suggesting that economic factors may exacerbate stress levels among medical students (Table 9). View this table: [Table 9.](http://medrxiv.org/content/early/2024/08/22/2024.08.22.24312406/T9) Table 9. Multivariate analysis of variance results for LF/HF stratified by marital status, living arrangement, family income, working and studying and family responsibilities. #### Consumption Habits The impact of consumption habits such as smoking, alcohol use, and caffeine intake on stress was less pronounced, with minimal influence on primary stress markers. However, subtle interactions were noted, indicating potential areas for further investigation. #### Medication Impact The combined use of anxiolytics and antidepressants showed complex interactions that significantly affected stress biomarkers, underscoring the need to consider medication history when analyzing stress responses in medical training settings. #### Academic and Perceived Stress Contrary to initial expectations, our study revealed no significant correlations between study habits, perceived academic stress, and physiological stress markers among last-year medical students. This finding suggests that these students, having progressed through rigorous academic training over several years, may have developed effective coping mechanisms to manage and mitigate the physiological manifestations of stress. This adaptation could explain why academic demands and self-perceived stress do not markedly alter physiological stress responses in this group. However, variations might exist with different student populations or in distinct educational settings, indicating a potential area for further research to explore how stress coping mechanisms evolve throughout medical education. ## DISCUSSION Our study suggests that virtual reality (VR) is comparable to traditional mannequin-based training in simulating stress responses during pericardiocentesis. While VR offers immersive and interactive environments, our findings cautiously support its equivalence rather than superiority over traditional methods. This aligns with research by Tsun-Ying and Marvin Mergen, who noted VR’s capability in procedural training without conclusively establishing its enhancement of clinical decision-making skills [3]. Consistent with Bolton’s findings, VR shows promise in enhancing access to medical training in low- and middle-income countries [9]. However, its broader application, including overcoming geographical and economic barriers, must be approached with a focus on sustainability and adaptability to local contexts. Concerns raised by Weissgalss regarding medical artificial intelligence underscore the need for robust policy frameworks to manage potential biases [10]. Our study also explored the integration of 3D modeling with VR, which enhances both the cost-effectiveness and educational accessibility of medical training. Ellen M. Hong and Tae Hoon Roh provided complementary evidence on the benefits of open-source software and photorealistic 3D models in surgical training [11] [12]. Such advancements support the broader application of these technologies in various medical disciplines, particularly in resource-limited environments. Moreover, our research delves into the psychological aspects of medical training. By incorporating 3D-printed mannequins and VR, we address both the technical and emotional challenges faced by healthcare professionals. This dual-focus approach is corroborated by studies from Ashley Towers and Lauryn R [13]. Rochlen noted significant improvements in student confidence and procedural accuracy with VR and 3D technologies [14]. The application of VR in critical care settings, as discussed by Bruno, highlights its potential to enhance educational outcomes while also presenting challenges related to technology integration, cost, and ethical considerations [15]. These insights underscore the necessity of navigating technical and human factors to fully exploit the benefits of immersive technologies in high-stakes environments. Our secondary analysis, which explored the relationships between stress parameters and demographic variables, revealed significant interactions between physiological biomarkers and age, sex, lifestyle, and medication intake. Tessa Helman’s work provides an important context for these findings, highlighting how stress responses can differ significantly between genders, potentially influencing the development of conditions such as coronary artery disease [16]. These insights emphasize the complexity of physiological responses to stress and underscore the importance of personalized training approaches that consider individual demographic and lifestyle factors. Moreover, the impact of lifestyle dynamics on physiological stress markers highlighted in our study emphasizes the complex interplay between socioeconomic factors and health. The findings suggest that family income and work-study balance significantly influence stress-related physiological responses, indicating the need for comprehensive lifestyle assessments in medical training. This finding aligns with Reinhold Kreutz’s observations on the effects of lifestyle changes induced by the COVID-19 pandemic on blood pressure and hypertension, further illustrating the interconnectedness of lifestyle factors and health outcomes [17]. Finally, our inclusion of heart rate variability (HRV) analysis in VR simulations provides a novel approach for assessing stress responses, which is critical for preparing medical professionals to perform under pressure. This aspect of our study aligns with findings from Sean L. Corrigan and David Narciso [18], who emphasized VR’s ability to replicate and manage real-world stress conditions effectively [19]. This study contributes to the understanding of VR’s potential in medical education, particularly for complex procedures. It emphasizes the need for careful integration of new technologies in training programs and highlights the importance of addressing infrastructural and ethical challenges to maximize their impact. This study’s limitations are notable and must be carefully considered when interpreting the findings. First, the nonrandomized design may introduce biases that could affect the robustness of the conclusions. The context of a high-income country also limits the generalizability of our results to settings with different economic conditions. Importantly, our participant group consisted exclusively of undergraduate medical students, which raises concerns about the applicability of our findings to practicing healthcare professionals. The experiences and stress responses of students may differ significantly from those of seasoned professionals who have developed more advanced coping mechanisms and clinical skills. Furthermore, focusing solely on a single training model restricts our understanding of how various VR platforms and technologies might perform under different educational scenarios Conversely, the study’s strengths are notable. The combination of traditional simulation with VR technologies fills a critical gap in medical education, making high-fidelity clinical simulation more accessible and reducing costs. Clinical validation by experienced physicians enhances the reliability of our findings and their practicality in real-world settings. The large sample size and objective methods used for assessing learning and stress parameters further solidify the study’s conclusions. Future research should aim to expand the geographic and demographic scope of studies, exploring VR’s utility and effectiveness across a broader spectrum of medical training environments. Randomized controlled trials and longitudinal studies are needed to provide more definitive evidence of VR’s impacts on learning outcomes and clinical performance. ## CONCLUSION This study highlights the effectiveness of virtual reality in simulating the stress-related aspects of complex medical procedures, such as pericardiocentesis. It demonstrates VR’s capability to replicate both the technical challenges and emotional pressures found in real-life medical scenarios, thereby significantly contributing to stress management training in medical education. VR’s immersive environments are instrumental in preparing medical students for the psychological demands they will face in their careers. The secondary analysis reveals that demographic and lifestyle factors significantly influence stress responses during training, with VR’s adaptability allowing for personalized scenarios to effectively address these varied responses. However, the study’s limitations, including its nonrandomized design and focus on undergraduate students, necessitate cautious interpretation of the findings, and suggest limited generalizability to practicing healthcare professionals. In conclusion, VR is evolving as a crucial tool in medical education, especially for stress management. Future research should aim to validate VR’s effectiveness in diverse and realistic settings and explore its potential across different medical training stages. Further studies are needed to integrate comprehensive stress analysis within VR training, tailoring it to the specific needs of medical trainees worldwide ## Data Availability The datasets used and/or analyzed during the current study are available in publicly available repositories. ## DECLARATIONS ### Ethics approval and consent to participate This study received ethical approval from the HM Montepríncipe University Hospital Research Ethics Committee (Code: 18.12.1339.GHM) and involved 132 consenting medical students. All participants provided written informed consent, and their data were anonymized to protect confidentiality ### Consent for publication Not applicable ### Availability of data and materials The datasets used and/or analyzed during the current study are available in publicly available repositories. ### Competing interests The authors declare that they have no competing interests. ### Authors’ contributions ARL designed the study, collected the data, performed the statistical analysis, and wrote the study as the principal investigator. RG participated in the development and correction of the virtual scenario. LZ performed the programming of the virtual scenario according to the created model. ARN participated in the performance tests and in the design and collection of biometric parameters AGP participated in the design and adaptation of the simulated scenarios PCF participated in proofreading and supervision of the study design, statistical analysis and writing of the paper. ## Acknowledgments Not applicable * Received August 22, 2024. * Revision received August 22, 2024. * Accepted August 22, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Xu M, Yong-He Luo, Zhang Y, Xia R, Qian H, Zou X: Game-based learning in medical education. Frontiers in Public Health 2023, 11. 2. 2.Mahajan A, Hawkins A: Current Implementation Outcomes of Digital Surgical Simulation in Low- and Middle-Income Countries: Scoping Review. JMIR medical education 2022, 9:e23287. 3. 3. Hui-Chen Lin, Gwo-Jen Hwang, Chou KR, Chia-Kuang Tsai: Fostering complex professional skills with interactive simulation technology: A virtual reality-based flipped learning approach. British Journal of Educational Technology 2022, 54(2):622–641. 4. 4.Javaid M, Haleem A, Ravi PS, Suman R: 3D printing applications for healthcare research and development. Global health journal 2022, 6(4):217–226. 5. 5.Urlings J, Guido dJ, Thomas J.J. Maal, Dylan JHAH: Views on Augmented Reality, Virtual Reality, and 3D Printing in Modern Medicine and Education: A Qualitative Exploration of Expert Opinion. J Digital Imaging 2023, 36:1930–1939. 6. 6.Brunzini A, Papetti A, Messi D, Germani M: A comprehensive method to design and assess mixed reality simulations. Virtual Reality 2022, 26(4):1257–1275. 7. 7.Jennifer LM, Lee E, Roberts T, Garrett R: Virtual Simulation to Increase Self-efficacy Through Deliberate Practice. Nurse Educ 2020, 45(4):202–205. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F22%2F2024.08.22.24312406.atom) 8. 8.Martínez-Pascual, B., Ramírez-Adrados, A., Fernández-Martínez, S. et al. Autonomic stress response of physiotherapy student in the different scenarios of an objective structured clinical examination. BMC Med Educ 22, 811 (2022). 9. 9.Bolton, W., Aruparayil, N., Please, H., Patel, R., Dalton, J., Agery, B., Mwachiro, M., Ndizeye, O., Daya, L., Bains, L., Mishra, A., Jani, P., Culmer, P., Gnanaraj, J., Jayne, D., & Dhanda, J. O018 A virtual reality training platform for surgeons, anaesthetists, and obstetricians in low and middle-income countries. British Journal of Surgery 2022, 109(Supplement_4). 10. 10.Weissglass, D. E Contextual bias, the democratization of healthcare, and medical artificial intelligence in low- and middle-income countries. Bioethics 2022, 36(2):201–209. 11. 11.Hong, E. M., Hakimi, A. A., Ho, D., Torkian, B. A., & Wong, B. J. F. Evaluating Open Source Software for 3D Imaging and Morphing in Cosmetic and Reconstructive Surgery. The Laryngoscope, 2021, 131(2):299–303. 12. 12.Roh, T. H., Oh, J. W., Jang, C. K., Choi, S., Kim, E. H., Hong, C.-K., & Kim, S.-H. Virtual dissection of the real brain: integration of photographic 3D models into virtual reality and its effect on neurosurgical resident education. Neurosurgical Focus, 2021, 51(2), E16. 13. 13.Towers, A., Dixon, J., Field, J., Martin, R., & Martin, N. Combining virtual reality and 3D-printed models to simulate patient-specific dental operative procedures-A study exploring student perceptions. European Journal of Dental Education, 2022, 26(2):393–403 14. 14.Rochlen, L. R., Putnam, E., Levine, R., & Tait, A. R. Mixed reality simulation for peripheral intravenous catheter placement training. BMC Medical Education, 2022, 22(1): 876. 15. 15.Bruno, R. R., Wolff, G., Wernly, B., Masyuk, M., Piayda, K., Leaver, S., Erkens, R., Oehler, D., Afzal, S., Heidari, H., Kelm, M., & Jung, C. Virtual and augmented reality in critical care medicine: the patient’s, clinician’s, and researcher’s perspective. Critical Care, 2022, 26(1): 326. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13054-022-04202-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36284350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F22%2F2024.08.22.24312406.atom) 16. 16.Helman, T. J., Headrick, J. P., Stapelberg, N. J. C., & Braidy, N. The sex-dependent response to psychosocial stress and ischaemic heart disease. Frontiers in Cardiovascular Medicine, 2023, 10. 17. 17.J. Rajzer, M., de Timary, P., Weber, T., Weisser, B., … Januszewicz, A. Lifestyle, psychological, socioeconomic, and environmental factors, and their impact on hypertension during the coronavirus disease 2019 pandemic. Journal of Hypertension, 2021, 39(6): 1077–1089. 18. 18.Corrigan, S. L., Roberts, S., Warmington, S., Drain, J., & Main, L. C. Monitoring stress and allostatic load in first responders and tactical operators using heart rate variability: a systematic review. BMC Public Health, 2021, 21(1): 1701. 19. 19.Narciso, D., Melo, M., Rodrigues, S., Cunha, J. P., Vasconcelos-Raposo, J., & Bessa, M. Using Heart Rate Variability for Comparing the Effectiveness of Virtual vs Real Training Environments for Firefighters. IEEE Transactions on Visualization and Computer Graphics, 2023, 29(7): 3238–3250.