作者muchmoa (...)
看板Perl
标题[问题] 问题请教
时间Fri Jan 11 23:09:18 2008
最近开始使用Perl来处理研究上的计算，
写了如下的程式:
--------------------------------------------------------------------------
use Cwd;
$dir=cwd;
open inf, "<$dir/RT.txt";
open out1, ">$dir/RT_output1.txt";
open out2, ">$dir/RT_output2.txt";
while(<inf>){
chomp;
@A=($mw, $seq)=split/\s+/;
print out1 "$seq\n";
foreach $seq(@A){
%h = ( W=> 11.0, F=> 10.5, L=> 9.6, I=> 8.4, M=> 5.8,
V=> 5.0, Y=> 4.0, A=> 0.8, T=> 0.4, P=> 0.2,
E=> 0.0, D=> -0.5, C=> -0.8, S=>-0.8, Q=> -0.9,
G=>-0.9, N=>-1.2, R=>-1.3, H=>-1.3, K=>-1.9 );
$sum+=$h{$_} for split//, $seq;
}
print out2 "$seq $sum\n";
}
close inf;
close out1;
close out2;
-------------------------------------------------------------------------
分成out1和out2只是为了检查，
在output2里的数值为何会是从上一列累计的啊？
问题是不是出在foreach那边呢？
另外我想要计算的只有英文字串的部份，
该如何略过RT档案里的第1,4,5列呢？
附上RT.txt:
http://www.badongo.com/file/7263659
烦请前辈们指点一下，
谢谢!!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 203.187.45.149
1F：→ Evenni:你的hash内容是固定,为什麽要放在loop里面呀.. 01/12 16:05
2F：推 Evenni:还有你的$seq在 @A=($mw,$seq)=split(..已经把文字拿出来 01/12 16:13
3F：→ Evenni:foreach $seq(@A) 就把原先的文字结果盖掉了 01/12 16:15
4F：→ muchmoa:我以为foreach $seq(@A)可以取出之前存在@A里的个别$seq 01/13 23:29
5F：→ muchmoa:那我该如何改写呢?谢谢 01/13 23:30
6F：→ Evenni:在foreach $seq(@A) 上面加一行 @A=split(//,$seq); 01/14 01:59
7F：→ muchmoa:这样还是不行耶，他会把每个$seq拆出来的字母混在一起 01/14 10:56
8F：→ Evenni:在@A=($mw, $seq)=split/\s+/;时,其实你已经取得字母部分了 01/14 13:06
9F：→ Evenni:而且字母已经存成$seq, @A=split(//,$seq)是把字母在切开 01/14 13:06
10F：→ Evenni:所以你其他部分都不用改,只要在foreach $seq(@A) 上面 01/14 13:08
11F：→ Evenni:补上@A=split(//,$seq)就好, 补看看,我试过是没问题 01/14 13:08
12F：→ muchmoa:这样算出来的一样是累计的耶，ACTPQK应该是-2.2而不是 01/14 13:22
13F：→ muchmoa:42.3(WFLIMS的44.5-2.2) 01/14 13:23
14F：→ Evenni:因为你印完$sum後没有把$sum归零,所以会累加上去是正常的 01/14 15:20
15F：→ Evenni:在print out2 "$seq $sum\n"; 下面补一行$sum = 0; 01/14 15:21
16F：→ muchmoa:谢谢喔，我以为每个foreach会有自己的$sum,那如何只处理有 01/14 21:11
17F：→ muchmoa:英文字母的列数呢？ 01/14 21:13
18F：→ Evenni:其实你不用烦恼这个问题，因为你是用hash，如果对应不到 01/14 21:52
19F：→ Evenni:字母，自然就不会作加法的动作， 01/14 21:53
20F：→ muchmoa:到也没错，谢谢喔 01/15 17:11