我知道在 R 中,它返回多重回归,它返回 βi=0 的假设检验,但是如果你想测试像 βi=1 这样的检验怎么办。是否有任何简单的命令,或者如果没有,您如何调用系数标准误差、系数值、回归自由度,以便我可以使用 t 分布 cdf 来计算 p 值。我想为通用程序执行此操作以运行多个数据
Kamster
问问题
2064 次
2 回答
3
R 中有几个包可以让您测试系数是否不同于 0 以外的值。例如,https://www.rforge.net/doc/packages/FSA/hoCoef.html。在这种情况下,您将在 hoCoef 函数中使用:指定 bo = 1。
于 2015-03-04T23:44:26.080 回答
2
虽然有包可以做到这一点,但它很简单,你可以编写一个小函数。
这将返回回归摘要的 p 值,sl
用于针对与 中的值相等的双尾检验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
于 2015-03-04T23:53:14.363 回答