Question: RNA-Seq and edgeR: remarkable differences between cpms and fitted values
0
gravatar for David Rengel
4 months ago by
David Rengel70
European Union
David Rengel70 wrote:

Hi, I am treating NovaSeq-produced RNA-Seq data with edgeR. In the analysis I look at what I call “strain” effect. There are DEGs for this factor that I find are biologically hard to explain. In the attached file, there are two examples. It is easy to see why they are declared as differentially expressed since fitted values vary among strains (right hand plots in the figure). What is hard for me to follow is how those fitted values might be obtained when the TMM-normalized cpms do not reveal that behavior (left-hand plots). In the latter case, cpms tend to 0 except for one repeat on the WT strain. I thank you in advance for any help/comment on the matter. Best, David

fitted vs cpm expression

ADD COMMENTlink modified 4 months ago by Gordon Smyth38k • written 4 months ago by David Rengel70

From the FAQ

How do I put images into my posts?

You can also put images and plots into your posts. All you need to do is create (or login to) a free account with http://imgur.com and upload any images or plots that you want to share. Then you can enter the imgur generated 'Direct Link' url into the 'URL' location that is available when you click on the 'image' button in the text editor. The plot will be embedded right into your post. To "attach" other files to your posts, upload them to Dropbox or another similar site and include the URL in your post. Make sure the URL is publicly accessible.

ADD REPLYlink written 4 months ago by James W. MacDonald51k

Hi James, Thanks, but is not it what I have done? I mean, there is a link to the figure. I've put a new imgur link just in case the previous one was not public. Hope it works. Best, David

ADD REPLYlink modified 4 months ago • written 4 months ago by David Rengel70

Edit I can't get the 'direct link' from your post, so it doesn't work for me.

ADD REPLYlink modified 4 months ago • written 4 months ago by James W. MacDonald51k

I can get a direct link to one of my own imgur posts, and it looks like this

enter image description here

ADD REPLYlink written 4 months ago by James W. MacDonald51k

I updated this post -- was the first time I used imgur, actually, so wasn't sure what URL hacking was necessary to get the image to show up.

You don't copy & paste the URL of the image page from imgur (which is https://imgur.com/a/oxlZBoo), but rather you have to right click on the image itself and get that URL (which is https://i.imgur.com/daGn4Q0.png).

I guess having to do this is obvious in retrospect, but still a bit of a PITA for most people the first time or so.

From the text in the FAQ, it's not clear where one gets the "Direct Link" that's written there ... I'm guessing that this is something that imgur shows when the image is first uploaded/created? Because I can't find it on the imgur webpage after the fact ...

... aaaaaaannnnnnyyyyyyyway ... there we have it.

ADD REPLYlink written 4 months ago by Steve Lianoglou12k

Yes, but if it's your imgur post, and you click on it, you are presented with a dialog box that looks like this:

enter image description here

And if you look to the right, there's a Direct Link that you can copy to post. Hence the part in the FAQ about a 'Direct Link'.

ADD REPLYlink written 4 months ago by James W. MacDonald51k

recursion

ADD REPLYlink written 4 months ago by Steve Lianoglou12k

perfect, thanks for the tips! :-)

ADD REPLYlink written 4 months ago by Steve Lianoglou12k
1

Keymaker, my friend

enter image description here

ADD REPLYlink written 4 months ago by James W. MacDonald51k

Can you add more details (code) to your question, please?

  1. Can you share a bit of your $samples data.frame and show us the code you used to generate the design and the fit?
  2. Where did you get these fitted values from? Is it off of some $fitted.values somewhere?
ADD REPLYlink written 4 months ago by Steve Lianoglou12k
Answer: RNA-Seq and edgeR: remarkable differences between cpms and fitted values
0
gravatar for Gordon Smyth
4 months ago by
Gordon Smyth38k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth38k wrote:

It doesn't make any sense to plot fitted values because they are proportional to the library sizes, which in turn vary independently of the groups or treatment conditions of interest. The reason why we recommend that you use CPMs or logCPMs for plotting purposes is that the library sizes have been divided out.

ADD COMMENTlink modified 4 months ago • written 4 months ago by Gordon Smyth38k

Thanks a lot Gordon, it is clearer now. Best, David

ADD REPLYlink written 4 months ago by David Rengel70
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 16.09
Traffic: 154 users visited in the last hour