Question: How to plot the batch corrected values in DESeq2
gravatar for dhwani.dholakia
11 months ago by
dhwani.dholakia0 wrote:

Hello, I have used DESeq2 for differential analysis of small RNA data. I wanted to plot jitter plot of top most differentialy expressed gene using the plotCounts function. When i plot them i can still see the batch effect in the image generated.I then refereed to your post There it is mentioned that it will model the batch effect in linear regression not remove them also in the following threads it is mentioned that we can use limma package removebatcheffect to visualize the batch effect removal. Based on this i have my question as follows:-

Will the batch effect modeled by limma will be same as that by DESeq2 and can that be used as exact proxy or is there any way to extract and plot batch effect removed modeled values

normalization deseq2 • 398 views
ADD COMMENTlink modified 11 months ago by Michael Love26k • written 11 months ago by dhwani.dholakia0
Answer: How to plot the batch corrected values in DESeq2
gravatar for Michael Love
11 months ago by
Michael Love26k
United States
Michael Love26k wrote:

The removeBatchEffect() function removes the shifts in mean values. This is not numerically identical to what happens when you include batch in the design of a GLM, but it’s very close, and close enough to be useful for visualization of the counts. You could plot VST data with removeBatchEffect applied and this would represent log2-scale normalized counts with shifts from batch effects removed.

ADD COMMENTlink written 11 months ago by Michael Love26k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 203 users visited in the last hour