How to get same results of Wilcoxon sign rank test

2019-08-02 07:29发布

R code:

x <- c(9, 5, 9 ,10, 13, 8, 8, 13, 18, 30)
y <- c(10, 6, 9, 8, 11, 4, 1, 3, 3, 10)
library(exactRankTests)
wilcox.exact(y,x, paired = TRUE, alternative = "two.sided")

The results: V = 3, p-value = 0.01562

SAS code:

data aaa;
  set aaa;
  diff=x-y;
run;
proc univariate;
  var diff;
run;

The results: S=19.5 Pr >= |S| 0.0156

How to get statistics S in R?

If n<=20 the exact P was same in SAS and R,but if n>20 the results were different.

x <- c(9, 5, 9 ,10, 13, 8, 8, 13, 18, 30,9, 5, 9 ,10, 13, 8, 8, 13, 18, 30,9,11,12,10)
y <- c(10, 6, 9, 8, 11, 4, 1, 3, 3, 10,10, 6, 9, 8, 11, 4, 1, 3, 3, 10,10,12,11,12)
wilcox.exact(y,x,paired=TRUE, alternative = "two.sided",exact = FALSE)

The results: V = 34, p-value = 0.002534

The SAS results:S=92.5 Pr >= |S| 0.0009

How to get the same statistics S and P value in SAS and R? Thank you!

标签: r sas
0条回答
登录 后发表回答