2009/02/17 (Tue) 単純明快でない単純主効果の多重比較

早口ことばみたいだなタイトル。

そしてチョコ祭りも終わったのでてきとーにテンプレ変える。
そのうちぺんぎん仕様を自ら作成すべし！！！

単純主効果の多重比較のしかたについて、明記されてる文献すくねええええ。
ってのが今回の話。
この分野ではバイブルのごとく崇め奉られている『心理学のためのデータ解析テクニカルブック』ですら、HSD（最小有意差）の求め方しか載ってねええええ。
多重比較の本を読んだら載ってるのだろうか。

とりあえず、狩野先生お手製の統計の教科書とテクニカルブックを併せて読みながら、今できる方法をめもっとく。
あ、今回はそれと某先輩のマクロも併用だった。

とりあえず先輩のマクロだけ先に宣伝。
今でも使ってます先輩マジでありがとう。
http://kiso.hus.osaka-u.ac.jp/matusita/index.html
先輩のこの飄々としたかんじがとてもすきだー。

それはさておき多重比較についてはつづき参照。

単純主効果の多重比較が適用されるのは、

2要因以上の分散分析
↓
やっべ交互作用有意だった
↓
単純主効果を見よう！
↓
しまった３つ以上水準があるとこで有意だ！←今ここ

ってケースですね。

とりあえず今回使った対応あり２要因の話だけで。
テクニカルブックにのってるHSDの求め方。
要因Bの水準における要因Aの単純主効果を見る場合。
Q1:要因Aと被験者Sの誤差の自由度と比較する平均の数から、有意になる臨界値をひっぱってくる（付録の表から）
＆
Q2:要因Aと要因Bと被験者Sの誤差の自由度と比較する平均の数から、有意になる臨界値を（ｒｙ
↓
HSD検定の臨界値＝（Q1*誤差A×S）＋（Q2*誤差A×B×S）*（Bの水準数－1)
↓
HSD＝臨界値＊（（Bでプールした誤差項÷被験者数）の平方根）
と、いうことは、
臨界値＝HSD÷（（Bでプールした誤差項÷被験者数）の平方根）
だよね！

だったら
見たい条件間の差÷（（Bでプールした誤差項÷被験者数）の平方根）＝それのスチューデントの値
ってことになるんじゃね！？
つーか実際記述するときってそれが必要になるんじゃね？

って思ったんだけど。
あってんのかな。どうなのかな。

対応なしとか対応なし・ありな計画法でも、細かい修正方法が違うだけで、だいたいの流れはおんなじだから、これを元に単純主効果の多重比較をするRプログラムとか書けるんじゃなかろうか。
Rだったら表なんか参照しなくてもスチューデントの値とか関数で返してくれるしな。
･･･いやまあ実際に書くかどうかはまた別問題ですが。

とはいえ、冒頭で紹介した先輩のマクロがあれば、多重比較に必要な誤差項とかあっさり出してくれるので、それを前提にして臨界値引いてくるだけのプログラムだったらそんなに苦労しない･･･かもしれない。
場合わけがめんどくさかったら、もう実験計画別にプログラム作ればいいわけだし。
うーむ。ありかも？

問題は未だに根強く残るRフォビアなわけですが。
だってCに慣れてなくて統計が不得意ならそらRも怖いわ。
へうー。