ANOVA
Tot nu toe heb je uitvoerobjecten gemaakt om elementen op te slaan en te extraheren voor een aangepaste en meer gedetailleerde presentatie van resultaten. Sommige uitvoerobjecten, zoals die van de functie aov() voor het uitvoeren van een variantieanalyse (ANOVA), kun je echter ook gebruiken als invoer voor andere functies zoals summary() en plot() om aangepaste rapportages en visualisaties te krijgen. De uitvoer van aov() kun je ook als invoer gebruiken voor TukeyHSD() om post-hoc paarsgewijze toetsen uit te voeren.
De gegevensset abaloneKeep en het pakket dplyr zijn voor je geladen.
Deze oefening maakt deel uit van de cursus
R voor SAS-gebruikers
Oefeninstructies
- Bereken
n, het gemiddelde, de standaarddeviatie en de variantie van age per sex. - Voer een ANOVA van
ageuit persex-groep en sla de uitvoer op alsabaov. - Voer
summary()uit opabaov. - Voer
TukeyHSD()post-hoc paarsgewijze toetsen uit opabaov.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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(___)