Question: DESeq2 log2 fold change
0
gravatar for bioinfonerd
7 months ago by
bioinfonerd30
bioinfonerd30 wrote:

Hi everyone,

I did the DESeq2 on my Nanostring gene expression data using the below code:

dds <- DESeqDataSetFromMatrix(countData = counts_set,colData = annotation,design = ~ pt_id+treatment)
sizeFactors(dds)  <- estimateSizeFactorsForMatrix(counts(dds))
dds <- DESeq(dds)
res <- results(dds,contrast=c("treatment","Treatment","Vehicle"))

But the log2FC I am getting are quite different from the log2FC I calculated manually. Could you let me know if that is normal and correct?

Thanks!

nanostring deseq2 • 585 views
ADD COMMENTlink written 7 months ago by bioinfonerd30

Every now and then this question appears here:

Discrepancy in DESeq2 fold change

How to recover treated/control count from DESeq2 output

Why can't I replicate DESeq's Log2FC calculation?

ADD REPLYlink written 7 months ago by h.mon24k

I guess that your FCs for the genes with smaller couts are much larger than the DESeq2-ones?

EDIT: You should have a look at this post from the DESeq2 developer towards DESeq2 on Nanostring data. He is concerned that the normalization might be a problem.

ADD REPLYlink modified 7 months ago • written 7 months ago by ATpoint14k
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: 1833 users visited in the last hour