SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
What is the meaning density in CummeRbund plots? Carlos Borroto Bioinformatics 16 12-29-2016 09:05 PM
volcano plots cummeRbund godzilla07 Bioinformatics 8 02-04-2015 12:55 PM
Help with how to interpret density plot (cummeRbund) DonDolowy RNA Sequencing 4 07-19-2013 02:33 AM
CummeRbund plots godzilla07 Bioinformatics 1 08-10-2012 04:00 AM
tag density plots bogdan RNA Sequencing 0 12-03-2010 10:44 PM

Reply
 
Thread Tools
Old 03-09-2013, 02:08 PM   #1
DrS
Member
 
Location: North Carolina

Join Date: Mar 2013
Posts: 17
Question Strange Dispersion and Density Plots in cummeRbund

Usually by the time I get to cummeRbund (having run TopHat, cufflinks, cuffmerge and cuffdiff) I see something in the dispersionPlot that looks at least vaguely diagonal, and in the csDensity plot that looks at least semi-normal.

I've just done a study comparing cells exposed to media only and to media + a chemical compound. The idea is to look for genes that are differentially expressed when exposed to the chemical compound.

There were 3 replicates for each (control and treatment).
I've run TopHat with the gtf file for hg19 (these are human cells), then cufflinks, cuffmerge, cuffdiff (in mostly the standard ways) and got plots that look like nothing I've seen before.

I re-rand cuffdiff with an upper-diagnol normalization, but honestly if these plots show fpkm I don't think that would have changed anything.

My csScatter is pretty well centered on the diagonal, so I think normalization worked well. csBoxplots, with and without replicates, also look like what I usually see, though the IQRs are a little longer than usual.

So what might be going on here, and should I be worried about my results?
What's causing that "tail" in the dispersion plot?

Any advice or insights would be greatly appreciated!

Hopefully these attached images work:


DrS is offline   Reply With Quote
Old 03-10-2013, 06:49 AM   #2
DrS
Member
 
Location: North Carolina

Join Date: Mar 2013
Posts: 17
Default

An Update:

I re-ran the analysis without using the gtf file for hg19 and got a dispersion plot that looks more like what I would usually see. My density plot is now more clearly bimodal as well, which isn't the approximately normal I like to see, but at least it's a more "regular" shape. (See below...)

The interesting (reassuring?) thing is: Results are identical! Same exact genes found to be differentially expressed in each.

I'm still wondering why this would happen - especially in light of the new information. Obviously, for publication reasons, I'd rather use the analysis with the gtf file included - but I still don't like the look of those plots, and would love to be able to explain why they look the way they do.







Anyone?
DrS is offline   Reply With Quote
Old 03-24-2013, 10:43 PM   #3
cookiemic
Member
 
Location: China

Join Date: Aug 2012
Posts: 11
Default

Hi DRS,
do you think it is normal for your density plot to have negative log10(fpkm) values?
cookiemic is offline   Reply With Quote
Old 03-25-2013, 04:35 AM   #4
DrS
Member
 
Location: North Carolina

Join Date: Mar 2013
Posts: 17
Default

Not at all - I noticed that later. The dispersion as well. I still have no idea what's going on with this data, either. If anyone has any clues why this might be happening, I'm still interested in some insights.
DrS is offline   Reply With Quote
Old 03-25-2013, 05:05 AM   #5
kmcarr
Senior Member
 
Location: USA, Midwest

Join Date: May 2008
Posts: 1,169
Default

Quote:
Originally Posted by cookiemic View Post
Hi DRS,
do you think it is normal for your density plot to have negative log10(fpkm) values?
Yes, it is certainly "normal", or at least understandable that you would see negative log10(fpkm) values. It simply means that FPKM was < 1.0. Let's imagine you have an mRNA that is 1,500 bp long with 20 fragments mapped to it out of a total of 100 million mapped reads.

FPKM = 20 1.5 100 = 0.133
log10(0.1333) = -0.875
kmcarr is offline   Reply With Quote
Old 03-25-2013, 05:24 AM   #6
DrS
Member
 
Location: North Carolina

Join Date: Mar 2013
Posts: 17
Default

Thank you - at least I can feel better about that. Still don't understand the "tail" I'm seeing on the one plot, but again, when I re-run the analysis with different parameters, I don't see that and get the exact same genes listed as DE.
*shrug*
DrS is offline   Reply With Quote
Reply

Tags
cummerbund plots, dispersion

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:14 PM.


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