Mercurial > repos > mir-bioinf > multi_join_left

changeset 4:46c880ae6db2 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author mir-bioinf
date Wed, 15 Apr 2015 17:51:39 -0400
parents 0aa0ebcd307c
children 1de2a8f041b3
files multi_join_left/multi_join_serial.xml multi_join_left/run-multi_join_serial.pl multi_join_left/test-data/multi_join_serial_in1.tab multi_join_left/test-data/multi_join_serial_in2.tab multi_join_left/test-data/multi_join_serial_in3.tab multi_join_serial.xml run-multi_join_serial.pl test-data/multi_join_serial_in1.tab test-data/multi_join_serial_in2.tab test-data/multi_join_serial_in3.tab
diffstat 10 files changed, 278 insertions(+), 278 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/multi_join_left/multi_join_serial.xml	Wed Apr 15 17:51:39 2015 -0400
@@ -0,0 +1,84 @@
+<tool id="Multi_Join_serial" name="Join multiple" version="0.0.1" force_history_refresh="True">
+  <description>tab delimited files serially</description>
+  <!-- cms commenting out to troubleshoot -->
+  <command interpreter="perl">
+	#for $j, $s in enumerate( $Files )
+		#silent	$j	
+	#end for
+
+	#for $i, $s in enumerate( $Files )
+	   run-multi_join_serial.pl --join_file $s.joinMe --join_col $s.joinCol --iteration $i --totalfiles $j --with_header $headerYes --resultsfile $Joined_all
+	  ##print "loop iteration $i.\n";
+          ;
+        #end for
+  </command> 
+  <inputs>
+	<repeat name="Files" title="Join file">
+		<param name="joinMe" type="data" checked="yes" format="tabular" label="Join" />
+		<param name="joinCol" label="using column" type="data_column" data_ref="joinMe" />
+	</repeat>
+        <param name="headerYes" type="select" label="Treat first line as header?" help="If header starts with #, it will NOT be read, so this field should be set to no. Otherwise it can be set to yes if first line is header for ALL FILES.">
+                <option value="yes" selected="true">Yes</option>
+                <option value="no">No</option>
+        </param>
+  </inputs>
+  <outputs>
+    	<data format="tabular" name="Joined_all" label="Multi-Join result"/>
+  </outputs>
+  <tests>
+     <test>
+	<param name="Files_0|joinMe" value="multi_join_serial_in1.tab" ftype="tabular"/>
+	<param name="Files_0|joinCol" value="1"/>
+	<param name="Files_1|joinMe" value="multi_join_serial_in2.tab" ftype="tabular"/>
+	<param name="Files_1|joinCol" value="1"/>
+	<param name="Files_2|joinMe" value="multi_join_serial_in3.tab" ftype="tabular"/>
+	<param name="Files_2joinCol" value="2"/>
+	<param name="headerYes" value="yes"/>
+	<output name="Joined_all" value="multi_join_serial_out.tab" ftype="tabular"/>
+     <test/>
+  <tests/>
+  <help>
+
+This tool performs a left-outer join on multiple (at least two) files using a perl script that Ron wrote (thanks, Ron!). The resulting joined file will have the same number of rows as the first file chosen and subsequent files' matches will be shown if present. Rows in the first file without matches in the other files will have empty cells. If none of the input files have a header present, a simple column number header will be added to the output file to denote the start of each set of matches (from each file, start denoted by "C1").
+
+
+.. class:: warningmark
+
+This tool may fail due to the system running out of memory depending on the number and size of input files and number of matching lines. The higher all of these are, the more likely the tool is to fail. A red output dataset saying "Job killed" typically means the system ran into an out of memory error and as a result the job was killed. This issue has yet to be addressed at the moment...
+
+
+**Steps:**
+	1. Click Add new File for each tab-delimited file you'd like to add and the column you want to join on.
+	2. After adding all files to join, select whether the headers should all be preserved (this should be Yes if all input datasets have headers).
+	3. Click Execute.
+
+-----
+
+**Example**
+
+Dataset1::
+
+  chr1 10 20 geneA 
+  chr1 50 80 geneB
+  chr5 10 40 geneL
+
+Dataset2::
+
+  geneA tumor-supressor
+  geneB Foxp2
+  geneC Gnas1
+  geneE INK4a
+
+Joining the 4th column of Dataset1 with the 1st column of Dataset2, no header, will yield::
+
+  C1   C2 C3 C4    C1    C2
+  chr1 10 20 geneA geneA tumor-suppressor
+  chr1 50 80 geneB geneB Foxp2
+  chr5 10 40 geneL
+
+
+</help>
+
+
+</tool>
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/multi_join_left/run-multi_join_serial.pl	Wed Apr 15 17:51:39 2015 -0400
@@ -0,0 +1,173 @@
+#!/usr/bin/perl 
+
+
+use Getopt::Long;
+use Pod::Usage;
+use IO::File;
+use Data::Dumper;
+
+#require '/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/multi_join_shell.pl';  ##comment this line out when finished testing
+#require '/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/addColumnsFromFile2ToFile1.pl';
+#require '/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/GetOptWC.pm';
+
+GetOptions(
+    "join_file=s"             => \$data_in,
+    "join_col=s"	     => \$coljoin,
+    "time"                   => \$mTime,
+    "q|quiet"                => \$quiet,
+    "iteration=i"	     => \$I,
+    "totalfiles=i"	     => \$N,
+    "with_header=s"	     => \$header_yes,
+    "input_name=s"	     => \$in_name,
+    "resultsfile=s"	     => \$out_file,
+#    "h|help"                 => \$help
+) or pod2usage( -exitval => 2, -verbose => 2 );
+
+
+#check parameters and options
+my $debug = scalar(@ARGV);
+
+$coljoin--;
+#pod2usage(-msg => "To troubleshoot. ARGV should be @ARGV with $debug arguments in it.");
+pod2usage(-msg => "Forward probability should be in [0, 1]!", -exitval => 2, -verbose => 2) if ($probF < 0 || $probF > 1);
+
+$N++;
+
+# #
+use IO::Handle;
+STDOUT->fdopen( \*OUTPUT, 'a' ) or die "cant open file $!\n";   #cms changing mode from 'w' to 'a' for multiple files in one run
+STDERR->fdopen( \*ERROR,  'a' ) or die "cant open file $!\n";   #cms changing mode from 'w' to 'a' for multiple files in one run
+# # #
+
+my @options;
+
+my $fileno = $I + 1; 
+
+##Keeping track of the input files (one per iteration of this script) in an external file:
+open $Filenames, '>>', "temp_filenames.txt" or die "cannot open the temporary file $!\n";
+print $Filenames "$data_in\t";
+print $Filenames "$coljoin\n";
+
+if (($I==$N-1)&&($N>=2)) {
+        ## At the end of the last iteration
+	close($Filenames);
+        
+
+	##Read in file temp_filenames.txt
+	open(my $tmpfile, "<", "temp_filenames.txt") or die "Cannot open temp file: $!";
+	my @fileArray = <$tmpfile>;
+	#unshift @fileArray,$conditions; ##don't need to do this since conditions aren't used here
+	close($tmpfile) or die "what is that??!!! $!";
+
+	
+	##Need to send output file name to shell script:
+	push @fileArray, $out_file;  ##adds out_file to the end of fileArray
+	##Also need to send yes/no for keeping header:
+	push @fileArray, $header_yes;
+
+
+	##@fileArray has one file per line,output,header_yes, so $N+1 rows
+	my $f=0;
+	my @first;
+	my @second;
+	do {
+
+		@first = split('\t',$fileArray[$f]);  ##was filename\tJoinCol
+
+		##CMS DEALING WITH HEADER OR NOT:
+                if ($header_yes eq "no") {
+			my $fh1;
+                	$fh1 = IO::File->new("<$first[0]");
+			my $line1file1 = $fh1->getline();
+			$line1file1 =~ s/\s+$//;
+			@cols = split "\t",$line1file1;
+			my $numcols1 = @cols;
+			my $head1;
+			for (my $i=1; $i<$numcols1; $i++) {
+				$head1.="C$i\t";
+			}
+			$head1.="C$numcols1\n";
+			open(my $fh_sub, '>', './header1.txt') or die "OOPIES: $!\n";
+			print $fh_sub $head1;
+			close $fh_sub;
+			system("cat $first[0] >> ./header1.txt");  ##put header in front of file
+			##now want to use ./header1.txt instead of what was in $first[0] earlier
+			$first[0] = "./header1.txt";
+		}
+		
+
+
+		@second = split('\t',$fileArray[$f+1]);
+
+		if ($header_yes eq "no") {
+                        my $fh2;
+                        $fh2 = IO::File->new("<$second[0]");
+                        my $line1file2 = $fh2->getline();
+                        $line1file2 =~ s/\s+$//;
+                        @cols = split "\t",$line1file2;
+                        my $numcols2 = @cols;
+                        my $head2;
+                        for (my $i=1; $i<$numcols2; $i++) {
+                                $head2.="C$i\t";
+                        }
+                        $head2.="C$numcols2\n";
+                        open(my $fh_sub, '>', './header2.txt') or die "ERROR: $!\n";
+			print $fh_sub $head2;
+                        close $fh_sub;
+			system("cat $second[0] >> ./header2.txt");
+			$second[0]="./header2.txt";
+                }
+
+		system("/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/addColumnsFromFile2ToFile1.pl", "-File1=$first[0]", "-File2=$second[0]", "-cola1=$first[1]", "-cola2=$second[1]", "-colb1=$first[1]", "-colb2=$second[1]");
+		$f+=2;
+		system("mv file1_file2.txt joined.txt");
+		if ($header_yes eq "no") {
+			system("rm ./header2.txt");
+			system("rm ./header1.txt");
+		}
+	} while ($f < 2);  ##FIRST TWO ONLY!!!
+	
+	for ($f; $f<$N; $f++) {
+                my @current = split('\t',$fileArray[$f]);  ##was filename\tJoinCol
+
+		if ($header_yes eq "no") {
+                        my $fh;
+                        $fh = IO::File->new("<$current[0]");
+                        my $line1file = $fh->getline();
+                        $line1file =~ s/\s+$//;
+                        @cols = split "\t",$line1file;
+                        my $numcols = @cols;
+                        my $head;
+                        for (my $i=1; $i<$numcols; $i++) {
+                                $head.="C$i\t";
+                        }
+                        $head.="C$numcols\n";
+                        open(my $fh_sub, '>', './header.txt') or die "OOPIES: $!\n";
+                        print $fh_sub $head;
+                        close $fh_sub;
+                        system("cat $current[0] >> ./header.txt");
+                        $current[0]="./header.txt";
+                }
+
+                system("/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/addColumnsFromFile2ToFile1.pl","-File1=joined.txt", "-File2=$current[0]", "-cola1=$first[1]", "-cola2=$current[1]", "-colb1=$first[1]", "-colb2=$current[1]");
+		system("mv file1_file2.txt joined.txt");
+		if ($header_yes eq "no") {
+			system("rm ./header.txt");
+		}
+        }
+	
+	system("mv joined.txt $fileArray[-2]");
+	
+	##NOT SURE WHAT TO DO WITH THIS FOR THE MULTI-JOIN TOOL:
+	##Now, make the EC files from the genes-results files (extract appropriate columns):
+	#my $condStr =  multi_join_shell(@fileArray);  ##RSEMgetTPMs needs to take care of carriage returns
+	##NEED TO MODIFY RSEMTOEBSEQ_SHELL SO IT TAKES THE OUTPUT FILENAME AS WELL	
+
+	system("rm temp_filenames.txt");
+
+}
+elsif ($N<2) {
+}
+
+
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/multi_join_left/test-data/multi_join_serial_in1.tab	Wed Apr 15 17:51:39 2015 -0400
@@ -0,0 +1,7 @@
+symbol	sample1	sample2	sample3
+CSNK2A1P	12345	1.2345	8.90
+APOM	0.0	0.0	0.0
+HIST1H2AJ	89.5	75.6	32.3
+ASPHD1	0.001	6.98	2.33
+STBD1	0.2545	600.7	0.05
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/multi_join_left/test-data/multi_join_serial_in2.tab	Wed Apr 15 17:51:39 2015 -0400
@@ -0,0 +1,7 @@
+Gene	Sample_4	Sample_5	Sample_6	Sample_7
+NOP10	11.1	22.2	33.3	0.0
+COX7B	0.05	0.0	6.75	0.0
+ASPHD1	10	0.075	3.57	1.14
+TSTD1	1.14	0.788	99.06	94.57
+APOM	123	123	123	123
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/multi_join_left/test-data/multi_join_serial_in3.tab	Wed Apr 15 17:51:39 2015 -0400
@@ -0,0 +1,7 @@
+Experiment	GENE	SampleX	SampleY	SampleZ
+1	HIST1H2AJ	0	0	0
+1	COX7B	1.1100	59.142	16602	
+1	ZYX	1.04	1.04	1.04
+1	APOM	132	85	97
+2	PFDN6	0.06	1.87	99.96
+2	NPIPL3	1234	80.14	0.075
--- a/multi_join_serial.xml	Wed Apr 15 16:31:04 2015 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,84 +0,0 @@
-<tool id="Multi_Join_serial" name="Join multiple" version="0.0.1" force_history_refresh="True">
-  <description>tab delimited files serially</description>
-  <!-- cms commenting out to troubleshoot -->
-  <command interpreter="perl">
-	#for $j, $s in enumerate( $Files )
-		#silent	$j	
-	#end for
-
-	#for $i, $s in enumerate( $Files )
-	   /opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/run-multi_join_serial.pl --join_file $s.joinMe --join_col $s.joinCol --iteration $i --totalfiles $j --with_header $headerYes --resultsfile $Joined_all
-	  ##print "loop iteration $i.\n";
-          ;
-        #end for
-  </command> 
-  <inputs>
-	<repeat name="Files" title="Join file">
-		<param name="joinMe" type="data" checked="yes" format="tabular" label="Join" />
-		<param name="joinCol" label="using column" type="data_column" data_ref="joinMe" />
-	</repeat>
-        <param name="headerYes" type="select" label="Treat first line as header?" help="If header starts with #, it will NOT be read, so this field should be set to no. Otherwise it can be set to yes if first line is header for ALL FILES.">
-                <option value="yes" selected="true">Yes</option>
-                <option value="no">No</option>
-        </param>
-  </inputs>
-  <outputs>
-    	<data format="tabular" name="Joined_all" label="Multi-Join result"/>
-  </outputs>
-  <tests>
-     <test>
-	<param name="Files_0|joinMe" value="multi_join_serial_in1.tab" ftype="tabular"/>
-	<param name="Files_0|joinCol" value="1"/>
-	<param name="Files_1|joinMe" value="multi_join_serial_in2.tab" ftype="tabular"/>
-	<param name="Files_1|joinCol" value="1"/>
-	<param name="Files_2|joinMe" value="multi_join_serial_in3.tab" ftype="tabular"/>
-	<param name="Files_2joinCol" value="2"/>
-	<param name="headerYes" value="yes"/>
-	<output name="Joined_all" value="multi_join_serial_out.tab" ftype="tabular"/>
-     <test/>
-  <tests/>
-  <help>
-
-This tool performs a left-outer join on multiple (at least two) files using a perl script that Ron wrote (thanks, Ron!). The resulting joined file will have the same number of rows as the first file chosen and subsequent files' matches will be shown if present. Rows in the first file without matches in the other files will have empty cells. If none of the input files have a header present, a simple column number header will be added to the output file to denote the start of each set of matches (from each file, start denoted by "C1").
-
-
-.. class:: warningmark
-
-This tool may fail due to the system running out of memory depending on the number and size of input files and number of matching lines. The higher all of these are, the more likely the tool is to fail. A red output dataset saying "Job killed" typically means the system ran into an out of memory error and as a result the job was killed. This issue has yet to be addressed at the moment...
-
-
-**Steps:**
-	1. Click Add new File for each tab-delimited file you'd like to add and the column you want to join on.
-	2. After adding all files to join, select whether the headers should all be preserved (this should be Yes if all input datasets have headers).
-	3. Click Execute.
-
------
-
-**Example**
-
-Dataset1::
-
-  chr1 10 20 geneA 
-  chr1 50 80 geneB
-  chr5 10 40 geneL
-
-Dataset2::
-
-  geneA tumor-supressor
-  geneB Foxp2
-  geneC Gnas1
-  geneE INK4a
-
-Joining the 4th column of Dataset1 with the 1st column of Dataset2, no header, will yield::
-
-  C1   C2 C3 C4    C1    C2
-  chr1 10 20 geneA geneA tumor-suppressor
-  chr1 50 80 geneB geneB Foxp2
-  chr5 10 40 geneL
-
-
-</help>
-
-
-</tool>
-
--- a/run-multi_join_serial.pl	Wed Apr 15 16:31:04 2015 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,173 +0,0 @@
-#!/usr/bin/perl 
-
-
-use Getopt::Long;
-use Pod::Usage;
-use IO::File;
-use Data::Dumper;
-
-#require '/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/multi_join_shell.pl';  ##comment this line out when finished testing
-#require '/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/addColumnsFromFile2ToFile1.pl';
-#require '/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/GetOptWC.pm';
-
-GetOptions(
-    "join_file=s"             => \$data_in,
-    "join_col=s"	     => \$coljoin,
-    "time"                   => \$mTime,
-    "q|quiet"                => \$quiet,
-    "iteration=i"	     => \$I,
-    "totalfiles=i"	     => \$N,
-    "with_header=s"	     => \$header_yes,
-    "input_name=s"	     => \$in_name,
-    "resultsfile=s"	     => \$out_file,
-#    "h|help"                 => \$help
-) or pod2usage( -exitval => 2, -verbose => 2 );
-
-
-#check parameters and options
-my $debug = scalar(@ARGV);
-
-$coljoin--;
-#pod2usage(-msg => "To troubleshoot. ARGV should be @ARGV with $debug arguments in it.");
-pod2usage(-msg => "Forward probability should be in [0, 1]!", -exitval => 2, -verbose => 2) if ($probF < 0 || $probF > 1);
-
-$N++;
-
-# #
-use IO::Handle;
-STDOUT->fdopen( \*OUTPUT, 'a' ) or die "cant open file $!\n";   #cms changing mode from 'w' to 'a' for multiple files in one run
-STDERR->fdopen( \*ERROR,  'a' ) or die "cant open file $!\n";   #cms changing mode from 'w' to 'a' for multiple files in one run
-# # #
-
-my @options;
-
-my $fileno = $I + 1; 
-
-##Keeping track of the input files (one per iteration of this script) in an external file:
-open $Filenames, '>>', "temp_filenames.txt" or die "cannot open the temporary file $!\n";
-print $Filenames "$data_in\t";
-print $Filenames "$coljoin\n";
-
-if (($I==$N-1)&&($N>=2)) {
-        ## At the end of the last iteration
-	close($Filenames);
-        
-
-	##Read in file temp_filenames.txt
-	open(my $tmpfile, "<", "temp_filenames.txt") or die "Cannot open temp file: $!";
-	my @fileArray = <$tmpfile>;
-	#unshift @fileArray,$conditions; ##don't need to do this since conditions aren't used here
-	close($tmpfile) or die "what is that??!!! $!";
-
-	
-	##Need to send output file name to shell script:
-	push @fileArray, $out_file;  ##adds out_file to the end of fileArray
-	##Also need to send yes/no for keeping header:
-	push @fileArray, $header_yes;
-
-
-	##@fileArray has one file per line,output,header_yes, so $N+1 rows
-	my $f=0;
-	my @first;
-	my @second;
-	do {
-
-		@first = split('\t',$fileArray[$f]);  ##was filename\tJoinCol
-
-		##CMS DEALING WITH HEADER OR NOT:
-                if ($header_yes eq "no") {
-			my $fh1;
-                	$fh1 = IO::File->new("<$first[0]");
-			my $line1file1 = $fh1->getline();
-			$line1file1 =~ s/\s+$//;
-			@cols = split "\t",$line1file1;
-			my $numcols1 = @cols;
-			my $head1;
-			for (my $i=1; $i<$numcols1; $i++) {
-				$head1.="C$i\t";
-			}
-			$head1.="C$numcols1\n";
-			open(my $fh_sub, '>', './header1.txt') or die "OOPIES: $!\n";
-			print $fh_sub $head1;
-			close $fh_sub;
-			system("cat $first[0] >> ./header1.txt");  ##put header in front of file
-			##now want to use ./header1.txt instead of what was in $first[0] earlier
-			$first[0] = "./header1.txt";
-		}
-		
-
-
-		@second = split('\t',$fileArray[$f+1]);
-
-		if ($header_yes eq "no") {
-                        my $fh2;
-                        $fh2 = IO::File->new("<$second[0]");
-                        my $line1file2 = $fh2->getline();
-                        $line1file2 =~ s/\s+$//;
-                        @cols = split "\t",$line1file2;
-                        my $numcols2 = @cols;
-                        my $head2;
-                        for (my $i=1; $i<$numcols2; $i++) {
-                                $head2.="C$i\t";
-                        }
-                        $head2.="C$numcols2\n";
-                        open(my $fh_sub, '>', './header2.txt') or die "ERROR: $!\n";
-			print $fh_sub $head2;
-                        close $fh_sub;
-			system("cat $second[0] >> ./header2.txt");
-			$second[0]="./header2.txt";
-                }
-
-		system("/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/addColumnsFromFile2ToFile1.pl", "-File1=$first[0]", "-File2=$second[0]", "-cola1=$first[1]", "-cola2=$second[1]", "-colb1=$first[1]", "-colb2=$second[1]");
-		$f+=2;
-		system("mv file1_file2.txt joined.txt");
-		if ($header_yes eq "no") {
-			system("rm ./header2.txt");
-			system("rm ./header1.txt");
-		}
-	} while ($f < 2);  ##FIRST TWO ONLY!!!
-	
-	for ($f; $f<$N; $f++) {
-                my @current = split('\t',$fileArray[$f]);  ##was filename\tJoinCol
-
-		if ($header_yes eq "no") {
-                        my $fh;
-                        $fh = IO::File->new("<$current[0]");
-                        my $line1file = $fh->getline();
-                        $line1file =~ s/\s+$//;
-                        @cols = split "\t",$line1file;
-                        my $numcols = @cols;
-                        my $head;
-                        for (my $i=1; $i<$numcols; $i++) {
-                                $head.="C$i\t";
-                        }
-                        $head.="C$numcols\n";
-                        open(my $fh_sub, '>', './header.txt') or die "OOPIES: $!\n";
-                        print $fh_sub $head;
-                        close $fh_sub;
-                        system("cat $current[0] >> ./header.txt");
-                        $current[0]="./header.txt";
-                }
-
-                system("/opt/galaxy/galaxy-dist/tools/ngs_rna/Unreleased/addColumnsFromFile2ToFile1.pl","-File1=joined.txt", "-File2=$current[0]", "-cola1=$first[1]", "-cola2=$current[1]", "-colb1=$first[1]", "-colb2=$current[1]");
-		system("mv file1_file2.txt joined.txt");
-		if ($header_yes eq "no") {
-			system("rm ./header.txt");
-		}
-        }
-	
-	system("mv joined.txt $fileArray[-2]");
-	
-	##NOT SURE WHAT TO DO WITH THIS FOR THE MULTI-JOIN TOOL:
-	##Now, make the EC files from the genes-results files (extract appropriate columns):
-	#my $condStr =  multi_join_shell(@fileArray);  ##RSEMgetTPMs needs to take care of carriage returns
-	##NEED TO MODIFY RSEMTOEBSEQ_SHELL SO IT TAKES THE OUTPUT FILENAME AS WELL	
-
-	system("rm temp_filenames.txt");
-
-}
-elsif ($N<2) {
-}
-
-
-
--- a/test-data/multi_join_serial_in1.tab	Wed Apr 15 16:31:04 2015 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,7 +0,0 @@
-symbol	sample1	sample2	sample3
-CSNK2A1P	12345	1.2345	8.90
-APOM	0.0	0.0	0.0
-HIST1H2AJ	89.5	75.6	32.3
-ASPHD1	0.001	6.98	2.33
-STBD1	0.2545	600.7	0.05
-
--- a/test-data/multi_join_serial_in2.tab	Wed Apr 15 16:31:04 2015 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,7 +0,0 @@
-Gene	Sample_4	Sample_5	Sample_6	Sample_7
-NOP10	11.1	22.2	33.3	0.0
-COX7B	0.05	0.0	6.75	0.0
-ASPHD1	10	0.075	3.57	1.14
-TSTD1	1.14	0.788	99.06	94.57
-APOM	123	123	123	123
-
--- a/test-data/multi_join_serial_in3.tab	Wed Apr 15 16:31:04 2015 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,7 +0,0 @@
-Experiment	GENE	SampleX	SampleY	SampleZ
-1	HIST1H2AJ	0	0	0
-1	COX7B	1.1100	59.142	16602	
-1	ZYX	1.04	1.04	1.04
-1	APOM	132	85	97
-2	PFDN6	0.06	1.87	99.96
-2	NPIPL3	1234	80.14	0.075