用david轉(zhuǎn)換成ENTREZID基因
gene<-read.csv(file="id.csv",header = T)
gene_diff<-gene$To#ENTREZID基因
kk.diff <- enrichKEGG(gene = gene_diff,
organism = 'mmu',
#可以更換物種雹有，這是小鼠的嬉探，人類是hsa
pvalueCutoff = 0.05)
head(kk.diff)[,1:6]
dotplot(kk.diff );ggsave('kk.diff.dotplot.png')#點(diǎn)圖
kegg_diff_dt <- as.data.frame(kk.diff)
diff_kegg<-kegg_diff_dt[kegg_diff_dt$pvalue<0.05,];diff_kegg$group=-1
對(duì)KEGG富集結(jié)果可視化
kegg_plot <- function(diff_kegg){
dat<-diff_kegg
colnames(dat)
dat$pvalue = -log10(dat$pvalue)
dat$pvalue=dat$pvalue*dat$group
dat=dat[order(dat$pvalue,decreasing =F),]
g_kegg<- ggplot(dat, aes(x=reorder(Description,order(pvalue, decreasing =T)), y=pvalue, fill=group)) +
geom_bar(stat="identity") +
scale_fill_gradient(low="blue",high="red",guide = FALSE) +
scale_x_discrete(name ="Pathway names") +
scale_y_continuous(name ="log10P-value") +
coord_flip() + theme_bw()+
theme(text = element_text(size=8),plot.title = element_text(hjust = 0.1))+
ggtitle("Pathway Enrichment")
}
g_kegg<-kegg_plot(diff_kegg)
print(g_kegg)
ggsave(g_kegg,filename = 'kegg_up_down.pdf')
畫出來選定某一通路的圖形
BiocManager::install("pathview")
require(pathview)
geneList<-gene#選定的基因
pathview(gene.data = "genelist", pathway.id = "hsa04110",
species="hsa",limit =list(gene=max(abs(geneList)),cpd=1))
畫出來選定某一通路的圖形奸汇，我們的基因集中在該通路出現(xiàn)的基因變成綠色顯示在該通路中
hsa04110.pathview.png
