始める無料で始める

年齢カテゴリと殻重量カテゴリ

先ほど第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カテゴリに分けた shellcatabaloneKeep に追加します。ラベルは "<= 0.235g""> 0.235g" にしてください。
  • agecatshellcat のクロス集計表 tableageshell を作成し、chisq.test()tableageshell に対して実行します。
  • agecatshellcat のモザイク図を作成します。色は "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 = ___)
コードを編集して実行