X-Git-Url: https://git.donarmstrong.com/?p=function2gene.git;a=blobdiff_plain;f=bin%2Fcombine_results;h=21b1c9276f32e8995dccb5bc52e9560643cc5cc1;hp=76a4e0c044e47f3080f9949c1ed43cb066576ddc;hb=4ec55f9073ff175b7f591c66eed2ef91a274268c;hpb=dd8ec1d4cecf282940831171ef0f796570d781fd diff --git a/bin/combine_results b/bin/combine_results index 76a4e0c..21b1c92 100755 --- a/bin/combine_results +++ b/bin/combine_results @@ -104,6 +104,7 @@ for my $file_name (@ARGV) { $genes{$gene[NAME]}{database}{$gene[DBNAME]}++; $genes{$gene[NAME]}{hits}++; $genes{$gene[NAME]}{terms}{$gene[KEYWORD]}++; + $genes{$gene[NAME]}{terms}{$gene[KEYWORD].'['.$gene[DBNAME].']'}++; add_unique_parts($genes{$gene[NAME]},'refseq',$gene[REFSEQ]); add_if_better($genes{$gene[NAME]},'description',$gene[DESCRIPTION]); add_if_better($genes{$gene[NAME]},'location',$gene[LOCATION]); @@ -114,7 +115,7 @@ for my $file_name (@ARGV) { print join(',',map {qq("$_")} @csv_fields),qq(\n); for my $gene (keys %genes) { - $genes{$gene}{rzscore} = scalar keys %{$genes{$gene}{terms}}; + $genes{$gene}{rzscore} = scalar grep {$_ !~ /\[/} keys %{$genes{$gene}{terms}}; next if $genes{$gene}{rzscore} == 1 and exists $genes{$gene}{terms}{antigen}; $genes{$gene}{rzscore} -= 1 if exists $genes{$gene}{terms}{antigen}; print STDOUT join (',',