[Bio] / FigKernelPackages / Observation.pm Repository:
ViewVC logotype

Diff of /FigKernelPackages/Observation.pm

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.57, Thu Apr 17 20:45:23 2008 UTC revision 1.58, Wed Apr 23 15:50:16 2008 UTC
# Line 5  Line 5 
5  use Data::Dumper;  use Data::Dumper;
6    
7  require Exporter;  require Exporter;
8  @EXPORT_OK = qw(get_objects);  @EXPORT_OK = qw(get_objects, get_sims_objects);
9    
10  use WebColors;  use WebColors;
11  use WebConfig;  use WebConfig;
# Line 374  Line 374 
374    
375  }  }
376    
377    =head3 get_sims_objects()
378    
379    This is the B<REAL WORKHORSE> method of this Package.
380    
381    =cut
382    
383    sub get_sims_objects {
384        my ($self,$fid,$fig,$parameters) = @_;
385    
386        my $objects = [];
387        my @matched_datasets=();
388    
389        # call function that fetches attribute based observations
390        # returns an array of arrays of hashes
391        get_sims_observations($fid,\@matched_datasets,$fig,$parameters);
392    
393        foreach my $dataset (@matched_datasets) {
394            my $object;
395            if ($dataset->{'class'} eq "SIM"){
396                $object = Observation::Sims->new($dataset);
397            }
398            push (@$objects, $object);
399        }
400        return $objects;
401    }
402    
403    
404  =head3 display_housekeeping  =head3 display_housekeeping
405  This method returns the housekeeping data for a given peg in a table format  This method returns the housekeeping data for a given peg in a table format
406    
# Line 666  Line 693 
693  =cut  =cut
694    
695  sub get_sims_observations{  sub get_sims_observations{
696        my ($fid,$datasets_ref,$fig,$parameters) = (@_);
697    
698      my ($fid,$datasets_ref,$fig) = (@_);      my ($max_sims, $max_expand, $max_eval, $sim_order, $db_filter);
699      #my $fig = new FIG;      if ($parameters->{flag}){
700      my @sims= $fig->sims($fid,500,10,"fig");        $max_sims = $parameters->{max_sims};
701          $max_expand = $parameters->{max_expand};
702          $max_eval = $parameters->{max_eval};
703          $db_filter = $parameters->{db_filter};
704          $sim_order = $parameters->{sim_order};
705          $group_by_genome = 1 if (defined ($parameters->{group_genome}));
706        }
707        else{
708          $max_sims = 50;
709          $max_expand = 5;
710          $max_eval = 1e-5;
711          $db_filter = "figx";
712          $sim_order = "id";
713        }
714    
715        my(@tmp, $id, $genome, @genomes, %sims);
716        my @tmp= $fig->sims($fid,$max_sims,$max_eval,$db_filter,$max_expand);
717        @tmp = grep { !($_->id2 =~ /^fig\|/ and $fig->is_deleted_fid($_->id2)) } @tmp;
718      my ($dataset);      my ($dataset);
719    
720      foreach my $sim (@sims){      if ($group_by_genome){
721          next if ($fig->is_deleted_fid($sim->[1]));        #  Collect all sims from genome with the first occurance of the genome:
722          foreach $sim ( @tmp ){
723            $id = $sim->id2;
724            $genome = ($id =~ /^fig\|(\d+\.\d+)\.peg\.\d+/) ? $1 : $id;
725            if (! defined( $sims{ $genome } ) ) { push @genomes, $genome }
726            push @{ $sims{ $genome } }, $sim;
727          }
728          @tmp = map { @{ $sims{$_} } } @genomes;
729        }
730    
731        foreach my $sim (@tmp){
732          my $hit = $sim->[1];          my $hit = $sim->[1];
733          my $percent = $sim->[2];          my $percent = $sim->[2];
734          my $evalue = $sim->[10];          my $evalue = $sim->[10];
# Line 718  Line 773 
773      my ($id) = (@_);      my ($id) = (@_);
774    
775      my ($db);      my ($db);
776      if ($id =~ /^fig\|/)              { $db = "FIG" }      if ($id =~ /^fig\|/)              { $db = "SEED" }
777      elsif ($id =~ /^gi\|/)            { $db = "NCBI" }      elsif ($id =~ /^gi\|/)            { $db = "NCBI" }
778        elsif ($id =~ /^gb\|/)            { $db = "GenBank" }
779      elsif ($id =~ /^^[NXYZA]P_/)      { $db = "RefSeq" }      elsif ($id =~ /^^[NXYZA]P_/)      { $db = "RefSeq" }
780        elsif ($id =~ /^ref\|/)           { $db = "RefSeq" }
781      elsif ($id =~ /^sp\|/)            { $db = "SwissProt" }      elsif ($id =~ /^sp\|/)            { $db = "SwissProt" }
782      elsif ($id =~ /^uni\|/)           { $db = "UniProt" }      elsif ($id =~ /^uni\|/)           { $db = "UniProt" }
783      elsif ($id =~ /^tigr\|/)          { $db = "TIGR" }      elsif ($id =~ /^tigr\|/)          { $db = "TIGR" }
# Line 729  Line 786 
786      elsif ($id =~ /^tr\|/)                          { $db = "TrEMBL" }      elsif ($id =~ /^tr\|/)                          { $db = "TrEMBL" }
787      elsif ($id =~ /^eric\|/)          { $db = "ASAP" }      elsif ($id =~ /^eric\|/)          { $db = "ASAP" }
788      elsif ($id =~ /^img\|/)           { $db = "JGI" }      elsif ($id =~ /^img\|/)           { $db = "JGI" }
789        elsif ($id =~ /^pdb\|/)           { $db = "PDB" }
790        elsif ($id =~ /^img\|/)           { $db = "IMG" }
791        elsif ($id =~ /^cmr\|/)           { $db = "CMR" }
792        elsif ($id =~ /^dbj\|/)           { $db = "DBJ" }
793    
794      return ($db);      return ($db);
795    
# Line 1432  Line 1493 
1493      #color is      #color is
1494      my $color = "6";      my $color = "6";
1495    
1496  =pod=  =head3
1497    
1498      if($cello_location){      if($cello_location){
1499          my $cello_descriptions = [];          my $cello_descriptions = [];
# Line 1675  Line 1736 
1736  =cut  =cut
1737    
1738  sub display {  sub display {
1739      my ($self,$gd,$array,$fig) = @_;      my ($self,$gd,$thing,$fig,$base_start,$in_subs,$cgi) = @_;
     #my $fig = new FIG;  
   
     my @ids;  
     foreach my $thing(@$array){  
         next if ($thing->class ne "SIM");  
         push (@ids, $thing->acc);  
     }  
   
     my %in_subs  = $fig->subsystems_for_pegs(\@ids);  
   
     my %sims_objects_evalue;  
     my $count = 0;  
     foreach my $thing (@$array){  
         if ($thing->class eq "SIM"){  
             $sims_objects_evalue{$count} = $thing->evalue;  
         }  
         $count++;  
     }  
1740    
1741      foreach my $index (sort {$sims_objects_evalue{$a}<=>$sims_objects_evalue{$b}} keys %sims_objects_evalue){      # declare variables
1742  #    foreach my $thing ( @$array){      my $window_size = $gd->window_size;
         my $thing = $array->[$index];  
         if ($thing->class eq "SIM"){  
1743              my $peg = $thing->acc;              my $peg = $thing->acc;
1744              my $query = $thing->query;      my $query_id = $thing->query;
   
1745              my $organism = $thing->organism;              my $organism = $thing->organism;
1746        my $abbrev_name = $fig->abbrev($organism);
1747        if (!$organism){
1748          $organism = $peg;
1749          $abbrev_name = $peg;
1750        }
1751              my $genome = $fig->genome_of($peg);              my $genome = $fig->genome_of($peg);
1752              my ($org_tax) = ($genome) =~ /(.*)\./;              my ($org_tax) = ($genome) =~ /(.*)\./;
1753              my $function = $thing->function;              my $function = $thing->function;
1754              my $abbrev_name = $fig->abbrev($organism);      my $query_start = $thing->qstart;
1755              my $align_start = $thing->qstart;      my $query_stop = $thing->qstop;
             my $align_stop = $thing->qstop;  
1756              my $hit_start = $thing->hstart;              my $hit_start = $thing->hstart;
1757              my $hit_stop = $thing->hstop;              my $hit_stop = $thing->hstop;
1758        my $ln_query = $thing->qlength;
1759        my $ln_hit = $thing->hlength;
1760        my $query_color = match_color($query_start, $query_stop, $ln_query, 1);
1761        my $hit_color = match_color($hit_start, $hit_stop, $ln_hit, 1);
1762    
1763              my $tax_link = "http://www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi?id=" . $org_tax;              my $tax_link = "http://www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi?id=" . $org_tax;
1764    
1765        # hit sequence title
1766              my $line_config = { 'title' => "$organism [$org_tax]",              my $line_config = { 'title' => "$organism [$org_tax]",
1767                                  'short_title' => "$abbrev_name",                                  'short_title' => "$abbrev_name",
1768                                  'title_link' => '$tax_link',                                  'title_link' => '$tax_link',
1769                                  'basepair_offset' => '0'                                  'basepair_offset' => '0'
1770                                  };                                  };
1771    
1772        # query sequence title
1773        my $replace_id = $peg;
1774        $replace_id =~ s/\|/_/ig;
1775        my $anchor_name = "anchor_". $replace_id;
1776        my $query_config = { 'title' => "Query",
1777                             'short_title' => "Query",
1778                             'title_link' => "changeSimsLocation('$replace_id', 1)",
1779                             'basepair_offset' => '0'
1780                             };
1781              my $line_data = [];              my $line_data = [];
1782        my $query_data = [];
1783    
1784              my $element_hash;              my $element_hash;
1785              my $links_list = [];      my $hit_links_list = [];
1786              my $descriptions = [];      my $hit_descriptions = [];
1787        my $query_descriptions = [];
1788              # get subsystem information  
1789              my $url_link = "?page=Annotation&feature=".$peg;      # get sequence information
1790              my $link;      # evidence link
1791              $link = {"link_title" => $peg,      my $evidence_link;
1792                       "link" => $url_link};      if ($peg =~ /^fig\|/){
1793              push(@$links_list,$link);        $evidence_link = "?page=Evidence&feature=".$peg;
1794        }
1795        else{
1796          my $db = &Observation::get_database($peg);
1797          my ($link_id) = ($peg) =~ /\|(.*)/;
1798          $evidence_link = &HTML::alias_url($link_id, $db);
1799          #print STDERR "LINK: $db    $evidence_link";
1800        }
1801        my $link = {"link_title" => $peg,
1802                    "link" => $evidence_link};
1803        push(@$hit_links_list,$link) if ($evidence_link);
1804    
1805              #my @subsystems = $fig->peg_to_subsystems($peg);      # subsystem link
1806              my @subs = @{$in_subs{$peg}} if (defined $in_subs{$peg});      my $subs = $in_subs->{$peg} if (defined $in_subs->{$peg});
1807              my @subsystems;              my @subsystems;
1808        foreach my $array (@$subs){
             foreach my $array (@subs){  
1809                  my $subsystem = $$array[0];                  my $subsystem = $$array[0];
1810                  push(@subsystems,$subsystem);                  push(@subsystems,$subsystem);
1811                  my $link;          my $link = {"link" => "?page=Subsystems&subsystem=$subsystem",
                 $link = {"link" => "?page=Subsystems&subsystem=$subsystem",  
1812                           "link_title" => $subsystem};                           "link_title" => $subsystem};
1813                  push(@$links_list,$link);          push(@$hit_links_list,$link);
1814              }              }
1815    
1816        # blast alignment
1817              $link = {"link_title" => "view blast alignment",              $link = {"link_title" => "view blast alignment",
1818                       "link" => "$FIG_Config::cgi_url/seedviewer.cgi?page=ToolResult&tool=bl2seq&peg1=$query&peg2=$peg"};               "link" => "$FIG_Config::cgi_url/seedviewer.cgi?page=ToolResult&tool=bl2seq&peg1=$query_id&peg2=$peg"};
1819              push (@$links_list,$link);      push (@$hit_links_list,$link) if ($peg =~ /^fig\|/);
1820    
1821        # description data
1822              my $description_function;              my $description_function;
1823              $description_function = {"title" => "function",              $description_function = {"title" => "function",
1824                                       "value" => $function};                                       "value" => $function};
1825              push(@$descriptions,$description_function);      push(@$hit_descriptions,$description_function);
1826    
1827              my ($description_ss, $ss_string);      # subsystem description
1828              $ss_string = join (",", @subsystems);      my $ss_string = join (",", @subsystems);
1829              $description_ss = {"title" => "subsystems",      $ss_string =~ s/_/ /ig;
1830        my $description_ss = {"title" => "subsystems",
1831                                 "value" => $ss_string};                                 "value" => $ss_string};
1832              push(@$descriptions,$description_ss);      push(@$hit_descriptions,$description_ss);
1833    
1834        # location description
1835        # hit
1836              my $description_loc;              my $description_loc;
1837              $description_loc = {"title" => "location start",      $description_loc = {"title" => "Hit Location",
1838                                  "value" => $hit_start};                          "value" => $hit_start . " - " . $hit_stop};
1839              push(@$descriptions, $description_loc);      push(@$hit_descriptions, $description_loc);
1840    
1841              $description_loc = {"title" => "location stop",      $description_loc = {"title" => "Sequence Length",
1842                                  "value" => $hit_stop};                          "value" => $ln_hit};
1843              push(@$descriptions, $description_loc);      push(@$hit_descriptions, $description_loc);
1844    
1845        # query
1846        $description_loc = {"title" => "Hit Location",
1847                            "value" => $query_start . " - " . $query_stop};
1848        push(@$query_descriptions, $description_loc);
1849    
1850        $description_loc = {"title" => "Sequence Length",
1851                            "value" => $ln_query};
1852        push(@$query_descriptions, $description_loc);
1853    
1854    
1855    
1856        # evalue score description
1857              my $evalue = $thing->evalue;              my $evalue = $thing->evalue;
1858              while ($evalue =~ /-0/)              while ($evalue =~ /-0/)
1859              {              {
# Line 1779  Line 1863 
1863              }              }
1864    
1865              my $color = &color($evalue);              my $color = &color($evalue);
   
1866              my $description_eval = {"title" => "E-Value",              my $description_eval = {"title" => "E-Value",
1867                                      "value" => $evalue};                                      "value" => $evalue};
1868              push(@$descriptions, $description_eval);      push(@$hit_descriptions, $description_eval);
1869        push(@$query_descriptions, $description_eval);
1870    
1871              my $identity = $self->identity;              my $identity = $self->identity;
1872              my $description_identity = {"title" => "Identity",              my $description_identity = {"title" => "Identity",
1873                                          "value" => $identity};                                          "value" => $identity};
1874              push(@$descriptions, $description_identity);      push(@$hit_descriptions, $description_identity);
1875        push(@$query_descriptions, $description_identity);
1876    
1877    
1878        my $number = $base_start + ($query_start-$hit_start);
1879        #print STDERR "START: $number";
1880        $element_hash = {
1881            "title" => $query_id,
1882            "start" => $base_start,
1883            "end" => $base_start+$ln_query,
1884            "type"=> 'box',
1885            "color"=> $color,
1886            "zlayer" => "2",
1887            "links_list" => $query_links_list,
1888            "description" => $query_descriptions
1889            };
1890        push(@$query_data,$element_hash);
1891    
1892        $element_hash = {
1893            "title" => $query_id . ': HIT AREA',
1894            "start" => $base_start + $query_start,
1895            "end" =>  $base_start + $query_stop,
1896            "type"=> 'smallbox',
1897            "color"=> $query_color,
1898            "zlayer" => "3",
1899            "links_list" => $query_links_list,
1900            "description" => $query_descriptions
1901            };
1902        push(@$query_data,$element_hash);
1903    
1904        $gd->add_line($query_data, $query_config);
1905    
1906    
1907              $element_hash = {              $element_hash = {
1908                  "title" => $peg,                  "title" => $peg,
1909                  "start" => $align_start,                  "start" => $base_start + ($query_start-$hit_start),
1910                  "end" =>  $align_stop,                  "end" => $base_start + (($query_start-$hit_start)+$ln_hit),
1911                  "type"=> 'box',                  "type"=> 'box',
1912                  "color"=> $color,                  "color"=> $color,
1913                  "zlayer" => "2",                  "zlayer" => "2",
1914                  "links_list" => $links_list,                  "links_list" => $hit_links_list,
1915                  "description" => $descriptions                  "description" => $hit_descriptions
1916                    };
1917        push(@$line_data,$element_hash);
1918    
1919        $element_hash = {
1920            "title" => $peg . ': HIT AREA',
1921            "start" => $base_start + $query_start,
1922            "end" =>  $base_start + $query_stop,
1923            "type"=> 'smallbox',
1924            "color"=> $hit_color,
1925            "zlayer" => "3",
1926            "links_list" => $hit_links_list,
1927            "description" => $hit_descriptions
1928                  };                  };
1929              push(@$line_data,$element_hash);              push(@$line_data,$element_hash);
1930    
1931              $gd->add_line($line_data, $line_config);              $gd->add_line($line_data, $line_config);
1932          }  
1933      }      my $breaker = [];
1934        my $breaker_hash = {};
1935        my $breaker_config = { 'no_middle_line' => "1" };
1936    
1937        push (@$breaker, $breaker_hash);
1938        $gd->add_line($breaker, $breaker_config);
1939    
1940      return ($gd);      return ($gd);
1941  }  }
1942    
# Line 1937  Line 2071 
2071      my $data = [];      my $data = [];
2072      my $count = 0;      my $count = 0;
2073      my $content;      my $content;
     #my $fig = new FIG;  
2074      my $cgi = new CGI;      my $cgi = new CGI;
2075      my @ids;      my @ids;
2076      $lineages = $fig->taxonomy_list();      $lineages = $fig->taxonomy_list();
# Line 1947  Line 2080 
2080          push (@ids, $thing->acc);          push (@ids, $thing->acc);
2081      }      }
2082    
2083      my (%box_column, %subsystems_column, %evidence_column, %e_identical);      my (%box_column, %subsystems_column, %evidence_column, %e_identical, $function_color);
2084      my @attributes = $fig->get_attributes(\@ids);      my @attributes = $fig->get_attributes(\@ids);
2085    
2086      # get the column for the subsystems      # get the column for the subsystems
# Line 1959  Line 2092 
2092      # get the column for pfam_domain      # get the column for pfam_domain
2093      %pfam_column = &get_pfam_column(\@ids, \@attributes,$fig);      %pfam_column = &get_pfam_column(\@ids, \@attributes,$fig);
2094    
2095        # get the colors for the function cell
2096        my $functions = $fig->function_of_bulk(\@ids,1);
2097        $function_color = &get_function_color_cell($functions, $fig);
2098        my $query_function = $fig->function_of($query_fid);
2099    
2100      my %e_identical = &get_essentially_identical($query_fid,$dataset,$fig);      my %e_identical = &get_essentially_identical($query_fid,$dataset,$fig);
2101      my $alias_col = &get_aliases(\@ids,$fig);      my $alias_col = &get_aliases(\@ids,$fig);
2102      #my $alias_col = {};      #my $alias_col = {};
2103    
2104      my $figfam_data = &FIG::get_figfams_data();      my $figfam_data = &FIG::get_figfams_data();
2105      my $figfams = new FFs($figfam_data);      my $figfams = new FFs($figfam_data);
 #    my $ff_hash = $figfams->families_containing_peg_bulk(\@ids);  
2106    
2107      my %sims_objects_evalue;      my $func_color_offset=0;
     my $simcount = 0;  
2108      foreach my $thing (@$dataset){      foreach my $thing (@$dataset){
         if ($thing->class eq "SIM"){  
             $sims_objects_evalue{$simcount} = $thing->evalue;  
         }  
         $simcount++;  
     }  
   
     foreach my $index (sort {$sims_objects_evalue{$a}<=>$sims_objects_evalue{$b}} keys %sims_objects_evalue){  
 #    foreach my $thing ( @$dataset){  
         my $thing = $dataset->[$index];  
   
2109          next if ($thing->class ne "SIM");          next if ($thing->class ne "SIM");
2110          my $single_domain = [];          my $single_domain = [];
2111          $count++;          $count++;
# Line 1995  Line 2121 
2121          my $e2      = $thing->hstop;          my $e2      = $thing->hstop;
2122          my $d1      = abs($e1 - $b1) + 1;          my $d1      = abs($e1 - $b1) + 1;
2123          my $d2      = abs($e2 - $b2) + 1;          my $d2      = abs($e2 - $b2) + 1;
2124          my $reg1    = "$b1-$e1 (<b>$d1/$ln1</b>)";          my $color1  = match_color( $b1, $e1, $ln1 );
2125          my $reg2    = "$b2-$e2 (<b>$d2/$ln2</b>)";          my $color2  = match_color( $b2, $e2, $ln2 );
2126            my $reg1    = {'data'=> "$b1-$e1 (<b>$d1/$ln1</b>)", 'highlight' => $color1};
2127            my $reg2    = {'data'=> "$b2-$e2 (<b>$d2/$ln2</b>)", 'highlight' => $color2};
2128    
2129            # organisms cell
2130            my ($org, $org_color) = $fig->org_and_color_of($id);
2131            my $org_cell = { 'data' =>  $thing->organism, 'highlight' => $org_color};
2132    
2133          # checkbox column          # checkbox cell
2134            my ($box_cell,$tax);
2135          my $field_name = "tables_" . $id;          my $field_name = "tables_" . $id;
2136          my $pair_name = "visual_" . $id;          my $pair_name = "visual_" . $id;
2137          my $box_col = qq(<input type=checkbox name=seq value="$id" id="$field_name" onClick="VisualCheckPair('$field_name', '$pair_name');">);          my $cell_name = "cell_". $id;
2138          my ($tax) = ($id) =~ /fig\|(.*?)\./;          my $replace_id = $id;
2139            $replace_id =~ s/\|/_/ig;
2140            my $anchor_name = "anchor_". $replace_id;
2141            if ($id =~ /^fig\|/){
2142              my $box = qq(<a name="$anchor_name"></a><input type="checkbox" name=seq value="$id" id="$field_name" onClick="VisualCheckPair('$field_name', '$pair_name','$cell_name');">);
2143              $box_cell = { 'data'=>$box, 'highlight'=>$org_color};
2144              ($tax) = ($id) =~ /fig\|(.*?)\./;
2145            }
2146            else{
2147              my $box = qq(<a name="$anchor_name"></a>);
2148              $box_cell = { 'data'=>$box, 'highlight'=>$org_color};
2149            }
2150    
2151          # get the linked fig id          # get the linked fig id
2152          my $fig_col;          my $anchor_link = "graph_" . $replace_id;
2153          if (defined ($e_identical{$id})){          my $fig_data =  "<table><tr><td>" . &HTML::set_prot_links($cgi,$id) . "</td>" . "&nbsp;" x 2;
2154              $fig_col = "<a href='?page=Annotation&feature=$id'>$id</a>";#&HTML::set_prot_links($cgi,$id) . "*";          $fig_data .= qq(<td><img height='10px' width='20px' src='./Html/anchor_alignment' alt='View Graphic View of Alignment' onClick='changeSimsLocation("$anchor_link", 0)'/></td></tr></table>);
2155            my $fig_col = {'data'=> $fig_data,
2156                           'highlight'=>"#ffffff"};
2157    
2158        my $replace_id = $peg;
2159        $replace_id =~ s/\|/_/ig;
2160        my $anchor_name = "anchor_". $replace_id;
2161        my $query_config = { 'title' => "Query",
2162                             'short_title' => "Query",
2163                             'title_link' => "changeSimsLocation('$replace_id')",
2164                             'basepair_offset' => '0'
2165                             };
2166            # function cell
2167            my $function_cell_colors = {0=>"#ffffff", 1=>"#eeccaa", 2=>"#ffaaaa",
2168                                        3=>"#ffcc66", 4=>"#ffff00", 5=>"#aaffaa",
2169                                        6=>"#bbbbff", 7=>"#ffaaff", 8=>"#dddddd"};
2170            my $current_function =  $thing->function;
2171            my $function_color = $function_cell_colors->{ $function_color->{$current_function} - $func_color_offset};
2172            my $function_cell;
2173            if ($current_function){
2174              if ($current_function eq $query_function){
2175                $function_cell = {'data'=>$current_function, 'highlight'=>$function_cell_colors->{0}};
2176                $func_color_offset=1;
2177              }
2178              else{
2179                $function_cell = {'data'=>$current_function,'highlight' => $function_color};
2180              }
2181          }          }
2182          else{          else{
2183              $fig_col = "<a href='?page=Annotation&feature=$id'>$id</a>";#&HTML::set_prot_links($cgi,$id);            $function_cell = {'data'=>$current_function,'highlight' => "#dddddd"};
2184          }          }
2185    
2186          push (@$single_domain, $box_col, $fig_col, $thing->evalue,          push (@$single_domain, $box_cell, $fig_col, {'data'=> $thing->evalue, 'highlight'=>"#ffffff"},
2187                "$iden\%", $reg1, $reg2, $thing->organism, $thing->function);   # permanent columns                 {'data'=>"$iden\%", 'highlight'=>"#ffffff"}, $reg1, $reg2, $org_cell, $function_cell);   # permanent columns
2188    
2189          my ($ff) = $figfams->families_containing_peg($id);          my ($ff) = $figfams->families_containing_peg($id);
2190    
2191          foreach my $col (sort keys %$scroll_list){          foreach my $col (sort keys %$scroll_list){
2192              if ($col =~ /associated_subsystem/)          {push(@$single_domain,$subsystems_column{$id});}              if ($col =~ /associated_subsystem/)          {push(@$single_domain,{'data'=>$subsystems_column{$id},'highlight'=>"#ffffff"});}
2193              elsif ($col =~ /evidence/)                   {push(@$single_domain,$evidence_column{$id});}              elsif ($col =~ /evidence/)                   {push(@$single_domain,{'data'=>$evidence_column{$id},'highlight'=>"#ffffff"});}
2194              elsif ($col =~ /pfam_domains/)               {push(@$single_domain,$pfam_column{$id});}              elsif ($col =~ /pfam_domains/)               {push(@$single_domain,{'data'=>$pfam_column{$id},'highlight'=>"#ffffff"});}
2195              elsif ($col =~ /ncbi_id/)                    {push(@$single_domain,$alias_col->{$id}->{"NCBI"});}              elsif ($col =~ /ncbi_id/)                    {push(@$single_domain,{'data'=>$alias_col->{$id}->{"NCBI"},'highlight'=>"#ffffff"});}
2196              elsif ($col =~ /refseq_id/)                  {push(@$single_domain,$alias_col->{$id}->{"RefSeq"});}              elsif ($col =~ /refseq_id/)                  {push(@$single_domain,{'data'=>$alias_col->{$id}->{"RefSeq"},'highlight'=>"#ffffff"});}
2197              elsif ($col =~ /swissprot_id/)               {push(@$single_domain,$alias_col->{$id}->{"SwissProt"});}              elsif ($col =~ /swissprot_id/)               {push(@$single_domain,{'data'=>$alias_col->{$id}->{"SwissProt"},'highlight'=>"#ffffff"});}
2198              elsif ($col =~ /uniprot_id/)                 {push(@$single_domain,$alias_col->{$id}->{"UniProt"});}              elsif ($col =~ /uniprot_id/)                 {push(@$single_domain,{'data'=>$alias_col->{$id}->{"UniProt"},'highlight'=>"#ffffff"});}
2199              elsif ($col =~ /tigr_id/)                    {push(@$single_domain,$alias_col->{$id}->{"TIGR"});}              elsif ($col =~ /tigr_id/)                    {push(@$single_domain,{'data'=>$alias_col->{$id}->{"TIGR"},'highlight'=>"#ffffff"});}
2200              elsif ($col =~ /pir_id/)                     {push(@$single_domain,$alias_col->{$id}->{"PIR"});}              elsif ($col =~ /pir_id/)                     {push(@$single_domain,{'data'=>$alias_col->{$id}->{"PIR"},'highlight'=>"#ffffff"});}
2201              elsif ($col =~ /kegg_id/)                    {push(@$single_domain,$alias_col->{$id}->{"KEGG"});}              elsif ($col =~ /kegg_id/)                    {push(@$single_domain,{'data'=>$alias_col->{$id}->{"KEGG"},'highlight'=>"#ffffff"});}
2202              #elsif ($col =~ /trembl_id/)                  {push(@$single_domain,$alias_col->{$id}->{"TrEMBL"});}              #elsif ($col =~ /trembl_id/)                  {push(@$single_domain,{'data'=>$alias_col->{$id}->{"TrEMBL"},'highlight'=>"#ffffff"});}
2203              elsif ($col =~ /asap_id/)                    {push(@$single_domain,$alias_col->{$id}->{"ASAP"});}              elsif ($col =~ /asap_id/)                    {push(@$single_domain,{'data'=>$alias_col->{$id}->{"ASAP"},'highlight'=>"#ffffff"});}
2204              elsif ($col =~ /jgi_id/)                     {push(@$single_domain,$alias_col->{$id}->{"JGI"});}              elsif ($col =~ /jgi_id/)                     {push(@$single_domain,{'data'=>$alias_col->{$id}->{"JGI"},'highlight'=>"#ffffff"});}
2205              elsif ($col =~ /taxonomy/)                   {push(@$single_domain,$lineages->{$tax});}              elsif ($col =~ /taxonomy/)                   {push(@$single_domain,{'data'=>$lineages->{$tax},'highlight'=>"#ffffff"});}
2206              #elsif ($col =~ /taxonomy/)                   {push(@$single_domain,$fig->taxonomy_of($taxid));}              #elsif ($col =~ /taxonomy/)                   {push(@$single_domain,$fig->taxonomy_of($taxid));}
2207              #elsif ($col =~ /figfam/)                     {push(@$single_domain,"<a href='?page=FigFamViewer&figfam=" . $ff_hash->{$id} . "' target='_new'>" . $ff_hash->{$id} . "</a>");}              #elsif ($col =~ /figfam/)                     {push(@$single_domain,"<a href='?page=FigFamViewer&figfam=" . $ff_hash->{$id} . "' target='_new'>" . $ff_hash->{$id} . "</a>");}
2208              elsif ($col =~ /figfam/)                     {push(@$single_domain,"<a href='?page=FigFamViewer&figfam=" . $ff . "' target='_new'>" . $ff . "</a>");}              elsif ($col =~ /figfam/)                     {push(@$single_domain,{'data'=>"<a href='?page=FigFamViewer&figfam=" . $ff . "' target='_new'>" . $ff . "</a>",'highlight'=>"#ffffff"});}
2209          }          }
2210          push(@$data,$single_domain);          push(@$data,$single_domain);
2211      }      }
# Line 2054  Line 2224 
2224      foreach my $id (@$ids){      foreach my $id (@$ids){
2225          my $field_name = "tables_" . $id;          my $field_name = "tables_" . $id;
2226          my $pair_name = "visual_" . $id;          my $pair_name = "visual_" . $id;
2227          $column{$id} = qq(<input type=checkbox name=seq value="$id" id="$field_name" onClick="VisualCheckPair('$field_name', '$pair_name');">);          my $cell_name = "cell_" . $id;
2228            $column{$id} = qq(<input type=checkbox name=seq value="$id" id="$field_name" onClick="VisualCheckPair('$field_name', '$pair_name', '$cell_name');">);
2229      }      }
2230      return (%column);      return (%column);
2231  }  }
# Line 2085  Line 2256 
2256      return (%column);      return (%column);
2257  }  }
2258    
2259    sub match_color {
2260        my ( $b, $e, $n , $rgb) = @_;
2261        my ( $l, $r ) = ( $e > $b ) ? ( $b, $e ) : ( $e, $b );
2262        my $hue = 5/6 * 0.5*($l+$r)/$n - 1/12;
2263        my $cov = ( $r - $l + 1 ) / $n;
2264        my $sat = 1 - 10 * $cov / 9;
2265        my $br  = 1;
2266        if ($rgb){
2267            return html2rgb( rgb2html( hsb2rgb( $hue, $sat, $br ) ) );
2268        }
2269        else{
2270            rgb2html( hsb2rgb( $hue, $sat, $br ) );
2271        }
2272    }
2273    
2274    sub hsb2rgb {
2275        my ( $h, $s, $br ) = @_;
2276        $h = 6 * ($h - floor($h));
2277        if ( $s  > 1 ) { $s  = 1 } elsif ( $s  < 0 ) { $s  = 0 }
2278        if ( $br > 1 ) { $br = 1 } elsif ( $br < 0 ) { $br = 0 }
2279        my ( $r, $g, $b ) = ( $h <= 3 ) ? ( ( $h <= 1 ) ? ( 1,      $h,     0      )
2280                                          : ( $h <= 2 ) ? ( 2 - $h, 1,      0      )
2281                                          :               ( 0,      1,      $h - 2 )
2282                                          )
2283                                        : ( ( $h <= 4 ) ? ( 0,      4 - $h, 1      )
2284                                          : ( $h <= 5 ) ? ( $h - 4, 0,      1      )
2285                                          :               ( 1,      0,      6 - $h )
2286                                          );
2287        ( ( $r * $s + 1 - $s ) * $br,
2288          ( $g * $s + 1 - $s ) * $br,
2289          ( $b * $s + 1 - $s ) * $br
2290        )
2291    }
2292    
2293    sub html2rgb {
2294        my ($hex) = @_;
2295        my ($r,$g,$b) = ($hex) =~ /^\#(\w\w)(\w\w)(\w\w)/;
2296        my $code = { 'A'=>10, 'B'=>11, 'C'=>12, 'D'=>13, 'E'=>14, 'F'=>15,
2297                     1=>1, 2=>2, 3=>3, 4=>4, 5=>5, 6=>6, 7=>7, 8=>8, 9=>9};
2298    
2299        my @R = split(//, $r);
2300        my @G = split(//, $g);
2301        my @B = split(//, $b);
2302    
2303        my $red = ($code->{uc($R[0])}*16)+$code->{uc($R[1])};
2304        my $green = ($code->{uc($G[0])}*16)+$code->{uc($G[1])};
2305        my $blue = ($code->{uc($B[0])}*16)+$code->{uc($B[1])};
2306    
2307        my $rgb = [$red, $green, $blue];
2308        return $rgb;
2309    
2310    }
2311    
2312    sub rgb2html {
2313        my ( $r, $g, $b ) = @_;
2314        if ( $r > 1 ) { $r = 1 } elsif ( $r < 0 ) { $r = 0 }
2315        if ( $g > 1 ) { $g = 1 } elsif ( $g < 0 ) { $g = 0 }
2316        if ( $b > 1 ) { $b = 1 } elsif ( $b < 0 ) { $b = 0 }
2317        sprintf("#%02x%02x%02x", int(255.999*$r), int(255.999*$g), int(255.999*$b) )
2318    }
2319    
2320    sub floor {
2321        my $x = $_[0];
2322        defined( $x ) || return undef;
2323        ( $x >= 0 ) || ( int($x) == $x ) ? int( $x ) : -1 - int( - $x )
2324    }
2325    
2326    sub get_function_color_cell{
2327      my ($functions, $fig) = @_;
2328    
2329      # figure out the quantity of each function
2330      my %hash;
2331      foreach my $key (keys %$functions){
2332        my $func = $functions->{$key};
2333        $hash{$func}++;
2334      }
2335    
2336      my %func_colors;
2337      my $count = 1;
2338      foreach my $key (sort {$hash{$b}<=>$hash{$a}} keys %hash){
2339        $func_colors{$key}=$count;
2340        $count++;
2341      }
2342    
2343      return \%func_colors;
2344    }
2345    
2346  sub get_essentially_identical{  sub get_essentially_identical{
2347      my ($fid,$dataset,$fig) = @_;      my ($fid,$dataset,$fig) = @_;
2348      #my $fig = new FIG;      #my $fig = new FIG;
# Line 2520  Line 2778 
2778                  $prev_stop = $stop;                  $prev_stop = $stop;
2779                  $prev_fig = $fid1;                  $prev_fig = $fid1;
2780    
2781                  if ((defined($reverse_flag{$region_genome})) && ($reverse_flag{$region_genome} eq $all_genes{$fid1})){                  if ((defined($reverse_flag{$region_genome})) && ($reverse_flag{$region_genome} eq $all_gnes{$fid1})){
2782                      $start = $gd_window_size - $start;                      $start = $gd_window_size - $start;
2783                      $stop = $gd_window_size - $stop;                      $stop = $gd_window_size - $stop;
2784                  }                  }
# Line 2760  Line 3018 
3018    
3019      return($content);      return($content);
3020  }  }
3021    

Legend:
Removed from v.1.57  
changed lines
  Added in v.1.58

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3