Normalized counts access in DESeq2 R: with and without batch
1
0
Entering edit mode
Matan G. ▴ 50
@matan-g-22483
Last seen 2.6 years ago

As described previously here: https://support.bioconductor.org/p/66067/ The way to access the normalized count table in DESeq2 can be done by the following script lines: dds <- estimateSizeFactors(dds) counts(dds, normalized=TRUE)

I'm trying to understand if the normalized count table that can extracted using the above is taking into consideration a batch effect that is modeled in the design of DESeq2.

For example if I have a batch in my design: design = ~ condition + Batch 1. Is the table accessed normalized for the batch? 2. If so, is there a way to get two separate count tables: a. normalized without batch. b. normalized with batch. 3. If not, is there a way to normalize for the batch and access it using the basic functions or arguments in DESEq2.

Thanks

deseq2 edger normalization • 624 views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 19 hours ago
United States

No, size factors dont use the design.

In order to get data corrected for batch, we recommend VST followed by removal of the batch effect with a linear regression, see the FAQ in the vignette.

ADD COMMENT

Login before adding your answer.

Traffic: 616 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