現(xiàn)在需要編一個(gè)數(shù)理統(tǒng)計(jì)的程序，弄不出來了，求大神幫忙，萬分萬分感謝�。�！

旋轉(zhuǎn)小馬

HWI-ST1106:755:H0V1MADXX:2:1101:1221:2186        Cluster_29235        XA:Z:Cluster_67972,+589,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1221:2186        Cluster_29235        XA:Z:Cluster_67972,-771,100M,1;
HWI-ST1106:755:H0V1MADXX:2:1101:1270:2135        Cluster_39664       
HWI-ST1106:755:H0V1MADXX:2:1101:1270:2135        Cluster_39664       
HWI-ST1106:755:H0V1MADXX:2:1101:1483:2147        Cluster_33745       
HWI-ST1106:755:H0V1MADXX:2:1101:1483:2147        Cluster_67963       
HWI-ST1106:755:H0V1MADXX:2:1101:1318:2150        Cluster_79602        XA:Z:Cluster_42584,-594,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1318:2150        Cluster_79602        XA:Z:Cluster_42584,+409,100M,1;
HWI-ST1106:755:H0V1MADXX:2:1101:1383:2209        Cluster_18421       
HWI-ST1106:755:H0V1MADXX:2:1101:1383:2209        Cluster_18421               
HWI-ST1106:755:H0V1MADXX:2:1101:1676:2131        Cluster_69157        XA:Z:Cluster_29478,+828,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1676:2131        Cluster_69157        XA:Z:Cluster_29478,-1069,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1543:2146        Cluster_41540       
HWI-ST1106:755:H0V1MADXX:2:1101:1543:2146        Cluster_41540       
HWI-ST1106:755:H0V1MADXX:2:1101:1725:2175        Cluster_40792       
HWI-ST1106:755:H0V1MADXX:2:1101:1725:2175        Cluster_40789        XA:Z:Cluster_40790,-372,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1516:2211        Cluster_6748        XA:Z:Cluster_54877,-968,100M,3;
HWI-ST1106:755:H0V1MADXX:2:1101:1516:2211        Cluster_6748        XA:Z:Cluster_54877,+711,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1902:2088        Cluster_6876       
HWI-ST1106:755:H0V1MADXX:2:1101:1902:2088        Cluster_6876       
HWI-ST1106:755:H0V1MADXX:2:1101:1883:2112        *       
HWI-ST1106:755:H0V1MADXX:2:1101:1883:2112        *       
HWI-ST1106:755:H0V1MADXX:2:1101:1921:2112        Cluster_42544        XA:Z:Cluster_78702,-3815,24M3D76M,3;Cluster_78704,-1041,2
4M3D76M,3;Cluster_78703,-3794,24M3D76M,3;
HWI-ST1106:755:H0V1MADXX:2:1101:1921:2112        Cluster_42544        XA:Z:Cluster_78702,+3629,100M,0;Cluster_78703,+3608,100M,
0;Cluster_78704,+855,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1758:2115        Cluster_66929        XA:Z:Cluster_24977,+1407,100M,1;
HWI-ST1106:755:H0V1MADXX:2:1101:1758:2115        Cluster_66929        XA:Z:Cluster_24977,-1585,100M,0;
文件是我經(jīng)過處理的，大體如上，先說明一下這個(gè)文件，第一列是序列，每?jī)尚惺且粚?duì)，第二列是比對(duì)上的染色體，第三列中有的有比對(duì)上的染色體，就是那些以“Cluster_*****”格式的，現(xiàn)在需要統(tǒng)計(jì)的是，兩行兩行的讀取文件，如果第二列中，每?jī)尚械娜旧�體是一樣的，就統(tǒng)計(jì)這條染色體的數(shù)目為一對(duì)，兩條，最后打印的結(jié)果為“染色體    對(duì)數(shù)    條數(shù)”格式的，而且還要考慮第三列中的情況，因?yàn)榈谌�列中有的也有染色體，同樣要統(tǒng)計(jì)其對(duì)數(shù)跟條數(shù)，注意：第二列中如果沒有比對(duì)上的染色體，其顯示的是“*”符號(hào)；程序?qū)ξ襾碇v比較繁瑣，求大神幫忙�。�！

旋轉(zhuǎn)小馬

大神們都哪里去了~~~

pitonas

郁悶。我怎么看來看去也看不懂{:2_170:}

pitonas

lz 應(yīng)該去學(xué)會(huì)如何提問

一個(gè)簡(jiǎn)單的例子說明如何提問

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

旋轉(zhuǎn)小馬

之前那個(gè)是輸入文件，輸出文件類似這種格式：
Cluster_61957   1647100 818718
Cluster_9141    1614933 802646
Cluster_30722   907721  451421
Cluster_69535   881255  438249
Cluster_26071   238105  118715
Cluster_40680   225737  110852
Cluster_42542   222271  110920
Cluster_42541   214644  107104
Cluster_6788    181676  90515
Cluster_38806   175830  74981
Cluster_42543   165019  81188
Cluster_40679   164901  81843
第一列就是染色體編號(hào)，第二列就是其出現(xiàn)的條數(shù)，第三列是其出現(xiàn)的對(duì)數(shù)，這樣講能明白了嗎？？？回復(fù) 3# pitonas


   

旋轉(zhuǎn)小馬

輸入文件為：
HWI-ST1106:755:H0V1MADXX:2:1101:1221:2186        Cluster_29235        XA:Z:Cluster_67972,+589,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1221:2186        Cluster_29235        XA:Z:Cluster_67972,-771,100M,1;
HWI-ST1106:755:H0V1MADXX:2:1101:1270:2135        Cluster_39664        
HWI-ST1106:755:H0V1MADXX:2:1101:1270:2135        Cluster_39664
HWI-ST1106:755:H0V1MADXX:2:1101:1883:2112        *        
HWI-ST1106:755:H0V1MADXX:2:1101:1883:2112        *        
HWI-ST1106:755:H0V1MADXX:2:1101:1483:2147        Cluster_29235        
HWI-ST1106:755:H0V1MADXX:2:1101:1483:2147        Cluster_67963               
HWI-ST1106:755:H0V1MADXX:2:1101:1921:2112        Cluster_42544        XA:Z:Cluster_78702,-3815,24M3D76M,3;Cluster_78704,-1041,24M3D76M,3;Cluster_78703,-3794,24M3D76M,3;
HWI-ST1106:755:H0V1MADXX:2:1101:1921:2112        Cluster_42544        XA:Z:Cluster_78702,+3629,100M,0;Cluster_78703,+3608,100M,0;Cluster_78704,+855,100M,0;
最后的輸出文件類似這樣格式：
Cluster_29235     3    1
Cluster_39664     2    1
Cluster_67972     2    1
Cluster_67963     1    0
Cluster_42544     2    1
Cluster_78702     2    1
Cluster_78703     2    1
Cluster_78704     2    1
這樣說的話，能明白了嗎？？？？？？回復(fù) 4# pitonas


   

pitonas

條數(shù)，對(duì)數(shù)
不明白{:2_176:}

旋轉(zhuǎn)小馬

好吧
HWI-ST1106:755:H0V1MADXX:2:1101:1221:2186        Cluster_29235        XA:Z:Cluster_67972,+589,100M,0;
HWI-ST1106:755:H0V1MADXX:2:1101:1221:2186        Cluster_29235        XA:Z:Cluster_67972,-771,100M,1;
HWI-ST1106:755:H0V1MADXX:2:1101:1483:2147        Cluster_29235        
HWI-ST1106:755:H0V1MADXX:2:1101:1483:2147        Cluster_67963  
在這四行文件中，每?jī)尚斜硎镜氖且粚?duì)序列，四行共兩對(duì)序列，第二行表示的是染色體， Cluster_29235在前兩行中出現(xiàn)的次數(shù)計(jì)算為1對(duì)，2條，后兩行中其又出現(xiàn)了1條，但是在后兩行中，他沒有成對(duì)出現(xiàn)，所以最后計(jì)算 Cluster_29235的次數(shù)的時(shí)候，計(jì)為出現(xiàn)1對(duì)，3條；而 Cluster_67963出現(xiàn)的次數(shù)為0對(duì)，1條； Cluster_67972出現(xiàn)的次數(shù)為1對(duì)，2條；這樣能明白了嗎，親？？？
回復(fù) 7# pitonas


   

pitonas

親親？
x Cluster_1 y
x Cluster_2 y
x Cluster_2 y
x Cluster_3 y
這樣輸出顯示？{:2_169:}

旋轉(zhuǎn)小馬

最好是這種格式顯示：
Cluster_29235     3    1
Cluster_39664     2    1
Cluster_67972     2    1
Cluster_67963     1    0
Cluster_42544     2    1
Cluster_78702     2    1
Cluster_78703     2    1
Cluster_78704     2    1
第一列表示的是每條染色體，第二列表示的是條數(shù)，第三列表示的是對(duì)數(shù)。。。。。。
回復(fù) 9# pitonas