The support.bioconductor.org editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: Mean-Variance Plot for Normalised Counts
0
gravatar for Dario Strbenac
8 months ago by
Dario Strbenac1.4k
Australia
Dario Strbenac1.4k wrote:

A SeqExpressionSet can store unnormalised counts and normalised counts, but it only seems to be possible to plot the unnormalised counts by meanVarPlot. It would be useful to be able to plot normalised counts to visualise the effect of applying a normalisation method (e.g. DESeq2's varianceStabilizingTransformation) to the count data set (i.e. is it visibly homoscedastic once it's transformed?). Also, if a SeqExpressionSet is created by newSeqExpressionSet, there is a table of NA created in set@assayData named normalizedCounts. If a table of normalised values is added later by using normCounts(set) <- measurementsVST, then an element named normCounts is added to set@assayData. So there are two tables for normalised values simultaneously in assayData. Can this inconsistency be resolved?

edaseq meanvarplot • 134 views
ADD COMMENTlink modified 8 months ago by davide risso810 • written 8 months ago by Dario Strbenac1.4k
Answer: Mean-Variance Plot for Normalised Counts
1
gravatar for davide risso
8 months ago by
davide risso810
Weill Cornell Medicine
davide risso810 wrote:

This is a bug and should be fixed. I'll take care of this ASAP.

normCounts(set) <- measurementsVST should store the normalized counts in the normalizedCounts slot.

The meanVarPlot function uses the normalized counts, if available. But of course it looks for it as normalizedCounts not normCounts. So fixing the bug will also make the behavior of the function what you expect.

ADD COMMENTlink modified 8 months ago • written 8 months ago by davide risso810
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 16.09
Traffic: 442 users visited in the last hour