DESeq2 vst function error
1
0
Entering edit mode
8 months ago

Hi everyone, I have a dataset with count reads of miRNA-Seq downloaded from TCGA: on the rows I have the miRNA names and on the columns the uuid barcodes(sample names). When I launch the vst function it gets me this error:

dd = DESeqDataSetFromMatrix(d, coldata, design= ~class) 
vst = vst(dd, blind = FALSE)

Error in vst(dd, blind = FALSE) : 
  less than 'nsub' rows with mean normalized count > 5, 
  it is recommended to use varianceStabilizingTransformation directly

This is the first attempt with this package, any suggestion? Thank you so much

DESeq2 RNA-Seq R normalization • 761 views
ADD COMMENT
1
Entering edit mode
8 months ago

Your dataset is small (< 1000 variables) and / or has a low number of expressed genes. Either way, the problem can be mitigated by using varianceStabilizingTransformation(), or just re-use vst() and set a lower value for nsub.

Type ?vst to see what is happening with the nsub parameter.

Kevin

ADD COMMENT
0
Entering edit mode

I set a lower number, and it works, thank you! My dataset has 1800 miRNAs. This does mean that only few of them are considered for the calculation? How does the algorithm operates in this cases?

ADD REPLY

Login before adding your answer.

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