相关系数矩阵

计算OTU间两两相关系数矩阵 ,数据量大时,可应用WGCNA中corAndPvalue, 但p值需要借助其他函数(the Benjamini and Hochberg false discovery rate

(FDR))矫正。具体如何实现呢,望大神赐教!!

请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

可以用R语言中的p.adjust() 函数完成矫正:

p.adjust(p, method = p.adjust.methods, n = length(p))
p.adjust.methods
# c("holm", "hochberg", "hommel", "bonferroni", "BH", "BY",
#   "fdr", "none")
请先 登录 后评论