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
微信打赏，为服务器增加50M流量 
支付宝打赏，为服务器增加50M流量 