Differential expression on Combat-adjusted FPKM
0
0
Entering edit mode
4.9 years ago
mare116 • 0

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?

RNA-Seq FPKM Combat differential expression • 2.3k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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 REPLY

Login before adding your answer.

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