求助，一个矩阵变换的问题！~

xiayuu119

首先献给大神们膝盖！~！涌泉相报~！

输入是这样：

ID1    ID2    a    b    c    d   ...
A     M       1    2    1    5
A     N       2    1    5    6
B     M       1    3    2    2
B     M       1    5    1    6
B     N       1    5    1    1
C     N       1    5    1    6
.
.
.


输出：

ID   a    A
M    1
N    2
ID   a    B
M    1
M    1
N    1
ID   a    C
N    1
.
.
.
ID   b    A
M    2
N    1
ID   b    B
M    3
M    5
N    5
ID   b    C
N    5

.
.
.

ID   c    A
M    1
N    5
ID   c    B
M    2
M    1
N    1
ID   c    C
N    1
.
.
.


ID   d    A
M    5
N    6
ID   d    B
M    2
M    6
N    1
ID   d    C
N    6
.
.
.



请各位帮忙~！~~~再次感谢~！

sunzhiguolu

回复 1# xiayuu119
你试下,

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my %hStat;
   
6. my @aHeader = grep {!m/\AI/} split ' ', scalar <DATA>;
   
7. map {push (@{$hStat{$_->[0]}}, $_)} map {[split]} <DATA>;
   
8. my $sOffset = 2;
   
9. foreach my $sVal (@aHeader){
   
10.     foreach my $sKey (sort keys %hStat){
    
11.         my @aData = @{$hStat{$sKey}};
    
12.         my $raFirst = splice (@aData, 0, 1);
    
13.         print "ID $sVal $raFirst->[0]\n";
    
14.         print "@$raFirst[1, $sOffset]\n";
    
15.         print "@{$_}[1, $sOffset]\n" for @aData;
    
16.     }
    
17.     $sOffset++;
    
18. }
    
19. 
    
20. __DATA__
    
21. ID1    ID2    a    b    c    d
    
22. A     M       1    2    1    5
    
23. A     N       2    1    5    6
    
24. B     M       1    3    2    2
    
25. B     M       1    5    1    6
    
26. B     N       1    5    1    1
    
27. C     N       1    5    1    6
    

复制代码

jason680

回复 1# xiayuu119

写完，才发现这是Perl板,自行改写就可以了..

$ awk '{for(n=3;n<=NF;++n){if(NR>1){if(!b[$1]++)b[++c]=$1;d[a[n],$1]=d[a[n],$1]$2" "$n"\n"}else a[n]=$n}}END{for(n=3;n<=NF;++n)for(m=1;m<=c;++m)printf("ID %s %s\n%s",a[n],b[m],d[a[n],b[m]])}' FILE
ID a A
M 1
N 2
ID a B
M 1
M 1
N 1
ID a C
N 1
ID b A
M 2
N 1
ID b B
M 3
M 5
N 5
ID b C
N 5
ID c A
M 1
N 5
ID c B
M 2
M 1
N 1
ID c C
N 1
ID d A
M 5
N 6
ID d B
M 2
M 6
N 1
ID d C
N 6

   

sunzhiguolu

输出结果:
------------------------------------------------------
ID a A
M 1
N 2
ID a B
M 1
M 1
N 1
ID a C
N 1
ID b A
M 2
N 1
ID b B
M 3
M 5
N 5
ID b C
N 5
ID c A
M 1
N 5
ID c B
M 2
M 1
N 1
ID c C
N 1
ID d A
M 5
N 6
ID d B
M 2
M 6
N 1
ID d C
N 6

sunzhiguolu

回复 3# jason680
大神, 能向您提个请求吗? (关于代码的)

   

jason680

回复 5# sunzhiguolu

本帖有关则直接发问（不必发两次...）
若与本帖无关，请发“消息”...
   

RE_HASH

$> perl -nae 'unless(@Hdr){@Hdr = @F}else{for my $i (2 .. $#F ){$O[$i] .= sprintf "%s$F[1] $F[$i]\n",($F[0] ne $Prev)?"ID $Hdr[$i] $F[0]\n":""} $Prev=$F[0]} END{print @O}' file
ID a A
M 1
N 2
ID a B
M 1
M 1
N 1
ID a C
N 1
ID b A
M 2
N 1
ID b B
M 3
M 5
N 5
ID b C
N 5
ID c A
M 1
N 5
ID c B
M 2
M 1
N 1
ID c C
N 1
ID d A
M 5
N 6
ID d B
M 2
M 6
N 1
ID d C
N 6

sunzhiguolu

向两位大神学习, 都太厉害了!!!

xiayuu119

向三位前辈学习！~~~太厉害了~！！~！再次感谢~！！！