Porównywanie modeli regresji wielorakiej mieszanej

Kiedy stosować:

Porównywanie modeli służy do wyboru takiego zestawu zmiennych niezależnych (efektów stałych), który opisuje dane najlepiej przy możliwie małej liczbie parametrów. Działa tak samo jak porównywanie modeli regresji wielorakiej: porównujemy model pełny z modelem zredukowanym (bez części zmiennych) albo pozwalamy programowi wybrać model automatycznie (dobór wsteczny, w przód, krokowy lub przegląd podzbiorów).

W modelach mieszanych obowiązuje jedna dodatkowa zasada: struktura losowa jest wspólna dla wszystkich porównywanych modeli. Dobór dotyczy wyłącznie efektów stałych, a losowy wyraz wolny, losowe nachylenia i typ kowariancji ustalone w oknie pozostają w każdym z modeli. Zmienna, której nadano losowe nachylenie, pozostaje w każdym porównywanym modelu jako efekt stały (nie można usunąć efektu stałego zmiennej, zachowując jej losowe nachylenie). O tym, czy same efekty losowe są potrzebne, decydują testy LR w raporcie regresji wielorakiej mieszanej.

Definicje miar

Modele porównuje się na podstawie kryteriów informacyjnych AIC, AICc i BIC (im mniejsza wartość, tym lepszy model) oraz testu ilorazu wiarygodności (LR) dla modeli zagnieżdżonych: statystyka \(2(\ell_P-\ell_Z)\) (\(\ell_P\), \(\ell_Z\) – logarytm wiarygodności modelu pełnego i zredukowanego) ma rozkład \(\chi^2\) o liczbie stopni swobody równej liczbie usuniętych zmiennych. Nieistotny wynik oznacza, że usunięte zmienne nie wnoszą informacji i należy wybrać model prostszy; wynik istotny – że model pełny jest istotnie lepszy.

Jak ustawić analizę

Statystyka \(\to\) Modele mieszane \(\to\) Regresja wieloraka mieszana – porównywanie modeli.

W oknie wskazujemy zmienną zależną \(Y\), zmienne modelu pełnego i zmienne modelu zredukowanego (lub kryterium doboru automatycznego), tak jak w porównywaniu modeli regresji wielorakiej. Opcje modelu mieszanego:

• Zmienna grupująca (ID) – kolumna identyfikująca grupę; wspólna dla wszystkich modeli,
• Losowe nachylenia (z X modelu pełnego) – zmienne, których efekt różni się między grupami; taka zmienna pozostaje w każdym porównywanym modelu,
• Kowariancja efektów losowych – niezależne lub pełna z korelacją, wspólna dla wszystkich modeli,
• Współczynniki modeli (dobór: ML) – REML lub ML dla tabel współczynników; sam dobór zawsze ML,
• Stopnie swobody efektów stałych – Satterthwaite lub test z Walda w tabelach współczynników.

Wyniki i interpretacja

Raport zawiera ramkę z opisem wspólnej struktury losowej (zmienna grupująca, liczba grup, składniki losowe, kowariancja), tabelę Kryteria informacyjne i testy LR z wartościami \(-2LL\) (ML), AIC, AICc, BIC, \(R^2\) brzegowym i warunkowym każdego modelu oraz testem LR między kolejnymi modelami, a także tabele współczynników porównywanych modeli (efekty stałe warunkowe, z wariancjami efektów losowych każdego modelu). Jeśli dla któregoś modelu wystąpiło dopasowanie osobliwe, jest to zaznaczone w tabeli.