How to export rld (rlog) or vsd values from DESeq2?
2
2
Entering edit mode
5.4 years ago
LuisNagano ▴ 90

Hi guys, i want to export the rlog values from DESeq2 or vsd, but i only got the dds counts, using the following script:

countsdds <- counts(dds, normalized = TRUE)
write.csv(countsdds, file="dds_norm_counts.csv")

Is there a way to export rlog values?

Thanks for your help!

DESeq2 rlog vsd RNA-seq • 8.2k views
ADD COMMENT
0
Entering edit mode
5.4 years ago

Googling rld DESeq2 already gives me the answer, e.g. here

ADD COMMENT
1
Entering edit mode

Thank you, WouterDeCoster, but i don't want to apply the functions "rlog" or "vst", I just did that. I want to export the output values after run the function, like .csv or .xlsx format table to my directory.

ADD REPLY
6
Entering edit mode
5.4 years ago

Right, what about

rld <- assay(rlog(dds, blind=FALSE))
ADD COMMENT
0
Entering edit mode

Worked it out! Thank you, WouterDeCoster

ADD REPLY
0
Entering edit mode

Thank you! I was trying to extract my variance stabilizing transformation and write it to a .csv. Got it to work with this:

vst_write <- assay(vst(dds, blind=FALSE))

write.csv(as.data.frame(vst_write), file="VST_results.csv")

ADD REPLY

Login before adding your answer.

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