#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