Ich muss zwei Regressionssteigungen vergleichen, bei denen:
$
y_1 ~ a + b_1x
y_2 ~ a + b_2x
$
Wie kann ich b1 und b2 vergleichen?
Oder in der Sprache meines speziellen Beispiels bei Nagetieren möchte ich vergleichen
antero-posterior diameter ~ a + b1 * humeral length
de naso-occipital length ~ a + b2 * humeral length
regression
regression-coefficients
Dra. Alejandra Echeverria
quelle
quelle
Antworten:
Okay, schauen wir uns Ihre Situation an. Sie haben grundsätzlich zwei Regressionen (APD = antero-posteriorer Durchmesser, NOL = naso-occipitale Länge, HL = Humeruslänge):
Schauen wir uns ein Beispiel mit erfundenen Daten (in
R
) an:dummy.var
x.new:dummy.var
Warnung: Dies funktioniert nur, wenn der antero-posteriore Durchmesser und die naso-okzipitale Länge (die beiden abhängigen Variablen) unabhängig sind. Sonst kann es sehr kompliziert werden.
BEARBEITEN
Diese beiden Beiträge auf der Website befassen sich mit derselben Frage: Erstens und zweitens .
quelle
dummy.var
, dass die Restvarianz auf beiden Ebenen , dh für beide DVs, gleich ist. Abhängig davon, was die DVs im ursprünglichen Kontext sind, ist es möglich, dass die Restvarianzen in den separaten Regressionen jedes DV radikal unterschiedlich sind. Ich frage mich, ob es besser wäre, den gleichen grundlegenden Ansatz zu verwenden, den Sie vorgeschlagen haben, aber mit einemgls
Modell, bei dem wir für jeden DV unterschiedliche Restvarianzen schätzen. Irgendwelche Gedanken dazu?library(nlme); mod4 <- gls(y.new~x.new*dummy.var, weights=varIdent(form= ~1 | dummy.var))