我寫的代碼最后的輸出結(jié)果有問題，不知道出在哪里，求大神幫忙指教一下，幫我改一下，

旋轉(zhuǎn)小馬

infile：
HWI-ST1106:755:H0V1MADXX:2:2204:4345:52798      99      Cluster_1   
HWI-ST1106:755:H0V1MADXX:2:2214:2460:68727      65      Cluster_1      
HWI-ST1106:755:H0V1MADXX:2:1205:4057:41509      153     Cluster_1   
HWI-ST1106:755:H0V1MADXX:2:2214:5236:82377      153     Cluster_1   
HWI-ST1106:755:H0V1MADXX:2:1202:16871:54180     89      Cluster_1   
HWI-ST1106:755:H0V1MADXX:2:1212:4117:41420      153     Cluster_1      
HWI-ST1106:755:H0V1MADXX:2:2204:4345:52798      147     Cluster_1     
HWI-ST1106:755:H0V1MADXX:2:2115:10141:17608     89      Cluster_10   
HWI-ST1106:755:H0V1MADXX:2:1212:17055:53562     163     Cluster_10   
HWI-ST1106:755:H0V1MADXX:2:2110:14814:19836     163     Cluster_10
HWI-ST1106:755:H0V1MADXX:2:2206:18216:74573     99      Cluster_10
HWI-ST1106:755:H0V1MADXX:2:2113:17626:93067     163     Cluster_10   
HWI-ST1106:755:H0V1MADXX:2:1207:3103:73263      99      Cluster_10   
HWI-ST1106:755:H0V1MADXX:2:1103:7010:88764      163     Cluster_10     
HWI-ST1106:755:H0V1MADXX:2:2111:3039:51639      99      Cluster_10   
HWI-ST1106:755:H0V1MADXX:2:1212:17055:53562     83      Cluster_10     
HWI-ST1106:755:H0V1MADXX:2:1107:5037:38527      163     Cluster_100  
HWI-ST1106:755:H0V1MADXX:2:1111:6627:73408      393     Cluster_100   
HWI-ST1106:755:H0V1MADXX:2:1215:16841:58015     73      Cluster_100  
HWI-ST1106:755:H0V1MADXX:2:1107:5037:38527      83      Cluster_100  
HWI-ST1106:755:H0V1MADXX:2:2203:2962:64393      393     Cluster_100  
HWI-ST1106:755:H0V1MADXX:2:1113:7244:73782      137     Cluster_100   
HWI-ST1106:755:H0V1MADXX:2:1104:14953:2121      163     Cluster_1000  
HWI-ST1106:755:H0V1MADXX:2:1212:12709:69650     153     Cluster_1000  
HWI-ST1106:755:H0V1MADXX:2:1108:11227:4178      99      Cluster_1000  
HWI-ST1106:755:H0V1MADXX:2:2112:18155:78922     99      Cluster_1000  
HWI-ST1106:755:H0V1MADXX:2:1203:2053:36493      99      Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:1112:20266:62035     163     Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:2104:3375:35332      89      Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:1105:17103:16844     153     Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:1215:2568:63103      89      Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:2205:3921:94905      129     Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:2107:17892:33129     89      Cluster_1000
HWI-ST1106:755:H0V1MADXX:2:2112:18155:78922     147     Cluster_1000       -
HWI-ST1106:755:H0V1MADXX:2:1104:14953:2121      83      Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:1108:11227:4178      147     Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:1112:20266:62035     83      Cluster_1000   
HWI-ST1106:755:H0V1MADXX:2:1203:2053:36493      147     Cluster_1000  
下面是我寫的代碼：
1 #!/usr/bin/perl
  2 use strict;
  3 use warnings;
  4 open IN,"hitsok.sam"or die"Can not open IN\n";
  5 open OUT,">oka.xls"or die"Can not open OUT\n";
  6 my (%hash,%chek,%pare);
  7 while(<IN>){
  8         chomp;
  9         my @line=split/\t/,$_;
10         if(exists $hash{$line[2]}){
11            $hash{$line[2]}++;
12            if(exists $chek{$line[0]}){
13               $pare{$line[2]}++;
14            }else{
15               $chek{$line[0]}=1;
16            }
17         }else{
18            $hash{$line[2]}=1;
19         }
20 }
21 foreach my $key (sort {$hash{$b}<=>$hash{$a}} keys %hash){
22         print OUT "$key\t$hash{$key}\t";
23         if(exists $pare{$key}){
24            print OUT "$pare{$key}\n";
25         }else{
26            print OUT "0\n";
27         }
28 }
29 close IN;
30 close OUT;  
正確的輸出結(jié)果應(yīng)該是這樣的：
Cluster_1000    16      4
Cluster_10      9       1
Cluster_1       7       1
Cluster_100     6       1
但是我寫的代碼的輸出結(jié)果是這樣的：
Cluster_1000    16      4
Cluster_10      9       1
Cluster_1       7       0
Cluster_100     6       0
求各位大神、各位高手幫忙找一下錯誤，萬分感謝！�。�

pitonas

lz 應(yīng)該去學(xué)會如何提問{:2_179:}

一個簡單的例子說明如何提問

http://www.72891.cn/thread-4111100-1-2.html

旋轉(zhuǎn)小馬

好吧，我現(xiàn)在嚴(yán)重懷疑自己的表達能力了，你哪里不明白，我單獨說那個地方？
回復(fù) 2# pitonas


   

jason680

回復(fù) 1# 旋轉(zhuǎn)小馬

Add some code for debug, and didn't find  Cluster_1 and Cluster_100

          if(exists $chek{$line[0]}){
              $pare{$line[2]}++;
              print "\$pare{\$line[2]}=$pare{$line[2]}, \$line[2]=$line[2]\n";
          }else{


$pare{$line[2]}=1, $line[2]=Cluster_10
$pare{$line[2]}=1, $line[2]=Cluster_1000
$pare{$line[2]}=2, $line[2]=Cluster_1000
$pare{$line[2]}=3, $line[2]=Cluster_1000
$pare{$line[2]}=4, $line[2]=Cluster_1000
   

pitonas

如何得到輸出結(jié)果：
Cluster_1000    16      4
Cluster_10      9       1
Cluster_1       7       1
Cluster_100     6       1

旋轉(zhuǎn)小馬

回復(fù) 4# jason680
你說的雖然不是全明白，但也明白了一點點，但是怎樣修改才能得到我想要的結(jié)果呢？？？

   

旋轉(zhuǎn)小馬

回復(fù) 5# pitonas
問題就在這嘛，我寫的代碼就是想要得到這個結(jié)果，但是為什么程序輸出的結(jié)果不對，現(xiàn)在就是不知道原因在哪啊？

   

jason680

回復(fù) 6# 旋轉(zhuǎn)小馬


would you help us to help you?

you didn't describe how to get result and I don't know how to modify.

q1208c

回復(fù) 7# 旋轉(zhuǎn)小馬

你至少得告訴我們, 你的這個結(jié)果是怎么來的吧? 靠大家來猜, 好象對解決問題沒什么幫助.

簡單的說吧, 你的 正確 輸出中, 第一字段, 我可以理解, 就是輸入的最后一個字段, 第二個字段, 我猜是 輸入最后一個字段的計數(shù), 但你輸出的第三個字段, 我實在猜不出來了, 也沒看出跟輸入有半毛錢關(guān)系. 那么, 你讓我們怎么幫你呢?

生成這樣的輸出很容易的, 一條 print就行了.  

1. 
   
2. print qq(
   
3. Cluster_1000    16      4
   
4. Cluster_10      9       1
   
5. Cluster_1       7       1
   
6. Cluster_100     6       1
   
7. );
   

復(fù)制代碼

只是一個紅薯

10         if(exists $hash{$line[2]}){
11            $hash{$line[2]}++;
12            if(exists $chek{$line[0]}){
13               $pare{$line[2]}++;
14            }else{
15               $chek{$line[0]}=1;
16            }
17         }
第二個if寫在第一個if里面是有問題的，導(dǎo)致第一個$line[0]統(tǒng)計不到