R - d3heatmap - 实施休息

Question

我正在尝试使用 d3heatmap 包绘制热图。 不幸的是，我还没有成功地使用breaks=... 选项实现某些中断，如heatmap 或heatmap.2。 这只会产生有趣的结果，我什至不确定我是否做错了什么或者该函数是否忽略了breaks。

例如，我试过：

breaks = c(seq(-10, -2), seq(-2, -1.65), seq(-1.65, 1.65), seq(1.65, 2), seq(2, 10)

和

breaks = c(-10, -2, -1.65, 1.65, 2, 10)

与

colors = c("red", "yellow", "green", "yellow", "red")

但似乎没有任何工作正常。

有什么建议吗？

这是我的数据的dput：

 > dput(mat)
 structure(c(-0.04, NA, 0.59, NA, 0.675, 0.96, 1.09, 0.445, NA, 
 0.545, NA, NA, 0.09, -1.11, NA, 0.99, 0.13, 0.215, 1.425, 0, 
 NA, 0.69, 0.805, NA, 0.69, 1.22, NA, 0.3, NA, 0.025, NA, 0.075, 
 0.36, -0.94, NA, -0.31, 0.26, 1.02, -1.19, NA, NA, -0.77, NA, 
 -1.48, 1.05, 0.48, NA, NA, NA, 1.49, -1.285, NA, 0.76, 1.14, 
 -0.62, NA, NA, NA, 0.95, NA, NA, -0.12, 0.49, NA, 2.31, NA, -0.33, 
 0.85, NA, -1.7, -1.63, NA, -1.12, 0.135, -0.18, NA, -0.245, NA, 
 -0.2, -0.2, 0.23, -0.11, NA, 0.3, -0.81, 0.04, 0.18, -0.7, 0.53, 
 0.44, -0.49, 0.28, 0.26, 0.06, 0.265, 0.21, 0.06, -0.175, 0.365, 
 0.255, 1.25, -0.35, 0.16, 0.125, 0.825, 0.08, 0.02, -0.02, 0.99, 
 0.79, -0.23, 0.06, NA, 0.36, -0.64, -0.195, 1.19, -0.29, 0.915, 
 NA, NA, NA, NA, 0.2, 0.1, NA, 0.04, 0.33, NA, 1.46, 2.36, NA, 
 -0.92, 1.295, NA, NA, 0.8, NA, 1.09, 1.45, 5.42, NA, NA, NA, 
 1.69, 3.43, NA, 0.55), .Dim = c(37L, 4L), .Dimnames = list(c("AT", 
 "BE", "BG", "CEE", "CH", "CN", "CZ", "DE", "DK", "EA", "EE", 
 "EMU", "ES", "EU", "FI", "FR", "GB", "GR", "HR", "HU", "IE", 
 "IT", "JP", "LU", "NL", "PL", "PT", "RO", "RS", "RU", "SE", "SI", 
 "SK", "TR", "UA", "UK", "US"), c("Credit Risk", "Funding and liquidity Risk", 
 "Macro Risk", "Market Risk")))

还有我正在运行的代码：

d3heatmap(abs(mat),
          dendrogram = "none",
          breaks = c(0,1.65,2,10),
          col = c("green", "yellow", "red"),
          na.rm = TRUE)

不过，使用 heatmap.2 的相同功能可以完美运行。

Answer 1

函数 d3heatmap 根本没有“breaks”参数。如果它作为参数传入，它会被默默地忽略。 （参见 ?d3heatmap。）

另一方面，gplots 包中的 heatmap.2 函数确实有一个“breaks”参数。这就解释了行为上的差异。

幸运的是，仍然可以通过将适当的“颜色”函数传递给 d3heatmap 来获得所需的行为。它的工作原理如下。

首先是示例数据：

mat <- structure(c(-0.04, NA, 0.59, NA, 0.675, 0.96, 1.09, 0.445, NA, 
                   0.545, NA, NA, 0.09, -1.11, NA, 0.99, 0.13, 0.215, 1.425, 0, 
                   NA, 0.69, 0.805, NA, 0.69, 1.22, NA, 0.3, NA, 0.025, NA, 0.075, 
                   0.36, -0.94, NA, -0.31, 0.26, 1.02, -1.19, NA, NA, -0.77, NA, 
                   -1.48, 1.05, 0.48, NA, NA, NA, 1.49, -1.285, NA, 0.76, 1.14, 
                   -0.62, NA, NA, NA, 0.95, NA, NA, -0.12, 0.49, NA, 2.31, NA, -0.33, 
                   0.85, NA, -1.7, -1.63, NA, -1.12, 0.135, -0.18, NA, -0.245, NA, 
                   -0.2, -0.2, 0.23, -0.11, NA, 0.3, -0.81, 0.04, 0.18, -0.7, 0.53, 
                   0.44, -0.49, 0.28, 0.26, 0.06, 0.265, 0.21, 0.06, -0.175, 0.365, 
                   0.255, 1.25, -0.35, 0.16, 0.125, 0.825, 0.08, 0.02, -0.02, 0.99, 
                   0.79, -0.23, 0.06, NA, 0.36, -0.64, -0.195, 1.19, -0.29, 0.915, 
                   NA, NA, NA, NA, 0.2, 0.1, NA, 0.04, 0.33, NA, 1.46, 2.36, NA, 
                   -0.92, 1.295, NA, NA, 0.8, NA, 1.09, 1.45, 5.42, NA, NA, NA, 
                   1.69, 3.43, NA, 0.55), .Dim = c(37L, 4L),
                   .Dimnames = list(c("AT", "BE", "BG", "CEE", "CH", "CN", "CZ", "DE", "DK", "EA", "EE", "EMU", "ES", "EU", "FI", "FR", "GB", "GR", "HR", "HU", "IE", "IT", "JP", "LU", "NL", "PL", "PT", "RO", "RS", "RU", "SE", "SI", "SK", "TR", "UA", "UK", "US"), c("Credit Risk", "Funding and liquidity Risk", "Macro Risk", "Market Risk")))

假设我们需要以下三个颜色箱：蓝色表示值 = 0 但 = 2。然后我们定义相应的颜色排序列表。

palette <- c("blue", "green", "red")

我们还定义了颜色箱的边界值。这些值必须包括域边界。

mi <- min(mat, na.rm = TRUE)
ma <- max(mat, na.rm = TRUE)
breaks <- c(mi, 0, 2, ma)

我们现在可以定义一个颜色插值函数，它将 [0,1] 中的值映射到颜色上，同时尊重我们的颜色箱。 'scales' 包可以在这里提供帮助。

install.package('scales') # if needed
library(scales)
colorFunc <- col_bin(palette, bins = rescale(breaks))

最初在我们的数据域中定义的中断需要重新调整为 [0,1]。 'scales' 包中的 'rescale' 函数处理了这个问题。

小细节：bin 的低边界包含，但高边界排除。所以值 0 将是绿色，介于 0 和 2 之间的任何值也将是绿色，但 2 将是红色。

我们现在可以绘制热图了。

d3heatmap(mat, dendrogram = "none", colors = colorFunc, na.rm = TRUE)

结果如下：