対応の無い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すると,データセット化する結果を選ぶことが出来ます.
0 件のコメント:
コメントを投稿