我用enricher函数做GO富集分析,代码如下 x <- enricher( query_gene, pvalueCutoff = 0.05, qvalueCutoff = 0.2, pAdjustMethod = "BH", universe = background, TERM2GENE = go2gene_bp, TERM2NAME = des ) x@ontology <- 'BP' GO_df <- x@result write.table(GO_df, '...
Plot Total Number of diffrentially Abundant Features (applying ratio cutoff) vs. qValue/pValue for all conditions