structure 分析

==============转化vcf为plink格式=============

vcftools --vcf all.var.flt.vcf --plink --out out

plink --file out --make-bed --out hapmap --noweb

======structure分类=====================

admixture hapmap.ped 5  //分成5类的情况 可以调整数字进行不同的类别分类

hapmap.5.Q

然后调整顺序并利用R画图

tbl_5=read.table("hapmap.5.Q")

barplot(t(as.matrix(tbl_5)), col=rainbow(5),ylab="Ancestry", border=NA,space=0)

K=3和K=4的结果
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容