2017年9月27日水曜日

wilcoxon順位和検定の話_npar1wayプロシジャ

対応の無い2群間の比較をする際,wilcoxon順位和検定を実行するときはあると思います.
この検定をSASで実行するには,npar1wayプロシジャのwilcoxonオプションを用います.
この検定は「2群間の比較」であることに注意してください.

univariateプロシジャでも2群の比較をすることはできますが,
univariateプロシジャではwilcoxonの「符号」順位和検定を行うことが出来ます.
この検定は今回のwilcoxon順位和検定と違い,対応のある2群間比較になるので割愛します.

/*---------- testdata ----------*/
data hoge ;
    cat = "A"; var1 = 1 ; output ;
    cat = "A"; var1 = 1 ; output ;
    cat = "A"; var1 = 1 ; output ;
    cat = "A"; var1 = 2 ; output ;
    cat = "A"; var1 = 2 ; output ;
    cat = "A"; var1 = 3 ; output ;
    cat = "A"; var1 = 3 ; output ;
    cat = "A"; var1 = 3 ; output ;
    cat = "A"; var1 = 3 ; output ;
    cat = "A"; var1 = 3 ; output ;
    cat = "A"; var1 = 4 ; output ;
    cat = "A"; var1 = 4 ; output ;
    cat = "A"; var1 = 4 ; output ;
    cat = "A"; var1 = 5 ; output ;
    cat = "A"; var1 = 5 ; output ;

    cat = "B"; var1 = 1 ; output ;
    cat = "B"; var1 = 1 ; output ;
    cat = "B"; var1 = 1 ; output ;
    cat = "B"; var1 = 1 ; output ;
    cat = "B"; var1 = 1 ; output ;
    cat = "B"; var1 = 1 ; output ;
    cat = "B"; var1 = 1 ; output ;
    cat = "B"; var1 = 1 ; output ;
    cat = "B"; var1 = 2 ; output ;
    cat = "B"; var1 = 2 ; output ;
    cat = "B"; var1 = 3 ; output ;
    cat = "B"; var1 = 4 ; output ;
    cat = "B"; var1 = 5 ; output ;
    cat = "B"; var1 = 6 ; output ;
    cat = "B"; var1 = 6 ; output ;
run ;

proc npar1way data=HOGE wilcoxon ;
      class CAT ;
      var VAR1 ;
      output out = NPAR ;
run;

正規分布で近似した時の片側P-値は「PR_WIL」
同じ時の両側P-値は「P2_WIL」の変数に格納されています.

結果のデータセットの全体像








おまけですが,npar1wayプロシジャを動かした時の結果の名前の一覧を以下に貼ります.
ods trace on した時に出てくるやつです.
上のプログラムでoutput outで指定すると全ての結果がデータセット化されますが,
結果の名前を指定してods outputすると,データセット化する結果を選ぶことが出来ます.