search

AnantharamanLab / METABOLIC

A scalable high-throughput metabolic and biogeochemical functional trait profiler
165 stars 42 forks source link

Calculation error Sequential_Transformation_input_2 #135

Open bioone opened 1 year ago

bioone commented 1 year ago

The code below in METABOLIC-C.pl may causing an error in the final calculation of the number of genomes involved in metabolism. my @tmp_enzymes_2 = split (/\;/,$CAZy_map{$GH_id}) will return a array， but the space will continue to have.

my %R_mh_02_summary = (); # step => gn => 1 or 0 foreach my $step (sort keys %R_mh_02){ foreach my $gn (sort keys %dbCANout){ $R_mh_02_summary{$step}{$gn} = 0; if ($step !~ /\+/){ my $enzymes = $R_mh_02{$step}; my @tmp_enzymes = split (/\;/, $enzymes); foreach my $enzyme_id (@tmp_enzymes){ foreach my $GH_id (sort keys %CAZy_map){ my @tmp_enzymes_2 = split (/\;/,$CAZy_map{$GH_id}); foreach my $enzyme_id_2 (@tmp_enzymes_2){ if ($enzyme_id_2 eq $enzyme_id and $dbCANout{$gn}{$GH_id}){ $R_mh_02_summary{$step}{$gn} = 1; } } } }

ChaoLab commented 1 year ago

Can you explain more about this issue? What do you mean by "but the space will continue to have"?