Analysing heart rate variability in preterm infants: the effect of temporal adjustment of NN peaks and missing data

The measurement of heart rate variability (HRV) in preterm infants provides important information on function to clinicians. Measuring the underlying electrocardiogram (ECG) in the neonatal intensive care unit is a challenge and there is a trade off between extracting accurate measurements of the HRV and the amount of ECG processed due to contamination. Knowledge on the effects of 1) quantization in the time domain and 2) missing data on the calculation of HRV features will inform clinical implementation. In this paper, we studied multiple 5 minute epochs from 148 ECG recordings on 56 extremely preterm infants. We found that temporal adjustment of NN peaks improves the estimate of the NN interval resulting in HRV features (m = 9) that are better correlated with age (median percentage increase in correlation of individual features: 0.2%, IQR: 0.0 to 5.6%; correlation with age predictor and age from 0.721 to 0.787). Improved (sub-sample) quantization of the NN intervals (via interpolation) reduced the overall value of HRV features (median percentage reduction in feature value: -1.3%, IQR: -18.8 to 0.0; m = 9), primarily through a reduction in the energy of high-frequency oscillations. HRV features were also robust to missing data, with measures such as mean NN, fractal dimension and the smoothed nonlinear energy operator (SNEO) less susceptible to missing data than features such as VLF, LF, and HF. Furthermore, age predictions derived from a combination of HRV measures were more robust to missing data than individual HRV measures.