Variance stabilized values and batch correction
1
0
Entering edit mode
le2336 ▴ 20
@le2336-10789
Last seen 2.1 years ago

I am interested in using variance stabilized values of my data (that are the outputs of DESeq2's varianceStabilizingTransformation function) for machine learning and other applications. I ran the following:

dds <- DESeqDataSetFromMatrix(countData = counts, colData = metadata, design = ~Institution + Condition)
vsd <- varianceStabilizingTransformation(dds, blind = FALSE)

Although I tried to include the batch (Institution) in the design argument, I still see batches in my PCA plot. This section of the vignette described a solution for this that I have pasted below:

mat <- assay(vsd)
mat <- limma::removeBatchEffect(mat, vsd$batch)
assay(vsd) <- mat
plotPCA(vsd)

Can the output of removeBatchEffect in this code now be used for all other downstream applications, and not just for visualization?

Thanks very much.

deseq2 vst • 526 views
ADD COMMENT
3
Entering edit mode
@mikelove
Last seen 5 days ago
United States

Yes, I have used this kind of workflow to remove unwanted variation from VST data.

ADD COMMENT
0
Entering edit mode

Thank you for the clarification!

ADD REPLY

Login before adding your answer.

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