|
0
|
1 #!usr/bin/perl
|
|
|
2 package IonFiltration;
|
|
|
3
|
|
|
4 ### Perl modules
|
|
|
5 use strict;
|
|
|
6 use warnings;
|
|
|
7
|
|
|
8
|
|
|
9
|
|
|
10
|
|
|
11
|
|
|
12
|
|
|
13 ########################################################################
|
|
|
14 ### Création of a hash containing all adduits and fragments possible ###
|
|
|
15 ########################################################################
|
|
|
16
|
|
|
17
|
|
|
18 sub MassCollecting{
|
|
|
19
|
|
|
20 my $mass_file = $_[0];
|
|
|
21 my %hmass;
|
|
|
22
|
|
|
23 open (F1, $mass_file);
|
|
|
24
|
|
|
25 while(my $line = <F1>){
|
|
|
26 chomp $line;
|
|
|
27 my @tline = split(/[\t;]/, $line);
|
|
|
28 if(defined($hmass{$tline[2]})){
|
|
|
29 print "The mass difference already exists : $tline[2] !\n";
|
|
|
30 }
|
|
|
31 $hmass{$tline[1]}{$tline[2]}=$tline[0];
|
|
|
32 }
|
|
|
33
|
|
|
34 close F1;
|
|
|
35 return %hmass;
|
|
|
36
|
|
|
37 }
|
|
|
38
|
|
|
39
|
|
|
40
|
|
|
41
|
|
|
42
|
|
|
43
|
|
|
44
|
|
|
45 ########################################################
|
|
|
46 ### Creation of a sif table + correlation filtration ###
|
|
|
47 ########################################################
|
|
|
48
|
|
|
49
|
|
|
50 sub sifTableCreation{
|
|
|
51
|
|
|
52 my $file = $_[0];
|
|
|
53 my $output_sif = $_[1];
|
|
|
54 # my $opt = $_[2];
|
|
|
55 # my $rt_threshold = $_[3];
|
|
|
56 # my $mass_threshold = $_[4];
|
|
|
57 my $correl_threshold = $_[5];
|
|
|
58 # my $dataMatrix = $_[6];
|
|
|
59 # my $output_tabular = $_[7];
|
|
|
60 my $combined_DMVM = $_[8];
|
|
|
61 # my $repres_opt = $_[9];
|
|
|
62 # my $intensity_threshold = $_[10];
|
|
|
63 # my $intensity_pourc = $_[11];
|
|
|
64 # my $refhmass = $_[12];
|
|
|
65
|
|
|
66
|
|
|
67
|
|
|
68
|
|
|
69 my %hheader_file;
|
|
|
70 my %hduplicate;
|
|
|
71
|
|
|
72 my %hcorrelgroup;
|
|
|
73 my $groupct=1;
|
|
|
74
|
|
|
75
|
|
|
76 my $linenb3=0;
|
|
|
77 my %hheader_line;
|
|
|
78 my %hrtmz;
|
|
|
79
|
|
3
|
80 open (F5, $combined_DMVM) or die "Impossible to open $file\n";
|
|
0
|
81 while(my $line = <F5>){
|
|
|
82 chomp $line;
|
|
3
|
83
|
|
0
|
84 my @tline = split(/\t/, $line);
|
|
|
85
|
|
|
86 if($linenb3 == 0){
|
|
|
87 for(my $i=0; $i<scalar(@tline);$i++){
|
|
|
88 my $a = $tline[$i];
|
|
|
89 $hheader_line{$a}=$i;
|
|
|
90 }
|
|
|
91 }
|
|
|
92 else{
|
|
|
93 if(defined($hheader_line{mzmed})){
|
|
|
94 my $b = $tline[$hheader_line{mzmed}];
|
|
|
95 $hrtmz{$tline[0]}{mz}=$b;
|
|
|
96 }
|
|
|
97 else{
|
|
|
98 my $b = $tline[$hheader_line{mz}];
|
|
|
99 $hrtmz{$tline[0]}{mz}=$b;
|
|
|
100 }
|
|
|
101 if(defined($hheader_line{rtmed})){
|
|
|
102 my $d = $tline[$hheader_line{rtmed}];
|
|
|
103 $hrtmz{$tline[0]}{rt}=$d;
|
|
|
104 }
|
|
|
105 else{
|
|
|
106 my $d = $tline[$hheader_line{rt}];
|
|
|
107 $hrtmz{$tline[0]}{rt}=$d;
|
|
|
108 }
|
|
|
109 }
|
|
|
110
|
|
|
111 $linenb3 ++;
|
|
|
112 }
|
|
|
113 close F5;
|
|
|
114
|
|
|
115
|
|
|
116 my $linenb=0;
|
|
|
117
|
|
|
118 open (F1, $file) or die "Impossible to open $file\n";
|
|
|
119 open(F2, ">$output_sif") or die "Impossible to open $output_sif\n";
|
|
|
120
|
|
|
121
|
|
|
122 while(my $line = <F1>){
|
|
|
123 chomp $line;
|
|
|
124 my @tline = split(/\t/, $line);
|
|
|
125
|
|
|
126 ###############################
|
|
|
127 ### Création of a sif table ###
|
|
|
128 ###############################
|
|
|
129
|
|
|
130 if($linenb == 0){
|
|
|
131 for(my $i=0; $i<scalar(@tline);$i++){
|
|
|
132 my $a = $tline[$i];
|
|
|
133 $hheader_file{$i}=$a;
|
|
|
134 }
|
|
|
135 }
|
|
|
136 else{
|
|
|
137 for(my $i=1; $i<scalar(@tline);$i++){
|
|
|
138 my $a=$tline[0];
|
|
|
139 my $b=$hheader_file{$i};
|
|
|
140 my $coef=$tline[$i];
|
|
|
141
|
|
|
142 if($a eq $b){
|
|
|
143 # print "This is a correlation between A ($a) and A ($b) !\n"
|
|
|
144 }
|
|
|
145 else{
|
|
|
146
|
|
|
147 #########################
|
|
|
148 ### Remove duplicates ###
|
|
|
149 #########################
|
|
|
150
|
|
|
151 my $y = $a."/".$b;
|
|
|
152 my $z = $b."/".$a;
|
|
|
153
|
|
|
154 if((!(defined($hduplicate{$y}))) && (!(defined($hduplicate{$z})))){
|
|
|
155
|
|
|
156 $hduplicate{$y}=1;
|
|
|
157 # my $abcoef=abs($coef); # Only when you want to consider negative correlations
|
|
|
158
|
|
|
159 # if($abcoef > $correl_threshold){ # Only when you want to consider negative correlations
|
|
|
160 if($coef > $correl_threshold){
|
|
|
161
|
|
|
162 print F2 "$a\t$coef\t$b\n";
|
|
|
163
|
|
|
164 my $count=0;
|
|
|
165
|
|
|
166 }
|
|
|
167 }
|
|
|
168 }
|
|
|
169 }
|
|
|
170 }
|
|
|
171 $linenb ++;
|
|
|
172 }
|
|
|
173 close F1;
|
|
|
174 close F2;
|
|
|
175 return ($output_sif, %hrtmz);
|
|
|
176 }
|
|
|
177
|
|
|
178
|
|
|
179
|
|
|
180
|
|
|
181
|
|
|
182 1; |
