(repeated-measures-ANOVA)を行う方法2第4章は統計を扱います。今回も「シロート統計学」のハルさんとコラボレーションすることとなりました。シロート統計学はEZRを使った統計分析をわかりやすく解説されています。第4章はシロート統計学で使われていたEZRをRで行うとどうなるのか?といった視点で進めていきます。今回使うデータもハルさんのサイトと同じものを使わせでいただく事になりました。それぞれ見比べることで参考にしてみてください!まずはハルさんのサイトをご覧ください。ここで行っている解析をRで実施する方法を紹介します。今回は前回のデータを使います。記事の1~4までをご準備ください。前回の記事を見るなんて面倒くさい!という方は以下のコードを実行してください。"B2:E32")View(grip_rep_anova)前回はある運動プログラムを行い、0週目・1週目・2週目で握力を測定した、ということを想定した仮想データとなっています。今回は握力の経時的検討に加えて男女の群分けも行います。2.被験者間要因と被験者内要因と分散分析のデザインについて分散分析を行う上で大切なのがvについてです。被験者間要因はいわゆる対応のないデータのことで被験者全体をA法、B法の2つに割り振るといった方法です。同じ人がA法とB法を行うことはありません。被験者内要因はいわゆる対応のあるデータのことで同じ人が繰り返し測定します。被験者間要因は列の左側に、被験者内要因は列の右側に並べると後々分析しやすくなります。被験者間要因は縦に、被験者内要因は横につなげるとEZRで分析を行うには都合がいいです。仮に下のデータがあったとします。この並びが分散分析を行う上での基本的な形になります。繰り返し)あるのでABC被験者間要因と被験者内要因の間にsを付ける上記の場合はABsCデザインと言います。被験者間要因と被験者内要因が両方あるデザインを混合計画とも呼んだりします。今回は以下のようになります。被験者間要因:sex(A)被験者内要因:W0,W1,W2の繰り返し(B)今回のExcelではsexの列が1番右にありますが被験者間要因を左、被験者内要因を右に書くのが習わしのようなのでAsBデザインということになります。被験者内要因(改善の傾向)に性別が影響しているかどうかを見たいということになります。行うには主に2つの方法があります。1.carパッケージのAnova関数を使う(EZRもcarパッケージを使ってる)2.ANOVA君を使うANOVA君は井関龍太さんが作成した分散分析に特化したスクリプトです。使い方に特徴がありますが使えるようになると慣れると簡単でどの方法よりも詳細に分析できます。行うには次の方法になります。20191carパッケージを読み込むここは前回の記事と全く同じです。carパッケージを読み込みます。carパッケージはEZRをインストールする時に一緒にインストールされています。インストールしていなければ先にインストールします。install.packages("car")その後library関数で呼び出します。library(car)2被験者要因を指定するdata.frameを準備するここは前回の内容と全く同じです。今回のデータではsexの列は「MとFの2つのデータがある」とパソコンは読み取れますがW0,W1,W2はwideデータとなっているため「この3列が繰り返しのデータである」ことがパソコンは読み取れません。それがわかるためのものを2つ作ります。繰り返しの順番がわかるtimeというベクトルを作ります。このデータはfactor型である必要があるためfactor関数を使います。timeのデータをdata.frame型に変更したidataを作ります。指定しないとfacorは五十音順になります。使わないとmid→post→preの順番になってしまいます。factor関数に関しては【1-10】Rでよく使われる型について説明します。3lm関数でモデルを作る次に分析するためのモデルを作ります。ここが前回の記事と唯一違うところです。モデルはlm関数を使います。model(名前は何でもいい)という変数名に結果を入れます。被験者内要因は3列になっているのでcbind関数でつなげます。被験者間要因は~の右に入れます。つなぎます(交互作用を見たい場合)。