| 查看: 307 | 回复: 1 | |||
[交流]
【求助】急！！！请教高手看一下perl程序里的循环问题？
|
|
为什么my $scafmatch =$allterm[7] - $allterm[6]不循环？之前的程序测试都对！急啊！！谢谢！ #!/usr/bin/perl -w 2 use strict; 3 unless (@ARGV == 3) { 4 die "please check your input files and output files"; 5 } 6 open IN,$ARGV[0] or die "$!\n"; 7 open FILE,$ARGV[1] or die "$!\n"; 8 open OUT,">$ARGV[2]" or die "$!\n"; 9 our @array; 10 $/ = ">"; 11 12 while( 13 chomp; 14 my @all = split /\n/,$_; 15 my $name = shift(@all); 16 my @term = split /\s+/,$name; 17 my $scafname = shift (@term); 18 my $scaf = join "",@all; 19 my @base = split //,$scaf; 20 my $num = 0; 21 for my $i(@base) { 22 if ($i =~/[ATGCN]/) { 23 $num++; 24 } 25 } 26 # print OUT $scafname,"\t",$num,"\n"; 27 my @quary = ($scafname,$num); 28 push @array,[@quary]; 29 } 30 close IN; 31 LABLE: while ( 32 chomp; 33 my @allterm = split /\s+/,$_; 34 my $a = ""; 35 my $rate = 0; 36 next LABLE if ($a eq $allterm[0]); 37 my $scafmatch = $allterm[7] - $allterm[6]; 38 for my $i(0 ..$#array) { 39 if ($array[$i][0] eq $allterm[0] ) { 40 $rate = ( $scafmatch / $array[$i][1]); 41 } 42 print OUT $array[$i][0],"\t",$array[$i][1],"\t","the rate 0f match:\t",$rate,"\n"; 43 $a = $allterm[0]; } 45 } 46 close FILE; 47 close OUT; [ Last edited by nono2009 on 2010-8-19 at 07:07 ] |
回复此楼» 猜你喜欢
332求调剂
已经有5人回复
-
085701求调剂初试286分
已经有3人回复
-
材料工程294求调剂
已经有4人回复
-
本科新能源科学与工程，一志愿华理能动285求调剂
已经有10人回复
-
070305高分子化学与物理 304分求调剂
已经有5人回复
-
083000学硕274求调剂
已经有8人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有10人回复
-
求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过
已经有11人回复
-
071000生物学求调剂，初试成绩343
已经有7人回复
-
0703化学调剂，求导师收
已经有8人回复
2楼2010-08-19 01:30:36
