VST in DESeq2 giving negative values

Hi, I used the varianceStabilizingTransformation function on my dds object. However, the VST output results in negative values for certain genes, is that supposed to happen? I know this can happen with the rlog transformation, but does it also happen with vst ?

deseq2 r bioconductor

• 1,780 views • link

updated 4.0 years ago by ATpoint

• written 4.0 years ago by sanatbhadsavle •

1 answer

https://support.bioconductor.org/p/85983/#85985

• 0 views • link

written 4.0 years ago by ATpoint

More posts like this

What is the right choice of normalization and preprocessing before WGCNA using RNAseq data in a coh…
written 1.8 years ago by sukeshinik5 •

Hello, I am performing WGCNA on RNA sequencing data for a cohort. We are struggling with the normalization method that will not interfere with the …
Weird correlation between RIP-seq replicates
written 1.9 years ago by Hiba •

Hello, I have a RIP-seq experiment with 3 replicates per condition. When following the classic RNA-seq analysis pipeline and after quantification with Salmon I've done …
should I use normalized counts or transformed of normalized counts for RNA-seq association analysis?
written 2.8 years ago by Alewa

Hi, I run differential gene expression with DeSeq2. ```r dds <- DESeq(dds) rld <- rlog(dds, blind=False) vsdata <- vst(dds, blind=FALSE) ``` If my goal is …
No differentially expressed genes
written 4.7 years ago by bart •

Hi, I have a dataset now with 36k lncRNAs and I'm using DESeq2 to find differentially expressed lncRNAs between a healthy group and a disease …
Precise order for RNA seq transformation and gene clustering.
written 5.9 years ago by TonyCN •

Hello all, I'm very much working alone on a particular body of work that requires RNA seq analysis. I've been using the tutorial: "RNA-seq workflow: …
How to get the expression of sample combined the biological replicates in DESeq2?
written 6.7 years ago by rysterzhu •

I have the data of 2 sample, and each sample has 3 replicates RNA-seq. design = ~ replicate + sample. i.e. ddsMatrix <- DESeqDataSetFromMatrix(countData = …
Is the heatmap matruix from vst expressed in fold change?
written 9.1 years ago by Lila M

Hi everybody, following the next code vsd <- varianceStabilizingTransformation(dds, blind = TRUE) vst <- assay(vsd) dds2 <- DESeq(dds) res2 <- results(dds2) mat <- assay(vsd)[head(order(res.05$padj), 300), …
Correct pipeline for DESeq analysis
written 9.2 years ago by Lila M

Hi everybody. I'm working with RNASeq data and I am relatively new in that field. I would like to know if my approach is correct …
PEER normalization after DESeq run
written 10.8 years ago by tonja.r

In PEER-paper they write: > When RNA-seq estimates are used for transcript abundance, we recommend using DESeq to estimate library sizes and variance-stabilize expression data …
what does "Error: is(cds, "CountDataSet") is not TRUE" mean?
written 11.7 years ago by Parham

Hi, I am practicing DESeq2 manual (Differential analysi s of count data - the DESeq2 package). When I do `varianceStabilizingTransformation `with dds created from my …