Mercurial > repos > xuebing > splicesite_max_entropy
diff splicesitescore/score3.pl~ @ 8:ec97befe722b
Uploaded
author | xuebing |
---|---|
date | Sun, 01 Apr 2012 08:14:45 -0400 |
parents | 3fd430270c0c |
children |
line wrap: on
line diff
--- a/splicesitescore/score3.pl~ Sun Apr 01 01:59:22 2012 -0400 +++ b/splicesitescore/score3.pl~ Sun Apr 01 08:14:45 2012 -0400 @@ -4,8 +4,8 @@ my $inputfile = $ARGV[0]; my $usemaxent = 1; - -my @metables = &makemaxentscores; +my $modelpath = "/home/wuxbl/rowley/tools/maxentscan/splicemodels/"; +my @metables = &makemaxentscores($modelpath); open (FILE,"<$inputfile") || die "can't open!\n"; @@ -51,7 +51,7 @@ } sub makemaxentscores{ - my $dir = "splicemodels/"; + my $dir = shift; my @list = ('me2x3acc1','me2x3acc2','me2x3acc3','me2x3acc4', 'me2x3acc5','me2x3acc6','me2x3acc7','me2x3acc8','me2x3acc9'); my @metables; @@ -71,48 +71,7 @@ } return @metables; } -sub makewmmscores{ - my $dir = "splicemodels/"; - my @list = ('me1s0acc1','me1s0acc2','me1s0acc3','me1s0acc4', - 'me1s0acc5','me1s0acc6','me1s0acc7','me1s0acc8','me1s0acc9'); - my @metables; - my $num = 0 ; - foreach my $file (@list) { - my $n = 0; - open (SCOREF,"<".$dir.$file) || die "Can't open $file!\n"; - while(<SCOREF>) { - chomp; - $_=~ s/\s//; - $metables[$num]{$n} = $_; - $n++; - } - close(SCOREF); - #print STDERR $file."\t".$num."\t".$n."\n"; - $num++; - } - return @metables; -} -sub makemmscores{ - my $dir = "splicemodels/"; - my @list = ('me2s0acc1','me2s0acc2','me2s0acc3','me2s0acc4', - 'me2s0acc5','me2s0acc6','me2s0acc7','me2s0acc8','me2s0acc9'); - my @metables; - my $num = 0 ; - foreach my $file (@list) { - my $n = 0; - open (SCOREF,"<".$dir.$file) || die "Can't open $file!\n"; - while(<SCOREF>) { - chomp; - $_=~ s/\s//; - $metables[$num]{$n} = $_; - $n++; - } - close(SCOREF); - #print STDERR $file."\t".$num."\t".$n."\n"; - $num++; - } - return @metables; -} + sub maxentscore{ my $seq = shift; my $table_ref = shift;