如何在R和SAS中获得相同的Wilcoxon符号秩检验结果

如何在R和SAS中获得相同的Wilcoxon符号秩检验结果,r,sas,R,Sas,R代码: 结果:S=19.5pr>=|S | 0.0156 如何获取R中的统计信息 如果是n20,结果是不同的 data aaa; set aaa; diff=x-y; run; proc univariate; var diff; run; x表示SAS没有对n>=20进行精确测试……这就是R计算“V”统计信息的方式(假设没有缺失值):aAlso,您似乎在程序之间有不同顺序的参数。将(y,x)传递给R,使其计算第一个减去第二个参数(或者在本例中为y-x),然后在SAS>中执行x-

R代码:

结果:S=19.5pr>=|S | 0.0156

如何获取R中的统计信息

如果是n20,结果是不同的

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

x表示SAS没有对n>=20进行精确测试……这就是R计算“V”统计信息的方式(假设没有缺失值):
aAlso,您似乎在程序之间有不同顺序的参数。将
(y,x)
传递给R,使其计算第一个减去第二个参数(或者在本例中为
y-x
),然后在SAS>中执行
x-y
,当n>=20时,SAS和R中的结果不同。@MrFlick在R中使用(y,x)和(y,x)时,P值相同。
data aaa;
  set aaa;
  diff=x-y;
run;
proc univariate;
  var diff;
run;
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)