Specifics of the HRV analysis have been supplied previously16 (link),17 (link) and were in accordance with the directives from the Task Force13 (link).
To study the parasympathetic regulation of HR we inspected the high frequency band of spectral analysis (HF: 0.15 Hz to 0.4 Hz) in absolute units in the frequency domain, RMSSD (root-mean square of differences between adjacent normal RR intervals) in the time domain and the SD1 Poincaré plot (standard deviation of the instantaneous variability of the beat-to-beat heart rate). We applied the Kubios® HRV v. 2.0 software to compute these indices18 .