ノンパラメトリック検定で多重比較したいとき
現在の（私が確認しうる限りの）SPSSでは3群以上のノンパラメトリック検定を行なった後に多重比較検定を行なう手順がありません．
SPSSはとても普及している統計ソフトですので，このソフトが持つ機能が統計処理の基準だと思われている節があります．
てっきり，「ノンパラメトリック検定は3群以上での多重比較ができない」
「3群以上の比較は分散分析をしないといけない」
と考えている学生や院生が多いようです．
今回はSPSSとExcelの両方を合わせて使って，ノンパラメトリック検定での多重比較をやってみます．
ノンパラメトリック検定（nonparametric test）というのは，
ｔ検定や分散分析といったパラメトリック検定（parametric test）とは違い ，
「母集団の正規性や等分散性を仮定しない検定」
のことです．
実験ではN数（対象数）が5とか，多くても20なんていうことがしばしばあります．
でも，こんな少数では母集団の正規性が認められることは少なく，最近では統計処理手法としてノンパラメトリック検定を薦める記述を目にすることが多いですね．
「母集団の正規性，等分散性」 なんていう意味不明なことは後で勉強しとくから，とりあえずノンパラメトリック検定での多重比較を教えてほしい，という人を対象にして以下に説明します．
例として取り扱うのは以下のデータです．
SPSS for Windows ver11.0を使用しています．
※後日，エクセルでもこの作業ができるよう，
■クラスカル・ウォリスの検定をエクセルでやる
■フリードマンの検定をエクセルでなんとかする
を取り上げましたので，ご参照ください．
クリックしたら画像が元の大きさで見れます．
これをノンパラメトリック検定の中にある「K個の独立したサンプルの検定」を選択すると，
以下のような画面がでるので，Kruskal-Wallisを選択して検定します．
すると，以下のような検定結果が出て，有意であることが認められました．
なので，多重比較を行なうことになります．
以前の記事ではt検定を用いた多重比較を紹介しました．
参考記事◆対応のあるt検定と多重比較
それと基本的には変わりません．
ノンパラメトリック検定であっても，２群比較を行ない，ボンフェローニ（Bonferroni）の方法でP値を修正をすればOKです．
SPSSにおけるノンパラメトリック検定での2群間の比較は以下のようにします．
※後日，エクセルでノンパラメトリック検定をする方法を記事にしました．
■マン・ホイットニーのU検定（エクセルでp値を出す）
■ウィルコクソンの符号付順位和検定（エクセルでp値を出す）
ご参照ください．
そしてOKを押します．
それで検定結果はこんなふうに出ます．正確有意確率というところを見るんだそうです．
0.690となっていますね．
有意ではありません．
SPSSの 「2個の独立サンプルの検定」 はグループ化変数をひとつずつ変えてやっていかなければなりません．
4群なので，全部で6通りの組み合わせです．
それを全部やって，結果をExcelに貼り付けるんですが．
そのExcelに貼り付ける方法はというと，まずは以下のように 「検定統計量」 のところをクリックして四角い枠をつけた状態にします．
ここでは最後の３*４の組み合わせをコピーし，あとの部分は邪魔なので消しました．
（１）ボンフェローニの方法
D列２行目に以下のような式を放り込みます．
=B2*6
多重比較する組み合わせの数をｐ値に乗じるだけです．
今回は4群なので６を乗じていますが，3群であれば３を，5群であれば10を乗じます．
これをオートフィルしてしまうのもいいのですが，今回はボンフェローニの検定以外にも，Excelで簡単にできるP値の多重比較補正を紹介します．
（２）サイダック（Sidak）の方法
次はSidakの検定です．サイダックの検定（もしくはシダックの検定）と読みます．
数式は以下のとおりです．E列２行目に入れました．
=1-(1-B2)^(6/1)
たしかにボンフェローニの検定よりはｐ値を甘く出してくるようです．
今回は4群なので，べき乗する値を「6/1」にしていますが，3群であれば「3/1」に，5群であれば「10/1」にします．
これをオートフィルします．
（３）ライアン（Ryan）の方法
次はライアン（Ryan）の方法を紹介します．
ライアンの方法のメリットは，ボンフェローニやサイダックといった修正方法よりも群数が多くても有意になりやすことです．
オートフィルで一気に，とはいけませんが，比較的簡単に計算できるので便利です．
この方法は平均値の大きさで多重比較する順番を決め，その順番ごとに計算する際の数値を少しずつ変える，というもの．
まず，以下をご覧ください．
グループ１，２，３，４の順で大きくなっています．
ライアンの方法では，平均値が最大・最小の群から多重比較を始めます．
「１*４」がそれに該当しますので，そのP値を修正する計算式を以下に示します．
=B4*6
次は最大・第２最小，最小と第２最大を比較します．
例では「２*４」，「１*３」に該当します．
「２*4」は，
=B6*4
=B3*4
最後に，残りの組み合わせを以下のように計算します．
画像は省略します．
「１*２」は，
=B2*2
「２*３」は，
=B5*2
「３*４」は，
=B7*2
ライアンの方法を使う上で重要なのは，ｐ値に乗ずる数値です．
それは，以下のようにして決めてください．
＝(群数×（（群数－それまでの検定回数）－１）÷２
（群数－それまでの検定回数） という数学や統計学が強い人からするとクールではない表現をしていますが，ここが肝です．
１回目の検定であれば０ということで計算してください．
具体例がないとイメージしにくいでしょう．
例えば４群であれば，１回目の検定は，
＝（４×（４－１）÷２
であり，P値に６をかけます．
２回目の検定では，
＝（４×（３－１）÷２
で，P値に４をかけます．以下，３回目は・・・，と続きます．
これが３群の検定であれば，１回目の検定では
＝（３×（３－１）÷２
でP値に３をかけ，２回目の検定では，
＝（３×（２－１）÷２
なので1.5をかける，ということです．
どうしてそんな計算式なのか？というのは，ライアンの方法を紹介した書籍やサイトをみてください．
ここではあくまで，目の前にあるデータを処理したい人用のサイトですので，統計学としての解説はすっ飛ばします．
全部を計算して，ボンフェローニとサイダック，そしてライアンを比べてみました．
ボンフェローニよりもサイダックのほうがやや甘め．
ライアンはさらに甘めです．
「２*３」のグループの組み合わせでは，ボンフェローニとサイダックでは有意差なしと判断されていますが，ライアンでは有意差ありということになっています．
一気にバーっと計算しちゃいたいならボンフェローニとサイダック．
一個ずつでも，有意差を出していきたいならライアン，ということになるでしょうか．
もちろん，今回紹介したP値の修正方法はノンパラメトリック検定にしか適用できないわけではありません．
t検定による多重比較にも適用できます．
対応のある／なしも関係ありません．
Excelだけで多重比較したい場合はこれらを利用するとよいでしょう．
後日,
■Excelで多重比較まとめ
に文字通りまとめましたので，こちらもご参照ください．
統計的有意にこだわらないのであれば，
■効果量（SE：effect size）をエクセルで算出する
次回は，Excelで膨大な量のデータを計算する上でのヒントを紹介します．
Amazon広告
※後日，こんな怪しいブログよりも信頼性が高いものに触れてもらうよう，
■独学で統計処理作業をスキルアップさせるための本
という記事を書きました．参照してください．■独学で統計処理作業をスキルアップさせるための本
外部サイトにも有益なリストがあります．こちらも参考にしてください．
■大学生が自力で「統計学」の勉強をするための良書１０選
■1ヶ月で統計学入門したので「良かった本」と「学んだこと」のまとめ
手計算で算出するのが面倒な人は，思い切ってエクセル統計の購入をオススメします．
エクセルや手計算で多重比較をしたい方は，以下の２冊がオススメです．
