DESEQ2 Batch effect
1
0
Entering edit mode
sherajilir • 0
@sherajilir-18567
Last seen 3.9 years ago

Hello everyone,

I am analyzing coad-tcga dataset and there are two sequencing platforms used, illuminaHiseq and IlluminaGA. Is there a way to remove batch effects in DESeq2 package prior to differential expression or does the GLM remove it automatically (which i dont hink happens).

Thanks!

deseq2 batch • 923 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 3 hours ago
United States

Take a look at our vignette. There we have a section on multi-factor design which shows how to control for covariates.

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#multi-factor-designs

ADD COMMENT
0
Entering edit mode

Thank you Michael for your prompt response. I did normalization and i'm getting an MA plot like this

https://drive.google.com/file/d/1QeCzKQcTR5A-Wx7kGU52MvhfEj-27eCH/view?usp=sharing

Do you think its acceptable? I have doubts because its not aligned to zero.

Thank you once again

ADD REPLY
0
Entering edit mode

Sorry for writting separately but just came to my mind to use shrinkage with apeglm package and got this

https://drive.google.com/file/d/1ygZ1n9J_3he0L5ecFRFRTIPJV6TRzmmE/view?usp=sharing

Now i think the problem is solved, though your feedback would be highly valuable.

Thanks again

ADD REPLY

Login before adding your answer.

Traffic: 974 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6