我正在比较两组大量观察结果，我有每个观察值的t统计量和自由度，但实际上我并没有原始数据。尽管如此，我想获取与t统计量相关的p值。我知道这是可能的，但我无法弄清楚如何让R报告回一个p值。请问有什么建议吗？
1个回答
8
如果您有类似的问题，首先应该查看分布任务视图
引用如下：
对于大多数经典分布，基本的R提供概率分布函数(p)、密度函数(d)、分位数函数(q)和随机数生成(r)。除了这个基本功能外，许多CRAN软件包提供了其他有用的分布。特别是，贡献包中提供了多元分布以及copulas。
您需要的是t分布，以及来自概率分布函数的值，所以您需要的函数是pt
pt(q, df)
其中q是您的分位数，在这种情况下是您观察到的测试统计量，而df则表示自由度。
我不确定您所说的每个观测值的自由度是什么意思，我想您可能是知道每个组中观测值的数量。
因此，如果您假设组内方差相等，则自由度将为n1+n2-2，否则您需要考虑不等方差并获得近似自由度（Welch-Satterthwaite修改）。
- mnel
网页内容由stack overflow 提供, 点击上面的可以查看英文原文，
原文链接
原文链接