Investigation of the Effect of eSport on HRV Signal by Using Poincaré Plot Analysis

Investigation of the Effect of eSport on HRV Signal by Using Poincaré Plot Analysis

Serdar Gündoğdu Ömer Halil Çolak Övünç Polat 

Bergama Vocational School, Department of Computer Technologies, Dokuz Eylül University, İzmir 35700, Turkey

Faculty of Engineering, Department of Electrical and Electronics Engineering, Akdeniz University, Antalya 07058, Turkey

Corresponding Author Email: 
ovuncpolat@akdeniz.edu.tr
Page: 
589-594
|
DOI: 
https://doi.org/10.18280/ts.390220
Received: 
10 April 2020
|
Accepted: 
12 November 2021
|
Published: 
30 April 2022
| Citation

OPEN ACCESS

Abstract: 

In this work, in order to determination of effect of a puzzle video game on HRV (Heart rate variability) signal, Poincare plot analysis is used. The linear and nonlinear dynamics of HRV were evaluated using HRV recordings obtained from each player at the time of play and rest. For this purpose, we have correlated the Poincaré plot descriptors with standard HRV measures derived from time- and frequency-domain methods. In addition, the correlations between the descriptors of the Poincaré plot along with the frequency and time domain measures of HRV have evaluated for the game stage. Statistically significant values were observed in eSport players in gaming and during resting. The obtained results show that play and rest can be separated by using the features obtained from the Poincare plot.

Keywords: 

heart rate variability, Poincaré plot analysis, eSport

1. Introduction

To date, competitive video gaming has not been widely studied in the literature [1]. Professional competitive video games, also known as electronic sports (eSports), consist of professional players who regularly compete in tournaments and are managed by the gaming community [2]. In recent years, the popularity of eSports events has grown tremendously [3], especially among youth. Video gamers may be happy with eSports events. But especially when teenagers play for a long time, the person may show pathological symptoms. Even situations such as depression, loneliness, and anxiety can be seen in addicted teens [4]. Due to the increase in popularity and its positive and negative effects, research interest in eSports has increased. Physical and psychological signals have been ever increasingly used for gaming studies. The most popular biomedical signals studied are electroencephalography, cardiac activity, facial electromyography, event-related potentials, skin conductance level, and functional magnetic resonance imaging [5].

In this work, physiological responses to the eSport game were measured using cardiac activity signals. The heart rate is the heartbeats number per minute, calculated by measuring the time interval between two R peaks in the electrocardiogram (ECG) signal. The heart rate is influenced by various facts, such as stress conditions and exercise [6, 7].

HRV is the phenomenon of the variation in the time intervals between sequential heartbeats [8]. The HRV analysis can be performed as long-term, short-term and ultra-short-term [8, 9]. While short recordings are generally more useful when the acute effect of any factor is wanted to be studied, long-term recordings are preferred for cardiac health evaluations [9].

There are many techniques including, frequency domain, time-domain and non-linear methods for the analysis of HRV. Time domain analysis derives HRV using either statistical or geometrical methods. The root mean square of successive differences (RMSSD) and the standard deviation of NN intervals (SDNN) values are obtained by using the time-domain technique. The method of frequency domain uses results at different frequency bandwidths by converting the NN interval using the FFT (fast Fourier transform). The used ranges in this technique are the VLF (0.003 Hz – 0.04 Hz), the LF (0.04 Hz –0.15 Hz), the HF (0.15 Hz –0.4 Hz), and the LF/HF ratio.

The HRV is closely associated with the autonomic nervous system and reflects the influences of sympathetic and parasympathetic nervous systems in the sinoatrial node. The RMSSD [10, 11], and the HF indicates the degree of the parasympathetic activity, while the LF reflects the degree of sympathetic activity [12]. The power of HF is highly correlated to the RMSSD [13, 14]. The ratio of LF/HF represents the balance among the sympathetic and parasympathetic nerves. The sympathetic nervous system activity conduces to LF power, and parasympathetic nervous system activity mainly conduces to HF power [15].

The Poincaré plot, a nonlinear technique used to visually search for patterns buried from within a time series, is a scatter graph by plotting each RR interval against the previous interval. This analysis is insensible to changes in trends in the R-R intervals differently from frequency-domain measurements [16]. Four nonlinear measurements as SD1, SD2, SD2 / SD1 and S were derived using the Poincaré plot technique in this study. The ellipse area (S) representing total HRV is correlated with BRS (baroreflex sensitivity), power of LF and HF, and RMSSD. SD1 measures short-term HRV (in ms) correlated with BRS, and the power of HF. The RMSSD is identical to the nonlinear parameter SD1, which reflects short-term HRV [17]. SD2 measures short-term and long-term HRV (in ms) correlated with the power of LF [18-21].

The Poincare plot is frequently used to study nonlinear dynamics. The plot simpler to study the dynamics of heart variability [22]. The visual analysis of the Poincare plot is referred to as a qualitative analysis technique [23]. In this study, the heart rate variability of eSports players was visualized and analyzed on the Poincare graph. It has been observed that the eSports game's Poincaré graphics effectively characterize the change on the players' cardiac function. The differences in HRV pattern for eSports game are shown from nonlinear analysis of Poincare parameters.

HRV analyzes can be used to evaluate the acute physiological effects of activity in the field of computer games, and the obtained results can be used to determine the effects of parameters such as gender, age, previous experience differences in the same game and the type of game.

Hamm et al. [24] found that importantly reduced HRV indices are primarily indicative of vagal withdrawal. Their results provided additional information regarding the physiological response of emotional stress to HRV measurements. Polli et al. [25] explored the relationship between oxidative stress and parasympathetic activity. They have found that the more the increase in parasympathetic activity during exercise, the more the reduction in oxidative stress levels (inversely related). Peçanha et al. [26] pointed out that RMSSD, SDNN, and HF parameters can be used to determine parasympathetic activity. Concerning LF, no consistent response to stress compared to previous studies. While some studies in the literature observed a rise of LF during stress, some other studies indicated an opposite trend [27]. In another study, Chen et al. [28] suggested that Poincaré plot (SD1) is a very good indicator candidate for stress level assessment.

Yeo et al. [12] reported a study evaluating the psychological effects of games on users by analyzing biomedical signals. As a result of their analysis on pulse transit time, skin temperature and, heart rate variability, they showed that sympathetic nerve activities during computer games were more active than parasympathetic activity. Shapiro et al. [29] aimed to examine the impact of the unfair game on different participants. During the game, ECGs were recorded for measuring HR. The results reflected an increase in stress while an unfair game playing. Haghshenas and Pooyan [30] showed that during playing computer games, the values of the nonlinear parameters (such as the SD1/SD2 ratio) were reduced. Some researchers in the literature [31-35] have asserted that in resting states, the LF band reflects activity of baroreflex, not sympathetic innervation of cardiac [36].

The aim of our work was to assess the linear and nonlinear dynamics of HRV in eSport players. We aimed to (1) evaluate four descriptors of the Poincaré plot, (2) to compare and (3) to contrast them with conventional parameters of HRV. To that purpose, we correlated Poincaré plot descriptors with standard HRV parameters obtained from time domain and frequency domain analysis.

2. Material and Methods

A study was carried out to investigate by using Poincare plot analysis of HRV during gaming and resting. In the experiment, nine adults (8 males and 1 female) with a mean of 33±7.6 years, ranging from 24 to 41 years took.

All volunteers provided written informed consent before participation. This study was approved by the Medical Ethical Committee of the Akdeniz University and the experiment was undertaken in compliance with national legislation and the Declaration of Helsinki.

The R-R intervals were records using a heart rate monitor (Polar V800) during the resting and gaming.

The Poincare plot analysis could be used in order to graphically define the beat to beat variability of R-R intervals. It is built by plotting every R-R interval as a function of the immediately previous one. R-R interval time series is indicated by the RRn.

The values of SD1 and SD2 represented the second moment measure of the distribution of points forming the Poincaré plot. The SD1 and SD2 can be described as follows [36]:

$S D 1=\frac{1}{N-1} \sum_{i=1}^{N-1}\left(D_{i}^{\prime}\right)^{2}$        (1)

$S D 2=\frac{1}{N-1} \sum_{i=1}^{N-1}\left(D_{i}^{\prime \prime}\right)^{2}$        (2)

where, D′ and D″ denote the distances of any point from the identity line (RRn=RRn+1) and the line perpendicular to the identity line (RRn+RRn+1=2RR). These operators can be expressed as follows [37]:

$D_{i}^{\prime}=\frac{1}{\sqrt{2}}\left|R R_{n}-R R_{n+1}\right|$        (3)

$D_{i}^{\prime \prime}=\frac{1}{\sqrt{2}}\left|R R_{n}+R R_{n+1}-2 \overline{R R}\right|$       (4)

S represents area of the ellipse and can be formulated as:

$S=\pi * S D 1 * S D 2$          (5)

SPSS Statistics 24.0 (IBM Corporation) software was used for the statistical analysis. The results are presented as mean. Wilcoxon signed-rank test in this study were used to determine if there was a significant statistical difference before and after the game. Correlations were evaluated with the Spearman's rank correlation method.

3. Results and Discussion

The nine volunteers who played the puzzle video game at sitting position. HRV data of the participants were recorded during the game and rest sessions. The linear and nonlinear dynamics of HRV were evaluated using HRV recordings obtained from each player. Figure 1 shows the instance of the graphic of change of the RR of the any volunteer during rest and game state. According to this figure, the RR values of the volunteer were decreased during the game. This means that the heart rate is accelerated during the game.

Figure 1. The change of the RR of the any volunteer during (a) rest and (b) game state

Figure 2 shows the Poincare plot with proper delay of HRV signals in states of rest and game for a participant. SD1, and S were importantly reduced (p≤ 0.05) in game state compared to rest state.

The resulting boxplots of the linear measures of HRV indices parameters are shown in Figure 3.

We used Wilcoxon signed‐rank tests to examine possible differences in linear measures of HRV parameters (used mean values) between rest and game states. By taking six HRV parameter variables (HR, SDNN, RMSSD, LF, HF and ratio of LF/HF) as independent and test scores as dependent variables, Wilcoxon signed‐rank has been conducted and obtained results are shown in Table 1.

Table 1. Comparison of linear measures of HRV indices parameters between rest and game status

Parameters

Rest

Game

p-value

HR

80.79

85.97

0.050*

SDNN

42.62

29.07

0.051

RMSSD

33.54

23.12

0.021*

LF

1646.11

172.82

0.028*

HF

785.08

270.73

0.008**

LF/HF

2.58

2.40

0.678

*significant at p≤0.05, **highly significant at p≤0.01.

Figure 2. Poincare plot of HRV signals in rest and game state

Figure 3. The box-plots of the linear measures of HRV indices parameters

Figure 4. The box-plots of the nonlinear measures of HRV indices parameters

As can be seen from the obtained results in Table 1, it was observed as other than SDNN and LF/HF, the other four variables (mean-HR, RMSSD, LF and HF) indicated significant difference.

Shapiro et al. [29] reported, there is an increase in stress while game playing (for unfair game on different participants).

The previous studies showed that RMSSD, SDNN as well as HF decreased during stress [27]. In our study, RMSSD and HF power also were decreased. Although games can evoke different emotions in individuals, in our study, considering the average of the volunteers, puzzle video game increased stress.

The resulting boxplots of the nonlinear measures of HRV indices parameters are shown in Figure 4.

By taking four mean of nonlinear measures of HRV parameters variables (the area of the Poincaré plots (S), SD1, SD2 and SD2/SD1) as independent and test scores as dependent variables, Wilcoxon signed‐rank has been applied and obtained results are shown in Table 2.

Table 2. Comparison of nonlinear measures of HRV indices parameters between rest and game status

Parameters

Rest

Game

p-value

S

6336.31

2466.58

0.028*

SD1

23.87

16.46

0.021*

SD2

54.76

36.96

0.051

SD2/SD1

2.54

2.81

0.767

*significant at p≤0.05, **highly significant at p≤0.01.

The use of the Poincaré plot allows us to visualize all points defined by consecutive RR intervals and easily identify points corresponding to outliers in heart rhythm fluctuations. The outliers do not affect the Poincaré plot as strongly as conventional time and frequency domain measurements of HRV data [38].

From Table 2 it was observed as other than SD2 and SD2/SD1, the other two variables (S and SD1) showed significant difference at 0.05 level. Haghshenas and Pooyan (2019) [30] and the current results showed that during playing computer games, the nonlinear parameters’ values were reduced. In our study, the reduction of the nonlinear parameter values of SD1 and S during the game was parallel to this study.

The acquired correlation results between the descriptors of the Poincaré plot of RR intervals (SD1, SD2, SD2/SD1 and S) and the mean RR interval in conjunction with the frequency domain and time domain measures of HRV are presented in Table 3. This correlation analysis was done by using Spearman's rank correlation method.

SD1 parameter correlated significantly with all HRV measurements studied (p≤0.05 and p≤0.01). There were high significantly correlations, all with p≤0.05 and p≤0.01, among SD2, S and all measures of HRV with the exception of LF/HF. No significant correlation between SD2/SD1 score and variables of interest was found.

There were significant negative correlations of HR with SD1 (p =0.007), SD2 (p = 0.010), S (p = 0.013). According to acquired correlation results, the RMSSD metric is almost fully correlated with the nonlinear metric SD1 (r=1, p≤0.001) at game status and this result is similar to previous study [17].

Whether SD1, SD2, SD2 / SD1, S and other parameters we determined are related to the short-term variability of heart rate and also whether they are related to the parasympathetic response can be associated with the results we have obtained.

Table 3. Correlation analysis of nonlinear measures with linear measures at game status

Parameters

SD1

SD2

SD2/SD1

S

r

p

r

p

r

p

r

p

RR

 0.845

0.004**

0.828

0.006**

-0.460

0.213

0.820

0.007**

HR

-0.817

0.007**

-0.800

0.010**

0.467

0.205

-0.783

0.013*

SDNN

 0.833

0.005**

0.983

0.000**

-0.117

0.765

0.933

0.000**

RMSSD

 1.000

0.000**

0.783

0.013*

-0.550

0.125

0.950

0.000**

LF

 0.783

0.013*

0.900

0.001**

-0.200

0.606

0.850

0.004**

HF

 0.950

0.000**

0.733

0.025*

-0.533

0.139

0.917

0.001**

LF/HF

-0.667

0.050*

-0.333

0.381

0.400

0.286

-0.617

0.077

**Correlation is highly significant at the 0.01 level (2-tailed). *Correlation is significant at the 0.05 level (2-tailed).

4. Conclusions

This work was designed to evaluate dynamics of HRV using Poincare plot in persons playing a puzzle video game. According to obtained results from linear analysis, mean-HR, RMSSD, LF and HF have significant in separation of rest and game state.

In nonlinear analysis, S and SD1 parameters are important in separation of rest and game status. SD1, SD2 and S, which are non-linear HRV measurements, are strongly correlated, except LF/HF from linear measurements. Obtained results show that linear and nonlinear dynamics of HRV has been a non-invasive marker in playing a puzzle video game. As a result, playing computer games, comparative values were obtained based on the reflection of heart rate. Using the proposed approach, the extent to which a Poincaré plot descriptors is related with low and high frequency power is determined, thus determining correlation among linear and non-linear analysis and frequency-domain spectral analysis techniques for puzzle video players.

With the nonlinear methods used in HRV analysis, a popular and non-invasive tool for evaluating cardiovascular activities, it is very useful for determining the effects and extent of esports games on the cardiovascular system. This proposed approach can be thought of as a practical way of evaluating the effects of esports play on autonomic control specifically or as a marker of sympathetic activity in terms of the cardiac effects of the related play activity.

In future studies, it can be examined whether the physiological responses that occur in the heart during playing a video game are individual and whether they are relatively independent from the performance of the game. In addition, the effects of playing video games in different body positions (as standing, lying positions etc.) on HRV and the process of returning these changes to normal in HRV can be examined.

Acknowledgment

We acknowledge the volunteers for their contribution to this study. This study was supported by the Research Project Department of Akdeniz University, Antalya, Turkey.

Ethic Statement

The study was approved by the Clinical Research Ethics Committee of the Akdeniz University.

  References

[1] Griffiths, M. (2017). Work stressed and play? A brief look at competitive gaming. The Neuropsychotherapist, 5(5): 49-54. https://doi.org/10.4103/ipj.ipj_195_20

[2] Griffiths, M. (2017). The psychosocial, impact of professional gambling, professional video gaming & eSports, Casino & Gaming International, 28: 59-63.

[3] Palanichamy, T., Sharma, M.K., Sahu, M., Kanchana, D.M. (2020). Influence of Esports on stress: A systematic review. Industrial Psychiatry Journal, 29(2): 191-199. https://doi.org/10.4103/ipj.ipj_195_20

[4] Sarda, E., Bègue, L., Bry, C., Gentile, D. (2016). Internet gaming disorder and well-being: A scale validation. Cyberpsychology, Behavior, and Social Networking, 19(11): 674-679. https://doi.org/10.1089/cyber.2016.0286

[5] Kivikangas, J.M., Chanel, G., Cowley, B., Ekman, I., Salminen, M., Järvelä, S., Ravaja, N. (2011). A review of the use of psychophysiological methods in game research. Journal of Gaming & Virtual Worlds, 3(3): 181-199. https://doi.org/10.1386/jgvw.3.3.181_1

[6] Freeman, J.V., Dewey, F.E., Hadley, D.M., Myers, J., Froelicher, V.F. (2006). Autonomic nervous system interaction with the cardiovascular system during exercise. Progress in Cardiovascular Diseases, 48(5): 342-362. https://doi.org/10.1016/j.pcad.2005.11.003

[7] Taelman, J., Vandeput, S., Spaepen, A., Huffel, S.V. (2009). Influence of mental stress on heart rate and heart rate variability. In 4th European Conference of the International Federation for Medical and Biological Engineering, pp. 1366-1369. https://doi.org/10.1007/978-3-540-89208-3_324

[8] Shaffer, F., Ginsberg, J.P. (2017). An overview of heart rate variability metrics and norms. Frontiers in public health, 258: 1-17. https://doi.org/10.3389/fpubh.2017.00258

[9] Dicle, A.R.A.S., Karakoç, B., Mitat, K.O.Z. (2014). Egzersiz yapan yetişkinlerde 1 saatlik koşu sonrasında KHD'deki 48 saatlik değişimin incelenmesi. SPORMETRE Beden Eğitimi ve Spor Bilimleri Dergisi, 12(1): 35-42. https://doi.org/10.1501/Sporm_0000000251

[10] Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology. (1996). Heart rate variability: Standards of measurement, physiological interpretation, and clinical use. Circulation, 93(5): 1043-1065. https://doi.org/10.1161/01.CIR.93.5.1043

[11] Nolan, J., Batin, P.D., Andrews, R., Lindsay, S.J., Brooksby, P., Mullen, M., Fox, K.A. (1998). Prospective study of heart rate variability and mortality in chronic heart failure: Results of the United Kingdom heart failure evaluation and assessment of risk trial (UK-heart). Circulation, 98(15): 1510-1516. https://doi.org/10.1161/01.cir.98.15.1510

[12] Yeo, M., Lim, S., Yoon, G. (2018). Analysis of biosignals during immersion in computer games. Journal of Medical Systems, 42(1): 1-9. https://doi.org/10.1007/s10916-017-0860-y

[13] Kleiger, R.E., Stein, P.K., Bigger Jr, J.T. (2005). Heart rate variability: Measurement and clinical utility. Annals of Noninvasive Electrocardiology, 10(1): 88-101. https://doi.org/10.1111/j.1542-474X.2005.10101.x

[14] Bigger Jr, J.T., Albrecht, P., Steinman, R. C., Rolnitzky, L.M., Fleiss, J.L., Cohen, R.J. (1989). Comparison of time-and frequency domain-based measures of cardiac parasympathetic activity in Holter recordings after myocardial infarction. The American Journal of Cardiology, 64(8): 536-538. https://doi.org/10.1016/0002-9149(89)90436-0

[15] Shaffer, F., McCraty, R., Zerr, C.L. (2014). A healthy heart is not a metronome: an integrative review of the heart's anatomy and heart rate variability. Frontiers in Psychology, 5: 1040. https://doi.org/10.3389/fpsyg.2014.01040

[16] Motie-Nasrabadi, A., Behbahani, S., Dabanloo, N.J. (2013). Ictal heart rate variability assessment with focus on secondary generalized and complex partial epileptic seizures. Advances in Bioresearch, 4(1): 50-58.

[17] Ciccone, A.B., Siedlik, J.A., Wecht, J.M., Deckert, J.A., Nguyen, N.D., Weir, J.P. (2017). Reminder: RMSSD and SD1 are identical heart rate variability metrics. Muscle & Nerve, 56(4): 674-678. https://doi.org/10.1002/mus.25573

[18] Brennan, M., Palaniswami, M., Kamen, P. (2001). Do existing measures of Poincare plot geometry reflect nonlinear features of heart rate variability? IEEE Transactions on Biomedical Engineering, 48(11): 1342-1347. https://doi.org/10.1109/10.959330

[19] Brennan, M., Palaniswami, M., Kamen, P. (2002). Poincare plot interpretation using a physiological model of HRV based on a network of oscillators. American Journal of Physiology-Heart and Circulatory Physiology, 283(5): H1873-H1886. https://doi.org/10.1152/ajpheart.00405.2000 

[20] Tulppo, M.P., Makikallio, T.H., Takala, T.E., Seppanen, T.H.H.V., Huikuri, H.V. (1996). Quantitative beat-to-beat analysis of heart rate dynamics during exercise. American Journal of Physiology-Heart and Circulatory Physiology, 271(1): H244-H252. https://doi.org/10.1152/ajpheart.1996.271.1.H244

[21] Tulppo, M.P., Makikallio, T.H., Seppänen, T., Laukkanen, R.T., Huikuri, H.V. (1998). Vagal modulation of heart rate during exercise: Effects of age and physical fitness. American Journal of Physiology-Heart and Circulatory Physiology, 274(2): H424-H429. https://doi.org/10.1152/ajpheart.1998.274.2.H424

[22] Ghatak, S.K., Aditya, S. (2018). Poincaré parameters and principal component analysis of Heart rate variability of subjects with health disorder. arXiv preprint arXiv:1802.10289. https://doi.org/10.48550/arXiv.1802.10289

[23] Khandoker, A.H., Karmakar, C., Brennan, M., Palaniswami, M., Voss, A. (2013). Poincaré plot methods for heart rate variability analysis. New York: Springer. https://doi.org/10.1007/978-1-4614-7375-6

[24] Hamm, W., Bogner-Flatz, V., Bauer, A., Brunner, S. (2020). FIFA World Cup 2018: Effect of emotional stress on conventional heart rate variability metrics. Clinical Research in Cardiology, 109(2): 266-270. https://doi.org/10.1007/s00392-019-01533-8

[25] Polli, A., Van Oosterwijck, J., Nijs, J., Marusic, U., De Wandele, I., Paul, L., Ickmans, K. (2019). Relationship between exercise-induced oxidative stress changes and parasympathetic activity in chronic fatigue syndrome: an observational study in patients and healthy subjects. Clinical Therapeutics, 41(4): 641-655. https://doi.org/10.1016/j.clinthera.2018.12.012

[26] Peçanha, T., Meireles, K., Pinto, A.J., Rezende, D.A.N., Iraha, A.Y., Mazzolani, B.C., Roschel, H. (2021). Increased sympathetic and haemodynamic responses to exercise and muscle metaboreflex activation in post-menopausal women with rheumatoid arthritis. The Journal of Physiology, 599(3): 927-941. https://doi.org/10.1113/JP280892

[27] Castaldo, R., Melillo, P., Bracale, U., Caserta, M., Triassi, M., Pecchia, L. (2015). Acute mental stress assessment via short term HRV analysis in healthy adults: A systematic review with meta-analysis. Biomedical Signal Processing and Control, 18: 370-377. https://doi.org/10.1016/j.bspc.2015.02.012

[28] Chen, Y., Zhang, L., Zhang, B., Chang’an, A.Z. (2020). Short-term HRV in young adults for momentary assessment of acute mental stress. Biomedical Signal Processing and Control, 57: 101746. https://doi.org/10.1016/j.bspc.2019.101746

[29] Shapiro, M.S., Rylant, R., de Lima, A., Vidaurri, A., van de Werfhorst, H. (2017). Playing a rigged game: Inequality's effect on physiological stress responses. Physiology & Behavior, 180: 60-69. https://doi.org/10.1016/j.physbeh.2017.08.006

[30] Haghshenas, L., Pooyan, M. (2019). Investigating the effect of computer games on the heart signal. In 2019 27th Iranian Conference on Electrical Engineering (ICEE), pp. 1775-1778. 10.1109/IranianCEE.2019.8786686

[31] Tiller, W.A., McCraty, R., Atkinson, M. (1996). Cardiac coherence: A new, noninvasive measure of autonomic nervous system order. Alternative Therapies in Health and Medicine, 2(1): 52-65. 

[32] Eckberg, D.L. (1997). Sympathovagal balance: A critical appraisal. Circulation, 96(9): 3224-3232. https://doi.org/10.1161/01.CIR.96.9.3224

[33] Porges, S.W. (2007). The polyvagal perspective. Biological Psychology, 74(2): 116-143. https://doi.org/10.1016/j.biopsycho.2006.06.009.

[34] Rahman, F., Pechnik, S., Gross, D., Sewell, L., Goldstein, D.S. (2011). Low frequency power of heart rate variability reflects baroreflex function, not cardiac sympathetic innervation. Clinical Autonomic Research, 21(3): 133-141. https://doi.org/10.1007/s10286-010-0098-y.

[35] Heathers, J.A. (2012). Sympathovagal balance from heart rate variability: An obituary. Experimental Physiology, 97(4): 556-556. https://doi.org/10.1113/expphysiol.2011.063867

[36] Shaffer, F., McCraty, R., Zerr, C. L. (2014). A healthy heart is not a metronome: An integrative review of the heart's anatomy and heart rate variability. Frontiers in Psychology, 5: 1040. https://doi.org/10.3389/fpsyg.2014.01040

[37] Karmakar, C.K., Gubbi, J., Khandoker, A.H., Palaniswami, M. (2010). Analyzing temporal variability of standard descriptors of Poincaré plots. Journal of Electrocardiology, 43(6): 719-724. https://doi.org/10.1016/j.jelectrocard.2010.09.001

[38] Guzik, P., Piskorski, J., Krauze, T., Schneider, R., Wesseling, K.H., Wykretowicz, A., Wysocki, H. (2007). Correlations between Poincaré plot and conventional heart rate variability parameters assessed during paced breathing. The Journal of Physiological Sciences, 57(1): 63. https://doi.org/10.2170/physiolsci.RP005506