DoaneAS / rseqc

Automatically exported from code.google.com/p/rseqc
0 stars 0 forks source link

bug in RPKM_saturation.py #4

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago
What steps will reproduce the problem?
1. execute RPKM_saturation.py on a bam file

What is the expected output? What do you see instead?
RPKM saturation table.  

I observe the output, but it is flawed (some lines have blank at 100%; and 
other lines exceed 100%):

for example (missing value at 100%):

chr start   end name    score   strand  5%  10% 15% 20% 25% 30% 35% 40% 45% 50% 55% 60% 6
5%  70% 75% 80% 85% 90% 95% 100%
chr18   46446222    46475703    SMAD7   0   -   1.228135854 4.912543414 6.550057886 6.75474244
1   6.140675811 5.731297961 5.438885164 6.140674946 6.004215973 5.772235262 5.6940
81498   5.526607452 5.762787569415.70205588067  5.895048778 5.910399636 5.996188725
    6.208905154 6.140675629 6.140674946 

and (values exceed 100%):

chr start   end name    score   strand  5%  10% 15% 20% 25% 30% 35% 40% 45% 50% 55% 60% 6
5%  70% 75% 80% 85% 90% 95% 100%                                                                                                                                                        
chr2    189156395   189460652   GULP1   0   +   5.632993606 4.4134327   5.615335319 4.828280234
    3.97208943  4.813144559 6.437706979 5.884576934 6.417526079 6.63888065  6.3994729
12  6.618069112 6.115818186526.00226509349  6.096646343 6.035347461 6.031688528 6.
016427876   6.897540418 6.430999349 6.875918034 6.73946975  6.289137172 6.718342886
    6.348290411 5.982649474 6.328389814 5.954875603 5.649190676 5.936208281035.9256
1361796 5.697337479 5.907038027 5.834167844 5.663901313 5.815878916 5.694890936 
5.431913641 5.677038613 5.805428684 5.485263048 5.787229848 5.686637963 5.413807
731 5.668811512 5.733578743 5.516786993 5.715605142 5.822334073 5.555729249 5.80
4082242 6.035346517 5.786497032 6.016426936 6.141229991 5.900059166 6.121978487 
6.316995527 6.134667137 6.297193033

Original issue reported on code.google.com by zumbo.p...@gmail.com on 24 Apr 2012 at 5:15

GoogleCodeExporter commented 9 years ago

This module is not fully tested currently. We will update it and get back to 
you shortly.

Original comment by wangligu...@gmail.com on 24 Apr 2012 at 6:11