Zrozumienie wariancji efektów losowych w modelach lmer ()

16

Mam problem ze zrozumieniem wyników mojego lmer()modelu. Jest to prosty model zmiennej wynikowej (Wsparcie) ze zmiennymi przechwytywaniami stanu / Losowymi efektami stanu:

mlm1 <- lmer(Support ~ (1 | State))

Wyniki summary(mlm1)są następujące:

Linear mixed model fit by REML 
Formula: Support ~ (1 | State) 
   AIC   BIC logLik deviance REMLdev
 12088 12107  -6041    12076   12082
Random effects:
 Groups   Name        Variance  Std.Dev.
 State    (Intercept) 0.0063695 0.079809
 Residual             1.1114756 1.054265
Number of obs: 4097, groups: State, 48

Fixed effects:
            Estimate Std. Error t value
(Intercept)  0.13218    0.02159   6.123

Rozumiem, że wariancja przechwytywania stanów zmiennych / efektów losowych jest 0.0063695. Ale kiedy wyciągam wektor tych losowych efektów i obliczam wariancję

var(ranef(mlm1)$State)

Wynik jest 0.001800869:, znacznie mniejszy niż wariancja zgłoszona przez summary().

O ile rozumiem, model, który podałem, można zapisać:

yi=α0+αs+ϵi, for i={1,2,...,4097}

αsN(0,σα2), for s={1,2,...,48}

Jeśli jest to poprawne, wówczas wariancja efektów losowych ( ) powinna wynosić . Jednak w rzeczywistości nie pasują one do mnie .σ 2 ααsσα2lmer()

nomad545
źródło
Czy masz wiedzę na temat sposobu szacowania parametrów lmer()? Wydaje się, że postulat ten jest szacowany przez wariancji empirycznej szacunkowej efektów losowych a s . Opis modelu nie jest jasne (perharps y i powinna być y i y ). Czy to zrównoważony projekt? σα2α^syiyis
Stéphane Laurent,
Oto bardzo podobne pytanie, z nieco inną odpowiedzią
Arne Jonas Warnke

Odpowiedzi:

11

To klasyczna anova w jedną stronę. Bardzo krótka odpowiedź na twoje pytanie brzmi: komponent wariancji składa się z dwóch terminów.

σ^α2=E[148s=148αs2]=148s=148α^s2+148s=148var(α^s)

Tak więc obliczony przez ciebie termin jest pierwszym terminem na rh (ponieważ losowe efekty oznaczają zero). Drugi termin zależy od tego, czy użyto REML ML, oraz od sumy kwadratowych błędów standardowych twoich losowych efektów.

prawdopodobieństwo prawdopodobieństwa
źródło
2
Ok, rozumiem! Tak więc suma kwadratów SE z RE - 1/48 * sum((se.ranef(mlm1)$State)^2)- jest 0.004557198. Wariancja punktowych oszacowań RE (uzyskanych, jak wyżej, przy użyciu var(ranef(mlm1)$State)) wynosi 0.001800869. Suma to 0.006358067, co jest zgłaszane przy użyciu wariancji summary()w lmer()modelu powyżej, 4 lub 5 cyfr przynajmniej. Wielkie dzięki @probability
nomad545
2
Dla tych, którzy szukają tej odpowiedzi i komentarza do pomocy, zauważ, że nomad545 również skorzystał z armpakietu R dla tej se.ranef()funkcji.
ndoogan
1
@probabilityislogic: Czy możesz podać więcej szczegółów na temat sposobu obliczania tego równania? W szczególności, w jaki sposób osiągnięto drugą równość? Czy nie powinno być też czapki na alfie po pierwszej równości?
użytkownik1357015
1
YN.ormzal(1nα0,Σ)Σ=Inσe2+σα2ZZTE(αs)=0var(αs)=E(αs2)