diff splicesitescore/score3.pl @ 8:ec97befe722b

Uploaded
author xuebing
date Sun, 01 Apr 2012 08:14:45 -0400
parents 85af77643e89
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;