5'/3' bias edge case bugs
Closed this issue · 1 comments
kachulis commented
A few bugs in calculation of 5'/3' bias, specifically for the edge case where only 1 gene is included in the calculation.
- In calculation of 25th and 75th percentiles, the index to read at is not checked against the size of the ratios vector. If the ratios vector is length 1, this leads to reading from past the end of the vector, giving undefined behavior. Also, the behaviors in both the if and else blocks are currently mathematically identical, I believe.
Lines 500 to 521 in 19967f3
- In the code to compute medians, the edge case when size is 1 hits a similar bug.
Lines 152 to 156 in 19967f3
agraubert commented
Thanks for pointing these out. Good catches!