群体结构：ld分析和画图

#PBS -N q338
#PBS -l nodes=1:ppn=8
#PBS -q middle

cd $PBS_O_WORKDIR


./lddecay.sh 338

lddecay.sh如下

#!/bin/bash
#lddecay,annovar
if [ ! -d "lddecay" ]
then 
    mkdir lddecay
    /public/home/caisl/lee/soft/PopLDdecay-3.41/bin/PopLDdecay -InVCF ${1}.filter.vcf -OutStat LD
    gunzip LD.stat*
    /public/home/caisl/lee/soft/PopLDdecay-3.41/bin/Plot_OnePop.pl -InFile LD.stat -output all.graph
    mv LD* lddecay
    mv all.graph* lddecay
else
    echo "lddecay has finished, next"
fi