The following example performs hierarchical clustering on the
transformed expression matrix subsetted by the DEGs identified in the above
differential expression analysis. It uses a Pearson correlation-based distance
measure and complete linkage for cluster joining.
library(pheatmap) geneids <- unique(as.character(unlist(DEG_list[]))) y <- assay(rlog(dds))[geneids, ] pdf("heatmap1.pdf") pheatmap(y, scale = "row", clustering_distance_rows = "correlation", clustering_distance_cols = "correlation") dev.off()