International Journal of Brain and Cognitive Sciences

p-ISSN: 2163-1840    e-ISSN: 2163-1867

2015;  4(1): 8-14

doi:10.5923/j.ijbcs.20150401.03

Performance in the Stroop Task and Simultaneously Recorded Heart Rate Variability before and after Meditation, Supine Rest and No - Intervention

Pailoor Subramanya 1, Shirley Telles 2

1Indian Council of Medical Research, Center for Advanced Research in Yoga and Neurophysiology, SVYASA, Bengaluru, India

2Patanjali Research Foundation, Patanjali Yogpeeth, Haridwar, Uttarakhand, India

Correspondence to: Shirley Telles , Patanjali Research Foundation, Patanjali Yogpeeth, Haridwar, Uttarakhand, India.

Email:

Copyright © 2015 Scientific & Academic Publishing. All Rights Reserved.

Abstract

In an earlier study heart rate variability (HRV) recordings after cyclic meditation (CM) suggested a shift towards vagal dominance. In a separate study, CM improved the performance in a task for attention. Generally, attention tasks result in sympathetic activation. The present study was intended to assess whether CM would change the autonomic status when practitioners performed an attention task, with simultaneous HRV recording Methods: Twenty-five male volunteers with group average age ± S.D., 23.9 ± 5.0 years were each assessed on three separate days in three sessions, (i) CM, (ii) supine rest (S R), and (iii) quiet sitting (QS). Each session was for 23 minutes. Simultaneous assessments were made before and after each of the three sessions and included the HRV, respiratory rate and performance in the Stroop task. Results: Following CM the LF power increased and HF power decreased (P < 0.05). There was an increase in the word scores, color scores and color-word scores of the Stroop task (P < 0.001, Repeated measures ANOVA, followed by Bonferroni adjusted post-hoc analyses). Conclusion: The results suggest that there was reduced parasympathetic activity after CM and possibly increased sympathetic activity, when participants simultaneously performed the Stroop task. Their improved performance in the Stroop task suggests better response inhibition, interference resolution, word reading and activation along with reduced physiological arousal based on HF power of the HRV.

Keywords: Stroop, HRV, Yoga

Cite this paper: Pailoor Subramanya , Shirley Telles , Performance in the Stroop Task and Simultaneously Recorded Heart Rate Variability before and after Meditation, Supine Rest and No - Intervention, International Journal of Brain and Cognitive Sciences, Vol. 4 No. 1, 2015, pp. 8-14. doi: 10.5923/j.ijbcs.20150401.03.

1. Background

Meditation has been described as a state of consciousness in which deep relaxation and increased internalized attention co-exist [1]. Meditation is not always easy for novices to practice, and in yoga texts meditation has been described as the seventh of eight stages required to reach ultimate mental transcendence, called Samadhi, in Sanskrit [2]. A meditation technique called CM, which has its’ origin in ancient yoga texts [3, 4], was derived especially for novices. CM has alternating cycles of practicing yoga postures slowly, with deep awareness, interspersed with periods of SR. CM differs from other meditation techniques as it includes the practice of yoga postures practiced slowly with heightened internal awareness. Tai-Chi-Qui-Gong is a comparable technique and these practices have been called ‘moving meditations’ [5].
Despite the fact that CM includes the practice of yoga postures, following twenty-three minutes of the practice there was a reduction in heart rate, as well as in the low frequency power of the HRV, while the HF power increased, suggesting a shift in cardiac autonomic balance towards vagal dominance after the practice (Cohen’s d = 0.87) [6]. The LF band of the HRV is believed to be chiefly but not entirely related to sympathetic activity when expressed in normalized units [7], while efferent vagal activity mainly contributes to the HF band [8].
Sympathetic activation is of interest in the present study since it is associated with increased vigilance and hence may be considered essential for performing an attentional task [9]. An early study showed improvement in different aspects of attention in practitioners of Transcendental Meditation [10].
More recently studies on practitioners of CM have shown improved performance in tasks requiring attention and other executive functions. In forty-two experienced meditators the P300 was measured before and after a session of CM on one day, and similarly before and after a comparable duration of SR on another day (Cohen’s d = 0.85) [11]. The P300 is an event related potential which reflects fundamental cognitive events requiring attentional and immediate memory- processes [12]. Following CM there was a decrease in the P300 peak latency and an increase in the P300 peak amplitude suggesting a greater improvement in cognitive processes after CM compared to SR [11].
Similarly the performance in a cancellation task was assessed in sixty-nine volunteers immediately before and after sessions of CM and SR on separate days [13]. Cancellation tasks assess attention, concentration, and visuo-spatial scanning abilities [14]. The net scores improved after both CM (by 24.9 percent) and SR (by 13.6 percent). These results also suggested that CM favorably impacts attention, concentration, visual scanning abilities and a repetitive motor response. It was speculated that the improvement in attention in these studies was related to a reduction in anxiety, as anxiety affects performance in tasks requiring attention [15].
However, as already mentioned attention is generally associated with increased sympathetic activity [16], which normally decreases when anxiety levels are low [17]. The association between HRV and trait anxiety was assessed in 92 healthy persons along with self-rating of trait anxiety and perceived emotional stress [18]. In an earlier study on healthy participants there was a shift towards vagal dominance in the HRV after CM [6]. Till now, the HRV has not been recorded in meditators in the same session as the attention task. Hence the present study was designed to record the HRV and respiratory rate before and after the Stroop task, in the same session [19]. The hypothesis being tested was that CM would improve Stroop task performance while reducing sympathetic activity.

2. Methods

Participants
There were twenty-five self-reported healthy male participants with ages ranging between 18 and 38 years (group mean ± S.D., 23.9 ± 5.0 years). The participants did not report any history of hospitalization or medication in the past six months. The temperature and blood pressure of the subjects were normal on the days of their participation in the study. The sample size was not determined prior to designing the study. However post-hoc analyses based on the effect size showed that the power of the tests used was 0.95 [20]. Participants were staying at a residential yoga c enter located in south India. For inclusion in the study participants had to meet the following criteria: (i) they had to be healthy based on a standard medical examination, (ii) literate and able to understand and perform the Stroop task, and (iii) have adequate experience of both CM and relaxation while supine i.e., 3 to 6 months. Exclusion criteria included: (i) color blindness based on the Ishihara charts, (ii) consumption of caffeine, nicotine containing substance s or any medicine which could influence cognition or autonomic function, and (iii) presence of extra systoles in the EKG. Participants were recruited by flyers on the notice board of the institution. The method of the study was explained to the participants, but not the hypothesis. No incentive was given to the participants for participation in the study. None of the participants was involved in any other ongoing yoga research at the center.
Design
Every participant was assessed in three types of practice sessions (CM, S R and QS sessions). The assessments were made on three different days for every participant. The sequence of practice was randomly allocated to balance the effect of the order of the sessions. The study design was explained to the participants and their signed informed consent was obtained. The study had the approval of the institutional ethics committee.
Time allocation within the sessions
The Stroop task, HRV and respiratory rate were measured ‘before’ and ‘after’ the interventions. The first 5 minutes of the ‘before’ period was in the sitting position, followed by 23 minutes of the ‘d uring’ period, where participants were given (i) CM practice using taped instructions, or (ii) SR in the corpse posture (= shavasana, in Sanskrit) without instructions or (iii) in the QS session participants were seated at ease. This was followed by 5 minutes of the ‘after’ period. In (i) ‘d uring’ periods of S R and QS session and (ii) ‘before’ and ‘after’ periods of all 3 sessions the participants were not given any specific instructions. The study design and assessment schedule are given schematically in Figure 1.
Figure 1. Schematic representation of the design of the study and assessment schedule
Assessments
The twenty-five participants were assessed before and after three types of sessions, (i) CM, (ii) SR, and (iii) QS. Assessments were made in a dimly lit, sound attenuated room (sound level 26 dB) with ambient temperature 24±2℃ while being observed throughout on a closed circuit television. The HRV and respiratory rate were recorded for 5 minutes using the Standard Limb Lead II configuration and volumetric pressure transducer respectively. The recordings were made immediately before the Stroop task prior to the intervention and after the Stroop task on completion of the intervention. There were no differences in results based on the order of the sessions, when the ‘after’ values were compared (P > 0.05, t-test for unpaired data).
Heart rate variability and Respiratory rate
HRV and respiratory rate were assessed using a MP100 BIOPAC data acquisition system (BIOPAC S ystems, Inc., U.S. A.). The EKG was recorded using Ag/AgC l pre-gelled electrodes. The recording was made with a standard limb Lead II configuration and an AC amplifier with 1.5 Hz low cut filter and 99 Hz high cut filter settings. The data were acquired at the sampling rate of 1024 Hz. Respiratory rate was recorded with a volumetric pressure transducer fixed around the trunk about 8 cm below the lower costal margin when the participants sat erect.
Stroop Task
The Stroop task consisted of a worksheet of 3 pages, each page having 100 items printed as 20 rows and 5 columns arranged randomly over a 8.5” 11” white sheet. The task used three colors, viz., b lue, green, and red. There were the standard Word page, Color page and Color- Word page. Instructions given were based on the Stroop task manual [21]. The participants were instructed to read the words, name the colors, and finally, name the color of the ink printed words which did not match, as quickly and as accurately as possible in the three tasks within 90 seconds for the complete task. All tasks were conducted and scored by a person who was unaware whether the participant was in a (i) CM or (ii) S R or (iii) QS session.
Inte rvention
Cyclic meditation session (CM)
During CM practice, participants kept their eyes closed and followed pre-recorded instructions. The instructions emphasized carrying out the practice slowly, with awareness and relaxation. The practice began by repeating a verse (0:40 min) from the yoga text, the Mandukya Upanishad [4]; followed by isometric contraction of the muscles of the body ending with supine rest (1:00 min.); slowly coming up from the supine position and standing at ease (called tadasana) and ‘balancing’ the weight on both feet, called centering (2:00 min.); then the first actual posture, bending to the right (ardhakaticakrasana, 1:20 min.); with 1:10 min. in tadasana for instructions about relaxation and awareness; bending to the left (ardhakaticakrasana, 1:20 min.); 1:10 min. as before; forward bending (padahastasana, 1:20 min.); another 1:10 min.; backward bending (ardhacakrasana, 1:20 min.); and slowly coming down to the supine position with instructions to relax different parts of the body in sequence (10:30 min.). The postures were practiced slowly, with awareness of all the sensations experienced. The total duration of the practice was 23 min. [6].
Supine rest session (SR)
SR was practiced as traditional shavasana (the corpse posture), which meant lying supine, with the legs apart, and arms away from the sides of the body, with the palms facing upwards, while the eyes were closed [22]. This practice lasted 23 min., so that the duration was the same as for CM. Participants were not given any instructions.
Quiet sitting session (QS)
Subjects were seated at ease and their thoughts were allowed to wander at random. They were not given any specific instructions except for the fact that they were told that they were to avoid meditating. The session was also at the same time of the day as the CM and SR sessions.
Thus in the present study, comparing the effects of practicing CM, SR and a QS session in the same subjects on separate days reduce d the possibility of bias. F urther, the order of the sessions (CM, S R and QS) was randomized to eliminate the effect of the order of sessions influencing the results.
Data extraction
Heart rate variability and Respiration
The frequency domain and time domain analys es of HRV data were carried out for 5 minute recordings for each of the following sessions CM, SR and QS. These 5 minutes epochs were recorded before and after each intervention. The data recorded were visually inspected off-line and only noise free data were included for analysis. None of the data had to be excluded due to noise. The heart rate was obtained based on R-R inter-beat interval analysis. The HRV power spectrum was obta ined using Fast Fourier transform analysis (FF T). The data were analyzed with an HRV analysis program developed by the Biosignal Analysis and Medical Imaging Group, University of Eastern F inland, Kuopio, F inland [23]. The energy in the HRV series in the following specific frequency bands was studied viz., the VLF band (0.0 - 0.04 Hz), LF band (0.04 - 0.15 Hz) and HF band (0.15 - 0.4 Hz). The LF and HF band values were expressed as normalized units. The following components of time domain HRV were analyzed: (i) mean RR interval (the mean of the intervals between adjacent QRS complexes or the instantaneous heart rate), (ii) RMSSD (root mean square of successive differences), (iii) NN50 (the number of interval differences of successive NN intervals greater than 50 ms), and (iv) pNN50 (the proportion of NN50 divided by total number of NN intervals). Mean respiratory rate was calculated before and after the sessions.
Stroop Task
As per the manual the correct responses were scored as ‘1’ depending on the response of the participant and if the responses were incorrect it was scored as ‘0’ [21]. Hence, the outcome measures were color scores, word scores, and color-word scores.
Data analysis
Statistical analysis was done using P ASW (Version 18.0). Repeated measures analyses of variance (ANOVA) were performed with two 'within subjects' factors, i.e., Factor 1: Sessions; CM, SR and QS and Factor 2: States; ‘before’, and ‘after’. This was followed by a post-hoc analysis with Bonferroni adjustment for multiple comparisons between the mean values of different states, and after- values were compared to the respective ‘before’ values. Data recorded in the before states of CM, S R, and QS sessions were compared for baseline differences with a one-way ANOVA.

3. Results

The group means and standard deviations for scores obtained in the HRV and respiratory rate are shown in Table 1, Stroop task [i.e., word task, color task, and color-word task scores] taken before and after CM, S R, and QS sessions are shown in Table 2.
Table 1. Group mean values ± S.D. of Frequency domain, Time domain analyses of the HRV and Respiratory rate recorded before and after, sessions of Cyclic Meditation (CM), Supine Rest (SR) and Quiet Sitting (QS)
Table 2. Group mean values ± S.D. for the Stroop Task recorded before and after sessions of Cyclic Meditation (CM), Supine Rest (SR), and Quiet Sitting (QS)
Repeated measures analyses of variance
Heart Rate Variability
Frequency domain analysis of the HRV spectrum
The repeated- measures analyses of variance (ANOVA) showed a significant difference between States (F = 5.75, df = 1, 24, Huynh-Feldt epsilon = 1, P < 0.05) for LF and HF. For VLF, LF, HF, and LF/HF, there were no significant differences between Sessions, or interaction between Sessions and S tates.
Time domain analysis of the HRV spectrum
There were no significant differences between Sessions, States or their interaction for mean RR, mean HR, RMSSD, NN50 and pNN50.
Respiratory Rate
For the respiratory rate there were no significant differences between Sessions, S tates or their interaction.
Stroop Task
The repeated- measures analyses of variance (ANOVA) of Stroop task scores showed a significant difference between (i) Sessions [for word scores, color scores and color-word scores] (F = 24.67, df = 2.26, 54.26, Huynh-Feldt epsilon = 0.28, P < 0.001) and (ii) S tates [Pre and Post] (F = 44.97, df = 1, 24, Huynh-Feldt epsilon = 1, P < 0.001). The interaction between Sessions and S tates was not significant for word scores, color scores or color-word scores.
Post-hoc analyses
Heart Rate Variability
Frequency domain analysis of the HRV spectrum
Post-hoc tests for multiple comparisons were performed with Bonferroni adjustment and all comparisons were made with the respective before states. After the Cyclic Meditation session, there was a significant increase in LF power (P < 0.05) and a decrease in HF power (P < 0.05). There were no significant differences after both supine rest and quiet sitting sessions.
Time domain analysis of the HRV spectrum
There was no significant difference between Sessions, States or their interaction for mean RR, mean HR, RMSSD, NN50 and pNN50, hence no further post-hoc analyses were performed.
Respiratory Rate
For the respiratory rate there were no significant differences between Sessions, S tates or their interaction, hence no further post-hoc analyses were performed.
Stroop Task
Post-hoc tests for multiple comparisons were performed with Bonferroni adjustment and all comparisons were made with the respective before states. After the CM session there was a significant increase in the word scores of the S troop task (P < 0.001), color scores (P < 0.001) and color-word scores (P < 0.001), compared to values of the before state. Also, after SR sessions there was a significant increase in word scores (P < 0.05) and after the QS session there was an increase in word scores (P < 0.05) and color scores (P < 0.01).

4. Discussion

Participants performed better in the Stroop task with increased color-word scores after CM but not after SR or QS. However the better performance was accompanied by an increased power in the LF component of the HRV, with a decrease in the HF component.
As already mentioned in the introduction the LF band of the HRV is believed to be chiefly but not exclusively due to sympathetic activity when expressed in normalized units [7], while efferent vagal activity contributes mainly to the HF band of the HRV [8]. The breath rate did not change after either of the interventions. This is relevant as it is known that there is an acute increase in LF, total spectrum HRV and in vagal baroreflex gain during slow breathing [24]. Training to increase the amplitude of respiratory sinus arrhythmia maximally increases the amplitude of heart rate oscillations at approximately 0.1 Hz [25]. However, considering that breath rate did not change in the present study, the breath frequency could not have influ enced the HRV. Hence the results suggest that after CM there was a shift in the autonomic balance towards parasympathetic decrease with a possible increase in sympathetic activity. As was already mentioned sympathetic activation is associated with increased vigilance considered necessary to perform an attentional task [6].
In an earlier study the LF band of the HRV reduced after CM, while the HF power increased [6]. However there were the following differences between that study and the present one: (i) in the earlier study participants were asked to practice CM and SR on two different days and the HRV and respiration were recorded. They were not given any task. In the present study participants were given the Stroop task to perform before and after the prac tice of CM or SR or QS, and the HRV and the respiration were recorded after the task. Hence, the earlier studies assessed practitioners of CM when they were doing the meditation without any task. The present results suggest that when given a task the physiologically relaxing effects of CM were not sufficient to reduce the sympathetic arousal associated with vigilance and performance of a task requiring attention in addition to other cognitive processes.
The absence of change in time domain measures which generally reflect parasympathetic tone, may be explained as follows: the sympathovagal balance is tonically and phasically modulated by at least three main factors, (i) central neural integration, (ii) peripheral inhibitory reflex mechanisms, and (iii) peripheral excitatory reflex mechanisms [26, 27]. While, variable phenomena such as the heart rate (recorded as the HRV) can be described as a function of time they may be better described as the sum of elementary oscillatory components [28]. Hence, the frequency domain may reflect those changes in the sympathovagal balance which are not detected by time domain analysis.
The S troop task requires the participant to suppress a habitual response and respond to an unusual one, in this case naming the color of the ink in which incongruously named color-words were printed [29]. The cognitive processes required to perform the Stroop task are complex. They include response inhibition, interference resolution, behavioral conflict resolution, word reading and production, and visual attention [30]. The present results suggest a possible benefit of practicing CM to improve abilities required to perform the Stroop task. Participants did not have any additional benefit of reduced sympathetic activation after performing CM and while doing the Stroop task. These results resemble those of a recent study in which Vipassana expert meditators showed greater P3b amplitudes to the target tones after meditation than they did both before meditation and after the no- meditation sessio n. They also simultaneously showed a larger LF/HF ratio during Vipassana meditation [31]. Hence Vipassana meditation was associated with better P3 potentials (suggestive of improved attention) along with increased sympathetic activity.
The limitations of the study are as follows: (i) in earlier studies [13, 32] the tasks used to assess performance before and after CM were different from the Stroop task used in the present study. The reason for this difference is that the participants in the present study had taken part in other studies in the residential yoga institution which used the tasks which were previously assessed. Hence, they would have approached the earlier tasks with a certain degree of familiarity and monotony. However in introducing a new task the results of the previous studies cannot be directly compared with the present results. (ii) All twenty- five participants were assessed in three types of sessions. Even though the order of the sessions was randomized the participants might have performed better after CM not merely because of the effects of CM per se, but because the CM session kept them busy and engaged while the other sessions did not. (iii) The sample size (n=25) in the present study was smaller than the sample sizes in the previous studies o n cancellation task performance (n=47), memory task performance (n=57) and HRV and respiration (n=47).
Despite these limitations the present results showed that when participants were given a task to perform following meditation, their perfor mance in the task which required focused attention was better, however there was also a shift in the autonomic balance towards decreased parasympathetic activity and possibly increased sympathetic activity.

ACKNOWLEDGEMENTS

The authors gratefully acknowledge the funding from the Department of Science and Technology (DS T), Government of India, as part of a grant (Project No. SR/CSI/22/2009) under Cognitive Science Research Initiative (CSI). The authors also acknowledge H.R. Nagendra P h.D. who derived the Cyclic Meditation technique from an ancient yoga text and Naveen K.V., P h.D., Co-P.I. of the project for his suggestions at different stages.

References

[1]  Murata, T., Takahashi, T., Hamada, T., Omori, M., Kosaka, H., Yoshida, H., & Wada, Y. (2004). Individual trait anxiety levels characterizing the properties of Zen meditation. Neuropsychobiology, 50(2), 189-194.
[2]  Taimini, I.K. (1986). The science of yoga. Madras, India: The Theosophical Publishing House.
[3]  Nagendra, H.R., & Nagarathna, R. (1997). New perspectives in stress management. Bangalore, India: Swami Vivekananda Yoga Prakashana.
[4]  Chinmayananda, S. (1984). Mandukya Upanishad. Bombay: Sachin P ublishers.
[5]  Chao, Y. F., Chen, S.Y., Lan, C., & Lai, J.S. (2002). The cardio respiratory responses and energy expenditure of Tai-C hi-Qui-Gong. The American Journal of Chinese Medicine, 30, 451-461.
[6]  Sarang, P., & Telles, S. (2006). Effects of two yoga based relaxation techniques on heart rate variability. International Journal of Stress Management, 13, 460-475.
[7]  Doğru, M.T., Başar, M.M., Yuvanç, E., Simşek, V., & Sahin, O. (2010). The relationship between serum sex steroid levels and heart rate variability parameters in males and the effect of age. Türk Kardiyoloji Derneği arşivi, 38, 459-465.
[8]  Malliani, A., Pagani, M., Lombardi, F., & Cerutti, S. (1991). Cardiovascular neural regulation explore in the frequency domain. Circulation, 84, 482-492.
[9]  Fredrikson, M., & Engel, B.T. (1985). Cardiovascular and electrodermal adjustments during a vigilance task in patients with borderline and established hypertension. Journal of Psychosomatic Research, 29(3), 235-246.
[10]  Peletier, K.R. (1972). Altered attention deployment in meditators. Berkeley, USA: Psychology Clinic, University of California.
[11]  Sarang, S. P., & Telles, S. (2006). Changes in P300 following two yoga-based relaxation techniques. International Journal of Neuroscience, 116(12), 1419-1430.
[12]  Polich, J. (1999). Electroencephalography: basic principles, clinical applications and related fields. In P300 in clinical applications. Edited by E. N iedermeyer, F. Lopes da S ilva. Baltimore- Munich: Urban and Schwarzenberg, 1073-1091.
[13]  Sarang, S.P., & Telles, S. (2007). Immediate effect of two yoga based relaxation techniques on performance in a letter-cancellation task. Perceptual and Motor Skills, 105(2), 379-385.
[14]  Uttl, B., & Pilkenton-Taylor, C. (2001). Letter cancellation performance across the adult life span. The Clinical Neuropsychologist, 15(4), 521-530.
[15]  Subramanya, P., & Telles, S. (2009). Performance on psychomotor tasks following two yoga-based relaxation techniques. Perceptual and Motor Skills, 109(2), 563-576.
[16]  Hamada, T., Murata, T., Takahashi, T., O htake, Y., Saitoh, M., K imura, H., Wada, Y., & Yoshida, H. (2006). Changes in autonomic function and EEG power during mental arithmetic task and their mutual relationship. Rinsho byori, 54(4), 329-334.
[17]  Van Honk, J., Hermans, E. J., d'Alfonso, A. A., Schutter, D. J., van Doornen, L., & de Haan, E. H. (2002). A left-prefrontal lateralized, sympathetic mechanism directs attention towards social threat in humans: evidence from repetitive transcranial magnetic stimulation. Neuroscience letters, 319(2), 99-102.
[18]  Dishman Rod, K., Nakamura, Yoshio., Garcia Melissa, E., Thompson Ray, W., Dunn Andrea, L., & Blair Steven, N. (2000). Heart rate variability, trait anxiety, and perceived stress among physically fit men and women. International Journal of Psychophysiology, 37(2), 121-133.
[19]  Hoshikawa, Y., & Yamamoto, Y. (1997). Effects of Stroop color-word conflict test on the autonomic nervous system responses. American Journal of Physiology, 272(3), H1113-H1121.
[20]  Faul, F., Erdfelder, E., Lang, A. G., & Buchner, A. (2007). G* Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behavior research methods, 39(2), 175-191.
[21]  Golden, C.J. (2003). Stroop Color and Word Test Adult’s Version: A manual for clinical and experimental uses. Wood Dale, IL: U.S.A. Stoelting Co.
[22]  Muktibodhananda, S. (2004). Hatha Yoga Pradipika. 2nd edition. Munger, Bihar, India: Yoga Publication Trust.
[23]  Tarvainen, M. P., Niskanen, J. P., Lipponen, J. A., Ranta-Aho, P. O., & Karjalainen, P. A. (2014). Kubios HRV–Heart rate variability analysis software. Computer methods and programs in biomedicine, 113(1), 210-220.
[24]  Lehrer, P. M., Vaschillo, E., Vaschillo, B., Lu, S. E., Eckberg, D. L., Edelberg, R., & Hamer, R. M. (2003). Heart rate variability biofeedback increases baroreflex gain and peak expiratory flow. Psychosomatic Medicine, 65(5), 796-805.
[25]  Lehrer, P.M., Vaschillo, E., & Vaschillo, B. (2000). Resonant frequency biofeedback training to increase cardiac variability: Rationale and manual for training. Applied Psychophysiology and Biofeedback, 25, 177-191.
[26]  Malliani, A. (1982). Cardiovascular sympathe tic afferent fibers. Reviews of Physiology, Biochemistry and Pharmacology, 94, 11-74.
[27]  Malliani, A., Pagani, M., & Lombardi, F. (1986). Positive feedback reflexes, in Zanchetti A, Tarazi RC (eds): Handbook of Hypertension: Volume 8. Pathophysiology of Hypertension. Amsterdam, Elsevier Science P ublishing Co, Inc, 69-81.
[28]  Malliani, A., Pagani, M., Lombardi, F., & Cerutti, S. (1991). Cardiovascular neural regulation explored in the frequency domain. Circulation, 84(2), 482-492.
[29]  Davidson, D. J., Zacks, R.T., & Williams, C.C. (2003). Stroop interference, practice, and aging. Neuropsychology, development, and cognition. Section B, Aging, neuropsychology and cognition, 10(2), 85-98.
[30]  Adleman, N. E., Menon, V., Blasey, C. M., White, C. D., Warsofsky, I. S., Glover, G. H., & Reiss, A. L. (2002). A developmental fMRI study of the Stroop color-word task. Neuroimage, 16(1), 61-75.
[31]  Delgado-Pastor, L. C., Perakakis, P., Subramanya, P., Telles, S., & Vila, J. (2013). Mindfulness (Vipassana) meditation: Effects on P3b event-related potential and heart rate variability. International Journal of Psychophysiology, 90(2), 207-214.
[32]  Subramanya, P., & Telles, S. (2009). Effect of two yoga-based relaxation techniques on memory scores and state anxiety. Biopsychosocial Medicine, 3(1), 8-12.