correlation - 与数千个变量的相关性 - FDR
问题描述
我想用一个变量获得我的 15k 基因的相关值,我创建了一个循环 - 下面的脚本,但是:
- 我对 FDR p 值感兴趣,我尝试添加 adjust="fdr",但它没有做任何事情
- 呈现它的最佳方式是什么,获取顶级变量并使用 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")
那是你会做的吗?
谢谢乔治
解决方案
推荐阅读
- sql - 仅连接该组的更新列
- csv - 将对象数组或 CSV 文件内容插入 Kusto 表
- python - 这行与 __str__(self) 相关的代码是什么意思?
- autohotkey - AHK 中的 #MaxHotkeysPerInterval 错误。如果陈述错误?
- jqgrid - 是什么导致 jqgrid 事件多次触发?
- c++ - const ref 的 I/O 操纵器错误或临时寿命延长?
- powershell - 在 Powershell 的计算字段中使用函数?
- java - 扫描仪无法从 zipinputstream 读取数据
- python - Python argparse 在调用函数时不会接受参数
- c# - 将代码从 Python 转换为 c# 时遇到问题