Day1 生信新手保护学习小组笔记20250121

## 1) run the GRNBoost2 w/ pysenic N times for this file
N=`expr $SLURM_ARRAY_TASK_ID % 12 + 1`
ITER=`expr $SLURM_ARRAY_TASK_ID / 12 + 1`
EXPRMAT=$(ls $DATADIR/loom/* | sed "${N}q;d")
LABEL=$(basename $EXPRMAT .loom | sed 's/.*\.//g')


OUTFILE=$DATADIR/grnboost2/OUD_Striatum_GRNBoost2.${LABEL}.${ITER}.tsv
echo "Working on Sample ID:${LABEL} for the ${ITER}th run."
if [[ ! -f $OUTFILE || $(cat $OUTFILE) == '' ]]; then
arboreto_with_multiprocessing.py \
-o ${OUTFILE} -m grnboost2 \
--seed ${SLURM_ARRAY_TASK_ID} \
--cell_id_attribute "CellID" \
--gene_attribute "Gene" \
--num_workers 15 \
${EXPRMAT} ${TF_LIST}
fi
exit