vcftools --vcf my.vcf --plink --out plink
plink--noweb--fileplink--geno0.1--maf0.05--hwe0.0001--make-bed--outout
forKin1 2 3 4 5;doadmixture --cv XXXXXX.bed $K| teelog${K}.out;done
grep -h CVlog*.out
取最小K值
tbl=read.table("hapmap3.3.Q")
barplot(t(as.matrix(tbl)), col=rainbow(3),xlab="Individual #", ylab="Ancestry", border=NA)
网友评论