python - Python での t 分布と自由度から両側 P 値を見つける

Question

自由度 n の t 分布の P 値を決定する方法を教えてください。

このテーマに関する調査は、このスタック交換の回答を示しています: https://stackoverflow.com/a/17604216

np.abs(tt) が T 値であると仮定しますが、自由度でどのように機能するのですか?それは n-1 ですか?

前もって感謝します

Answer 1

はい、n-1その例の自由度です。

t 値と自由度が与えられると、(補完 CDF とも呼ばれる) の「生存関数」sfを使用してscipy.stats.t、片側 p 値を計算できます。最初の引数は T 値で、2 番目の引数は自由度です。

たとえば、このページの表の最初のエントリは、自由度 1 の場合、p=0.1 の臨界 T 値が 3.078 であることを示しています。でそれを確認する方法は次のt.sfとおりです。

In [7]: from scipy.stats import t

In [8]: t.sf(3.078, 1)
Out[8]: 0.09999038172554342   # Approximately 0.1, as expected.

両側 p 値の場合は、片側 p 値を 2 倍にします。