14

我有两组数据，我想使用带有 viridis 色标的热图来呈现它们。对于第一个数据集，我的值范围从 0 到 1.2，我可以很容易地看到我想看到的差异。然而，我的第二个数据集有一些异常值，导致范围从 0 到 2。现在更难看到 0 和 1 之间有趣范围内的差异，并且更难以直接比较这两个图像。是否有可能使用 viridis 色阶显示从 0 到 1.2 的数据，同时以黄色显示更高的值（viridis 色阶的“最高”颜色）？这是一个例子：

library(viridis)

#Create Data
DataSet1 <- expand.grid(x = 0:5, y = 0:5)
DataSet1$z <- runif(36, 0, 1.2)

DataSet2 <- expand.grid(x = 0:5, y = 0:5)
DataSet2$z <- runif(36, 0, 2)

#Plot Data
ggplot(DataSet1, aes(x, y, fill = z)) + 
  geom_tile() +
  scale_fill_viridis() +
  geom_text(aes(label = round(z, 2)), size = 2)

DataSet1：0.5和0.7的区别很容易看出

在此处输入图像描述

ggplot(DataSet2, aes(x, y, fill = z)) + 
  geom_tile() +
  scale_fill_viridis() +
  geom_text(aes(label = round(z, 2)), size = 2)

DataSet2：0.5 和 0.7 之间的差异很难看出

在此处输入图像描述

4

4 回答 4

14

@ClausWilke 的解决方案更好，因为它显示在图例中，但有时只需要一个快速解决方案而无需编写太多特定代码。这个也依赖scales包

ggplot(DataSet2, aes(x, y, fill = z)) + 
  geom_tile() +
  scale_fill_viridis(limits = c(0.2, 1), oob = scales::squish) +
  geom_text(aes(label = round(z, 2)), size = 2)

在此处输入图像描述

于 2019-02-19T10:01:53.487 回答
12

您可以定义任意重新缩放功能。不确定这看起来有多棒，可能需要对图例进行一些处理，但原则上，这种机制允许您以任何您想要的方式将数据值映射到比例上。

ggplot(DataSet2, aes(x, y, fill = z)) + 
  geom_tile() +
  scale_fill_viridis(rescaler = function(x, to = c(0, 1), from = NULL) {
    ifelse(x<1.2, 
           scales::rescale(x,
                           to = to,
                           from = c(min(x, na.rm = TRUE), 1.2)),
           1)}) +
  geom_text(aes(label = round(z, 2)), size = 2)

在此处输入图像描述

于 2018-01-24T18:50:32.977 回答
5

你在寻找这样的东西吗？

ggplot(DataSet2, aes(x, y, fill = z)) + 
  geom_tile() +
  scale_fill_gradient(low="green", high="red", limits=c(0, 1.2), 
                      na.value = "yellow") +
  geom_text(aes(label = round(z, 2)), size = 2)

在此处输入图像描述

使用viridis颜色，asper jazzurro 推荐。

ggplot(DataSet2, aes(x, y, fill = z)) + 
  geom_tile() + 
  scale_fill_gradientn(colors = viridis_pal()(9), limits=c(0, 1.2), 
                       na.value = "#FDE725FF") + 
  geom_text(aes(label = round(z, 2)), size = 2)

在此处输入图像描述

于 2018-01-24T16:06:29.300 回答
2

这不一定是一种改进，但您可以执行以下操作以黄色显示较高的值：

DataSet2A <- DataSet2 %>% filter(z <= 1.2)
DataSet2B <- DataSet2 %>% filter(z > 1.2)

ggplot(DataSet2A, aes(x, y, fill = z)) +
  geom_tile() +
  scale_fill_viridis(begin = 0, end = .75) +
  geom_text(aes(label = round(z, 2)), size = 2) +
  geom_tile(data = DataSet2B, aes(x, y), fill = "yellow")

在此处输入图像描述

也许如果您使用截止以及比例中的begin=和end=参数来控制viridis您正在使用的比例部分，您可以获得您想要的结果。（请注意，每个绘图只能有一个填充比例，但您可以设置额外的恒定填充，就像我在这里用黄色所做的那样。）

于 2018-01-24T15:56:54.593 回答