r - 从 at stat 获取 ap 值

Question

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

Answer 1

Distributions 上的任务视图应该是您回答此类问题的第一个停靠点

去引用

对于大多数经典分布，基数 R 提供概率分布函数 (p)、密度函数 (d)、分位数函数 (q) 和随机数生成 (r)。除了这个基本功能之外，许多 CRAN 包还提供了其他有用的发行版。特别是，多变量分布以及 copula 在贡献包中可用。

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

pt(q, df)

q你在哪里quantile，在这种情况下是你观察到的检验统计量，以及df自由度

我不确定您的意思是什么the degrees of freedom for *each* observation，我想您的意思是您知道每组中的观察次数。

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