【群体结构】structure结构图绘制R语⾔绘制structure图
临渊而立
上⼀篇⽂章:
前⾯讲到利⽤CLUMPP合并多个K值的结果,这边主要叙述R语⾔绘制structure图。
代码相对简单,我主要利⽤R语⾔的低级函数rect()完成。
美白有效
输⼊⽂件准备
包含四列,第1列名字(也可以去掉),第2-4列为K的成分⽐例。
R语⾔代码besure
导⼊⽂件名为,格式如上。通过循环每⼀⾏,将不同的列的值进⾏绘制⽅块,达到效果。
plotData<-read.table("",header=T)
plot(c(0,241),c(0,1),axes=F,xlab="",ylab="",type="n")
col1<-c('blue','green','red')
for(i in 1:241){五月桃花>名声的近义词
#利⽤累加数值来改变y轴绘图的位置,来达到堆积的效果
rect(i-1,plotData[i,2],i,0,col=col1[1],border=NA)
rect(i-1,plotData[i,2]+plotData[i,3],i,plotData[i,2],col=col1[2],border=NA)
rect(i-1,1,i,plotData[i,2]+plotData[i,3],col=col1[3],border=NA)
# K值较⼤,需要增加rect,这边适⽤于K=3
姜汉
}
成品图如下:可爱短发发型
>11月2号是什么星座