Ploting of normalized counts after batch removal
2
0
Entering edit mode
altuda ▴ 10
@91f63ac0
Last seen 10 days ago
Czechia

Hello, I tried to search for a correct answer, but I don't think I've found one. I have an experiment with paired samples, I've performed DESeq2, and now I would like to plot DESeq2 normcounts. All the discussions and tutorials I found are regarding plotting heatmaps and PCA plots, where you remove the batch effect of the paired samples from vst. Now, how do you think I should proceed regarding the normcounts? Since I am introducing negative values, when I apply limma::removeBatchEffect it doesn't seem correct.

Thank you for your time.

DESeq2 • 228 views
ADD COMMENT
0
Entering edit mode

You can use ComBat-seq from the sva package to regress covariates such as the pairing information. Advantage is that it returns "raw corrected" counts that are still integers, so you can put them into vst again and plot that.

ADD REPLY
0
Entering edit mode
@mikelove
Last seen 12 hours ago
United States

You can't really plot counts with variance from nuisance covariates removed. I would go the route of using VST data instead (see FAQ).

ADD COMMENT
0
Entering edit mode
@gordon-smyth
Last seen 1 hour ago
WEHI, Melbourne, Australia

I am introducing negative values, when I apply limma::removeBatchEffect it doesn't seem correct.

The values input to removeBatchEffect() are always on the log-scale. Getting negative log-scale values is no problem. The unlogged expression values that they correspond to are always positive.

ADD COMMENT

Login before adding your answer.

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