ANOVA
Do tej pory tworzyłeś/-aś obiekty wynikowe, by zapisywać i wyodrębniać elementy potrzebne do spersonalizowanej i bardziej szczegółowej prezentacji rezultatów. Niektóre obiekty wynikowe – takie jak te zwracane przez funkcję aov() służącą do przeprowadzania analizy wariancji (ANOVA) – mogą jednak pełnić rolę danych wejściowych dla innych funkcji, np. summary() i plot(), umożliwiając elastyczne raportowanie wyników i tworzenie wizualizacji. Obiekt zwrócony przez aov() można też przekazać do funkcji TukeyHSD(), aby wykonać parami testy post hoc.
Zbiór danych abaloneKeep oraz pakiet dplyr są już wczytane.
To ćwiczenie jest częścią kursu
R dla użytkowników SAS
Instrukcje do ćwiczenia
- Oblicz
n, średnią, odchylenie standardowe i wariancję wieku (age) według płci (sex). - Przeprowadź analizę ANOVA zmiennej
agew grupach wyznaczonych przezsexi zapisz wynik jakoabaov. - Uruchom funkcję
summary()na obiekcieabaov. - Wykonaj parami testy post hoc za pomocą funkcji
TukeyHSD()na obiekcieabaov.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# Compute n(), mean(), sd() and var() of age by sex
abaloneKeep %>% group_by(___) %>% select(___, ___) %>%
summarise(across(everything(), list(mean = ___,
___ = ~sd(.x),
var = ___)),
N = ___)
# Run aov() of age by sex, save as abaov
abaov <- aov(___ ~ ___, data = ___)
# Run summary() of abaov
___
# Perform TukeyHSD posthoc pairwise tests on abaov
TukeyHSD(___)