第4章では統計を扱っています。その中で何度か分散分析のを紹介してます。(repeated-measures-ANOVA)を行う方法2Rで分散分析を行う関数はいくつかあり、今までaov関数、lm関数→anova関数、Anova関数など紹介してきました。RにはANOVA君という分散分析に特化したスクリプトがあります。ANOVA君を使えばどんな形式の分散分析でも行え、多重比較も一気にできます。今回はANOVA君を紹介します。今回は前回のデータを使います。記事の1~4までをご準備ください。前回の記事を見るなんて面倒くさい!という方は以下のコードを実行してください。"B2:E32")View(grip_rep_anova)前回はある運動プログラムを行い、0週目・1週目・2週目で握力を測定した、ということを想定した仮想データとなっています。男女の群分けも行います。2.被験者間要因と被験者内要因と分散分析のデザインについて分散分析を行う上で大切なのが被験者間要因と被験者内要因についてです。被験者間要因はいわゆる対応のないデータのことで被験者全体をA法、B法の2つに割り振るといった方法です。同じ人がA法とB法を行うことはありません。被験者内要因はいわゆる対応のあるデータのことで同じ人が繰り返し測定します。被験者間要因は列の左側に、被験者内要因は列の右側に並べると後々分析しやすくなります。被験者間要因は縦に、被験者内要因は横につなげると分析を行うには都合がいいです。仮に下のデータがあったとします。この並びが分散分析を行う上での基本的な形になります。繰り返し)あるのでABC被験者間要因と被験者内要因の間にsを付ける上記の場合はABsCデザインと言います。被験者間要因と被験者内要因が両方あるデザインを混合計画とも呼んだりします。ANOVA君を使うためにはこのデザインの名付け方とデータの並べ方が鍵になります。3.ANOVA君の導入ANOVA君は通常のインストールでは行なえません。ANOVA君の作者である井関龍太さんのホームページでダウンロードする必要があります。上記記事でtxtファイルを保存します。保存場所は作業フォルダに入れておくとすぐに使えますが、どこでも大丈夫です。作業フォルダはgetwd関数を使えばわかります。()内は何も入れません。Fileでダウンロードしたファイルを指定すればOKです。4.ANOVA君を使うためのデータの準備ANOVA君を使うには次の情報が大切です。1何デザインか?2デザインに併せたデータの並び順にする1何デザインか?今回は以下の列があります。(repeated-measures-ANOVA)を行う方法1ではBのみで分散分析を行いました。この場合は被験者間要因がありませんのでsAデザインとなります。(repeated-measures-ANOVA)を行う方法2では被験者間要因がA被験者内要因がBなのでAsBデザインとなります。分散分析には一元配置とか二元配置とか反復測定とか混合モデルとかいろいろあるのですが、なれない方はまず上記のデザインが付けられることを優先してもいいと思います。その後書籍で二元配置とか出てきた時にそのデータをみながら「あぁ、これは○○デザインのことだな」と理解できればいいと思います。2デザインに併せたデータの並び順にするそしてANOVA君はデザインの順にデータを並べることが必須です。これはtidyverseパッケージのselect関数を使えばできます。今はW0,W1,W2,sexの順に並んでいます。sAデザインならsexを消し、W0,W1,W2のみの列にします。select関数は【2-4】Rで指定した列や行だけを取り出すselect関数、slice関数、filter関数を紹介しますで紹介しています。今回はgrip_anovakun_1_wideという変数名を付けました。tidyverseパッケージを読み込んでから行います。W2)もしAsBデザインならsex,W0,W1,W2の並び順にする必要があります。W2)これでデータができました。5.ANOVA君を使うANOVA君で分散分析を行うにはanovakun関数を使います。各要素のグループの数)ポイントはAやBの要素の数を入れることです。今回の結果をそれぞれ表示します。3)sAデザインと書いています。AがW0,W1,W2だったので、a1,a2,a3と表記されます。それぞれの平均と標準偏差が書いてあります。つぎに反復測定分散分析に必要な球面性の検定を行っています。(repeated-measures-ANOVA)を行う方法1をご参照ください。sAデザインは要素が1つしかないので交互作用はありません。有意差があれば多重比較を行います。多重比較をは何も指定しないとShafferの方法になります。TRUEを付け加えます。ひつような球面仮定性についてみています。次に分散分析表がでています。timeで有意差がでているので時間経過の影響はあるようですが、その傾向は男女で違うようです。timeで有意差が出ていたので男女合わせたtimeで多重比較の結果が出ています。timeの交互作用についてです。女性・男性それぞれ時系列に差があるようです。W2はW0においてFとMに差があるか?という意味ですが今回知りたいことではありません。男女ともに時系列で有意差が見られたのでそれぞれの多重比較の結果が出ています。W2で大きいようです。7.まとめ今回はANOVA君について紹介しました。