Photoplethysmography sampling frequency: pilot assessment of how low can we go to analyze pulse rate variability with reliability?

Pulse rate variability (PRV) analysis appears as the first alternative to heart rate variability analysis for wearable devices; however, there is a constraint on computational load and energy consumption for the limited system resources available to the devices. Considering that adjustment of the sampling frequency is one of the strategies for reducing computational load and power consumption, this study aimed to investigate the influence of sampling frequency (f(s)) on PRV analysis and to find the minimum sampling frequency while maintaining reliability. We generated 5000, 2500, 1000, 500, 250, 100, 50, 25, 20, 15, 10, 5 Hz down-sampled photoplethysmograms from 10 kHz-sampled PPGs and derived time-and frequency-domain variables of the PRV. These included AVNN, SDNN, SDSD, RMSSD, NN50, pNN50, total power, VLF, LF, HF, LF/HF, nLF and nHF for each down-sampled signal. Derived variables were compared with heart rate variability of the 10 kHz-sampled electrocardiograms, and then statistically investigated using one-way ANOVA test and Bland-Altman analysis. As a result, significant differences (P < 0.05) were found for SDNN, SDSD, RMSSD, NN50, pNN50, TP, HF, LF/HF, nLF and nHF, but not for AVNN, VLF and LF. Based on the post hoc tests, it was found that the NN50 and pNN50, SDSD and RMSSD, LF/HF and nHF, SDNN, TP and nLF analysis had significant differences at f(s) <= 20 Hz, f(s) <= 15 Hz, f(s) <= 10 Hz; f(s) = 5 Hz, respectively. In other words, a significant difference was not seen for any variable if the f(s) was greater than 25 Hz. Consequently, our pilot study suggests that analysis of variability in the time and frequency domain from pulse rate obtained through PPG may be potentially as reliable as that derived from the analysis of the electrocardiogram, provided that f(s) >= 25 Hz sampling frequency is used.