| 查看: 1914 | 回复: 3 | |||
csdyg木虫 (正式写手)
|
[求助]
Perl 生物信息学编程求助，运行出现问题，烦请帮忙查找一下错误，谢谢！
|
|
下面是我用Perl 写的一个有关统计成对比对氨基酸序列中氨基酸取代的一个程序。运行时从第一个序列文件中提取出双序列进行比对后统计氨基酸取代，然后进行读取第二个文件中的氨基酸序列进行比对再统计总的氨基酸取代，这样依次进行。 但是目前的程序在读取第二个文件中的序列后得到的统计结果不是加上第一个文件中的结果，而是出现一些莫名其妙的增减。 烦请高手帮忙检查错误，非常感谢！ #!/usr/bin/perl -w use SOAP::Lite; my($service) = SOAP::Lite->proxy('http://www.ebi.ac.uk/soap$ if (@ARGV!=1){ print "The file is missing\n"; exit; } $filename=$ARGV[0]; if (open(E,$filename)){ while($file= chomp $file; if (open(F,$file)){ $seq=""; while($line= $seq=$seq.$line; } } @a=split(/\n\n/,$seq ); $seq1=$a[0]; $seq2=$a[1]; my($soap_response)= $service->runAndWaitFor( SOAP:: Data->type(map=>{ asequence_direct_data => $seq1, bsequence_direct_data => $seq2 })); my(@res) = $soap_response->paramsout; my($res) = $soap_response->result; my($outfile) = %$res->{'outfile'}; @b=split(/#[=-]+/,$outfile); $align=$b[2]; @c=split(/ [0-9]+/,$align); $i=0; $n=@c; $s1=""; $s2=""; while ($i<$n){ if ($i%4==1){ $s1=$s1.$c[$i]; } if ($i%4==3){ $s2=$s2.$c[$i]; } $i++; } $s1=~ s/\s*//g; $s2=~ s/\s*//g; $n=length $s1; for($i=0; $i<$n; $i++){ $a1=substr($s1,$i,1); $a2=substr($s2,$i,1); if (($a1 ne "-" ) and ($a2 ne "-" )){ $p="$a1,$a2"; $count{$p}++; } } close F; } } close E; foreach $key (keys %count){ print "$key $count{$key}\n"; } [ Last edited by jjdg on 2012-3-6 at 02:22 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 perl语言专栏 |
» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 该贴已过期5
已经有11人回复
- 求助高人帮忙分析一下污水厂运行情况
已经有15人回复
- 求助：word2003公式编辑器
已经有8人回复
- 2012整合生物信息学国际论坛 浙江大学
已经有5人回复
- 小白一枚 请各位大侠帮帮忙吧~~
已经有10人回复
- 请教Windows系统下perl脚本的ARGV应用问题，谢谢
已经有5人回复
- c#高级编程全本下载！
已经有84人回复
- 生物信息学 找工作
已经有20人回复
- 求助！win7下运行vc6.0时出现msdev错误怎么办？？？
已经有1人回复
- 中国科学院青岛生物能源与过程研究所招聘全职/兼职生物信息学研究人员
已经有10人回复
- 【求助】请教生物信息学方面的专家，一个进化树构建的问题？
已经有10人回复
- 【求助】程序运行已经无错误，但是输入文件有问题，求教。
已经有23人回复
csdyg
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 3737.5
- 散金: 277
- 红花: 9
- 帖子: 422
- 在线: 529.6小时
- 虫号: 94520
- 注册: 2005-09-27
- 专业: 生物大分子结构与功能
2楼2012-03-05 20:45:43
sharon0812
金虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 1506.4
- 红花: 1
- 帖子: 50
- 在线: 21.4小时
- 虫号: 1157787
- 注册: 2010-11-28
- 专业: 药物分析
3楼2012-04-10 11:32:16
4楼2012-04-26 17:41:16
