forked from tanghaibao/goatools
-
Notifications
You must be signed in to change notification settings - Fork 0
/
run.sh
executable file
·48 lines (35 loc) · 1.04 KB
/
run.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
36
37
38
39
40
41
42
43
44
45
46
47
48
#!/bin/bash
# configuration for Gene Ontology files
GO_OBO_FILE=go-basic.obo
GOSLIM_OBO_FILE=goslim_generic.obo
GO_OBO_DOWNLOAD=http://geneontology.org/ontology/go-basic.obo
GOSLIM_OBO_DOWNLOAD=http://www.geneontology.org/ontology/subsets/goslim_generic.obo
# if the gene ontology files don't exist, download them
if [ ! -f $GO_OBO_FILE ]
then
echo "downloading GO file: $GO_OBO_FILE"
wget -O $GO_OBO_FILE $GO_OBO_DOWNLOAD
fi
if [ ! -f $GOSLIM_OBO_FILE ]
then
echo "downloading GOslim file: $GOSLIM_OBO_FILE"
wget -O $GOSLIM_OBO_FILE $GOSLIM_OBO_DOWNLOAD
fi
echo "select the test from below"
select TEST in 'Test enrichment' 'Test plotting go terms' 'Test mapslim' 'Test mapslim on assocation'
do
case $REPLY in
1)
python scripts/find_enrichment.py --pval=0.05 --indent data/study data/population data/association
;;
2)
python scripts/plot_go_term.py --term=GO:0008135
;;
3)
python 'tests/test_mapslim.py'
;;
4)
python scripts/map_to_slim.py --association_file=data/association --slim_out=direct $GO_OBO_FILE $GOSLIM_OBO_FILE
;;
esac
done