Age categories by shellweight categories
Earlier, in chapter 2, you created agecat for abalones with ages younger than 10.5 years (coded as "< 10.5") versus those 10.5 years and older (coded as "10.5 and older"). Let's also break the abalone shell weights into categories equal to or lower than 0.235 grams versus those that weigh more than 0.235 grams.
The abaloneKeep dataset with agecat has been loaded for you plus the dplyr package.
Deze oefening maakt deel uit van de cursus
R For SAS Users
Oefeninstructies
- Use
ifelse()function to addshellcattoabaloneKeepfor shell weight split into 2 categories above and below 0.235g, labeled as"<= 0.235g"and"> 0.235g". - Create contingency table called
tableageshellofchisq.test()foragecatbyshellcat. Runchisq.test()ontableageshell. - Make a mosaic plot of
agecatbyshellcat. Use colors"light blue"and"grey". Set title to"Mosaicplot of age by shell categories".
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Add shellcat to abaloneKeep
___ <- ___ %>%
mutate(shellcat = ifelse(___, ___, ___))
# Run chisq.test of agecat by shellcat
tableageshell <- ___ %>%
with(___)
chisq.test(___)
# Make mosaicplot of agecat by shellcat
mosaicplot(___ ~ ___, data = ___,
color = c(___, ___),
main = ___)