I am currently working on a project with RNA-Seq data. I need to compute a ratio of the number of reads at each position for two different protocols so I have log(protocol1value/protocol2value). The problem is that its quite often I will have something like log(25/0), which is impossible to compute. Is there a way to normalize the values so that I can make this sort of computation accurately?
I want to add that, based on advice by some colleagues, I have considered and dismissed the possibilities of adding 1 to all counts before computation or replacing all 0's with very small numbers like .000001.
I'm sure there must be a way, I am just having trouble finding it. Thanks in advance.
I want to add that, based on advice by some colleagues, I have considered and dismissed the possibilities of adding 1 to all counts before computation or replacing all 0's with very small numbers like .000001.
I'm sure there must be a way, I am just having trouble finding it. Thanks in advance.
Comment