powered by
Rolling estimate of expression fold change phi
phi_hat_roll(Y_obs, lambda_ref, d_obs, mu, sig, h)