ウィルコクソンの符号付順位和検定     Last modified: Aug 25, 2015

例題

 「10 人の被検者について,ある測定値を得た。同じ被検者に対して,1 年後にもう一度測定した。その結果を表 1 に示す。1 年間で母代表値に差があったかどうか検定しなさい。」

表 1.対応のあるデータ
被験者 1 2 3 4 5 6 7 8 9 10
$X_i$ 269 230 365 282 295 212 346 207 308 257
$Y_i$ 273 213 383 282 297 213 351 208 294 238


R による解析:

> x <- c(269, 230, 365, 282, 295, 212, 346, 207, 308, 257)
> y <- c(273, 213, 383, 282, 297, 213, 351, 208, 294, 238)
> wilcox.test(x, y, paired=TRUE)

	Wilcoxon signed rank test with continuity correction

data:  x and y
V = 22, p-value = 1
alternative hypothesis: true location shift is not equal to 0

 警告メッセージ: 
1:  wilcox.test.default(x, y, paired = TRUE) で: 
   タイがあるため、正確な p 値を計算することができません 
2:  wilcox.test.default(x, y, paired = TRUE) で: 
   ゼロ値のため、正確な p 値を計算することができません 

このようなエラーが出たときは,exactRankTests パッケージの wilcox.exact 関数を使えばよい

> library(exactRankTests)
> wilcox.exact(x, y, paired=TRUE)

	Exact Wilcoxon signed rank test

data:  x and y
V = 22, p-value = 0.9766
alternative hypothesis: true mu is not equal to 0


・ 手法の解説ページ
・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI