年齢カテゴリと殻重量カテゴリ
先ほど第2章で、アワビの年齢を 10.5 年未満("< 10.5")と、10.5 年以上("10.5 and older")の2カテゴリに分けた agecat を作成しました。ここでは、殻の重量も 0.235 g 以下と 0.235 g より大きいの2カテゴリに分けてみましょう。
agecat を含む abaloneKeep データセットと、dplyr パッケージは読み込まれています。
この演習はコースの一部です
SAS ユーザーのための R
演習の手順
ifelse()関数を使って、殻重量を 0.235 g を境に上下の2カテゴリに分けたshellcatをabaloneKeepに追加します。ラベルは"<= 0.235g"と"> 0.235g"にしてください。agecatとshellcatのクロス集計表tableageshellを作成し、chisq.test()をtableageshellに対して実行します。agecatとshellcatのモザイク図を作成します。色は"light blue"と"grey"を使い、タイトルは"Mosaicplot of age by shell categories"としてください。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# 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 = ___)