head(dat,3) dat$p<- pnorm(q=dat$z_score, lower.tail=F)## 是否大于某个值,使用右侧检验head(dat,3) 。 003、对P值进行FDR矫正 dat <- read.table("fst.fst", header =T) head(dat,3) dat$z_score<- (dat$FST - mean(dat$FST))/sd(dat$FST) head(dat,3) dat$p<- pnorm(q=dat...