diff gomwu_a.pl @ 1:f7287f82602f draft

"planemo upload commit 486235d6560c9e95bd42152ad19bf7c3941cdc1b"

--- a/gomwu_a.pl	Thu Apr 14 13:28:05 2022 +0000
+++ b/gomwu_a.pl	Tue Apr 19 08:28:43 2022 +0000
@@ -13,6 +13,8 @@
 
 ";
 
+use File::Basename;
+
 print "@ARGV";
 
 my $onto=$ARGV[0] or die $usage;
@@ -47,10 +49,13 @@
 elsif ($div eq "CC") { $division="cellular_component";}
 else { die "unrecognized division: $div\n";}
 
-my $inname2=$measure.".".$div.".tmp";
-my $inname3=$div."_".$measure;
-my $inname31="dissim0_".$div."_".$gen2go;
-my $inname4="dissim_".$div."_".$measure."_".$gen2go;
+($mname,$mdir,$mext) = fileparse($measure,'\..*');
+($aname,$adir,$aext) = fileparse($gen2go,'\..*');
+print "$mname - $mdir - $mext\n";
+my $inname2=$mdir.$mname.".".$div.".tmp";
+my $inname3=$mdir.$mname."_".$div.".tsv";
+my $inname31=$mdir."dissim0_".$div."_".$aname.$aext;
+my $inname4=$mdir."dissim_".$div."_".$mname."_".$aname.$aext;
 
 my @donealready=();
 
@@ -131,7 +136,7 @@
 		next;
 	}
 	chomp;
-	($seq,$ns)=split(/,/, $_);
+	($seq,$ns)=split(/\t/, $_);
 	if ($seq=~/SEQ/) { $seq.="_s";}
 	$dnds{$seq}=$ns;
 }
@@ -217,7 +222,7 @@
 my %value={};
 my $des;
 my $go; 
-my $l;
+my $ll;
 my $gn;
 my $val;
 my @gos=();
@@ -228,7 +233,7 @@
 
 while (<TAB>){
 	chomp;
-	($des,$go,$l,$val,$gn)=split(/\t/,$_);
+	($des,$go,$ll,$val,$gn)=split(/\t/,$_);
 	$value{$gn}=$val;
 	$desc{$go}=$des;
 	push @{$genes{$go}},$gn;
@@ -316,7 +321,7 @@
 		}
 	}
 
-	open OUT, ">$inname31" or die "gomwu_b: cannot create output $inname31\n";
+	open OUT, ">$inname31" or die "gomwu_a: cannot create output $inname31\n";
 	
 	print {OUT} join("\t",@gos),"\n";