transformed read counts in deseq2

Question: transformed read counts in deseq2

6 months ago by

anna • 0

anna • 0 wrote:

How can I extract the normalized and transformed read counts in deseq2 which are used to calculate the log fold change?

read counts deseq2 data • 294 views

ADD COMMENT • link •

modified 6 months ago by Prakash • 1.8k • written 6 months ago by anna • 0

See answer here: retrieve normalised count data from DESeq2

ADD REPLY • link written 6 months ago by genomax ♦ 78k

i've tried that function but the data are not transformed. Also i've read that the function (rlog) retrieve the transformed counts but i can't get any result by doing this: write.csv(as.data.frame(rlog(dds, blind=TRUE, fitType = "parametric")), file="transformed.csv")

ADD REPLY • link modified 6 months ago • written 6 months ago by anna • 0

What function exactly did you try? rlog values are not used for calculating log fold changes.

ADD REPLY • link written 6 months ago by swbarnes2 ♦ 7.5k

6 months ago by

ATpoint ♦ 30k

Germany

ATpoint ♦ 30k wrote:

rld <- rlog(dds)
assay(rld)

How to export rld (rlog) or vsd values from DESeq2?

As far as I know essentially all DESeq2 object types are SummarizedExperiment containers and therefore can be accessed with assay to select counts. data.frame(rlog) is not possible as you first have to access the counts inside the SummarizedExperiments container.

ADD COMMENT • link modified 6 months ago • written 6 months ago by ATpoint ♦ 30k

6 months ago by

Prakash • 1.8k

India

Prakash • 1.8k wrote:

you can also extract normalized read count from dds object

res <- as.data.frame(counts(dds, normalized=TRUE))

ADD COMMENT • link written 6 months ago by Prakash • 1.8k

Please log in to add an answer.

Similar posts • Search »