Evaluating removing batch effect tools with PCA or MDS?
0
1
Entering edit mode
8.3 years ago
Lluís R. ★ 1.2k

It is the first time doing a meta-analysis on microarrys. I have two studies from two different xips I want to use. I used MergeMaid to combine them and now I am looking to correct the batch effect. I have done it with ComBat and removeBatchEffect. To asses the best batch removing tool I performe a PCA and a MDS plot. In both cases the PCA and the MDS show an improvement against not using it.

The PCAs are quite different, when calculated with ComBat and removeBatchEffect, due to the underlaying the different system. However the MDS are quite similar. Does this imply something about the batch effect?

PCA of raw data:

image: PCA of raw data

MDS of raw data:

image: MDS of raw data

PCA of combat correction:

image: combat PCA

MDS of combat correction:

image: combat MDS

PCA of removeBatchEffect:

image: removeBatchEffect PCA

MDS of removeBatchEffect:

image: removeBatchEffect MDS

bioconductor limma sva r batch-effect • 4.3k views
ADD COMMENT
0
Entering edit mode

I have worked with combat on RNA-seq data and I have definitely seen improvement in the overall results.the data can be compared after batch removing batch effects.

ADD REPLY
0
Entering edit mode

Thanks Ron, do you know if it is better with MDS or PCA to asses the changes? I mean, both function ComBat and removeBatchEffect work find according to MDS but on a PCA only Combat performs better. And ComBat is quite more "aggressive" with data than other methods, that's why it is more effective in big batch effects.

ADD REPLY
0
Entering edit mode

I have not worked with MDS yet,so I cant comment on that.But,PCA gives good results.Whether be it library effects or just batches from different times,after doing combat.

ADD REPLY

Login before adding your answer.

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