Rでノンパラメトリックのwilcox.testを行うと数値に同値があるばあい”タイがあるため,正確な p 値を計算することができません”と言ってくる.どうやら正規近似したものを出してくるらしい.
正確なP値を計算するにはexactRankTestsパッケージかcoinパッケージを使うといいらしい.
ちなみにcoinパッケージはexactRankTestsの後継になるようだ.
でもexactRankTestsの方が単純な関数ですむので,プラグラムに慣れていない人には楽.
exactRankTestsパッケージでやってみると,,,,
wilcox.exact(x,y)
これだけ.
で,対応があるような場合は,
wilcox.exact(x,y,paired=TRUE)
ちなみに前者が順位和検定で後者が符号順位和検定となる.
コメント