#!/usr/bin/perl -w 

#
#set xlabel "diff of averages"
#set output "vplot_diffave_wty2h2_y2h2y2h2pb_lab.png"
#plot "wilcoxout6.txt" u 15:14:1 w labels title 'wt y2h2'  , "wilcoxout6.txt" u 29:28:1 w labels title 'y2h2 y2h2pb' 

#set xlabel "diff of averages"
#set output "vplot_diffave_wtwtpb_wtpby2h2pb.png"
#plot "wilcoxout6.txt" u 8:7:1 w labels title 'wt wtpb'  , "wilcoxout6.txt" u 22:21:1 w labels title 'wtpb y2h2pb' 


open(IF,"wilcoxout6.txt");
#gene info0_info1 z01 prob01 log_prob01 diff01 info0_info2 z02 prob02 log_prob02 diff02 info1_info3  z13 prob13 log_prob13 diff13 info2_info3 z23 prob23 log_prob23 diff23
#YAL003W wt_wtpb z -1.554420 p 0.120084 0.920514854167563 -0.00607395687598244 wt_y2h2 z -0.806515 p 0.41994 0.376812756102941 -0.00633111378107565 wtpb_y2h2pb z -2.800758 p 0.0050982 2.29258313136329 0.0118069917354435 y2h2_y2h2pb z -3.585686 p 0.0003362 3.47340229089655 0.0120641486405367
#YBL040C wt_wtpb z 1.706976 p 0.087826 1.05637689656572 0.00172670234562688 wt_y2h2 z 0.441275 p 0.65902 0.18110140519683 0.0130150276572483 wtpb_y2h2pb z -20.221901 p 6.2814e-91 90.2019435498117 0.492986968792204 y2h2_y2h2pb z -20.337395 p 6.0036e-92 91.221588251069 0.481698643480583
$head = <IF>;
open(OFWWPS,">signif_wo6_wtwtpb.txt");
open(OFWWPN,">notsignif_wo6_wtwtpb.txt");

open(OFWPYPS,">signif_wo6_wtpby2h2pb.txt");
open(OFWPYPN,">notsignif_wo6_wtpby2h2pb.txt");

open(OFWYS,">signif_wo6_wty2h2.txt");
open(OFWYN,">notsignif_wo6_wty2h2.txt");

open(OFYYPS,">signif_wo6_y2h2y2h2pb.txt");
open(OFYYPN,">notsignif_wo6_y2h2y2h2pb.txt");


print OFWWPS $head;
print OFWWPN $head;
print OFWPYPS $head;
print OFWPYPN $head;

print OFWYS $head;
print OFWYN $head;
print OFYYPS $head;
print OFYYPN $head;



while($a=<IF>){
@b = split(/ /,$a);
$p1 = $b[5];
$lp1 = $b[6];
$p2 = $b[12];
$lp2 = $b[13];
$p3 = $b[19];
$lp3 = $b[20];
$p4 = $b[26];
$lp4 = $b[27];

if($lp1>10 ){print OFWWPS $a;}
else{print OFWWPN $a;}

if($lp2>10 ){print OFWYS $a;}
else{print OFWYN $a;}

if($lp3>10 ){print OFWPYPS $a;}
else{print OFWPYPN $a;}

if($lp4>10 ){print OFYYPS $a;}
else{print OFYYPN $a;}


}