I have RNA Seq (RPKM) data. I am trying to look for fold change in the gene expression for some genes.
Say the RPKM value in condition1 is A and for condition 2 is B
For few genes I have A as zero values. So when I calculate fold change(B/A). For the genes which has A as zero the value would become infinite. Question1: How can the infinity value be taken care of. Is Adding 1 to all the data a good way to handle it, specially when I am interested in fold change. Question2: Is the fold change as B/A correct of doing it.