2019-01-06

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)