Bayesian denoising algorithm dealing with colored, non-stationary noise in continuous glucose monitoring timeseries

Introduction: The retrospective analysis of continuous glucose monitoring (CGM) timeseries can be hampered by colored and non-stationary measurement noise. Here, we introduce a Bayesian denoising (BD) algorithm to address both autocorrelation of measurement noise and temporal variability of its variance. Methods: BD utilizes adaptive, a-priori models of signal and noise, whose unknown variances are derived on partially-overlapped CGM windows, via smoothing approach based on linear mean square estimation. The CGM signal and noise variability profiles are then reconstructed using a kernel smoother. BD is first assessed on two simulated datasets, DS1 and DS2. On DS1, the effectiveness of accounting for colored noise is evaluated by comparison against a literature algorithm; on DS2, the effectiveness of accounting for the noise variance temporal variability is evaluated by comparison against a Butterworth filter. BD is then evaluated on 15 CGM timeseries measured by the Dexcom G6 (DR). Results: On DS1, BD allows reducing the root-mean-square-error (RMSE) from 8.10 [6.79–9.24] mg/dL to 6.28 [5.47–7.27] mg/dL (median [IQR]); on DS2, RMSE decreases from 6.85 [5.50–8.72] mg/dL to 5.35 [4.48–6.49] mg/dL. On DR, BD performs a reasonable tracking of noise variance variability and a satisfactory denoising. Discussion: The new algorithm effectively addresses the nature of CGM measurement error, outperforming existing denoising algorithms.