Je dois comparer deux pentes de régression où:
$
y_1 ~ a + b_1x
y_2 ~ a + b_2x
$
Comment comparer b1 et b2?
Ou dans la langue de mon exemple spécifique chez les rongeurs, je veux comparer
antero-posterior diameter ~ a + b1 * humeral length
de naso-occipital length ~ a + b2 * humeral length
regression
regression-coefficients
Dra. Alejandra Echeverria
la source
la source
Réponses:
D'accord, regardons votre situation. Vous avez essentiellement deux régressions (APD = diamètre antéro-postérieur, NOL = longueur naso-occipitale, HL = longueur humérale):
Pour tester l'hypothèse , vous pouvez procéder comme suit:β1 , 1= β1 , 2
Regardons un exemple avec des données composées (en
R
):dummy.var
x.new:dummy.var
Attention: cela ne fonctionne que si le diamètre antéro-postérieur et la longueur naso-occipitale (les deux variables dépendantes) sont indépendants. Sinon, cela peut devenir très compliqué.
ÉDITER
Ces deux articles sur le site traitent de la même question: Premier et deuxième .
la source
dummy.var
, c'est-à-dire pour les deux DV. Selon ce que sont les DV dans le contexte d'origine, il est possible que les variances résiduelles soient radicalement différentes dans les régressions distinctes de chaque DV. Je me demande s'il serait préférable d'utiliser la même approche de base que vous avez proposée, mais avec ungls
modèle où nous estimons différentes variances résiduelles pour chaque DV. Des pensées à ce sujet?library(nlme); mod4 <- gls(y.new~x.new*dummy.var, weights=varIdent(form= ~1 | dummy.var))