How to use the Combat-seq for batch effect removal?
2
1
Entering edit mode
15 months ago
Yoomi ▴ 20

I try to remove the batch effect using Combat-seq in RNA-seq count data. Should I normalize RNA-seq raw count data before or after using the Combat-seq? Show me the example code of how to use Combat-seq and Normalization.

batch-effect normalization Combat-seq • 4.3k views
ADD COMMENT
0
Entering edit mode

Show me the example code of how to use Combat-seq and Normalization.

This might be a language barrier, but please do not use the imperative tone without any "please" or "could you" or such polite phrases. It comes across as demanding and on forums, could lead to antagonizing people unintentionally.

ADD REPLY
0
Entering edit mode

Thank you so much~!!:)

ADD REPLY
2
Entering edit mode
15 months ago
ATpoint 85k

Combat-seq accepts raw counts and returns batch-corrected raw counts. Once you have the corrected counts treat them as any raw counts with any RNA-seq package of choice.

ADD COMMENT
0
Entering edit mode

Thank you so much~!^^

ADD REPLY
1
Entering edit mode
15 months ago
bk11 ★ 3.0k

should i normalized the RNA-seq raw countdata before or after using the Combat-seq?

Please see in this post below which had been discussed earlier.

In Which Order Use Normalization And Batch Effects Removal?

Also remember that standard methods like DESeq2, edgeR etc. can handle batch effect very well. You do not have to perform batch correction using outside softwares if you are planning to perform differential expression analysis using these methods.

ADD COMMENT
1
Entering edit mode

Thank you~!:)

Can I handle the batch effect in this manner?

dds.dex <- DESeqDataSetFromMatrix(countData = round(ga.ta.dex.merge), colData = coldata, design = ~ condition)
dds.dex <- DESeq(dds.dex)
vsd <- vst(dds.dex, blind=FALSE)
mat <- assay(vsd)
mm <- model.matrix(~condition, colData(vsd))
mat <- limma::removeBatchEffect(mat, batch=vsd$batch , design=mm)
assay(vsd) <- mat
ADD REPLY

Login before adding your answer.

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