| 查看: 417 | 回复: 2 | ||
[求助]
用perl 统计dna三核苷酸频数
|
|
我的文件是一行序列信息，一行序列。信息以>开头，接下来是核苷酸序列。程序如下，可计算结果前几条都没有问题，后面的计算结果就错了，请问这个程序问题在哪？谢谢！ #!/usr/bin/perl-w use strict; my @n =qw(a c g t); my ($a, $frag, @frag, $f, $i, $j, $k, $r, $trint, %count); open (MSET, "F:\\m.txt"  ;open (OUT, ">F:\\64m.txt" ;while ( chomp; if(/^>/){next;} for ($a=0; $a<=(length($_)-3); $a++){ $frag[$a]= substr($_, $a, 3); # print "$frag[$a]\t"; } #print "\n"; for ($i=0; $i<4; $i++){ for ($j=0; $j<4; $j++){ for ($k=0; $k<4; $k++){ $trint = $n[$i].$n[$j].$n[$k]; #print "$trint\n"; $count{$trint}=0; foreach $frag(@frag){ if ($trint =~/$frag/){ $count{$trint} +=1; } # print "$trint is $count{$trint} \n"; } } } } foreach $f(sort keys %count){ print OUT "$f\t"; } print OUT "\n"; foreach $f(sort keys %count){ print OUT "$count{$f}\t"; } print OUT "\n"; } close MSET; close OUT; |
回复此楼» 猜你喜欢
315分求调剂
已经有3人回复
-
材料与化工304求B区调剂
已经有3人回复
-
一志愿天津大学339材料与化工求调剂
已经有3人回复
-
297求调剂
已经有3人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有3人回复
-
打过很多竞赛，085406控制工程300分，求调剂
已经有3人回复
-
北京林业大学硕导招生广告
已经有5人回复
-
一志愿北京化工大学材料与化工 264分各科过A区国家线
已经有4人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有4人回复
-
274求调剂
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 国家自然科学基金资助项目统计资料公布
已经有24人回复
- DNA浓度
已经有13人回复
- 求助高手，我用手提的土壤DNA，用什么方法纯化
已经有11人回复
- DNA 条带 和 OD值
已经有13人回复
- 国家自然科学基金资助项目统计资料
已经有16人回复
- 知道氨基酸序列，如何得知核苷酸序列
已经有5人回复
- 定点突变最多能够突变多少相连的核苷酸残基?
已经有6人回复
- 控制理论中多智能体系统自然科学基金项目统计及分析
已经有145人回复
- dna各种提取方法简介
已经有104人回复
- DNAssist 很好用的DNA分析软件
已经有7人回复
- 有木有人知道“分离相似分子量的DNA分子的方法＂呢？
已经有5人回复
- DNA的连接反应
已经有5人回复
- 【求助/交流】植物DNA入境
已经有3人回复
2楼2011-12-26 15:37:09
3楼2011-12-26 15:40:54
