python里怎么用F检验来比较两个样本的方差?
1个回答
F test的原理非常简单,所以不妨自己写
先求出两个样本的方差的比值,再写出两个样本的自由度,然后就去查F分布的概率累计函数,就可以得到p value了。
from scipy.stats import f
F = np.var(a) / np.var(b)
df1 = len(a) - 1
df2 = len(b) - 1
p_value = 1 - 2 * abs(0.5 - f.cdf(F, df1, df2))
谢啦,给你点个赞!
-
wlk1993
2017-12-01 15:07