ScipyにはExcelのftestやRのvar.testのように単純に等分散かどうかのp値を返す関数がない。
python用ftestを作ってみた。因みに両側検定。
import numpy as np
import scipy.stats as stats
def ftest(x1, x2):
dfn = len(x1)-1
dfd = len(x2)-1
F = np.var(x1) / np.var(x2)
f_pf = stats.f.cdf(F, dfn, dfd)
return min(f_pf, 1 - f_pf)*2
import scipy.stats as stats
def ftest(x1, x2):
dfn = len(x1)-1
dfd = len(x2)-1
F = np.var(x1) / np.var(x2)
f_pf = stats.f.cdf(F, dfn, dfd)
return min(f_pf, 1 - f_pf)*2
0 件のコメント:
コメントを投稿