Hello all,

I have 2 scripts for visualizing the TPM values.

1. In one script the variance of each gene is computed across the samples.
2. In another, the counts are log2 transformed, followed by a heatmap.

Which is the correct way of visualizing TPM values.

1.

 V <- apply(countdata, 1, var)

selectedGenes <- names(V[order(V, decreasing = T)][1:1000])

pheatmap(newdata[selectedGenes,], scale = 'row', 
         show_rownames = T, clustering_distance_rows = "correlation",
         cluster_rows = TRUE, cluster_cols = TRUE,  fontsize_row = 5,
         fontsize = 5, 
       )  

2.

 logtransformed <- log2(countdata + 1 )

my_palette <- colorRampPalette(c("green", "black", "red"))(n = 1000)

z.mat = t(scale(t(logtransformed), scale = T, center = T))

heatmap.2(z.mat, dendrogram="both", scale="none", trace="none", col = my_palette, cexCol =0.5