如何测试多个回归系数是否没有统计学差异?

机器算法验证 回归 假设检验 多重回归
2022-03-24 19:17:58

假设我估计以下多元线性回归 我如何测试

y=β0+β1x1+β2x2+β3x3+β4x4+ϵ
β1=β2=β3

我知道要测试是否你可以简单地构造一个测试 β1=β2Z

Z=β1β2seβ12+seβ22

是否有多个系数估计的类比?

1个回答

您可以使用检验来测试系数上的FL

让你的零假设是和你的设计矩阵秩为那么统计量将是:H0:Lβ=cXkF

F=(Lβ^c)(σ^2L(XX)1L)1(Lβ^c)q

其中是您正在测试的限制数量。在 null 下,这将具有自由度为分布。qFqnk

您可以使用包R的功能轻松做到这一点例如:linearHypothesiscar

library(car) 
lm.model <- lm(mtcars)
linearHypothesis(lm.model, c("cyl = 0", "disp = 0", "hp = 0")) # all 3 zero
linearHypothesis(lm.model, c("cyl = disp", "disp = hp")) # all 3 equal