対応の無いt検定
対応の無いt検定
import scipy as sp b = [7,8,10,5,8,7,9,5,6,9,10,6,7,8,7,9,10,10] a = [9,9,6,10,9,8,10,7,9,10,6,8,9,9,10,7,8,8,10,9] t_value = (sp.mean(a)-sp.mean(b))/sp.sqrt((sp.var(b,ddof=1)/len(b) + sp.var(a,ddof=1)/len(a))) t_value
1.463060101595841
一行で書くなら
t_value,p_value =sp.stats.ttest_ind(a,b,equal_var=False)
t_value
1.463060101595841
sp.stats.ttest_ind(a,b,equal_var=False)
だとp値も表示してくれる。
Ttest_indResult(statistic=1.463060101595841, pvalue=0.15335375977552246)
対応の無いt検定は、サンプル数が異なっても大丈夫