Muszę porównać dwa stoki regresji, gdzie:
$
y_1 ~ a + b_1x
y_2 ~ a + b_2x
$
Jak mogę porównać b1 i b2?
Lub w języku mojego konkretnego przykładu u gryzoni, chcę porównać
antero-posterior diameter ~ a + b1 * humeral length
de naso-occipital length ~ a + b2 * humeral length
regression
regression-coefficients
Dra. Alejandra Echeverria
źródło
źródło
Odpowiedzi:
Ok, spójrzmy na twoją sytuację. Masz w zasadzie dwie regresje (APD = średnica przednio-tylna, NOL = długość nosowo-potyliczna, HL = długość ramienia):
Aby przetestować hipotezę , możesz wykonać następujące czynności:β1 , 1= β1 , 2
Spójrzmy na przykład z gotowymi danymi (w
R
):dummy.var
x.new:dummy.var
Ostrzeżenie: Działa to tylko wtedy, gdy średnica przednio-tylna i długość nosowo-potyliczna (dwie zmienne zależne) są niezależne. W przeciwnym razie może się to bardzo skomplikować.
EDYTOWAĆ
Te dwa posty na stronie dotyczą tego samego pytania: pierwszego i drugiego .
źródło
dummy.var
, to znaczy dla obu DV. W zależności od tego, jakie DV są w oryginalnym kontekście, możliwe jest, że rezydualne wariancje są radykalnie różne w oddzielnych regresjach każdego DV. Zastanawiam się, czy lepiej byłoby zastosować to samo zaproponowane przez ciebie podstawowe podejście, ale zgls
modelem, w którym szacujemy różne odchylenia rezydualne dla każdego DV. Masz jakieś przemyślenia na ten temat?library(nlme); mod4 <- gls(y.new~x.new*dummy.var, weights=varIdent(form= ~1 | dummy.var))