require(DOSE)
require(clusterProfiler)
data(geneList)
gene = names(geneList)[abs(geneList) > 2]
david = enrichDAVID(gene = gene, idType=“ENTREZ_GENE_ID”,
listType=“Gene”, annotation=“KEGG_PATHWAY”)
barplot(david)
比较不同cluster基因的富集
data(gcSample)
x=compareCluster(gcSample, fun=“enrichDAVID”, annotation=“KEGG_PATHWAY”)
plot(x)
如若转载,请注明出处:https://www.ouq.net/1783.html