Python p-wartość z t-statystyka

Question

Python p-wartość z t-statystyka

Mam kilka wartości t i stopni swobody i chcę znaleźć wartości p z nich (jest dwuogniskowy). W prawdziwym świecie używałbym tabeli t-test z tyłu podręcznika Statystyki; jak zrobić odpowiednik w Pythonie?

Tzn.

t-lookup(5, 7) = 0.04156 albo coś w tym stylu(wymyślone liczby).

Wiem, że w SciPy gdybym miał tablice, mógłbym zrobić scipy.stats.ttest_ind, ale nie mam. mam tylko statystyki t i stopnie swobody.

21

python statistics

Author: Andrew Latham, 2013-07-10

Source

1 answers

score 25 · Accepted Answer

Od http://docs.scipy.org/doc/scipy/reference/tutorial/stats.html

Jako ćwiczenie możemy obliczyć nasze ttest również bezpośrednio bez użycia podanej funkcji, która powinna dać nam taką samą odpowiedź, a więc czy:

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