Batch effect and DE using TCGA and my data
1
0
Entering edit mode
5.2 years ago

Hello,

I'm trying to perform a differential gene expression analysis using some TCGA data that I reprocessed and my data. However I'm worried about the batch effect.

Normally to perform deseq I use this:

dds <- DESeqDataSetFromHTSeqCount(sampleTable = d, directory = directory, design= ~ batch)
dds <- estimateSizeFactors(dds)
dds <- DESeq(dds)

and then to plot an remove the batch effect I use this:

vsd <- vst(dds)
plotPCA(vsd, "batch")
assay(vsd) <- limma::removeBatchEffect(assay(vsd), vsd$batch)
plotPCA(vsd, "batch")

my doubt is about to perform differential gene expression after removing this batch effect with limma, is there any way? or DeSeq itself do this when I put batch as design?

Thanks

Thaise

RNA-Seq R • 1.5k views
ADD COMMENT
0
Entering edit mode

Please use the formatting bar (especially the code option) to present your post better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.
code_formatting

ADD REPLY
2
Entering edit mode
5.2 years ago

or DeSeq itself do this when I put batch as design?

Yes. Don't try to correct for it, just let DESeq take it into account when making its models. You can also incorporate the design when calculating the PCA plot.

ADD COMMENT

Login before adding your answer.

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