[解決済み] サンプルデータから信頼区間を計算する

質問

正規分布と仮定して信頼区間を計算したいサンプルデータがあります。

numpyとscipyのパッケージを見つけてインストールし、numpyが平均と標準偏差を返すようにしました（numpy.mean(data)、dataはリスト）。サンプルの信頼区間を取得するための任意のアドバイスは非常に感謝されます。

どのように解決するのですか？

import numpy as np
import scipy.stats


def mean_confidence_interval(data, confidence=0.95):
    a = 1.0 * np.array(data)
    n = len(a)
    m, se = np.mean(a), scipy.stats.sem(a)
    h = se * scipy.stats.t.ppf((1 + confidence) / 2., n-1)
    return m, m-h, m+h

というように計算することができます。