Question: Taking out , vsd <- varianceStabilizingTransformation(dds) from deseq2
1
gravatar for krushnach80
2.8 years ago by
krushnach80810
krushnach80810 wrote:

I have five different cell type i m doing both rld and vsd transformation .I want the vsd data to perform other analyiss such as Self organising maps as well as WGCNA . How can i take out the vsd transformed data

Im using this code

vsd <- varianceStabilizingTransformation(dds, blind=FALSE)

To get all my normalised data set im using this code

resOrdered <- res[order(res$padj),]
summary(resOrdered)
resdata <- merge(as.data.frame(res), as.data.frame(counts(dds, normalized=TRUE)), by="row.names", sort=FALSE)
names(resdata)[1] <- "Gene"

head(resdata)
dim(resdata)
## Write results
write.csv(resdata, file="Primary_cell_line__normalised_data.csv",row.names = FALSE,quote = F)

But this is for the rlog transformation what i saw in the manual if im not wrong , how can i do the same for vsd transformed data ?Any suggestion or help would be highly appreciated

rna-seq R • 2.7k views
ADD COMMENTlink modified 2.8 years ago by andrew.j.skelton736.0k • written 2.8 years ago by krushnach80810
3
gravatar for andrew.j.skelton73
2.8 years ago by
London
andrew.j.skelton736.0k wrote:

Your second chunk of code is not for rlog transformed data but for the normalised counts. If you want the data in the same shape as the normalised counts but for rlog or vst transformed data then you need something along the lines of:

vsd <- assay(varianceStabilizingTransformation(dds, blind=FALSE))
rlog <- assay(rlogTransformation(dds, blind=FALSE))
ADD COMMENTlink written 2.8 years ago by andrew.j.skelton736.0k

so i have to do use your code then for rlog or vst ,but is it after normalisation ?

ADD REPLYlink written 2.8 years ago by krushnach80810

Have a look at ?varianceStabilizingTransformation

This function calculates a variance stabilizing transformation (VST) from the fitted dispersion-mean relation(s) and then transforms the count data (normalized by division by the size factors or normalization factors), yielding a matrix of values which are now approximately homoskedastic (having constant variance along the range of mean values). The transformation also normalizes with respect to library size. The rlog is less sensitive to size factors, which can be an issue when size factors vary widely. These transformations are useful when checking for outliers or as input for machine learning techniques such as clustering or linear discriminant analysis.

ADD REPLYlink written 2.8 years ago by andrew.j.skelton736.0k

yes i looked into the function my question is the way im getting a normalised data set for all the sample and gene can i get the same for vst or rlog

i using this code

rld <- rlog(dds, blind=FALSE)
vsd <- varianceStabilizingTransformation(dds, blind=FALSE)

am i doing it correctly or there is some error?

ADD REPLYlink written 2.8 years ago by krushnach80810
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: 1402 users visited in the last hour