From what I can find in papers, heatmaps using RNA seq data are created in several ways: using log-fold changes, z-scores, etc.

The edgeR vignette states:

> Inputing RNA-seq counts to clustering or heatmap routines designed for microarray data is
not straight-forward, and the best way to do this is still a matter of research. To draw a
heatmap of individual RNA-seq samples, we suggest using moderated log-counts-per-million.
This can be calculated by cpm with positive values for prior.count, for example :

    > logcpm <- cpm(y, log=TRUE)

Just out of curiosity, I was wondering, how would it differ from calculating z-scores using the fitted.values (derived from the glmQLFit step) in the RNA seq analysis pipeline. Would the heat maps created using z-scores calculated from fitted.values turn out all that different?


The purpose of making heatmap of logCPMs is to display sample to sample variability, which allows you to see variability both between groups and between replicates.

Plotting fitted values instead would be pointless because fitted values do not show variability between replicates, and also incorrect because fitted values are not normalized by library size.