annotate PostAnalysis @ 5:2ebca9da5e42 draft default tip

planemo upload
author bioitcore
date Thu, 07 Sep 2017 17:39:24 -0400
parents adc0f7765d85
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
1 #!/usr/bin/perl
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
2 my $SrcFolder="/home/galaxy/galaxy-dist/tools/SpliceTrap.0.90.1/bin";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
3 # this script is a wrapup for Post analysis based on the ratio file output
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
4
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
5 use strict;
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
6 use Getopt::Long;
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
7 my $RatioFile="";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
8 my $OutputFile = "";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
9 my $JunctionCut=5;
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
10 my $CutoffLevel="M";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
11 my $noIRM = 0;
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
12 my $noIRMstr="";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
13
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
14 GetOptions (
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
15 "i:s"=>\$RatioFile,
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
16 "o:s"=>\$OutputFile,
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
17 "c:s"=>\$CutoffLevel,
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
18 "noIRM|noirm"=>\$noIRM,
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
19 "j:i"=>\$JunctionCut
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
20 );
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
21
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
22 my $InputParaDes=" Usage of the script:
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
23 -i input file (.ratio file)
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
24 -o output file
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
25 -c Cutoff Level:H/[M]/L
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
26 Means High, Middle or Low
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
27 -j Junction reads per junction requirement for each exon-isoform [5]
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
28 --noIRM Use the unadjusted inclusion ratios (before IRM correction)
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
29 ";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
30
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
31 if($RatioFile eq "")
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
32 {
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
33 print $InputParaDes;
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
34 exit;
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
35 }
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
36
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
37 if($CutoffLevel ne "H" and $CutoffLevel ne "M" and $CutoffLevel ne "L")
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
38 {
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
39 print $InputParaDes;
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
40 exit;
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
41 }
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
42 if($noIRM)
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
43 {
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
44 $noIRMstr= "noirm";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
45 }
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
46
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
47
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
48 system("perl $SrcFolder/ApplyCutoff.jie.pl $RatioFile $CutoffLevel $JunctionCut $noIRMstr >$OutputFile.raw");
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
49
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
50 open(rawfile, "$OutputFile.raw");
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
51 open(outfile, ">$OutputFile");
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
52 while(my $line=<rawfile>)
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
53 {
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
54 chomp($line);
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
55 my @a=split("\t",$line);
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
56 if($noIRM)
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
57 {
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
58 print outfile join("\t",$a[21],$a[1],$a[3],$a[4],$a[5],$a[6],$a[7],$a[11],$a[12],$a[13],$a[14]),"\n";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
59 }
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
60 else
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
61 {
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
62 print outfile join("\t",$a[21],$a[2],$a[3],$a[4],$a[5],$a[6],$a[7],$a[11],$a[12],$a[13],$a[14]),"\n";
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
63 }
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
64 }
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
65 close(outfile);
adc0f7765d85 planemo upload
bioitcore
parents:
diff changeset
66 close(rawfile);