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