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 值，第二个参数是自由度。

例如，本页表格的第一个条目表示，对于 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 值加倍。