1.vcf2ped vcftools --vcf ps.single.210.obs0.8.12135.recode.vcf --plink --out ps.single.210.obs0.8.12135 if Error: Could not open temporary file. terminal: ulimit -n 3000 2.export bed file plink --noweb --file rank --geno 0.05 --maf 0.05 --hwe 0.0001 --make-bed --out rank 3.admixture.sh for K in 11; \ do ./admixture -j20 --cv ./rank.bed $K > ./log/log${K}.out; done 4.plot.sh for i in 11; \ do python2 distruct2.2_admixture.py -K $i --input=./R/rank --output=ps.single.210.12135.203.filter_$i.pdf --title=K=$i --popfile=popmap.txt --poporder=poporder.txt; done