Rコマンダーで分割プロットANOVA 「理学療法」Vol28(8)のデータ http://www.hs.hirosaki-u.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検定と同義)