Here, we used the Welch’s method to estimate the power spectrum. The number of points used for FFT was chosen to be a power of 2 and larger than 4 times the sampling frequency. This corresponds to frequency resolution below 0.5 Hz when Hamming window was used (Full-width half maximum for this window is 2 bins).