用python進行F測試，找到臨界值

Question

使用python，是否可以計算具有x和y自由度的F分布的臨界值？ 換句話說，我需要計算給定 ax 自由度和 5% 置信度的臨界值，但我沒有看到統計書籍中的表格，是否可以使用 python 中的任何函數獲取它？

例如，我想找到具有 3 和 39 自由度的 F 分布的臨界值，置信水平為 5%。 答案應該是：2.85

Answer 1

IIUC，您可以使用scipy.stats.f.ppf ，它給出了 cdf 的倒數：

>>> import scipy.stats
>>> scipy.stats.f.ppf(q=1-0.05, dfn=3, dfd=39)
2.8450678052793514
>>> crit = _
>>> scipy.stats.f.cdf(crit, dfn=3, dfd=39)
0.95000000000000007

Answer 2

只是為了讓 dfn 和 dfd 更清楚地了解 DSM 的答案/scipy.stats：

dfn表示分子中使用的方差估計的自由度數。

dfd是分母中使用的方差估計的自由度數。

dfn = a-1
dfd = N-a

其中a是組數， N是實驗中的受試者總數。
F 分布的形狀取決於dfn和dfd 。 自由度越低，需要顯着的 F 值就越大。