SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Log fold change Erwan Bioinformatics 3 03-25-2014 05:05 AM
Standardized Fold Change dbaud Introductions 5 01-27-2014 10:47 AM
fold change value-cuffdiff madsaan Bioinformatics 4 02-10-2011 06:51 AM
fold change for genes with 0 reads rnaseq RNA Sequencing 15 11-05-2010 01:23 AM
fold change for genes with 0 reads rnaseq General 2 11-04-2010 07:44 PM

Reply
 
Thread Tools
Old 04-22-2014, 01:47 PM   #1
ppm
Junior Member
 
Location: USA

Join Date: Mar 2014
Posts: 3
Default RPKM vs Fold Change

Hi,

I have been trying to make TSS plots. I see that the RPKM values for both the input and enriched ChIPseq samples are close while when I plot them for fold change, the global enrichment in one of the samples around TSS is more than the other. Also, I have two different cell lines, therefore, two different inputs as well.

THank you,
ppm is offline   Reply With Quote
Old 04-28-2014, 06:50 AM   #2
ppm
Junior Member
 
Location: USA

Join Date: Mar 2014
Posts: 3
Default

Would anyone know what could be the reason for this?
ppm is offline   Reply With Quote
Old 04-28-2014, 09:48 AM   #3
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

My guess would be that it has something to do with how you're calculating RPKM. Is the "M" from the total mapped reads or just the total mapped reads in the regions that you're looking at?
dpryan is offline   Reply With Quote
Old 04-28-2014, 02:48 PM   #4
jparsons
Member
 
Location: SF Bay Area

Join Date: Feb 2012
Posts: 62
Default

Why are you using RPKM here?

I'm not too familiar with your study, but I'd think that you want to be looking at windows of very specific width (around the TSSs, +/- N nt) and therefore will be comparing sequences of identical length. No need to length normalize when everything is the same length.
jparsons is offline   Reply With Quote
Old 05-06-2014, 07:49 AM   #5
ppm
Junior Member
 
Location: USA

Join Date: Mar 2014
Posts: 3
Default

So, I have ChIPseq data for two cell-lines, therefore, two inputs.

I am plotting TSS plots +/-1500 bp. I am using ngsplot to do this.

When I just use enriched bam file or total input bam file, it gives me RPKM plotted on the Y-axis. These look very similar in case of both the cell lines.

Next, when I plot enrichedvstotalInput for these cell lines. I see that one of the cell line is more enriched than the other. I am not sure why would this happen. Views will be much appreciated.

Thanks
ppm is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:20 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO