Entering edit mode
6.0 years ago
rice.researcher
▴
220
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)
i'd rather suggest BatchQC for batch effect removal + visualization + statistics https://bioconductor.org/packages/release/bioc/html/BatchQC.html