scipy.stats

scipy.statsモジュールには統計には乱数に限らず統計学に関する関数がたくさん。

from scipy import *
# データはあんまり意味なし
x = arange(0.0, 2*pi, 0.01)
y = sin(x)

# 幾何平均
print stats.mean(y)      # -1.72551106802e-06
# stats.varが分散 (sigma^2)
print sqrt(stats.var(y)) # 0.707286084386

どっちともaxisというoptional argumentを取れる。
多次元の場合はaxisの方向にそって平均や分散を求めることができる。
大したこと無い関数だけど多次元の場合にも使えるし、
イチイチ考える必要がないからちょっと便利。