R言語で一元配置分散分析(oneway.test関数)と交互作用図(interaction.plot関数)をしてみる。www.crosshyou.infowww.crosshyou.info政府統計の総合窓口のデータや、OECDのデータを使って、Rの練習をしています。ときどき、読書記録も載せています。R言語で一元配置分散分析(oneway.test関数)と交互作用図(interaction.plot関数)をしてみる。データ分析www.crosshyou.infoの続きです。都道府県知事の所轄の団体数と文部科学大臣の管轄の団体数の比率が宗派によって違うのかを調べてみます。mutate関数で比率を表す変数を作成しましょう。キリスト教系は13.8ぐらいですね。tapply関数とmean関数で宗派ごとのhiritsuを計算します。いうこのchidaというデータフレーム内には存在しないファクタ水準が残っていますね。as.character関数で文字列に変換し、as.factor関数でファクタに変換します。した出力になりましたね。仏教系は177.6と宗派によって全然違いますね。分散分析で宗派によって比率が違うかどかを検定します。この場合は一元配置分散分析でしょうか。oneway.test関数を使います。なっていますから、各宗派の比率はすべて等しいという帰無仮説は棄却されました。宗派によって比率はことなるということです。同じようにYearによってhiritsuが異なるか見てみましょう。2005年は186、2017年は174と違うようです。グラフにしてみましょうか。barplot関数を使いました。比率が低下しているようですね。それではoneway.test関数で一元配置分散分析をします。Yearの違いでhiritsuに違いがあることがわかりました。hiritsuはType(宗派)にもYear(年)にも影響されているということですね。それでは次に、TypeとYearの交互作用図を作成してみましょう。intaraction.plot関数を使います。4本の線は、ほぼ同じような形状なので年と宗派の交互作用は無いようですね。見るの「レシピ8.20一元配置分散分析を行う」と「レシピ8.21交互作用図を作成する」を参考にしました。次回はwww.crosshyou.infoです。最後に今回のR言語のコマンドを記載しておきます。睡眠、身の回りの用事、食事の時間は女性のほうが長い。R言語で法人と個人の事業所の割合を比較する。はてなブログPro主にクロス表(分割表)分析をしようかなと思います。5引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。