W książce statystyk, którą czytam, zaleca się omega kwadrat do zmierzenia efektów moich eksperymentów. Udowodniłem już, że stosując projekt podzielonego wykresu (połączenie projektu między podmiotami i projektu między podmiotami), moje czynniki w obrębie podmiotów są statystycznie istotne przy p <0,001 i F = 17.
Teraz chcę zobaczyć, jak duża jest różnica ... czy istnieje gdzieś implementacja omegi kwadratowej dla R (lub python? Wiem ... można marzyć;) Wyszukiwanie w Internecie rzeczy związanych z R jest ból * , nie wiem, jak udaje mi się znaleźć rzeczy z C.
dzięki!
r
anova
effect-size
split-plot
levesque
źródło
źródło
Odpowiedzi:
Funkcja obliczania kwadratu omega jest łatwa do napisania. Ta funkcja pobiera obiekt zwrócony przez test aov i oblicza i zwraca i omega do kwadratu:
edycja: zaktualizowana funkcja dla modeli n-way aov:
źródło
Niedawno musiałem zgłosić .ω2
Jest to brudna funkcja, którą można łatwo wyczyścić. Oblicza częściowy i prawdopodobnie powinien być stosowany tylko w układach silnych między podmiotami.ω2
źródło
Sugerowałbym, że uogólniony kwadrat eta jest uważany ( ref , ref ) za bardziej odpowiednią miarę wielkości efektu. Jest zawarty w wyjściu ANOVA w pakiecie ez dla R.
źródło
Znalazłem funkcję kwadratu omega w czyimś profilu .Rprofile, który udostępnili online:
http://www.estudiosfonicos.cchs.csic.es/metodolo/1/.Rprofile
źródło
Pakiet Daniela „wzmacniejacke” Lüdecke
sjstats
nie może wykonywać kwadratów omega-kwadrat, częściowo-omega-kwadrat itp. Dla modeli ANOVA. Sprawdź to.Oto winieta, która pokazuje, że:
https://cran.r-project.org/web/packages/sjstats/vignettes/anova-statistics.html
źródło