aild-cojo.sb

#!/bin/bash

#SBATCH --ntasks=1
#SBATCH --job-name=_cojo
#SBATCH --time=12:00:00
#SBATCH --partition=medium
#SBATCH --array=1-233%5
#SBATCH --mem=128800
#SBATCH --output=aild/cojo/_cojo_%A_%a.out
#SBATCH --error=aild/cojo/_cojo_%A_%a.err
#SBATCH --export ALL

. /etc/profile.d/modules.sh
module load default-cardio
module load slurm
module load use.own

export TMPDIR=/scratch/jhz22/tmp
export p=$(awk 'NR==ENVIRON["SLURM_ARRAY_TASK_ID"]{print $1}' work/INF1.aild)
export r=$(awk 'NR==ENVIRON["SLURM_ARRAY_TASK_ID"]{print $2}' work/INF1.aild)
if [ -f aild/cojo/${p}-${r}.jma.cojo ]; then
   rm aild/cojo/${p}-${r}.jma.cojo aild/cojo/${p}-${r}.ldr.cojo
fi
gcta-1.9 --bfile /scratch/jhz22/data/INTERVAL/INTERVAL \
         --extract aild/AILD.${r}.snplist \
         --cojo-file work/${p}.ma \
         --cojo-slct \
         --cojo-p 5e-10 \
         --cojo-wind 30000 \
         --cojo-collinear 0.9 \
         --maf 0.01 \
         --out aild/cojo/${p}-${r}