计算OTU间两两相关系数矩阵 ,数据量大时,可应用WGCNA中corAndPvalue, 但p值需要借助其他函数(the Benjamini and Hochberg false discovery rate
(FDR))矫正。具体如何实现呢,望大神赐教!!
可以用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")
如果觉得我的回答对您有用,请随意打赏。你的支持将鼓励我继续创作!