How to eliminate batch effect using removeBatchEffect function for tissue expression profiles?
0
1
Entering edit mode
6.1 years ago

For meta-analysis work, I would like to combine two rice GEO series( GSE27988 and GSE27856 ) where each series analyze multiple anatomical tissue related expression profiles.To eliminate batch effect, I'm not sure whether I can use limma removeBatchEffect function as given below. Similar posts I see related to batch effect are generally control vs treatment.

libraray(limma)

batch <- c("A","A","A",... "B","B","B")  # 24 samples in batch A and 14 samples in batch B
corrected<- removeBatchEffect(data, batch)
Microarray batch-effect R limma • 2.4k views
ADD COMMENT
1
Entering edit mode

i'd rather suggest BatchQC for batch effect removal + visualization + statistics https://bioconductor.org/packages/release/bioc/html/BatchQC.html

ADD REPLY

Login before adding your answer.

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