Jak interpretować wariancję i korelację efektów losowych w modelu efektów mieszanych?


28

Mam nadzieję, że wszystkim wam to nie przeszkadza, ale potrzebuję pomocy w interpretacji wyników dla liniowego modelu efektów mieszanych, o których starałem się nauczyć w R. Jestem nowy w analizie danych podłużnych i regresji liniowych efektów mieszanych. Mam model, który dopasowałem do tygodni jako predyktor czasu, a moim wynikiem jest kurs pracy. Modelowałem wynik z tygodniami (czasem) i kilkoma ustalonymi efektami, płcią i rasą. Mój model zawiera efekty losowe. Potrzebuję pomocy w zrozumieniu, co oznacza wariancja i korelacja. Dane wyjściowe są następujące:

Random effects  
Group   Name    Variance  
EmpId intercept 680.236  
weeks           13.562  
Residual 774.256  

Korelacja wynosi .231.

Potrafię zinterpretować tę korelację, ponieważ istnieje pozytywny związek między tygodniami a wynikiem, ale chcę móc to powiedzieć w kategoriach „23% ...”.

Naprawdę doceniam pomoc.


Dzięki „gość” i Makro za odpowiedź. Przepraszam, że nie odpowiedziałem, byłem na konferencji i teraz nadrabiam zaległości. Oto wynik i kontekst.

Oto podsumowanie modelu LMER, który prowadziłem.

>summary(LMER.EduA)  
Linear mixed model fit by maximum likelihood  
Formula: Score ~ Weeks + (1 + Weeks | EmpID)   
   Data: emp.LMER4 

  AIC     BIC   logLik   deviance   REMLdev   
 1815     1834  -732.6     1693    1685

Random effects:    
 Groups   Name       Variance Std.Dev. Corr  
 EmpID   (Intercept)  680.236  26.08133        
          Weeks         13.562 3.682662  0.231   
 Residual             774.256  27.82546        
Number of obs: 174, groups: EmpID, 18


Fixed effects:    
            Estimate Std. Error  t value  
(Intercept)  261.171      6.23     37.25    
Weeks          11.151      1.780    6.93

Correlation of Fixed Effects:  
     (Intr)  
Days -0.101

Nie rozumiem, jak interpretować wariancję i resztkę dla efektów losowych i wyjaśnić to komuś innemu. Nie wiem też, jak interpretować korelację, poza tym, że jest dodatnia, co wskazuje, że osoby z wyższymi punktami przechwytywania mają wyższe nachylenia, a te z tymi z niższymi punktami przechwytywania mają niższe nachylenia, ale nie wiem, jak wyjaśnić korelację w kategoriach z 23% z. . . . (Nie wiem, jak dokończyć zdanie, a nawet czy ma to sens). Jest to dla nas analiza innego typu, ponieważ my (ja) próbujemy przejść do analiz podłużnych.

Mam nadzieję, że to pomoże.

Dziękuję za pomoc do tej pory.

Zeda


1
Zeda, przydałoby się zobaczyć więcej wyników R, w tym podsumowanie efektów stałych
gość

1
Jedno, co widzę, to to, że szacunkowa korelacja wewnątrzklasowa dla EmpID wynosi . Oznacza to, że szacunkowa korelacja między dwiema osobami na tym samym poziomie EmpID to . Zgadzam się z @guest, że więcej danych wyjściowych (i pewien kontekst) byłoby pomocne. ρρ^=680,236/(680,236+13,562+774,256)ρ^
Makro

Zeda, przekonwertowałem twoją odpowiedź jako edycję i połączyłem twoje dwa niezarejestrowane konta. Zarejestruj go, aby móc samodzielnie śledzić i aktualizować swój post.
chl

Odpowiedzi:


40

Twój dopasowany model lme()można wyrazić jako

yjajot=α0+α1xjot+δ0ja+δ1jaxjot+ϵjajot

gdzie jest wynikiem i- tego pracownika w x j tygodniach, α 0 i α 1 oznaczają odpowiednio stały punkt przecięcia i nachylenie, δ 0 i i δ 1 i są losowym punktem przecięcia i nachylenia, a ϵ i j jest wartością resztkową . Założenia dla efektów losowych δ 0 i , δ 1 i resztkowych ϵ i jyjajotjaxjotα0α1δ0jaδ1jaϵjajotδ0jaδ1jaϵjajot

(δ0ja,δ1ja)T.reN.((0,0)T.,sol)ϵjajotreN.(0,σ2))

sol

(sol12)sol122)sol122)sol2)2))

Możesz uzyskać macierz wariancji między warunkami efektów losowych od VarCorr(LMER.EduA)$ID.

Twój wynik w zasadzie to mówi

α0α1

sol12)sol2)2)σ2)

sol122)VarCorr(LMER.EduA)0,23×sol12)sol2)2)

sol12)sol2)2)


2
L.ZAT.miX

@chl: Naprawdę doceniam twoją odpowiedź w tak ładnym formacie (nic nie wiem o LaTex). Co ważniejsze, poprawiłeś moją niechlujną odpowiedź dotyczącą części kowariancji. Jeszcze raz dziękuję, chl!
bluepole

Kredyty należy kierować do @GGeco, który podał szczegółowe informacje na temat matrycy VC; tak jak powiedziałem, tekstowałem tylko część twojej odpowiedzi (i +1).
chl

2
Jak by to działało, gdybyś miał wiele efektów losowych?
user124123
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.