which, while preserving the streaming nature of the algorithm (for
     parallel computation purposes) is numerically much more stable. Note
     that it does *not* require two passes as a common fix to numerical
     instability does.
ENH: removed the static function that calculates centered moments from
     raw moments since the latter are no longer calculated.
ENH: added G1 skewness estimator.
ENH: updated tests accordingly.