annotate generate_matrix_for_pca_lda.pl @ 0:263a3b92368b draft default tip

Imported from capsule None
author devteam
date Mon, 28 Jul 2014 11:30:18 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
1 #!/usr/bin/perl -w
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
2
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
3 use strict;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
4 use warnings;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
5
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
6 my $Input_Matrix = $ARGV[0];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
7 my $Input_Label = $ARGV[1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
8
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
9 my %Hash_X = ();
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
10 my %Hash_Y = ();
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
11 my $My_Num_X = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
12 my $My_Num_Y = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
13
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
14 open (OUT, "> $ARGV[2]");
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
15
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
16 open (LABEL, "< $Input_Label") ||
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
17 die "Sorry, I couldn't open the escape.txt for clone: $!\n";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
18
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
19 my $Label_Index = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
20 my $X_Label;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
21 my $input_Label;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
22 while (defined($input_Label = <LABEL>)){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
23 chomp($input_Label);
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
24 my @cArray_Label = $input_Label =~ /(\S+)\s*/g;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
25 if ($input_Label =~ /\w/){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
26 if ($Label_Index == 0){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
27 $Hash_X{$cArray_Label[0]} = $cArray_Label[1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
28 $X_Label = $cArray_Label[1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
29 $Label_Index = 1;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
30 }else{
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
31 if ($cArray_Label[1] eq $X_Label){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
32 $Hash_X{$cArray_Label[0]} = $cArray_Label[1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
33 }else{
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
34 $Hash_Y{$cArray_Label[0]} = $cArray_Label[1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
35 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
36 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
37 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
38 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
39 close(LABEL);
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
40
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
41 open (MATRIX, "< $Input_Matrix") ||
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
42 die "Sorry, I couldn't open the escape.txt for clone: $!\n";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
43
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
44 my %Hash_Matrix = ();
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
45 my %Hash_Features = ();
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
46 my @cArray_Features = ();
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
47
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
48 my %Hash_Sum = ();
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
49 my $Matrix_Index = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
50 my $input_Matrix;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
51 while (defined($input_Matrix = <MATRIX>)){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
52 chomp($input_Matrix);
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
53 my @cArray_Matrix = $input_Matrix =~ /(\S+)\s*/g;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
54 if ($input_Matrix =~ /\w/){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
55 if ($Matrix_Index == 0){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
56 @cArray_Features = @cArray_Matrix;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
57 my $Temp_Num_Array = scalar(@cArray_Matrix);
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
58 my $Temp_Index = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
59 for(;$Temp_Index < $Temp_Num_Array; $Temp_Index++){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
60 $Hash_Features{$cArray_Matrix[$Temp_Index]} = "BOL";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
61 $Hash_Sum{$cArray_Matrix[$Temp_Index]} = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
62 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
63 $Matrix_Index = 1;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
64 }else{
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
65 $Hash_Matrix{$cArray_Matrix[0]} = $input_Matrix;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
66 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
67 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
68 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
69 close(MATRIX);
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
70
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
71 my $Trace_Key;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
72
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
73 foreach $Trace_Key (sort {$a cmp $b} keys %Hash_X){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
74 my @cArray_Trace_X = $Hash_Matrix{$Trace_Key} =~ /(\S+)\s*/g;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
75 my $Num_Array_Feature_X = scalar(@cArray_Features);
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
76 my $Index_Feature_X = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
77 for(;$Index_Feature_X < $Num_Array_Feature_X; $Index_Feature_X++){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
78 if ($Hash_Features{$cArray_Features[$Index_Feature_X]} eq "BOL"){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
79 $Hash_Features{$cArray_Features[$Index_Feature_X]} = $cArray_Trace_X[$Index_Feature_X + 1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
80 }else{
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
81 $Hash_Features{$cArray_Features[$Index_Feature_X]} = $Hash_Features{$cArray_Features[$Index_Feature_X]} . "\t" . $cArray_Trace_X[$Index_Feature_X + 1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
82 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
83
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
84 $Hash_Sum{$cArray_Features[$Index_Feature_X]} += $cArray_Trace_X[$Index_Feature_X + 1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
85 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
86 $My_Num_X ++;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
87 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
88
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
89 my $Append_Key;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
90 foreach $Append_Key (keys %Hash_Features){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
91 $Hash_Features{$Append_Key} = $Hash_Features{$Append_Key} . "\t" . $Hash_Sum{$Append_Key};
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
92 $Hash_Sum{$Append_Key} = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
93 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
94
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
95 foreach $Trace_Key (sort {$a cmp $b} keys %Hash_Y){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
96 my @cArray_Trace_Y = $Hash_Matrix{$Trace_Key} =~ /(\S+)\s*/g;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
97 my $Num_Array_Feature_Y = scalar(@cArray_Features);
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
98 my $Index_Feature_Y = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
99 for(;$Index_Feature_Y < $Num_Array_Feature_Y; $Index_Feature_Y++){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
100 if ($Hash_Features{$cArray_Features[$Index_Feature_Y]} eq "BOL"){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
101 $Hash_Features{$cArray_Features[$Index_Feature_Y]} = $cArray_Trace_Y[$Index_Feature_Y + 1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
102 }else{
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
103 $Hash_Features{$cArray_Features[$Index_Feature_Y]} = $Hash_Features{$cArray_Features[$Index_Feature_Y]} . "\t" . $cArray_Trace_Y[$Index_Feature_Y + 1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
104 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
105
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
106 $Hash_Sum{$cArray_Features[$Index_Feature_Y]} += $cArray_Trace_Y[$Index_Feature_Y + 1];
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
107 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
108 $My_Num_Y ++;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
109 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
110
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
111 foreach $Append_Key (keys %Hash_Features){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
112 $Hash_Features{$Append_Key} = $Hash_Features{$Append_Key} . "\t" . $Hash_Sum{$Append_Key} . "\t" . "EOL";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
113 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
114
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
115 my $Prt_Key;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
116 print OUT " \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
117 foreach $Prt_Key (sort {$a cmp $b} keys %Hash_X){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
118 print OUT "$Prt_Key \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
119 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
120 print OUT "X(SUM) \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
121
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
122 foreach $Prt_Key (sort {$a cmp $b} keys %Hash_Y){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
123 print OUT "$Prt_Key \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
124 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
125 print OUT "Y(SUM) \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
126 print OUT "\n";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
127
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
128 my $Prt_Index = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
129 my $Prt_Array_Num = scalar (@cArray_Features);
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
130 for(;$Prt_Index < $Prt_Array_Num; $Prt_Index++){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
131 print OUT "$cArray_Features[$Prt_Index] \t$Hash_Features{$cArray_Features[$Prt_Index]}\n";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
132 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
133
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
134 print OUT " \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
135 my $My_Label_Index = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
136 for(;$My_Label_Index < $My_Num_X; $My_Label_Index++){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
137 print OUT "X \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
138 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
139 print OUT " \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
140
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
141 $My_Label_Index = 0;
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
142 for(;$My_Label_Index < $My_Num_Y; $My_Label_Index++){
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
143 print OUT "Y \t";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
144 }
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
145 print OUT " \t\n";
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
146
263a3b92368b Imported from capsule None
devteam
parents:
diff changeset
147 close(OUT);