where d_y1 and d_y2 are the running 2nd and 4th moments of the input, respectively.
The 2.0* is not present in the estimate _expression_ from the reference (Pauluzzi and Beaulieu eq. 39).
The 3 dB offset can also be seen by generating a signal (say random BPSK) with AWGN, measuring the signal and noise powers separately while estimating the SNR via the m2m4 on the sum.
It's an easy fix, I just wanted to get verification and notify other users.