cd /mnt/d/peng/climate/paper-climate/climate-1/NJtree
plink --file 72breeds --indep-pairwise 100 25 0.05 --chr-set 29 --out goatld0.2
plink --file 72breeds --indep-pairwise 100 50 0.05 --chr-set 29 --out goatld50
plink --file 72breeds --extract goatld50.prune.in --chr-set 29 --recode --out goatforpca
用spider 软件map,ped转换成Arlequin格式
计算距离，转换成矩阵，在splitstree中构建进化树
格式为：
3
ABR 0 0.08517 0.0922
ANK 0.08517 0 0.07226
ARG 0.0922 0.07226 0
