r - 从 t 统计量中获取 p 值

Question

很难说出这里问的是什么。这个问题是模棱两可的、模糊的、不完整的、过于宽泛的或修辞的，无法以目前的形式得到合理的回答。为了帮助澄清这个问题以便可以重新打开它， visit the help center .

关闭 9年前。

我正在比较两组的一长串观察结果，我有每个观察结果的 t 统计量和自由度，但实际上我没有原始数据。尽管如此，我还是想获得与 t 统计量相关的 p 值。我知道这是可能的，但我不知道如何让 R 报告 p 值。有什么建议吗？

Answer 1

task view on Distributions应该是您回答此类问题的第一站

引用

For most of the classical distributions, base R provides probability distribution functions (p), density functions (d), quantile functions (q), and random number generation (r). Beyond this basic functionality, many CRAN packages provide additional useful distributions. In particular, multivariate distributions as well as copulas are available in contributed packages.

你在 t 分布之后，并且是概率分布函数中的一个值，所以你所追求的函数是 pt

pt(q, df)

其中 q 是您的 quantile，在这种情况下是您观察到的测试统计数据，而 df 是自由度

我不确定您所说的*每个*观察的自由度是什么意思，我想您的意思是您知道每个组中的观察数量。

因此，如果您假设组内方差相等，则自由度将为 n1+n2 - 2，否则您需要考虑不等方差并获得近似度自由的(Welch-Satterthwaite modification)