R语言ggplot2画图Pathway富集分析气泡图热文聚热

　　转载于http：www。omicshare。comforumthread14611。htmlOmicShareForum专业的生物信息学论坛安装包install。packages（ggplot2）导入包帮助文档链接：library（ggplot2）设置工作路径到数据存放的文件夹下读数据pathwayread。table（R0vsR3。path。richFactor。head20。tsv，headerT，sept）画图ppggplot（pathway，aes（richFactor，Pathway））ppgeompoint（）改变点的大小ppgeompoint（aes（sizeR0vsR3））四维数据的展示pbubbleppgeompoint（aes（sizeR0vsR3，color1log10（Qvalue）））自定义渐变颜色pbubblescalecolourgradient（lowgreen，highred）绘制pathway富集散点图prpbubblescalecolourgradient（lowgreen，highred）labs（colorexpression（log〔10〕（Qvalue）），sizeGenenumber，xRichfactor，yPathwayname，titleTop20ofpathwayenrichment）改变图片的样式（主题）prthemebw（）保存图片ggsave（out。pdf）保存为pdf格式，支持pdf，png，svg多重格式ggsave（out。png）保存为png格式ggsave（out2。png，width4，height4）设定图片大小备注，绘图数据的说明：（1）绘图数据来自我们公司KEGG富集分析的结果，相应文件是结题报告中存在的，略作调整即可；（2）绘图数据每一列的意思：1）Pathway：通路的名称2）R0vsR3：差异表达基因中，属于这个通路的基因的数量3）AllUnigene：所有基因中属于这个通路的基因的数量4）Pvalue：富集分析p值5）Qvalue：富集分析的Q值6）richFactor：在我们分析报告中，没有提供这一列，但很容易计算。是第二列除以第三列得到；7）PathwayID：通路ID8）Genes：通路中基因的ID9）KOs：通路中基因的KO号补充一点：绘图仅仅用到4类，分布是第1，2，5，6列再补充一点（来源25，26，28楼的讨论）（1）在pathway名称中如有重名，这会导致错误。在表中每个pathway只能出现一次；（2）文本中，出现了引号会导致错误。例如，Alzheimersdisease，Huntingtonsdisease这样的名称。这两个pathway名称中的引号需要删除。引号的出现，会导致R无法识别两个引号间的其他符号（退格符，换行符等），导致文件读取错误。如果一定要保留引号，则引号的内容再用引号囊括起来，例如：Alzheimersdisease，Huntingtonsdisease，从而避免单个出现的引号对其他字符的影响。