R语言-秩和检验--根据W值求标准化的Z值

本文介绍了秩和检验的基础知识,包括其原理和R语言中的应用。通过实例展示了如何使用`wilcox.test`函数进行秩和检验,并讨论了在不同样本量下W值的标准化过程。标准化公式为Z=(|W-n*(n+1)/4|-0.5)/sqrt(n*(n+1)*(2*n+1)/24),以适应不同数据向量长度的比较需求。
摘要由CSDN通过智能技术生成

秩和检验原理解释戳这里:
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的公式

  • 1
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值