DESEq2 VST error
Last seen 4.1 years ago
I'd like to perform variance stabilization transformation in DESeq2 but got this error message. Any input would be much appreciated. Thank you.
vsd <- vst(dds, blind = TRUE)
Error in vst(dds, blind = TRUE) : less than 'nsub' rows,
it is recommended to use varianceStabilizingTransformation directly
Last seen 4 hours ago
To put more meat on my comment above, take a look at the code in the
vst function (look at
?vst and just type
vst at the console and hit return, the code will reveal itself)
You'll find that the
vst function does its magic on a subset of rows (ie. the size of the subset is equal to
nsub) then can apply it to the entire dataset.
If you call
varianceStabilizingTransform directly on your expression object
dds, it doesn't try to work on a subset of the genes in your expression matrix, and just does its thing on the whole thing.
Your problem is that your expression matrix doesn't have more than 1000 genes (rows) in it.
Traffic: 323 users visited in the last hour