Comment tester si les coefficients de régression multiples ne sont pas statistiquement différents?


8

Supposons que régression linéaire multivariée suivante Comment puis-je tester cette ?

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

Je sais que pour tester si vous pouvez simplement construire un test avec β1=β2Z

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

Existe-t-il un analogue pour les estimations de coefficients multiples?


2
Le test d'égalité de et suppose implicitement que les estimations de sont pas corrélées. En général, ce sera incorrect; le dénominateur doit inclure un terme pour leur covariance. β1β2βi
whuber

Si vos variables X sont dans des unités différentes, les coefficients bêta sont également dans des unités différentes. Dans ce cas, je ne vois pas en quoi il serait logique de les comparer.
Harvey Motulsky

Réponses:


9

Vous pouvez utiliser le test pour tester toute restriction linéaire sur vos coefficients.FL

Soit votre hypothèse nulle et votre matrice de conception de rang . La statistique sera alors:H0:Lβ=cXkF

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

où est le nombre de restrictions que vous testez. Sous le nul, cela aura une distribution avec des degrés de liberté et .qFqnk

En Rvous pouvez facilement le faire avec la fonction linearHypothesisdu carpackage. Par exemple:

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
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.