Dear all,
I think this is a pretty trivial question, but how should one normalize count RNAseq data for clustering/mds/etc?
Will his suffice?
library (DESeq2)
cnts <- matrix(rnbinom(n=1000, mu=100, size=1/0.5), ncol=10)
normed= varianceStabilizingTransformation(cnts)
pheatmap::pheatmap(cor(normed))