R语言实现定性资料的秩和检验

介绍

定性资料好比等级，毒性，应答等，能够以具备分级的因子的形式表示，好比(+ ++， +++)，分别对应因子的1,2,3种水平，这样不一样组样本只要看这些数据的等级的排列是否一致就能够判断这两个群体的分布是否有差别(秩和检验)。

数据

好比数据以下
咱们这里想要知道A、B处理毒性是否有差别

代码

library(data.table) 
setwd("/个人文档/project/test")

dt <- fread("example.tsv") #读入文件，这里用fread依赖data.table包，也能够用read.table
dt$分组 <- factor(dt$分组)  # 分组转变成因子
dt$毒性 <- as.numeric(factor(dt$毒性))  #定性资料转变成因后再转成能够计算的数字

# 进行检验计算p值和zscore
pvalue <- wilcox.test(毒性 ~ 分组, data = dt)$p.value
zscore <-  qnorm(pvalue/2)   # 默认wilcox双尾检验因此要除以2，若是你wilcox作的是单尾检验，这里不用除以2
print(zscore)

参考资料

https://stats.stackexchange.com/questions/330129/how-to-get-the-z-score-in-wilcox-test-in-r