有没有更简单的方法？

静女，其姝

#!usr/bin/perl -w
use strict;

open (IA,"$ARGV[0]") || die "input file can't open $!";
open (IB,"$ARGV[1]") || die "input file can't open $!";
my (@a,@b,@site,$s);
while (<IA>){
        chomp;
        @a=split;
        if ($a[1] ne 1){
        push (@site,$a[0]);}

}

while (<IB>){
        chomp;
        @b=split;
        foreach $s(@site){
                if ($s eq $b[0])
                {print "$_\n";}
        }
}


IA:
1        1
2        2
3        1
4        2
5        1


IB:
1        a        b
2        b        b
3        c        e
4        a        c
5        c        d

OUT:
2        b        b
4        a        c

sunzhiguolu

回复 1# 静女，其姝

if ($s eq $b[0])

确定上面这个判断逻辑可以得到下面的结果?

1. 
   
2. OUT:
   
3. 2        b        b
   
4. 4        a        c
   

复制代码

非常抱歉, 我没有仔细看.

if ($a[1] ne 1)

进行数值比较, 最好还是使用 (== | !=) 比较运算符较好, 不过在您的例子中没有问题.

sunzhiguolu

试下:

1. 
   
2. #!/usr/bin/perl -w
   
3. use strict;
   
4. 
   
5. my %hKeys;
   
6. while (<>){
   
7.     my ($sId, $sVal) = split;
   
8.     if ($sVal !~ m{\A[a-z]\z}){
   
9.         $hKeys{$sId} = '' if ($sVal != 1);
   
10.     }else{
    
11.         print if (exists $hKeys{$sId});
    
12.     }
    
13. }
    

复制代码

1. 
   
2. perl iab.pl ia ib
   

复制代码

静女，其姝

非常感谢，我先自己看看每一行的含义，有不懂之处再请教
另外，我用ne比较数值时确实碰壁过，多谢指点

回复 3# sunzhiguolu


   

静女，其姝

这个程序并没有简便太多，而且那个匹配只是适用example的数据，
请问，有没有更简单的方法（一行code）？grep 可以达到吗？


回复 3# sunzhiguolu


   

sunzhiguolu

回复 5# 静女，其姝

这个程序并没有简便太多，而且那个匹配只是适用example的数据，
请问，有没有更简单的方法（一行code）？grep 可以达到吗？

1.> 真实的数据啥样?
2.> 您能否提供一个一行 code 代码 show 一下?

我并没有冒犯您的意思, 如果言语不周请您见谅. 我只是针对问题而言.

   

Okelani

1. perl -ane '$h{$F[0]}?print:$F[1]!=1&&$h{$F[0]}++' file1 file2
   

复制代码

静女，其姝

是我太抱歉，索要太多

谢谢你不吝赐教
IA:
chr1        1        1
chr1        2        1
chr1        3        2
chr1        4        2
chr1        5        1       
chr1        6        2

IB：
chr1        1        A        A        A        A
chr1        2        T        T        T        T
chr1        3        C        G        C        G
chr1        4        G        G        G        A
chr1        5        A        A        A        A
chr1        6        A        A        T        A

output：
chr1        3        C        G        C        G
chr1        4        G        G        G        A
chr1        6        A        A        T        A

解释：两个文件的第一列都是染色体编号（忽略），第二列都是坐标（忽略），IB的三列到最后一列是

碱基，如果一行只有一种碱基，IA的第三列就是1，如果有2种碱基，IA的第三列就是2，现在要将IB中有

两种碱基的输出。虽然不用IA也可以得出结果，但是有其他比较复杂的情况，总之一定要根据IA来判断

。注：chr1有时也可能是Chr1或者包含T、C、G字母的其他string


回复 6# sunzhiguolu


   

sunzhiguolu

回复 8# 静女，其姝
一行 code 7 楼的前辈已经给您回复了. 参考 (他 | 她) 的代码即可.


   

j_cle

回复 3# sunzhiguolu


    很简单的做法，支持一个