# step1: calculate the LDAK weights

ldak=./ldak5.linux
bfile=xx

$ldak --cut-weights sections --bfile $bfile --no-thin YES
num_sections=$(cat sections/section.number)
qsub -t 1-$num_sections ldak_calc_weights.sh

ldak_calc_weights.sh
```shell
#$ -cwd
#$ -m a
#$ -l h_data=8G,h_rt=0:30:00
#$ -j y
#$ -o ./job_out

ldak=./ldak5.linux
bfile=xx

$ldak --calc-weights sections --bfile $bfile --section $SGE_TASK_ID
```
## After the above finishes
$ldak --join-weights sections --bfile $bfile

# step2: calculate tagging based on the weights

qsub -t 1-22 ldak_calc_tagging.sh

```shell
#$ -cwd
#$ -m a
#$ -l h_data=10G,h_rt=6:30:00
#$ -j y
#$ -o ./job_out


ldak=./ldak5.linux
bfile=xxx
$ldak --calc-tagging tagging/sumldak$SGE_TASK_ID --bfile $bfile --weights sections/weights.all --power -0.25  --window-kb 1000 --chr $SGE_TASK_ID --bit-size 2200
```

Then it's done!