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.
Cet exercice fait partie du cours
R For SAS Users
Instructions
- Use
ifelse()
function to addshellcat
toabaloneKeep
for shell weight split into 2 categories above and below 0.235g, labeled as"<= 0.235g"
and"> 0.235g"
. - Create contingency table called
tableageshell
ofchisq.test()
foragecat
byshellcat
. Runchisq.test()
ontableageshell
. - Make a mosaic plot of
agecat
byshellcat
. Use colors"light blue"
and"grey"
. Set title to"Mosaicplot of age by shell categories"
.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# 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 = ___)