求助一程序

bioinfor

有两个文件，分别是list.txt与input.txt。

list.txt放的文件如下（部分）：
CLOCK        CLOCK        9575
PPARD        PPARD        5467
DSG4        DSG4        147409
PTGER3        PTGER3        5733
NUP35        NUP35        129401

文件共三列，第一和二列相同，第三列对应相应的数字


input.txt放的文件如下（部分）：
147409
9575
129401


input.txt只有一列，如果input.txt里的数字与list.txt里的数字相同则输出list.txt 文件里相应的值并保存到result.txt文件中


如上面的列子应该输出
DSG4
CLOCK
NUP35



我把list.txt，input.txt，以及两者应该得到的结果以附件形式列上去。

lkk2003rty

发错地方了 sorry  编辑掉
不水贴  写个awk
awk 'NR==FNR{a[$3]=$1;next}{if($0 in a)printa[$0]}' list.txt input.txt

x9x9

1. use strict;
   
2. use warnings;
   
3. 
   
4. open LIS,"<list.txt" or die "can't open file:$!";
   
5. open INP,"<input.txt" or die "can't open file:$!";
   
6. open RES,">result.txt" or die "can't creat file:$!";
   
7. 
   
8. chomp(my @arr=<INP>);
   
9. close INP;
   
10. 
    
11. 
    
12. while (<LIS>){
    
13. chomp;
    
14. my ($x,$y,$z) =split /\t/,$_;
    
15. foreach my $list (@arr){
    
16. print RES $x."\n" if ($list == $z);
    
17. }
    
18. }
    
19. 
    
20. close LIS;
    
21. close RES;
    

复制代码

bioinfor

发错地方了 sorry  编辑掉
不水贴  写个awk
awk 'NR==FNR{a[$3]=$1;next}{if($0 in a)printa[$0]}' list.txt input.txt

awk不太会用，运行了一下没有输出内容，不知是怎么回事？


x9x9战友的程序运行了一下附件的test，出来的结果只有一条，应该是55条，不知是什么问题

lkk2003rty

回复 4# bioinfor


    貌似俺俺错键了 漏了个空格

1. awk 'NR==FNR{a[$3]=$1;next}{if($0 in a)print a[$0]}' list.txt input.txt

复制代码

toniz

1. perl -ne '$tmp=$_;if($ARGV eq "list.txt"){push @tmp,$tmp}else{for(@tmp){print $1."\n" if(/(\S+)\s+$tmp/)}} ;' list.txt input.txt
   

复制代码

x9x9

回复 4# bioinfor


    修改过了，应该可以了。你的rusult.txt 里是按字母排序的，我的没做排序。

bioinfor

运行了一下，程序都很好用，谢谢！学习了

liyangole

回复 6# toniz


    这个不错就是初学者很难理解。