我知道在R中它返回一个多元回归它返回βi= 0的假设检验但是如果你想测试像βi= 1这样的测试怎么办。有没有任何简单的命令,如果不是如何调用系数标准误,系数值,回归自由度,所以我可以使用t分布cdf来计算p值。我想为一般程序执行此操作以运行多个数据
答案 0 :(得分:3)
R中有几个包可以测试系数是否与0以外的值不同。例如,https://www.rforge.net/doc/packages/FSA/hoCoef.html。在这种情况下,您使用: 在hoCoef函数中指定,bo = 1。
答案 1 :(得分:2)
虽然有这样做的软件包,但是你可以编写一些函数这么简单。
这将返回sl
中回归摘要的p值,以进行双尾检验,以确保与b0
中的值相等:
testb0=function(sl,b0) {
slm=sl$coefficients #$
t0=(slm[,1]-b0)/slm[,2]
pt(abs(t0),sl$df[2],lower.tail=FALSE)
}
对该功能的测试:
testb0( summary(lm(dist~speed+I(speed^2),cars)), b0=c(0,1,0) )
返回三个p值
(Intercept) speed I(speed^2)
0.43415754 0.48308979 0.06820122