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

Diff of /FigKernelScripts/svr_corr_by_exp.pl

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

revision 1.1, Tue Feb 8 19:58:53 2011 UTC revision 1.2, Tue Feb 8 20:04:44 2011 UTC
# Line 7  Line 7 
7  #  #
8    
9    
10  =head1 svr_corr_by_exp  =head1 svr_corr_by_exp [-m MinPCC]
11    
12  Get genes that have similar expression profiles.  Get genes that have similar expression profiles.
13    
# Line 44  Line 44 
44    
45  This is used only if the column containing PEGs is not the last.  This is used only if the column containing PEGs is not the last.
46    
47    =item -m Minimum value for the Pearson correlation coefficient
48    
49  =back  =back
50    
51  =head2 Output Format  =head2 Output Format
# Line 64  Line 66 
66  my $usage = "usage: svr_corr_by_exp [-c column]";  my $usage = "usage: svr_corr_by_exp [-c column]";
67    
68  my $column;  my $column;
69  my $rc  = GetOptions('c=i' => \$column);  my $min_pcc = 0;
70    my $rc  = GetOptions('c=i' => \$column,
71                         'm=i' => \$min_pcc);
72    
73  if (! $rc) { print STDERR $usage; exit }  if (! $rc) { print STDERR $usage; exit }
74    
75  my @lines = map { chomp; [split(/\t/,$_)] } <STDIN>;  my @lines = map { chomp; [split(/\t/,$_)] } <STDIN>;
# Line 81  Line 86 
86          foreach my $peg2 (@pegs2)          foreach my $peg2 (@pegs2)
87          {          {
88              my $pcc = sprintf("%0.3f",$x->{$peg2});              my $pcc = sprintf("%0.3f",$x->{$peg2});
89                if ($pcc >= $min_pcc)
90                {
91              print join("\t",(@$_,$pcc,$peg2)),"\n";              print join("\t",(@$_,$pcc,$peg2)),"\n";
92          }          }
93      }      }
94  }  }
95    }

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3