Python valeur p à partir de la statistique t
J'ai des valeurs t et des degrés de liberté et je veux en trouver les valeurs p (c'est bilatéral). Dans le monde réel, j'utiliserais un tableau de test t à l'arrière d'un manuel de statistiques; comment faire l'équivalent en Python?
par exemple.
t-lookup(5, 7) = 0.00245 ou quelque chose comme ça.
Je sais dans SciPy si j'avais des tableaux que je pouvais faire scipy.stats.ttest_ind, mais pas moi. J'ai juste des statistiques t et des degrés de liberté.
De http://docs.scipy.org/doc/scipy/reference/tutorial/stats.html
En tant qu'exercice, nous pouvons également calculer directement notre test sans utiliser la fonction fournie, ce qui devrait nous donner la même réponse, et c'est donc le cas:
tt = (sm-m)/np.sqrt(sv/float(n)) # t-statistic for mean
pval = 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