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