Question: Differential expression on Combat-adjusted FPKM
0
gravatar for mare116
12 months ago by
mare1160
mare1160 wrote:

I have RNA-seq for 62 samples from 4 groups sequenced in 2 known batches, with a clear batch-effect in the PCA. I've derived expression using FPKM, and I've used Combat on log2(FPKM + 1) values. This removes the batch effect successfully, but now I would like to analyze differential gene expression between the 4 groups.

My question/concern: I've read a few people saying that you shouldn't use these Combat-adjusted values for differential expression analysis, but I don't understand why. Can someone please comment on this or explain?

ADD COMMENTlink written 12 months ago by mare1160

For differential expression,you should use counts. Refer to this post: A: RNA sequencing data batch effect removal

ADD REPLYlink written 12 months ago by Ron660

An alternative to account for a batch effect problems is to model the batch effect in the design formula of DESeq2/edgeR/limma-voom.

ADD REPLYlink written 12 months ago by WouterDeCoster23k
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 2.3.0
Traffic: 942 users visited in the last hour