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(X′X)−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 .qFqn−k

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.