ComBat-seq adjustment of batch effects
0
1
Entering edit mode
12 months ago

Hi, If I have done my RNA seq with 4 samples X 3 replicates each in two batches (see below for the design matrix): 1) Is it correct to adjust batch effects on the raw count file first by ComBat-seq (Zhang et.al., 2020) and then process the output by edgeR?

CONDITIONS  SAMPLES BATCH
T1            R1      1
T1            R2      1
T1            R3      1
T2            R1      1
T2            R2      2
T2            R3      2
T3            R1      1
T3            R2      1
T3            R3      1
T4            R1      1
T4            R2      2
T4            R3      2

I am asking this question since the only R2 and R3 for conditions T2 and T4 were processed in batch 2 and I could not include such a batch effect in the edgeR model matrix (see below) which needs batch 2 to be present in each of T1, T2, T3 and T4.

design = model.matrix(~ 0 + conditions + batch, data=y$samples)

Ref: Zhang et.al., 2020: https://academic.oup.com/nargab/article/2/3/lqaa078/5909519 Kindly help. Thanks. Arumoy.

RNA-Seq ComBat-seq edgeR batch effects SVA • 1.1k views
ADD COMMENT

Login before adding your answer.

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