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

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

Adding 1 would overcome the issues.

ADD REPLYlink written 16 months ago by Hussain Ather910

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 16 months ago by genomax65k
0
gravatar for Hussain Ather
16 months ago by
Hussain Ather910
National Institutes of Health, Bethesda, MD
Hussain Ather910 wrote:

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

ADD COMMENTlink written 16 months ago by Hussain Ather910
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: 777 users visited in the last hour