|
楼主 |
发表于 2014-5-31 10:56:01
|
显示全部楼层
#双样本秩和检验
rank=0;
a=x-y;
b=0;
k=1;
n=length(a)
a=a[order(abs(a))];
for(i in 1:n){
if(a>0){
rank=rank+i;
b[k]=a;
k=k+1;
}
}
#统计与w+绝对值相同数出现的次数
s=rep(0,length(b))
for(i in 1:length(b)){
for(j in 1:length(a))
if(b==abs(a[j])){
s=s+1; }
}
for(i in 1:length(s))
rank=rank-(s-1)/s;
p=2*min(psignrank(rank,n),1-psignrank(rank,n))
cat("\nP值为:",p)
整个儿的程序是这样的 |
|