HRV4Training

Blog post by Marco Altini.
In the last post I covered the importance of context when interpreting HRV data, and pointed out that the optimal routine involves a morning HRV measurement, since this is the easiest way to isolate context and limit the impact of other stressors.
Other important aspects, beyond when to take the measurement, are the measurement duration, type (lying, standing or orthostatic), paced breathing, and frequency (how often?). In this post, I will cover measurement duration. In particular, I will reference to a recent study by Michael Esco and Andrew Flatt, and try to replicate and extend some of their work. Esco and Flatt (see the full paper here) were interested in determining if a $60$ seconds measurement is sufficient for valid assessment of HRV, in particular of $ln(rMSSD)$, a metric often used as a proxy to parasympathetic activity when monitoring training load. Their reference was the same metric ($ln(rMSSD)$) over a longer (standard) window of $5$ minutes. In their study population ($23$ collegiate athletes), $60$ seconds seemed sufficient to obtain HRV measures as good as the reference, while shorter measurements were a bit more problematic, especially the $10$ seconds one. However, no statistically significant difference was found between any time window. In my analysis I will look at windows of $10$ s, $30$ s, $60$ s and $5$ minutes, as done by Esco and Flatt, as well as for $3$ minutes recordings, since HRV4Trainings provides a choice between $60$ s, $3$ minutes and $5$ minutes, and in general, the shorter the more practical. For this analysis I will also look at rMSSD and not only $ln(rMSSD)$. ## Dataset
For this analysis I used a dataset I collected for another study I designed during my PhD. More details can be found here. Long story short, in my research I was measuring energy expenditure, and therefore I also recorded most of the data in the morning in laboratory settings. Participants were required to restrain from drinking (except for water), eating and smoking in the 2 hours preceding the experiment, similarly to what was required by Esco and Flatt. Data was collected using the ECG Necklace, a device that measures ECG and acceleration.
My participants were however mainly university employees or students, and therefore less fit than the collegiate athletes of the original study. For my protocol, participants went through a series of physical activities (lying, sedentary, household, walking at different speeds, running, biking - $VO_2$max test). Here I isolated the first 20 minutes of data for each participant and highlighted the "lying" part in the RR intervals time series, which I detected using accelerometer data (since the sensor is on the chest, this step is rather easy):
The only preprocessing here was the removal of consecutive RR intervals differing more than $20\%$ and discarding participants where the overall noise was more than $1\%$ of the data. This step was to avoid HRV features to be affected by ectopic beats or other motion artifacts or noise in general. I ended up with $25$ participants.
## Analysis of different time windows
As a next step I computed time and frequency domain HRV features over the RR intervals marked as "lying", so discarding all the remaining data. This is typically a $10$-$15$ minutes segment (see data highlighted in red above). Features are computed over windows of $10$ seconds to $5$ minutes (non-overlapping). Therefore we will end up with many $10$ s windows and just a few $5$ minutes windows (our reference). If I understood it correctly, in the original paper Esco and Flatt randomly picked one of these time window as a comparison. Here I decided to keep all of them so that I could look at their distribution. For the Bland-Altman plots, I computed the mean of the values over a certain window length and participant, and used that value as the comparison one, instead of selecting one randomly.
Let's look at rMSSD for different time windows and participants. Each color is a participant, and for each participant we see boxplots for the different conditions (i.e. time windows). By looking at boxplots we can easily see how the values are distributed and visually compare across time windows. Time windows are ordered by length, from $10$ s to $5$ minutes:
What can we derive from the plot above?
In general, we can see positive results, with rMSSD values being quite stable over different time windows. For each participant, the median of the boxplot (black line) doesn't differ much across different time windows. However, the variability (spread of the distribution) is quite high for $10$ s windows, while it seems to get already much better for $30$ s. It is important to note that HRV is never constant and even the reference measure (i.e. recording over 5 minutes) can change quite a bit over a $10$-$15$ minutes period. We can spot these differences in our $5$ minute reference window for the few participants were "lying" was longer than $10$ minutes, for example participant $107$. Similarly, plotting $ln(rMSSD)$:
Finally, we can look at bias and limits of agreement using Bland-Altman plots. First rMSSD:
and then the logarithmic version:
As we could expect, we have reduced bias with increased window length. An analysis of variance (ANOVA) shows no significant differences between conditions (i.e. between time windows), when using the mean over the different windows to represent the HRV value. The same was reported by Esco and Flatt, who also computed the intra-class correlation (ICC), and showed a decrease in ICC as the window size decreased. From the spread that we can see in the previous boxplots, I would also suggest longer windows to be preferable.
## Summary
The dataset used for this analysis was not collected with the intention of measuring HRV and the protocol was not exactly the same as the one of the original study by Esco and Flatt, however, I could replicate fairly well their findings. $60$ seconds seem to be a sufficient time window for morning HRV measurements even though, according to my data, the best results are shown when we use a $3$ minutes window.
While $60$ seconds showed limited spread and good agreement with the reference $5$ minutes, shorter measurements are probably not a good idea. No statistically significant differences were found in either the study by Esco and Flatt and my analysis, however, $10$ seconds is a very short window even for a rather stable measure such as rMSSD. Especially at rest conditions and in trained participants, $10$ seconds could be as short as one single breathing cycle. Clearly, there are many trade-offs to consider. Especially when using the camera version of HRV4Training, I strongly advise to use the short $60$ seconds time window, since performing the PPG measurement for more than $60$ seconds will increase the risk of moving the finger, and therefore affecting negatively the measurement due to motion artifacts.
0 Comments
Your comment will be posted after it is approved.
## Leave a Reply. |
This blog is curated byMarco Altini, founder of HRV4Training Blog Index The Ultimate Guide to HRV1: Measurement setup 2: Interpreting your data 3: Case studies and practical examples How To1. Intro to HRV 2. How to use HRV, the basics 3. HRV guided training 4. HRV and training load 5. HRV, strength & power 6. Overview in HRV4Training Pro 7. HRV in team sports HRV MeasurementsBest Practices 1. Context & Time of the Day 2. Duration 3. Paced breathing 4. Orthostatic Test 5. Slides HRV overview 6. Normal values and historical data Data Analysis1a. Acute Changes in HRV (individual level) 1b. Acute Changes in HRV (population level) 1c. Acute Changes in HRV & measurement consistency 1d. Acute Changes in HRV in endurance and power sports 2a. Interpreting HRV Trends 2b. HRV Baseline Trends & CV 3. Tags & Correlations 4. Ectopic beats & motion artifacts 5. HRV4Training Insights 6. HRV4Training & Sports Science 7. HRV & fitness / training load 8. HRV & performance 9. VO2max models 10. Repeated HRV measurements 11. VO2max and performance 12. HR, HRV and performance 13. Training intensity & performance 14. Publication: VO2max & running performance 15. Estimating running performance 16. Coefficient of Variation 17. More on CV and the big picture 18. Case study marathon training 19. Case study injury and lifestyle stress 20. HRV and menstrual cycle 21. Cardiac decoupling 22. FTP, lactate threshold, half and full marathon time estimates 23. Training Monotony Camera & Sensors1. ECG vs Polar & Mio Alpha 2a. Camera vs Polar 2b. Camera vs Polar iOS10 2c. iPhone 7+ vs Polar 2d. Comparison of PPG sensors 3. Camera measurement guidelines 4. Validation paper 5. Android camera vs Chest strap 6. Scosche Rhythm24 7. Apple Watch 8. CorSense 9. Samsung Galaxy App Features1. Features and Recovery Points 2. Daily advice 3. HRV4Training insights 4. Sleep tracking 5. Training load analysis 6a. Integration with Strava 6b. Integration with TrainingPeaks 6c. Integration with SportTracks 6d. Integration with Genetrainer 6e. Integration with Apple Health 6f. Integration with Todays Plan 7. Acute HRV changes by sport 8. Remote tags in HRV4T Coach 9. VO2max Estimation 10. Acute stressors analysis 11. Training Polarization 12. Lactate Threshold Estimation 13. Functional Threshold Power(FTP) Estimation for cyclists 14. Aerobic Endurance analysis 15. Intervals Analysis 16. Training Planning 17. Integration with Oura 18. Aerobic efficiency and cardiac decoupling Other1. HRV normal values 2. HRV normalization by HR 3. HRV 101 |