#!/usr/bin/perl
use CGI qw(param);
use GD;
$id = param("id");
$id =~ tr/a-z/A-Z/;
print "Content-type:text/html\n\n";
print "\n";
print "The correlated genes with $id
\n";
$im = new GD::Image(300, 300);
$white = $im->colorAllocate(255,255,255);
$black = $im->colorAllocate(0,0,0);
$grey = $im->colorAllocate(200,200,200);
$red = $im->colorAllocate(250,0,0);
$green = $im->colorAllocate(0,255,0);
$blue = $im->colorAllocate(0,0,255);
$greenblue = $im->colorAllocate(0,255,255);
print "Using method of local pattern
";
open(RAW_DATA, "pseu_dp_worm_6.txt");
$ii = 0;
while($buffer= ){
chomp($buffer);
($gene1, $gene2, $score, $x, $y, $neg_pos) = split(/ +/, $buffer);
if($gene1 eq $id){
$level1{$ii} = $gene2;
if($x != $y){
$anno{$ii} = "t";
}
else{
$anno{$ii} = "i";
}
$ii++;
}
if($gene2 eq $id){
$level1{$ii} = $gene1;
if($x != $y){
$anno{$ii} = "t";
}
else{
$anno{$ii} = "i";
}
$ii++;
}
}
close(RAW_DATA);
$num_level1 = $ii;
if($num_level1 == 0){
print "No strong correlated genes found
";
}
if($num_level1 == 1){
$centerx = 120;
$centery = 150;
$im->arc($centerx,$centery,12,12,0,360,$red);
$im->fill($centerx, $centery, $red);
# $gn = $gene_ORF{$id};
$gn = $id;
if(length($gn) == 0){
$gn = $id;
}
$im->string(gdMediumBoldFont,$centerx-16,$centery+8,$gn, $red);
}
else{
$centerx = 150;
$centery = 150;
$im->arc($centerx,$centery,12,12,0,360,$red);
$im->fill($centerx, $centery, $red);
# $gn = $gene_ORF{$id};
$gn = $id;
if(length($gn) == 0){
$gn = $id;
}
$im->string(gdMediumBoldFont,$centerx-16,$centery+8,$gn, $red);
}
$dis = 80;
print "";
$imgname = $id . $num_method . ".png";
open(DISPLAY,">./tmp/$imgname");
flock(DISPLAY, LOCK_EX);
binmode DISPLAY;
print DISPLAY $im->png;
close(DISPLAY);
print "![](\"./tmp/$imgname\")
";
}
#print "
";
print "
black dots: time-shifted relationship
";
print "blue dots: inverted relationship
";
print "";
sub draw_node
{
my($xx, $yy, $nn, $gn, $cc) = @_;
$im->arc($xx,$yy,12,12,0,360,$cc);
$im->fill($xx, $yy, $cc);
$im->string(gdMediumBoldFont,$xx-16,$yy+8,$gn, $black);
$x1 = $xx-4;
$y1 = $yy-4;
$x2 = $xx+4;
$y2 = $yy+4;
print <
imgarea
}
sub draw_edge
{
my($x1, $y1, $x2, $y2) = @_;
$im->line($x1, $y1, $x2, $y2, $black);
}