ｔ分布のPythonコードを紹介します。
そして、ｔ分布と正規分布を同時に描画して、差異を観察するコードを紹介します。
ｔ分布のコード【Python】
ｔ分布のPythonコード
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import norm #normは正規分布
from scipy.stats import t #t分布
n = 10 # サンプルの個数
deg = n-1 # 自由度
#ｔ分布を表示（自由度n-1）
X = np.arange(-3, 3, 0.1)
Y = t.pdf(X, deg)
plt.plot(X, Y, color="red", label = "t of n-1")
#グラフの表示
plt.legend() #凡例表示
plt.show()ｔ分布の表示結果
ｔ分布と正規分布の比較【Python】
分布を比較するコード
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import norm #normは正規分布
from scipy.stats import t #t分布
n = 2 # サンプルの個数
mu = 50 # 母平均
sigma = 10 # 母標準偏差
deg = n-1 # 自由度
#正規分布を表示
X = np.arange(-3, 3, 0.1)
Y = norm.pdf(X, 0, 1)
plt.plot(X, Y, color="orange", label = "N(0,1)")
#ｔ分布を表示（自由度n-1）
W = np.arange(-3, 3, 0.1)
Z = t.pdf(W, deg)
plt.plot(W, Z, color="red", label = "t of n-1")
#ｔ分布を表示（自由度n）
WW = np.arange(-3, 3, 0.1)
ZZ = t.pdf(WW, deg+1)
plt.plot(WW, ZZ, color="black", label = "t of n")
#グラフの表示
plt.legend() #凡例表示
plt.show()分布の表示結果
自由度 $2$, $1$ のｔ分布と、標準正規分布を同時に出力しました。
自由度を変えて、いくつも観察してみましょう。
