Variance stabilized values and batch correction
1
0
Entering edit mode
le2336 ▴ 20
@le2336-10789
Last seen 12 months 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 • 190 views
ADD COMMENT
3
Entering edit mode
@mikelove
Last seen 4 hours 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: 258 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