Question: RNASeq data RPKM
0
gravatar for tanyabioinfo
2.6 years ago by
tanyabioinfo0 wrote:

Hi

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.

Best

Tanya

rna-seq • 1.2k views
ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by tanyabioinfo0

Yes adding 1 to FPKMS is a good way to take care of this. The fold change of B/A is fine but you have to take log of the fold change.

ADD REPLYlink written 2.6 years ago by Ron990

Hi Hussain

Thanks for your reply.

I have this data at three time point with the ligand binding at the second and third time point. I want to see the effect of the ligand on gene expression. I am looking at overall pattern of the gene so I need to consider these cases. Do you think adding 1 to the entire dataset a good way to overcome such issues.

Best

ADD REPLYlink written 2.6 years ago by tanyabioinfo0

Adding 1 would overcome the issues.

ADD REPLYlink written 2.6 years ago by Hussain Ather940

tanyabioinfo : Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized.

This comment should have gone under @Hussain's answer.

ADD REPLYlink written 2.6 years ago by genomax85k
0
gravatar for Hussain Ather
2.6 years ago by
Hussain Ather940
National Institutes of Health, Bethesda, MD
Hussain Ather940 wrote:

You can't compute fold change if the gene wasn't detected. Ignore them. The fold change is B/A.

ADD COMMENTlink written 2.6 years ago by Hussain Ather940
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1060 users visited in the last hour