I would add all stats, including m, s, something like
``assert_allclose((m, v, s, k), (0, m2, 0, m4/m2**2 -3.), atol=1e-14, rtol=1e-??)``