R：計算給定卡方和自由度的 p 值

如何計算給定卡方和自由度的 p 值？例如，當 df = 2 時，卡方 = 15 的精確 p 值是多少？

在應用統計中，卡方檢驗統計作為殘差平方和，或來自平方效應之和或來自對數似然差。在所有這些應用程序中，目的是測試某個向量參數是否為零，以及是否為非零，並且卡方統計量與觀察到的效果的平方大小有關。所需的 p 值是卡方值的右尾概率，在 R 中為您的示例：

> pchisq(15, df=2, lower.tail=FALSE)
[1] 0.0005530844

對於其他 df 或統計值，您顯然只需將它們替換為上述代碼。

默認情況下，R 中的所有累積概率函數都會計算左尾概率。然而他們也有一個lower.tail論點，你總是可以設置它FALSE來獲得正確的尾部概率。這樣做而不是計算是一種很好的做法 $ 1-p $ 正如你可能在一些初級教科書中看到的那樣。

該函數qchisq進行反向計算，找到與任何給定尾部概率相對應的卡方統計量的值（“q”代表分位數）。例如，對應於 p 值為 0.05 的卡方統計量由下式給出

> qchisq(0.05, df=2, lower.tail=FALSE)
[1] 5.991465

引用自：https://stats.stackexchange.com/questions/250819