 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 312 | 回復(fù): 1 | |||
[交流]
【求助】急�。�！請教高手看一下perl程序里的循環(huán)問題？
|
|
為什么my $scafmatch =$allterm[7] - $allterm[6]不循環(huán)？之前的程序測試都對！急啊��！謝謝！ #!/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 ] |
回復(fù)此樓» 猜你喜歡
309分085801求調(diào)劑
已經(jīng)有7人回復(fù)
-
生物學(xué)308分求調(diào)劑（一志愿華東師大）
已經(jīng)有3人回復(fù)
-
求調(diào)劑：一志愿：南京大學(xué) 專業(yè)：0705 總分320 ，本科985，四六級已過
已經(jīng)有3人回復(fù)
-
環(huán)境工程 085701，267求調(diào)劑
已經(jīng)有14人回復(fù)
-
308求調(diào)劑
已經(jīng)有12人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有21人回復(fù)
-
311求調(diào)劑
已經(jīng)有8人回復(fù)
-
294分080500材料科學(xué)與工程求調(diào)劑
已經(jīng)有13人回復(fù)
-
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有19人回復(fù)
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有23人回復(fù)
★
nono2009(金幣+1):鼓勵應(yīng)助。 2010-08-19 07:07:26
shuangyang(金幣+2): 2010-08-19 09:09:48
nono2009(金幣+1):鼓勵應(yīng)助。 2010-08-19 07:07:26
shuangyang(金幣+2): 2010-08-19 09:09:48
|
有個錯誤蠻明顯的。 你在while（ 不知這個和你說的那個錯誤有沒有關(guān)系。 ps：你貌似是先處理一個fasta文件，那個計數(shù)的循環(huán)是不是可以用length（）函數(shù)代替，可以的話速度會提高很多。 |
2樓2010-08-19 01:30:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 shuangyang 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)308分求調(diào)劑（一志愿華東師大） +3 | 相信必會光芒萬?/a> 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[電化學(xué)] 070300化學(xué)調(diào)劑 +7 | 山頂見α 2026-03-25 | 7/350 |
|
|
|
[考研] 349求調(diào)劑 +4 | zwjjjjjj 2026-03-31 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-29 | 10/500 |
|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +13 | 深V宿舍吧 2026-03-29 | 13/650 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +9 | 手機(jī)用戶 2026-03-26 | 9/450 |
|
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +6 | 鐵憨憨123425 2026-03-31 | 6/300 |
|
|
|
[考研] 08開頭看過來�。。� +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
|
[考研] 340求調(diào)劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +7 | 朵朵- 2026-03-26 | 9/450 |
|
|
|
[考研] 英一數(shù)一總分334求調(diào)劑 +4 | 陳陽坤 2026-03-31 | 4/200 |
|
|
|
[考研] 調(diào)劑310 +13 | 溫柔的晚安 2026-03-25 | 14/700 |
|
|
|
[考研] 一志愿大連理工大學(xué)，機(jī)械工程學(xué)碩，341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
|
[考研] 085600材料與化工329分求調(diào)劑 +5 | 這是有太陽哇 2026-03-25 | 5/250 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 312，生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|