如何使用 R 计算临界 t 值?

机器算法验证 r 统计学意义 多重回归
2022-01-25 01:50:59

对不起,如果这是一个新手问题;我正在尝试第一次自学统计数据。我想我已经完成了基本程序,但我正在努力用 R 执行它。

所以,我试图在形式的多元线性回归中评估回归系数的重要性

y^=Xβ^

我认为用于测试的 t 统计量H0:β^j=0,Ha:β^j0是(谁)给的

t0=β^j0se(β^j)=β^jσ^2Cjj=β^jCjjSSRes/(np)
在哪里Cjj是个jth对角线进入(XX)1.

到现在为止还挺好。我知道如何使用 R 中的矩阵运算来计算所有这些值。但是为了拒绝空值,这本书说我需要

|t0|>tα/2,np

我如何计算这个临界值tα/2,np使用 R? 现在,我知道如何找到这些值的唯一方法是查看书后的表格。肯定有更好的办法。

1个回答

欢迎来到 Stats Stackexchange。您可以计算tα/2,np通过做 R 中的临界值qt(1-alpha/2, n-p)

在以下示例中,我要求 R 给我95%临界值df=1,2,,10. 结果是给定置信水平下 t 分布的前十个临界值的列表:
> qt(.975, 1:10)
[1] 12.706205 4.302653 3.182446 2.776445 2.570582 2.446912 2.364624 2.306004 2.262157 2.228139