例題:
「10 人の被検者について,ある測定値を得た。同じ被検者に対して,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