transformed read counts in deseq2
2
0
Entering edit mode
3.3 years ago
anna ▴ 10

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

deseq2 data read counts • 2.3k views
ADD COMMENT
1
Entering edit mode
ADD REPLY
0
Entering edit mode

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
0
Entering edit mode

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

ADD REPLY
2
Entering edit mode
3.3 years ago
ATpoint 66k
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
1
Entering edit mode
3.3 years ago
Prakash ★ 2.2k

you can also extract normalized read count from dds object

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

Login before adding your answer.

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