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検定と同義）