写在前面——书接上回,通过绘制差异基因火山图,能够看出上下调基因的分布情况。这次我们通过对差异基因的GO富集分析,可以看到涉及到的具体通路,更进一步的了解实验变化。本文使用的数据集为私有数据集,不过绘图并不难,弄懂原理之后,套用在自己的数据上即可。
library(readxl)
library(dplyr)
library(ggplot2)
# 整理数据格式
bp_up_data <- read_xlsx("go/up_bp.xlsx")
bp_down_data <- read_xlsx("go/down_bp.xlsx")
bp_up <- bp_up_data[,c(2,5,11)]
bp_up$change <- c("up")
bp_down <- bp_down_data[,c(2,5,11)]
bp_down$change <- c("down")
bp <- rbind(bp_up, bp_down)
colnames(<