我正在比较两组的一长串观察,我有每个观察的t统计量和自由度,但我实际上没有原始数据。尽管如此,我想得到与t统计量相关的p值。我知道这是可能的,但我无法弄清楚如何让R报告一个p值。有什么建议吗?
答案 0 :(得分:7)
task view on Distributions应该是您回答此类问题的第一个停靠点
引用
对于大多数经典分布,基数R提供概率 分布函数(p),密度函数(d),分位数函数 (q)和随机数生成(r)。超越这个基础 功能,许多CRAN包提供额外的有用 分布。特别是,多变量分布以及 copulas以提供的包装形式提供。
你在t
分布之后,以及概率分布函数中的值,所以你所追求的函数是pt
pt(q, df)
其中q
是您的quantile
,在这种情况下是您观察到的测试统计信息,df
是自由度
我不确定the degrees of freedom for *each* observation
你的意思,我猜你的意思是你知道每组的观察次数。
所以,如果你假设群内方差相等,则自由度为n1+n2 - 2
,否则你需要考虑不等方差并获得近似自由度({{3} })