Zrozumienie efektu ciągłego czynnika losowego w modelu efektów mieszanych

10

Rozumiem wpływ kategorycznego efektu losowego na model efektów mieszanych, ponieważ wykonuje on częściowe łączenie obserwacji według poziomu efektu losowego, skutecznie zakładając, że obserwacje same w sobie nie są niezależne, a jedynie ich częściowe pule. Również dla mojego zrozumienia, w takim modelu obserwacje dzielące ten sam poziom efektu losowego, ale różniące się poziomem ustalonego efektu, przeważą nad obserwacjami różniącymi się zarówno poziomem efektu losowego, jak i poziomu efektu stałego.

Jaki jest zatem wpływ ciągłego czynnika losowego? Biorąc pod uwagę, że model bez efektu losowego wykazał, że efekt stały miał wielkość efektu X. Czy powinienem oczekiwać, że jeśli obserwacje na różnych poziomach efektu stałego pochodzą z odległych krańców kontinuum efektu losowego, rozmiar efektu będzie mniejszy w model, który obejmował czynnik losowy, a jeśli obserwacje na różnych stałych poziomach czynników miały podobne wartości efektu losowego, to rozmiar efektu by się zwiększył?

Roey Angel
źródło
1
Czy możesz podać formuły i / lub kod R / Stata jako przykład swojego myślenia? Używasz nieco niezwykłego języka ... przynajmniej dla mnie niezwykłego. Myślę, że twój „ciągły czynnik losowy” jest tym, co nazwałbym „losowym nachyleniem”, ale najpierw chciałem to sprawdzić.
StasK
@StasK W kategoriach R: jeśli czynnik losowy jest kategoryczny (czynnik w R), wówczas obserwacje są częściowo połączone, tj. Średnie grupy (poziomy współczynników losowych) są ważonymi średnimi średniej populacji, a niepodzielona grupa średnich z wagami proporcjonalnymi do wielkości próby i odwrotności wariancji. Moje pytanie brzmi: co się dzieje, gdy czynnik losowy jest ciągły (liczbowo w kategoriach R). Jak to wpływa na model?
Roey Angel,
1
@RoeyAngel: prawdopodobnie nie wpływa to w żaden rozsądny sposób. Specjalnie dla R„y lmer, na przykład model, w którym efekt losowy ma znaczącej wartości dla każdego punktu danych, nie zdołają nawet obliczeniowej. Pomyśl o tym w kategoriach czysto pojęciowych: jeśli macierz jest kwadratowa, to wektor γ trzymający realizację efektów losowych będzie miał rozmiar N ( N : # punktów próbnych), a zatem będziesz miał niemożliwą do zidentyfikowania strukturę błędów. Czy na pewno o to pytasz? Jako StasK trudno mi również odpowiedzieć na twoje pytanie. ZγNN
usεr11852
@ user11852 hmmm Szczerze mówiąc, nigdy nie próbowałem tego z przypadkowym efektem, w którym każdy punkt ma unikalną wartość. Więc w zasadzie mówisz, że efekt losowy jest zawsze traktowany jako czynnik kategoryczny (tj. Nie ma analogii do tego, jak na przykład ciągłe zmienne są traktowane w ANCOVA).
Roey Angel,
@RoeyAngle: Nie wiem konkretnie o ANCOVA, ale z pewnością to, co powiedziałem o stojakach niemożliwych do zidentyfikowania. Nie możesz oszacować jeśli γ jest równe rozmiarowi twoich danych. Zostało to potraktowane tak kategorycznie, jak Z odzwierciedla strukturę (tj. Kategoryzację) samych danych (np. Partia, grupa, lokalizacja itp.). Pomyśl o tym w kontekście modeli hierarchicznych (podzbioru modeli mieszanych): jeśli hierarchia zdefiniuje na pewnym poziomie tyle potomków, ile punktów danych, wówczas będzie zbędna. γγZ
usεr11852

Odpowiedzi:

1

Musiałem mocno przemyśleć to, o co prosiłeś. Na początku pomyślałem @ user11852, że chciałeś, aby każda obserwacja miała swój własny wyjątkowy efekt losowy. To sprawiłoby, że model byłby beznadziejnie niezidentyfikowany, ponieważ nie byłoby możliwego sposobu na odróżnienie losowej zmienności efektu od błędu modelu.

Uważam jednak, że w zakresie zamierzonego pytania wszystkie efekty losowe są w rzeczywistości ciągłe i prawdopodobnie normalnie rozłożone. Jednak twoja aluzja do „kategorialnego” nie jest zła, ponieważ macierz projektowa dla losowego przechwytywania (zwykle nazywana Z) wyglądałaby jak macierz projektowa dla zmiennej jakościowej.

(α¯+αi)+(β¯+βi)xij,
α¯β¯αiβiiβiαii

Pomyślmy teraz o twojej proponowanej sytuacji:

różne poziomy ustalonego efektu pochodziły z odległych krańców kontinuum efektów losowych

β¯xijxijβiixijxijβi

β

Ben Ogorek
źródło