Jestem autorem pakietu ez dla R i pracuję nad aktualizacją w celu włączenia automatycznego obliczania ilorazów wiarygodności (LR) w danych wyjściowych ANOVA. Chodzi o to, aby zapewnić LR dla każdego efektu, który jest analogiczny do testu tego efektu, który osiąga ANOVA. Na przykład LR dla efektu głównego reprezentuje porównanie modelu zerowego z modelem, który zawiera efekt główny, LR dla interakcji reprezentuje porównanie modelu, który obejmuje oba składowe efekty główne w porównaniu z modelem, który obejmuje zarówno efekty główne, jak i ich interakcja itp.
Teraz rozumiem obliczenia LR na podstawie Glovera i Dixona ( PDF ), który obejmuje podstawowe obliczenia, a także poprawki na złożoność, oraz dodatku do Bortolussi & Dixon ( załącznik PDF ), który obejmuje obliczenia obejmujące zmienne wielokrotnego pomiaru. Aby przetestować moje zrozumienie, opracowałem ten arkusz kalkulacyjny , który bierze dfs i SS z przykładowej ANOVA (wygenerowanej z projektu 2 * 2 * 3 * 4 przy użyciu fałszywych danych) i krok po kroku obliczania LR dla każdego efektu.
Byłbym bardzo wdzięczny, gdyby ktoś z większą pewnością siebie przy takim obliczeniu mógł rzucić okiem i upewnić się, że wszystko zrobiłem poprawnie. Dla tych, którzy wolą kod abstrakcyjny, oto kod R implementujący aktualizację do ezANOVA () (patrz zwłaszcza linie 15-95).
źródło