0

我想用一个变量获得我的 15k 基因的相关值,我创建了一个循环 - 下面的脚本,但是:

  1. 我对 FDR p 值感兴趣,我尝试添加 adjust="fdr",但它没有做任何事情
  2. 呈现它的最佳方式是什么,获取顶级变量并使用 corplot?

这是我正在使用的代码:

spearman <- pval <- NULL
for(i in 2:ncol(data)){
   rho[i-1]  <- cor.test(data[,1],data[,i], method = "spearman")$estimate
   pval[i-1] <- cor.test(for_cor[,1],for_cor[,i], method = "spearman")$p.value
 }
res <- cbind(spearman,pval)
rownames(res) <- paste(colnames(data)[1],"vs",colnames(data)[2:ncol(data)])
head(res)

我可以手动计算,

p.adjust(mypvalues,method="bonferroni")

那是你会做的吗?

谢谢乔治

4

0 回答 0