DEXSeq cluster heat map like DESeq is desirable

[katzman]

We really like the cluster heat map that we get when using DESeq.

Is it possible to get a similar plot from DEXSeq?

This is the code to do it in DESeq:

# need a "full blind" version for variance-stabilized-dispersion heatmap. get it before modifying cds
fullBlind.method      <- "blind"
fullBlind.sharingMode <- "fit-only"
fullBlind.fitType     <- "local"
cdsFullBlind <- estimateDispersions(cds, method=fullBlind.method, sharingMode=fullBlind.sharingMode, fitType=fullBlind.fitType  )
# per the vignette, should use "blind" dispersion estimates to avoid bias
vsd <- getVarianceStabilizedData(cdsFullBlind)
dists <- dist( t(vsd) )
heatmap (as.matrix(dists), symm=TRUE)