24小时热门版块排行榜    

查看: 433  |  回复: 5

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

[交流] 【求助】perl问题求解!【已完结】

各位大牛们,小弟会一点入门级perl语言,有一个问题想请教大家:
我想知道,perl对变量的模式匹配是如何操作的,例如,我有一个文件,里面数据的格式如下

我有些个数据里面的某一些的数据名称(>后面的那些),在另一个文件中,我想提取出那些匹配的数据名称下面的字符(ttacg之类的),自己弄了个难看的脚本如下
#!/usr/bin/perl

if (@ARGV < 2)
               {
               die "usage: outputseq See script for details.\n";
               }
open(LOCI, "$ARGV[0]" || die "usage: outputseq location file \nCan not open the seq_file $ARGV[0].\nSee script for details.\n";
open(DNASEQ, "$ARGV[1]" || die "usage: genome seq file \nCan not open the seq_file $ARGV[1].\nSee script for details.\n";            

while($seq=
{   @seq=split(/\s+/,$seq);
   $name=$seq[1];
for($info=
   {  
   if($info=~/$name/)
      {
    print "$info";
    $info=;
    print $info;
    last;
      }
   }
}
close LOCI;
close DNASEQ;

结果总是弄不出来,不得已请教各位,希望大家能帮我看看,非常感谢!

[ Last edited by 余泽成 on 2010-8-28 at 22:14 ]
回复此楼

» 猜你喜欢

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)


xp198766(金币+2):谢谢回答! 2010-08-28 12:41:23
余泽成(金币+1):谢谢参与应助! 2010-08-28 21:22:27
$name=$seq[0] ??
2楼2010-08-28 12:31:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

谢谢!不好意思,我的是在第二列,是seq[1]的……我还在那下面print $name过,是对的……
3楼2010-08-28 12:41:03
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

请高手帮忙看下啊,谢谢了!
4楼2010-08-28 15:04:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

高手解答一下吧!
5楼2010-08-28 17:26:35
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

结果已找到,请求关贴~~
6楼2010-08-28 21:33:12
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xp198766 的主题更新
信息提示
请填处理意见