秩和检验原理解释戳这里:
https://zhuanlan.zhihu.com/p/170354683
对于W值得解释戳上面链接(也可直接看下图):
R语言进行秩和检验的函数
wilcox.test(c(1,3,5),c(5,3,1))#举个例子
最后得到的结果:这里看到有W值和p-value
但是一些时候要把W标准化(比如不同分析中(数据)向量的长度不同----e.g.不同簇中细胞数量不同,不能将W直接比较,这时候就需要标准化)
R语言得到的只是W值,这个时候可以根据公式转化,查阅书籍
Z值根据公式得出:(n值表示样本量)
编写成代码
Z=(abs(W-n*(n+1)/4)-0.5)/sqrt(n*(n+1)*(2*n+1)/24)#秩和检验根据W求Z的公式