上一篇文章:【群体结构】CLUMPP软件使用
前面讲到利用CLUMPP合并多个K值的结果,这边主要叙述R语言绘制structure图。
代码相对简单,我主要利用R语言的低级函数rect()完成。
输入文件准备
- 包含四列,第1列名字(也可以去掉),第2-4列为K的成分比例。
R语言代码
导入文件名为plot.clumpp.txt
,格式如上。通过循环每一行,将不同的列的值进行绘制方块,达到效果。
plotData<-read.table("plot
上一篇文章:【群体结构】CLUMPP软件使用
前面讲到利用CLUMPP合并多个K值的结果,这边主要叙述R语言绘制structure图。
代码相对简单,我主要利用R语言的低级函数rect()完成。
导入文件名为plot.clumpp.txt
,格式如上。通过循环每一行,将不同的列的值进行绘制方块,达到效果。
plotData<-read.table("plot