文件批量输出问题

最爱在深秋

有如下格式文件(每列用\t隔开的)：

ID_REF  GSM717702       GSM717703       GSM717704       GSM717705       GSM717706
AFFX-BioB-3_at  7.081480562     7.562581761     7.470186073     7.109637558     7.925162211
AFFX-BioB-5_at  7.129958916     7.45944498      7.640895702     7.389891485     8.02738222
AFFX-BioB-M_at  7.271660603     7.932285748     7.903917024     7.751718444     8.232053431


我想得到5个文件，文件名分别是第一行第二列到最后一列的字符，即GSM717702.xls、GSM717703.xls、 GSM717704.xls、 GSM717705.xls、  GSM717706.xls

对应的内容分别是第一列和该文件名对应的列，例如对于文件GSM717702.xls的内容为：
ID_REF  GSM717702
AFFX-BioB-3_at  7.081480562
AFFX-BioB-5_at  7.129958916
AFFX-BioB-M_at  7.271660603

对于文件GSM717703.xls的内容为：
ID_REF  GSM717703
AFFX-BioB-3_at  7.562581761
AFFX-BioB-5_at  7.45944498
AFFX-BioB-M_at  7.932285748


我的文件有很多列，不可能写cut命令一个个把文件cut出来，请教各位大侠，有什么高招。

最爱在深秋

想到一个比较笨的方法：


#!/usr/bin/perl
use strict;
use warnings;

my $in=shift;
open IN,$in;
my $head =<IN>; chomp $in;
my @a=split /\t+/,$head,2;
my @b=split /\t+/,$a[1];
my $n=2;
foreach my $i (@b){
        system "cut -f1,$n $in > $i.exp";
        $n++;
}

最爱在深秋

不对，最后一个文件的输出文件名有问题。。。

回复 2# 最爱在深秋


   

rubyish

biru:

1. #!/usr/bin/perl
   
2. 
   
3. my ( undef, @file ) = split /\s+/, <DATA>;
   
4. my @cache;
   
5. 
   
6. while (<DATA>) {
   
7.     my ( $affx, @data ) = split;
   
8.     for my $i ( 0 .. $#data ) {
   
9.         push @{ $cache[$i] }, "$affx\t$data[$i]\n";
   
10.     }
    
11. 
    
12. }
    
13. 
    
14. for my $i ( 0 .. $#file ) {
    
15.     open my $fh, '>', "$file[$i].xls";
    
16.     print $fh @{ $cache[$i] };
    
17. }
    
18. 
    
19. __DATA__
    
20. ID_REF  GSM717702       GSM717703       GSM717704       GSM717705       GSM717706
    
21. AFFX-BioB-3_at  7.081480562     7.562581761     7.470186073     7.109637558     7.925162211
    
22. AFFX-BioB-5_at  7.129958916     7.45944498      7.640895702     7.389891485     8.02738222
    
23. AFFX-BioB-M_at  7.271660603     7.932285748     7.903917024     7.751718444     8.232053431
    

复制代码

最爱在深秋

貌似不行噢~~
回复 4# rubyish


   

jason680

回复 2# 最爱在深秋

why to use chomp $in ??
   

my $head =<IN>; chomp $in;
   

最爱在深秋

It's a mistake.
Finally I modify it.

#!/usr/bin/perl

use strict;
use warnings;
my $in=shift;

open IN,$in;
my $head =<IN>; chomp $head;
my @a=split /\t+/,$head,2;
my @b=split /\t+/,$a[1];
my $n=2;
foreach my $i (@b){
        next if ($n >145);   ## total number of files
        system "cut -f1,$n $in > $i.exp";
        $n++;
}

回复 6# jason680


   

q1208c

回复 7# 最爱在深秋

我很好奇, 你已经使用了 perl了, 为啥又调用 system 命令 使用 shell 命令来切分文件? 这样做有啥好处么?

直接把你的每一行 split 成 @row, 然后 把 $row[$i] print 到 $FILE[$i] 里不就行了?
   

rubyish

sorry~ {:3_188:}:
update.

1. #!/usr/bin/perl
   
2. 
   
3. my ( $id, @file ) = split /\s+/, <DATA>;
   
4. my @cache;
   
5. 
   
6. while (<DATA>) {
   
7.     my ( $affx, @data ) = split;
   
8.     for my $i ( 0 .. $#data ) {
   
9.         push @{ $cache[$i] }, "$affx\t$data[$i]$/";
   
10.     }
    
11. }
    
12. 
    
13. for my $i ( 0 .. $#file ) {
    
14.     open my $fh, '>', "$file[$i].xls";
    
15.     print $fh $id, "\t", $file[$i], $/;
    
16.     print $fh @{ $cache[$i] };
    
17. }
    
18. 
    
19. __DATA__
    
20. ID_REF  GSM717702       GSM717703       GSM717704       GSM717705       GSM717706
    
21. AFFX-BioB-3_at  7.081480562     7.562581761     7.470186073     7.109637558     7.925162211
    
22. AFFX-BioB-5_at  7.129958916     7.45944498      7.640895702     7.389891485     8.02738222
    
23. AFFX-BioB-M_at  7.271660603     7.932285748     7.903917024     7.751718444     8.232053431
    

复制代码

output:
GSM717706.xls

1. ID_REF        GSM717706
   
2. AFFX-BioB-3_at        7.925162211
   
3. AFFX-BioB-5_at        8.02738222
   
4. AFFX-BioB-M_at        8.232053431
   

复制代码