ComBat-seq adjustment of batch effects
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, 2020) and then process the output by edgeR?

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, 2020: Kindly help. Thanks. Arumoy.

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

Login before adding your answer.

Traffic: 1635 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6