Rコマンダーで分割プロットANOVA 「理学療法」Vol28（8）のデータ

Presentation on theme: "Rコマンダーで分割プロットANOVA 「理学療法」Vol28（8）のデータ"— Presentation transcript:

1 Rコマンダーで分割プロットANOVA 「理学療法」Vol28（8）のデータ

2 データの読み込み ③Rコマンダーで[データ]-[データのインポート]-[テキスト…]を選び ①範囲指定する （変数名も込みで）
②右クリック，コピーする ②右クリック，コピーする ④ここではデータ名を「data3」とする．「ファイル内に…」にチェックを入れ，「クリップボード」に変更してOK

3 読み込み完了 データが読み込まれたら，ここにデータ名が出る 出ないようなら前スライドの方法を繰り返します

4 群分けのデータを因子変数へ変更する必要はない
このデータでは高齢群が0，若年群が1 いままでは，群を分ける変数を因子変数に指定しないと，解析できなかったが，分割プロットに関しては不要

5 分割プロット分散分析の選択 ① ② ③ メニューから選ぶ

6 注意！ この時点で，対応のない要因（つまり群を分ける要因）が，一番上に表示されないとならない この場合は，「X1ヵ月後」なので修正する
変数名の先頭に　．　ドットを付け加えると良い 一旦キャンセルする

7 変数名の修正① 以上のように選ぶ

8 変数名の修正② 群を選ぶ→OKボタンをクリック

9 変数名の修正③ .群と付ける

10 解析に戻って… .群が，先頭に来る この時点で，全ての変数を選んでOKする 「X1ヵ月後」と「開始時」の順番は問わない

11 結果① 平均と標準偏差を見て，どれがどの変数かを押さえておく
メンドーザMendozaの多標本球面性検定が行われる ここでは有意ではないので，ANOVAの修正が不要である ※そもそも，2要因とも2水準なので，この検定は不要である． 平均と標準偏差を見て，どれがどの変数かを押さえておく ここでは，Aが群（a1＝高齢群，a2＝若年群），Bが治療前後の水準（b1=1ヵ月後，b2=開始時）である．

12 結果② メンドーザの検定結果から，自動的にグリーンハウスカイザーのε修正が行われる
結果では，全ての要因で有意となっている（結果の表示がずれます）． 2要因とも2水準なので，多重比較法は行うまでもない．しかし交互作用が有意なので，全ての水準で差の検定が必要→次スライド

13 結果③ b1におけるAの差，b2におけるAの差，…，と全てのケースでANOVAを行っている 1元配置分散分析である（2標本t検定と同義）