-
Notifications
You must be signed in to change notification settings - Fork 0
/
onePAGE.sh
35 lines (27 loc) · 655 Bytes
/
onePAGE.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
inDir=$1
exp=$2
geneset=$3
if [ -z ${TEISERDIR+x} ]; then
echo "Error: TEISERDIR is not set. Please set TEISERDIR to the TEISER directory.";
return;
fi
cd $inDir
outDir=${exp/.txt/};
outDir=${outDir}_onePAGE_${geneset}
echo $exp $geneset
base=`basename $exp`
base=${base/.txt/}
perl ${TEISERDIR}/run_mi_gene_list.pl \
--expfile=$exp \
--genefile=${geneset}.txt \
--exptype=continuous \
--ebins=11 \
--draw_min=-3 \
--draw_max=3 \
--species=human \
--doremovedups=0 \
--doremoveextra=0
rm -r $outDir
mv ${exp}_GENESET $outDir
mv ${outDir}/${base}.txt.summary.pdf ${outDir}.summary.pdf
echo 'done!'