来自
http://docs.scipy.org/doc/scipy/reference/tutorial/stats.html
作为练习,我们也可以不使用提供的函数而直接计算ttest,该函数应该给我们相同的答案,所以它做到了:
tt = (sm-m)/np.sqrt(sv/float(n)) # t-statistic for meanpval = stats.t.sf(np.abs(tt), n-1)*2 # two-sided pvalue = Prob(abs(t)>tt)print 't-statistic = %6.3f pvalue = %6.4f' % (tt, pval)t-statistic = 0.391 pvalue = 0.6955