pythonで信頼区間を計算するコードちょうだい

2024年6月8日

Pythonで信頼区間を計算するためには、scipyパッケージの t または norm モジュールを使用できます。以下は、標本標準偏差が与えられた場合に95％信頼区間を計算する例です。

python
import numpy as np
from scipy.stats import t

# データセット
data = [10, 12, 14, 16, 18]

# 標本標準偏差と標本平均の計算
s = np.std(data, ddof=1)  # 不偏標準偏差を計算
mean = np.mean(data)

# データのサンプル数
n = len(data)

# 95%信頼区間の計算
alpha = 0.95
t_value = t.ppf((1 + alpha) / 2, df=n-1)  # 自由度 n-1 の t 分布の上側 2.5 パーセント点
margin_of_error = t_value * s / np.sqrt(n)
confidence_interval = (mean - margin_of_error, mean + margin_of_error)

print("95%信頼区間:", confidence_interval)

このコードでは、与えられたデータセットから標本標準偏差 s と標本平均 mean を計算し、95％信頼区間を計算しています。t.ppf() 関数は、t分布の累積分布関数の逆関数であり、95％信頼水準に対応するt値を計算するために使用されます。