Question: correcting the batch effects in Limma and SVA
1
gravatar for Bogdan
4 weeks ago by
Bogdan570
Palo Alto, CA, USA
Bogdan570 wrote:

Dear all,

about correcting the batch effects in LIMMA and SVA, 'd appreciate having your comments :

assuming that we have a set of RNA-seq data (no treatment, + treatment) in many distinct BATCHES, would the linear model in LIMMA :

design <- model.matrix(~0 + BATCH + Treatment, data=RNA)

suffice to correct for batch effects, or shall we additionally use COMBAT (as it is described in the SVA package : https://bioconductor.org/packages/release/bioc/vignettes/sva/inst/doc/sva.pdf ) in order to obtain the corrected counts ?

thanks a lot,

-- bogdan

limma sva combat • 141 views
ADD COMMENTlink modified 4 weeks ago by Gordon Smyth38k • written 4 weeks ago by Bogdan570
Answer: correcting the batch effects in Limma and SVA
4
gravatar for Gordon Smyth
4 weeks ago by
Gordon Smyth38k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth38k wrote:

I would myself just put BATCH in the limma linear model and omit Combat, and that is also recommended by Vegard Nygaard's group:

https://support.bioconductor.org/p/72926/

The Combat author (Dr Johnson) has agreed with this advice unless there are large difference in batch variances, in which case he has recommended using Combat but also including BATCH in the limma linear model:

https://support.bioconductor.org/p/72815/

I'm not sure how that would be done for RNA-seq though. I guess Combat would be run on logCPM values, in which case you would have to follow it up with limma-trend rather than limma-voom.

ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by Gordon Smyth38k

Thanks a lot, Gordon. Very precious advise !

ADD REPLYlink written 29 days ago by Bogdan570

Dear Gordon, if I may add and ask please (on a related topic):

would COMBAT and LIMMA be suited for assessing the differential expression between clusters of single cells (that were generated by using 10X Genomics protocols/kits) ?

've noted a tutorial from UC Davis that propose SVA/COMBAT and LIMMA for batch correction and differential expression :

https://ucdavis-bioinformatics-training.github.io/2019-single-cell-RNA-sequencing-Workshop-UCDUCSF/scrnaseqanalysis/scRNA_Workshop-PART6.html

thanks,

bogdan

ADD REPLYlink written 29 days ago by Bogdan570

Dear Gordon, if I may add and ask please (on a related topic regarding single cell RNA-seq):

would COMBAT and LIMMA be suited for assessing the differential expression between clusters of single cells (that were generated by using 10X Genomics protocols/kits) ?

've noted a tutorial from UC Davis that propose SVA/COMBAT and LIMMA for batch correction and differential expression :

https://ucdavis-bioinformatics-training.github.io/2019-single-cell-RNA-sequencing-Workshop-UCDUCSF/scrnaseqanalysis/scRNA_Workshop-PART6.html

thanks,

bogdan

ADD REPLYlink modified 29 days ago • written 29 days ago by Bogdan570
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 278 users visited in the last hour