[Bio] / FigKernelScripts / exp_pc_compute.pl Repository:
ViewVC logotype

Diff of /FigKernelScripts/exp_pc_compute.pl

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

revision 1.7, Fri Jul 22 19:12:45 2011 UTC revision 1.8, Tue Aug 9 18:49:15 2011 UTC
# Line 196  Line 196 
196      require Statistics::Descriptive;      require Statistics::Descriptive;
197      for (my $i = 0; $i < @$gene_ids-1; $i++) {      for (my $i = 0; $i < @$gene_ids-1; $i++) {
198          my $stat = Statistics::Descriptive::Full->new();          my $stat = Statistics::Descriptive::Full->new();
199          $stat->add_data(@{$gxp_hash->{$gene_ids->[$i]}});          my $iData = $gxp_hash->{$gene_ids->[$i]};
200            $stat->add_data(@$iData);
201    
202          for (my $j = $i+1; $j < @$gene_ids; $j++)          for (my $j = $i+1; $j < @$gene_ids; $j++)
203          {          {
204              my ($q, $m, $r, $err) = $stat->least_squares_fit(@{$gxp_hash->{$gene_ids->[$j]}});              my $jData = $gxp_hash->{$gene_ids->[$j]};
205                my ($q, $m, $r, $err) = $stat->least_squares_fit(@$jData);
206              $values{$gene_ids->[$i]}->{$gene_ids->[$j]} = $r;              $values{$gene_ids->[$i]}->{$gene_ids->[$j]} = $r;
207              $values{$gene_ids->[$j]}->{$gene_ids->[$i]} = $r;              $values{$gene_ids->[$j]}->{$gene_ids->[$i]} = $r;
208          }          }

Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3