[code]
set.seed(1234) #设置种子点
p <- runif(10,min=0,max=0.1)
p <- sort(p)
p_adj0 <- p.adjust(p,method = 'bonferroni')
p_adj1 <- p.adjust(p,method = 'holm')
p_adj2 <- p.adjust(p,method = 'BH')
p0 = length(p)*p
p1 = (length(p) - rank(p) + 1) * p
p2 = length(p) * p / rank(p)
pbind = cbind(p,p0,p_adj0,p1,p_adj1,p2,p_adj2)
pbind
|