R言語のaov関数とTukeyHSD関数を使ってみる。前回の分析が間違っていたことに気づく。crosshyouスマートフォン用の表示で見るcrosshyou主にクロス表(分割表)分析をしようかなと思いはじめましたが、あまりクロス表の分析はできず。R言語の練習ブログになっています。R言語のaov関数とTukeyHSD関数を使ってみる。前回の分析が間違っていたことに気づく。データ分析www.crosshyou.infoの続きです。比率は全て同じではなくて、違う年がある、ということがわかりました。今回は、R言語のaov関数とTukeyHD関数を使って、どの年とどの年に違いがあるかを調べようと思います。aov関数でANOVAオブジェクトを作成し、TukeyHSD関数を実行します。あれ?エラーになってしまいました。「因子でないものは無視されました」とエラーメッセージがあります。Yearは数値型のままでした。as.factor関数で因子(ファクタ)に変換してからもう一回やってみましょう。p値なので、上の図ではすべて違いがない年の組み合わせですね。年ごとの違いは無い、っていうことですね。oneway.test関数で確認してみます。前回と違うのはなんででしょう?。あ、前回はYearを数値型のままで分析していましたね。今回はYearをファクタにしたので年ごとに有意な差は無いという結果になりました。ねんのため、Yearを数値型にしたYear2という変数を作成してoneway.test関数で分析してみます。とここまで考えて、前回のコマンドを見たら、なんと、コマンドがhiritsuとTypeのままでした。年によって比率の違いは無い、ということだったんですね。前回の分析が間違っていました。それでは気をとりなおして、宗派でのANOVAオブジェクトでTukeyHSD関数を使ってみます。adjがすべての行で0ですから、全てのペアで比率には有意な違いがある、ということです。plot関数でグラフで表しましょう。どの組み合わせも0をまたいでいないので、有意な違いだとわかります。最後に今回のR言語のコマンドを記載しておきます。R言語のaov関数とTukeyHSD関数を使ってみる。前回の分析が間違っていたことに気づく。