Transcript 分割プロットデザインのANOVAの場合

Rコマンダーで分割プロット
ANOVA
「理学療法」Vol28（8）のデータ
http://www.hs.hirosakiu.ac.jp/~pteiki/data/data3.xls
データの読み込み
③Rコマンダーで[データ]-[データのインポート][テキスト…]を選び
①範囲指定する
（変数名も込みで）
②右クリック，コピーする
②右クリック，コピーする
④ここではデータ名を「data3」とする．「ファイル内に
…」にチェックを入れ，「クリップボード」に変更してOK
読み込み完了
• データが読み込ま
れたら，ここにデー
タ名が出る
– 出ないようなら前
スライドの方法を
繰り返します
群分けのデータを因子変数へ変更する必要はない
• このデータでは高齢群が0，若年群が1
• いままでは，群を分ける変数を因子変数に指定しないと，解
析できなかったが，分割プロットに関しては不要
分割プロット分散分析の選択
①
②
• メニューから選ぶ
③
注意！
• この時点で，対応のない要因（つまり群を分ける要因）が，一番上に表示
されないとならない
• この場合は，「X1ヵ月後」なので修正する
• 変数名の先頭に ． ドットを付け加えると良い
• 一旦キャンセルする
変数名の修正①
• 以上のように選ぶ
変数名の修正②
• 群を選ぶ→OKボタンをクリック
変数名の修正③
• .群と付ける
解析に戻って…
• .群が，先頭に来る
• この時点で，全ての変数を選んでOKする
• 「X1ヵ月後」と「開始時」の順番は問わない
結果①
•
•
•
メンドーザMendozaの多標本球面性検定が行われる
ここでは有意ではないので，ANOVAの修正が不要である
※そもそも，2要因とも2水準なので，この検定は不要である．
平均と標準偏差を見て，どれ
がどの変数かを押さえておく
ここでは，Aが群（a1＝高齢
群，a2＝若年群），Bが治療
前後の水準（b1=1ヵ月後，
b2=開始時）である．
結果②
•
•
•
メンドーザの検定結果から，自動的にグリーンハウスカイザーのε修正が行われる
結果では，全ての要因で有意となっている（結果の表示がずれます）．
2要因とも2水準なので，多重比較法は行うまでもない．しかし交互作用が有意なので，全ての水準で差の検
定が必要→次スライド
結果③
•
b1におけるAの差，b2におけるAの差，…，と全てのケースでANOVAを行っている
– 1元配置分散分析である（2標本t検定と同義）