# HG changeset patch
# User xuebing
# Date 1333282485 14400
# Node ID ec97befe722b32092f4b7e94492b9e13435f078b
# Parent 3fd430270c0c855f4ec49016c4afd35007f29530
Uploaded
diff -r 3fd430270c0c -r ec97befe722b splicesite.xml
--- a/splicesite.xml Sun Apr 01 01:59:22 2012 -0400
+++ b/splicesite.xml Sun Apr 01 08:14:45 2012 -0400
@@ -7,6 +7,13 @@
+
+
+
+
+
+
+
diff -r 3fd430270c0c -r ec97befe722b splicesitescore/score3.pl
--- 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() {
- 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() {
- 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;
diff -r 3fd430270c0c -r ec97befe722b splicesitescore/score3.pl~
--- 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() {
- 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() {
- 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;
diff -r 3fd430270c0c -r ec97befe722b splicesitescore/score5.pl
--- a/splicesitescore/score5.pl Sun Apr 01 01:59:22 2012 -0400
+++ b/splicesitescore/score5.pl Sun Apr 01 08:14:45 2012 -0400
@@ -4,7 +4,7 @@
my $inputfile = $ARGV[0];
my $usemaxent = 1;
-my $modelpath = "splicesitescore/";
+my $modelpath = "/home/wuxbl/rowley/tools/maxentscan/";
my %me2x5 = &makescorematrix($modelpath.'me2x5');
my %seq = &makesequencematrix($modelpath.'splicemodels/splice5sequences');
diff -r 3fd430270c0c -r ec97befe722b splicesitescore/score5.pl~
--- a/splicesitescore/score5.pl~ Sun Apr 01 01:59:22 2012 -0400
+++ b/splicesitescore/score5.pl~ Sun Apr 01 08:14:45 2012 -0400
@@ -4,7 +4,7 @@
my $inputfile = $ARGV[0];
my $usemaxent = 1;
-my $modelpath = "/Users/xuebing/galaxy-dist/tools/mytools/splicesitescore/";
+my $modelpath = "/home/wuxbl/rowley/tools/maxentscan/splicemodels/";
my %me2x5 = &makescorematrix($modelpath.'me2x5');
my %seq = &makesequencematrix($modelpath.'splicemodels/splice5sequences');