参考来自 https://mp.weixin.qq.com/s/2eNKtz4j0STKcANOqSEHaQ
library(GOplot)
circ <- circle_dat(ECgenelist)
chord <- chord_dat(circ, ECprocess)
head(chord)
423c71d153ff99ce0c3a8e710e1e37a2.png
d46bda39e0982b112beb03bf756f8262.png
作图的数据是这样的,我们模仿这个数据准备自己的数据
我们需要用到的是Term和genes两列
df<-read.csv("D:/Jupyter/GOplot/example.txt",
header = T,
sep="\t",
row.names = 1,
check.names = F)
df<-df[,1:7]
df$logFC<-sample(c(-1,1),dim(df)[1],replace = T)
GOChord(as.matrix(df), space = 0.02, gene.order = 'logFC', gene.space = 0.25, gene.size = 5)
f0ff283cc87f579f8c0d08b10431fc12.png