R语言绘制弦状图

1.数据

2.程序代码

library(statnet)
library(circlize)

data<-read.csv("Bacteria otu.csv",header = T,row=1)

my.data<-as.matrix(data)

row.names(my.data)<-c("A1","B1","C1","D1")

colnames<-c("Proteobacteria","Actinobacteria","Bacteroidetes","Cyanobacteria","Firmicutes","Planctomycetes","Verrucomicrobia","Chloroflexi","Chlorobi","Thermomicrobia","Others")

grid.col=NULL

grid.col[c("A1","B1","C1","D1")]=c("blue","black","orenge","chocolate")

grid.col[colnames(my.data)]=c("lavender","khaki","mistyrose","siennal","skyblue","brownl","gold","maroon","scalmon","moccasin","wheat")

circos.par(gap.degree=c(rep(2,nrow(my.data)-1),10,rep(2,ncol(my.data)-1),10),start.degree=180)

chordDiagram(my.data,directional=TRUE,diffHeight=0.06,grid.col=grid.col,transparency=0.5)

legend("right",pch = 20,legend = colname(my.data),col = grid.col[colnames(my.data)],bty = "n",cex=1,pt.cex=3,border="black")

pdf(file="circlize.pdf",width=8,height=5,pointsize=8)

chordDiagram(my.data,directional=TRUE,diffHeight=0.06,grid.col=grid.col,transparency=0.5)

legend("right",pch=20,legend=colnames(my.data),col=grid.col[colnames(my.data)],bty="n",cex=1,pt.cex=3,border="black")

3.图形效果