使用的get113.csv文件格式:
Description,GeneNumber,GO_Function
biologl_pro,8,biological_process
…
biologl_pro,8,biological_process
画图前先对所用数据进行排序处理
横版的数据为MF,CC,BP;每组数据都是GeneNumber升序,代码如下:
dat <- read.csv("gst11 vs wt go.csv",header=TRUE)
bp=dat[which(dat$GO_Function %in% "biological_process"),]
library(dplyr)
newbp=arrange(bp,bp[,2])
cc=dat[which(dat$GO_Function %in% "cellular_component")