Calculates signal-to-noise. Axis can equal None (ravel array first), an integer (the axis over which to operate).
scipy.stats.samplestd
scipy.stats.bayes_mvs