 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 310 | 回復(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ù)此樓» 猜你喜歡
一志愿西電085401數(shù)一英一299求調(diào)劑 六級521
已經(jīng)有3人回復(fù)
-
生物考研337分求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有14人回復(fù)
-
285求調(diào)劑
已經(jīng)有10人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有29人回復(fù)
-
理學(xué)07化學(xué) 303求調(diào)劑
已經(jīng)有5人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
土木304求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國科學(xué)院大學(xué)265求調(diào)劑
已經(jīng)有7人回復(fù)
-
286求調(diào)劑
已經(jīng)有8人回復(fù)
caucast
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 21.1
- 帖子: 4
- 在線: 2.5小時
- 蟲號: 1068313
- 注冊: 2010-08-02
- 性別: GG
- 專業(yè): 生物信息
★
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ā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +4 | zzll03 2026-03-31 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085600材料與化工329分求調(diào)劑 +5 | 這是有太陽哇 2026-03-25 | 5/250 |
|
|
|
[考研] 323分 食品與營養(yǎng)調(diào)劑 +3 | 嘿ooo 2026-03-31 | 3/150 |
|
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 12/600 |
|
|
|
[考研]
 304材料求調(diào)劑
+5
|
鐘llll 2026-03-26 | 5/250 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7 | nhj_ 2026-03-25 | 7/350 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
|
[考研] 0856求調(diào)劑 +8 | 楒桉 2026-03-28 | 8/400 |
|
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 機械學(xué)碩總分317求調(diào)劑！�。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-25 | 7/350 |
|