miRNA芯片数据挖掘-GO圈图

install.packages("digest")
install.packages("GOplot")

library(GOplot)
setwd("F:\\")              #设置工作目录

ego=read.table("GO.txt", header = T,sep="\t",check.names=F)   #输入文件  
go=data.frame(Category = "All",ID = ego$ID,Term = ego$Description, Genes = gsub("/", ", ", ego$geneID), adj_pval = ego$p.adjust)


id.fc <- read.table("id.txt", header = T,sep="\t",check.names=F) #输入文件
genelist <- data.frame(ID = id.fc$id, logFC = id.fc$logFC)
row.names(genelist)=genelist[,1]

circ <- circle_dat(go, genelist)
termNum = 5                                     
geneNum = nrow(genelist)                       

chord <- chord_dat(circ, genelist[1:geneNum,], go$Term[1:termNum])
pdf(file="circ.pdf",width = 13,height = 11)
GOChord(chord, 
        space = 0.001,           
        gene.order