1

我知道在 R 中,它返回多重回归,它返回 βi=0 的假设检验,但是如果你想测试像 βi=1 这样的检验怎么办。是否有任何简单的命令,或者如果没有,您如何调用系数标准误差、系数值、回归自由度,以便我可以使用 t 分布 cdf 来计算 p 值。我想为通用程序执行此操作以运行多个数据

4

2 回答 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 回答