现在需要编一个数理统计的程序，弄不出来了，求大神帮忙，万分万分感谢！！！

旋转小马

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;
文件是我经过处理的，大体如上，先说明一下这个文件，第一列是序列，每两行是一对，第二列是比对上的染色体，第三列中有的有比对上的染色体，就是那些以“Cluster_*****”格式的，现在需要统计的是，两行两行的读取文件，如果第二列中，每两行的染色体是一样的，就统计这条染色体的数目为一对，两条，最后打印的结果为“染色体    对数    条数”格式的，而且还要考虑第三列中的情况，因为第三列中有的也有染色体，同样要统计其对数跟条数，注意：第二列中如果没有比对上的染色体，其显示的是“*”符号；程序对我来讲比较繁琐，求大神帮忙！！！

旋转小马

大神们都哪里去了~~~

pitonas

郁闷。我怎么看来看去也看不懂{:2_170:}

pitonas

lz 应该去学会如何提问

一个简单的例子说明如何提问

http://bbs.chinaunix.net/thread-4111100-1-2.html

旋转小马

之前那个是输入文件，输出文件类似这种格式：
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
第一列就是染色体编号，第二列就是其出现的条数，第三列是其出现的对数，这样讲能明白了吗？？？回复 3# pitonas


   

旋转小马

输入文件为：
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
这样说的话，能明白了吗？？？？？？回复 4# pitonas


   

pitonas

条数，对数
不明白{:2_176:}

旋转小马

好吧
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  
在这四行文件中，每两行表示的是一对序列，四行共两对序列，第二行表示的是染色体， Cluster_29235在前两行中出现的次数计算为1对，2条，后两行中其又出现了1条，但是在后两行中，他没有成对出现，所以最后计算 Cluster_29235的次数的时候，计为出现1对，3条；而 Cluster_67963出现的次数为0对，1条； Cluster_67972出现的次数为1对，2条；这样能明白了吗，亲？？？
回复 7# pitonas


   

pitonas

亲亲？
x Cluster_1 y
x Cluster_2 y
x Cluster_2 y
x Cluster_3 y
这样输出显示？{:2_169:}

旋转小马

最好是这种格式显示：
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
第一列表示的是每条染色体，第二列表示的是条数，第三列表示的是对数。。。。。。
回复 9# pitonas