再问一个文本处理的问题！谢谢几位啦！

yn521yn

#@        2492        NA20349        0        0        1        0        ASW
#@        2494        NA20357        0        0        2        0        ASW
#@        2494        NA20356        0        0        1        0        ASW
#@        2495        NA20359        0        0        2        0        ASW
#@        2496        NA20363        0        0        2        0        ASW
rs# alleles chrom pos strand assembly# center protLSID assayLSID panelLSID QCcode NA19625 NA19700 NA19701
rs7235612 A/G chr18 1543 + ncbi_b36 broad urn:LSID:affymetrix.orgrotocol:GenomeWideSNP_6.0:3 urn:LSID
rs10853286 A/G chr18 2842 + ncbi_b36 sanger urn:LSID:illumina.orgrotocol:Human_1M_BeadChip:3 urn:LSID
。。。。
。。。。
。。。。
再问一个弱弱的问题，我有如上的一个数据，其中#@开头和rs#开头的行可看做表头，需要对表头以下的数据进行分割，具体分割不用管，想生成多个分割的小文件，其中表头都与大文件的表头一样，perl语句应该怎么写啊？想了几天了。。。
谢谢了！

[ 本帖最后由 yn521yn 于 2009-9-26 14:34 编辑 ]

冰雪兰心

你的描述感觉不太清楚哦
你指的数据分隔是按行分隔，还是按列分隔？
“其中表头都与大文件的表头一样”哪个大文件，没看懂

cxfcxf

很简单啊 思路就是 吧表头的文件全部读出来 放到一个数组里面去

如果是每一行的话
然在读文件的时候 不是表头的行 每一行 就print一次这个数组+本身这一行 到一个文件

你也没说具体需求....反正思路是这样

cxfcxf

1. 
   
2. #!/usr/bin/perl
   
3. 
   
4. use strict;
   
5. use warnings;
   
6. 
   
7. my @data = <DATA>;
   
8. my @array;
   
9. foreach (@data) {
   
10.     chomp;
    
11.     push @array, $_ if (/^\#\@/);
    
12.     push @array, $_ if (/^rs\#/);
    
13. }
    
14. 
    
15. foreach (@data) {
    
16.     chomp;
    
17.     if (!/^\#\@/) {
    
18.         if (!/^rs\#/) {
    
19.             foreach my $head (@array) {
    
20.                 print "$head\n";
    
21.             }
    
22.             print "$_\n\n\n\n";
    
23.         }
    
24.     }
    
25. }
    
26. 
    
27. 
    
28. __DATA__
    
29. #@        2492        NA20349        0        0        1        0        ASW
    
30. #@        2494        NA20357        0        0        2        0        ASW
    
31. #@        2494        NA20356        0        0        1        0        ASW
    
32. #@        2495        NA20359        0        0        2        0        ASW
    
33. #@        2496        NA20363        0        0        2        0        ASW
    
34. rs# alleles chrom pos strand assembly# center protLSID assayLSID panelLSID QCcode NA19625 NA19700 NA19701
    
35. rs7235612 A/G chr18 1543 + ncbi_b36 broad urn:LSID:affymetrix.orgrotocol:GenomeWideSNP_6.0:3 urn:LSID
    
36. rs10853286 A/G chr18 2842 + ncbi_b36 sanger urn:LSID:illumina.orgrotocol:Human_1M_BeadChip:3 urn:LSID
    

复制代码

结果

1. 
   
2. Press ENTER or type command to continue
   
3. #@        2492        NA20349        0        0        1        0        ASW
   
4. #@        2494        NA20357        0        0        2        0        ASW
   
5. #@        2494        NA20356        0        0        1        0        ASW
   
6. #@        2495        NA20359        0        0        2        0        ASW
   
7. #@        2496        NA20363        0        0        2        0        ASW
   
8. rs# alleles chrom pos strand assembly# center protLSID assayLSID panelLSID QCcode NA19625 NA19700 NA19701
   
9. rs7235612 A/G chr18 1543 + ncbi_b36 broad urn:LSID:affymetrix.orgrotocol:GenomeWideSNP_6.0:3 urn:LSID
   
10. 
    
11. 
    
12. 
    
13. #@        2492        NA20349        0        0        1        0        ASW
    
14. #@        2494        NA20357        0        0        2        0        ASW
    
15. #@        2494        NA20356        0        0        1        0        ASW
    
16. #@        2495        NA20359        0        0        2        0        ASW
    
17. #@        2496        NA20363        0        0        2        0        ASW
    
18. rs# alleles chrom pos strand assembly# center protLSID assayLSID panelLSID QCcode NA19625 NA19700 NA19701
    
19. rs10853286 A/G chr18 2842 + ncbi_b36 sanger urn:LSID:illumina.orgrotocol:Human_1M_BeadChip:3 urn:LSID
    

复制代码

改成写文件就可以了

yn521yn

按行分割的，大文件就是原文件。

HongLian3

这样子应该可以吧
#!/usr/local/bin/perl -w
use warnings;
use strict;

open(IN,"infile";

my $i = 0;
while(<IN>
{
  if( $_ =~ /^#@/ || $_ =~ /^rs#/ )
  {
    open(OUT,">outfile$i";
    print $_ OUT;
    $i++;
  }
  else
  {
    print $_ OUT;
  }
}