heatmap3如何从价值转变为色彩？
我使用R的heatmap3软件包来生成一些基因表达数据热图。
我的问题是：基因表达的价值如何被“映射”为颜色？
作为相同（可重现）问题的示例，让我们使用mtcars数据集上的小插图中的示例代码：
library(heatmap3)
heatmap3(mtcars,scale="col",margins=c(2,10),RowSideColors=RowSideColors, balanceColor = TRUE)
图例显示颜色范围从大约-2到加3.如果我们使用玛莎拉蒂宝来作为例子
mpg cyl disp hp drat wt qsec vs am gear carb 15.0 8 301.0 335 3.54 3.570 14.60 0 1 5 8
我们可以在热图中看到它具有鲜红的颜色，即对于碳水化合物，分数为〜3，对于qsec，蓝色分数为约-1.5。
看看数据集中的其他车，他们有以下碳水化合物：
4 4 1 1 2 1 4 2 2 4 4 3 3 3 4 4 4 1 2 1 1 2 2 4 2 1 2 2 4 6 8 2
所以玛莎拉蒂宝来的碳排放量绝对是异常的，法拉利迪诺（carb = 6）也是如此，这在热图中也是预期的红色。
同样，看着qseq
16.46 17.02 18.61 19.44 17.02 20.22 15.84 20.00 22.90 18.30 18.90 17.40 17.60 18.00 17.98 17.82 17.42 19.47 18.52 19.90 20.01 16.87 17.30 15.41 17.05 18.90 16.70 16.90 14.50 15.50 14.60 18.60
玛莎拉蒂的14.60接近14.5的最低点，在我们预期的福特Pantera L中也是蓝色的。
但是，我们如何从价值到颜色？ 计算这个z分数的公式是什么？ 它是平均每车或每参数（英里，cyl等）的基础？
好的，我想我只是想出了它 - 但是在这里发帖以防其他人有同样的问题：
它使用基本的Z分数公式，所以
(mtcars$qsec - mean(mtcars$qsec))/sd(mtcars$qsec)
(value - pop_mean(value))/ sd(value)
上一篇: How does heatmap3 go from value to colour?
下一篇: bar chart