Question: RNASeq data RPKM
0
gravatar for tanyabioinfo
22 months 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 • 976 views
ADD COMMENTlink modified 22 months ago • written 22 months 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 22 months ago by Ron970

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 22 months ago by tanyabioinfo0

Adding 1 would overcome the issues.

ADD REPLYlink written 22 months 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 22 months ago by genomax73k
0
gravatar for Hussain Ather
22 months 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 22 months 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: 1655 users visited in the last hour