R3.6.0版
参考： Gene_ID转换，gokeggreactome|clusterprofilerreactomepa
转录组学习(功能丰富分析) ) ) ) ) )。
下载biocmanager :3360安装(annotation hub )安装包，所有缺少的包都是这样#### 用于库(分析中心)库导入的包library ) org.hs.eg.db (人类注释数据库library )集群配置文件库) DP 库(读取XL )、 read_excel(old_excel.xls ) ) f(c )1) )取所需列) ENSEMBL编号为entrezideg2ensembl=totable ) org.hs.egensensensensense me(ensembl_id=f ) results=merge ) genelists，EG2Ensembl， 转换为by='ensembl_id '，all.x=T ) id=NA.omit(results$gene_id )非na的ENTREZID#GO富集分析EgO-Enrichgo ) org DD showCategory=10，title='Enrichment GO Top10 ' )气泡图barplot ) ego， show category=title=' enrichment go ' (条形图plotgograph ) ego ) go图表，如果看不清楚，左上角显示pdf#KEGG分析#ekk-enrichkegg ) gene=id 请尝试另存为organism='hsa ) (qvalueCutoff=0.05 ) #KEGG富集分析dotplot(ekk，font.size=8) #绘制气泡图browsekegg ) ekk，