#!/bin/bash

# model directory
export model='/oak/stanford/groups/akundaje/avsec/basepair/data/processed/comparison/output/nexus,peaks,OSNK,0,10,1,FALSE,same,0.5,64,25,0.004,9,FALSE,[1,50],TRUE'

export imp_file='mm10-tRNAs-1000bp_resized.imp_score.h5'
# generate mm10-tRNAs-1000bp_resized.imp_score.h5
basepair imp-score-seqmodel \
	 $model \
	 $imp_file \
	 --dataspec ${model}/dataspec.yaml  \
	 --regions mm10-tRNAs-1000bp_resized.bed3 \
	 --gpu=0 \
	 --method=deeplift \
	 --peak-width=1000 \
	 --seq-width=1000 \
	 --batch-size=16 \
	 --memfrac=1 \
	 --num-workers=5


for task in Oct4 Sox2 Nanog Klf4; do
    echo $task
    mkdir -p $task
    basepair modisco-score2 $model/deeplift/$task/out/profile/wn $task/instances.parq --imp-scores=$imp_file --ignore-filter --trim-frac=0.08 --n-jobs=20
done
