Mercurial > repos > cristian > rbgoa
diff gomwu_a.pl @ 1:f7287f82602f draft
"planemo upload commit 486235d6560c9e95bd42152ad19bf7c3941cdc1b"
author | cristian |
---|---|
date | Tue, 19 Apr 2022 08:28:43 +0000 |
parents | 91261b42c07e |
children | 5acf9dfdfa27 |
line wrap: on
line diff
--- 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";