大侠们看看这个问题该怎么处理

newfinder

有一个文件，例如下：
VS1  A B C D A B
VS2  B A D C C D
VS3  A B A C B A
EM1  C D D C C D
EM2  A B B A D C
EM3  B A C D A B
EM4  A B A B B A
DK1  A B D C D C
DK2  B A C D C D

每两列为一个单位，我想统计一下分别含VS、EM的所有个体中第一个两列中A的频率，B的频率，C的频率，D的频率；第二个两列中A的频率，B的频率，C的频率，D的频率；第三个两列中A的频率，B的频率，C的频率，D的频率……（第四个两列……）

即最终想要的结果为：

                               VS                                                                     EM                                                  DK
first col        A: 0.5  B: 0.5  C: 0  D:0                             A: 0.375  B: 0.375  C: 0.125  D: 0.125       A: 0.5  B: 0.5  C: 0  D: 0
second col    A: 0.167  B: 0 C: 0.5 D: 0.333                    A: 0.25  B:0.25  C: 0.25  D: 0.25                A: 0  B: 0  C: 0.5  D: 0.5
third col       A: 0.333  B: 0.333  C: 0.167  D: 0.167        A:0.25  B: 0.25  C: 0.25  D: 0.25                A: 0  B: 0  C: 0.5  D: 0.5

跪求各位大侠，这样的问题用perl怎么处理呢。

bikong0411

用hash搞起啊

qq27898

我的方法比较笨。。。。。求大侠指导，我也好好学习一下

qq27898

回复 2# bikong0411


    用hash知道，那你准备定义几个hash？说说你的思路

xiumu2280

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my %hash_all_data;
   
6. while (<DATA>) {
   
7.         chomp;
   
8.         my @line = split;
   
9.         my ($title) = $line[0]=~/\D+/gs;
   
10.         shift @line;
    
11.         my $m=1;
    
12.         for (my $n = 0;$n<@line ;$n=$n+2) {
    
13.                 my @d = @line;
    
14.                 my @data = splice (@d,$n,2);
    
15.                 push @{$hash_all_data{$title}{$m}},@data;
    
16.                 $m++;
    
17.         }
    
18. }
    
19. foreach my $k (keys %hash_all_data) {
    
20.         foreach my $k1 (sort {$a <=> $b} keys %{$hash_all_data{$k}}) {
    
21.                 print "$k\t$k1 col\t";
    
22.                 my $all_n = @{$hash_all_data{$k}{$k1}};
    
23.                 my @A = grep {$_ eq "A"} @{$hash_all_data{$k}{$k1}};
    
24.                 my @B = grep {$_ eq "B"} @{$hash_all_data{$k}{$k1}};
    
25.                 my @C = grep {$_ eq "C"} @{$hash_all_data{$k}{$k1}};
    
26.                 my @D = grep {$_ eq "D"} @{$hash_all_data{$k}{$k1}};
    
27.                 my $a = @A/$all_n;
    
28.                 my $b = @B/$all_n;
    
29.                 my $c = @C/$all_n;
    
30.                 my $d = @D/$all_n;
    
31.                 print "A:$a\tB:$b\tC:$c\tD:$d\n";
    
32.         }
    
33. }
    
34. 
    
35. __DATA__
    
36. VS1         A B  C D  A B
    
37. VS2  B A  D C  C D
    
38. VS3  A B  A C  B A
    
39. EM1  C D  D C  C D
    
40. EM2  A B  B A  D C
    
41. EM3  B A  C D  A B
    
42. EM4  A B  A B  B A
    
43. DK1  A B  D C  D C
    
44. DK2  B A  C D  C D
    

复制代码

、
小数保留什么的自己调一下

newfinder

回复 5# xiumu2280
这个太牛了，必须得赞一个！


   

mcshell

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. my ( $n, $element_num );
   
5. while (<DATA>) {
   
6.     chomp;
   
7.     my ( $_, @line ) = split /\s{2}/;
   
8.     $element_num = $#line;
   
9.     /([A-Z]+)/;
   
10.     push @{ $n->{$1} }, [@line];
    
11. }
    
12. for my $key ( keys $n ) {
    
13.     print $key, "\n";
    
14.     for my $num ( 0 .. $element_num ) {
    
15.         my $s;
    
16.         my $m = { A => 0, B => 0, C => 0, D => 0 };
    
17.         for my $key2 ( 0 .. $#{ $n->{$key} } ) {
    
18.             $s += split " ", $n->{$key}[$key2]->[$num];
    
19.             $m->{$_}++ for split " ", $n->{$key}[$key2]->[$num];
    
20.         }
    
21.         for ( sort keys $m ) {
    
22.             my $p = sprintf "%0.3f", $m->{$_} / $s;
    
23.             print " $_ $p ";
    
24.         }
    
25.         print "\n";
    
26.     }
    
27. }
    
28. 
    
29. __DATA__
    
30. VS1  A B  C D  A B
    
31. VS2  B A  D C  C D
    
32. VS3  A B  A C  B A
    
33. EM1  C D  D C  C D
    
34. EM2  A B  B A  D C
    
35. EM3  B A  C D  A B
    
36. EM4  A B  A B  B A
    
37. DK1  A B  D C  D C
    
38. DK2  B A  C D  C D
    

复制代码

xiumu2280

大神膜拜一个~~~！！
回复 7# mcshell


   

yestreenstars

1. #!/usr/bin/perl
   
2. 
   
3. @a = qw/A B C D/;
   
4. 
   
5. while (<DATA>) {
   
6.         $len=split;
   
7.         @_[0] =~ s/\d+//;
   
8.         $h1{@_[0]}++;
   
9.         for ($i=2;$i<=@_;$i+=2) {
   
10.                 $h2{@_[0]}{$i/2}{@_[$i-1]}++;
    
11.                 $h2{@_[0]}{$i/2}{@_[$i]}++;
    
12.         }
    
13. }
    
14. 
    
15. END {
    
16.         foreach $k1 (keys %h1) {
    
17.                 print "$k1\n";
    
18.                 for ($i=1;$i<=($len-1)/2;$i++) {
    
19.                         print "$i ";
    
20.                         for ($j=0;$j<@a;$j++) {
    
21.                                 printf "%s: %.3f ",$a[$j],$h2{$k1}{$i}{$a[$j]}/($h1{$k1}*2);
    
22.                         }
    
23.                         print "\n";
    
24.                 }
    
25.         }
    
26. }
    
27. 
    
28. __DATA__
    
29. VS1  A B  C D  A B
    
30. VS2  B A  D C  C D
    
31. VS3  A B  A C  B A
    
32. EM1  C D  D C  C D
    
33. EM2  A B  B A  D C
    
34. EM3  B A  C D  A B
    
35. EM4  A B  A B  B A
    
36. DK1  A B  D C  D C
    
37. DK2  B A  C D  C D

复制代码

newfinder

回复 7# mcshell
VS1  A B C D A B
VS2  B A D C C D
VS3  A B A C B A
EM1  C D D C C D
EM2  A B B A D C
EM3  B A C D A B
EM4  A B A B B A
DK1  A B D C D C
DK2  B A C D C D
大侠，如果每列之间都是用空格隔开，只在个体号与后面的字母是tab键，其他列之间全是空格键，如上面的数据，程序要怎么改就可以了呢