| 查看: 2141 | 回复: 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语言专栏 |
» 猜你喜欢
filecode，4个jtjc了
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有117人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有4人回复
-
今天放榜吗？
已经有13人回复
-
今天基金会出结果吗？20260819
已经有8人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
-
时间戳变了，能看出什么问题？
已经有14人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有30人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 该贴已过期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
