all-igi.sh 1.26 KB
Newer Older
Philip Lijnzaad's avatar
Philip Lijnzaad committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
#!/bin/sh -x
# -*- mode: sh; -*-
# $Id$
# 
# run the complete IGI show. May not be all that useful (since often I have
# to repair things by hand ... in fact, a Makefile would be more convenient),
# but at least can serve as some kind of documentation of what's supposed to
# happen. See also README in this directory, of course.

# prepare the sources into a form that can be merged; results go to
# subdir out/collated
all-collates.sh > all-collates.out 2>all-collates.log

# do the merge; results go to subdir out/merged
all-merges.sh > all-merges.out 2>all-merges.log

# run statistics on the merges, and also produce the summary and final gtf
# file(s). Results go to  out/{stats,mapping,summary,final}
all-stats.sh > all-stats.out 2>all-stats.log

Philip Lijnzaad's avatar
Philip Lijnzaad committed
21
22
23
24
25
26
27
28
29
30
31
32
33
34

# do the peptide business; first, collate them:
pep-collate.sh 2> pep-collate.log

# Now pull out the longest peptide of each original igi
# peptide file, using the 'valid' igi's from the summary files (which are
# proudced by all-stats.sh.
summary=out/summary/ens_affy_fgenesh.summary
peptidefiles="ensembl/ensembl.pep affymetrix/affymetrix.pep fgenesh/fgenesh.pep"
outdir=pep
gtfsummary2pep.pl $summary $peptidefiles > $outdir/pep.out 2> $outdir/pep.log
gzip -c $outdir/pep.out > $outdir/pep.out.gz
# ftp stuff, and we're done.