5

我正在使用 R 来求解具有 R 中的两个变量的方程。我通过序列改变这两个变量，然后使用带有插值的 geom_raster 来生成解图。

但是，当我绘制数据时，我注意到绘图区域超出了我作为输入提供的 y 轴值的范围。

如果我使用 ylim 对该区域施加限制，则插值图会缩小并且不再限制值。如果我使用 vjust = 0，则图的顶部限定了该区域，但底部没有。

我的可重现示例代码如下：

library(ggplot2)
library(reshape2)


x_range <- seq(0.001,10, by=0.001)
y_range <- (seq(1*10^-10, 1*10^-9, by = 1*10^-10))

FUN <- function(x, y) log10((1)/(x*y))

data <- outer(x_range, y_range, FUN)

colnames(data) <- y_range
rownames(data) <- x_range

melted_data <- melt(data)

p <- ggplot(data=melted_data)

# basic plot
p + geom_raster(aes(x=Var1, y=Var2, fill=value), interpolate = TRUE) +
  geom_hline(yintercept = 1*10^-10) +
  geom_hline(yintercept = 1*10^-9)

在此处输入图像描述 

# with imposed ylim
p + geom_raster(aes(x=Var1, y=Var2, fill=value), interpolate = TRUE) +
  geom_hline(yintercept = 1*10^-10) +
  geom_hline(yintercept = 1*10^-9) + 
  ylim(1*10^-9, 1*10^-10)

在此处输入图像描述

# with modified vjust
p + geom_raster(aes(x=Var1, y=Var2, fill=value), interpolate = TRUE, vjust = 0) +
  geom_hline(yintercept = 1*10^-10) +
  geom_hline(yintercept = 1*10^-9)

在此处输入图像描述

我希望最终的绘图能够限制给定的实际数据（由 y = 1*10^-9 和 y = 1*10^-10 处的两条水平线指定。

我认为我观察到的是由于插值像素大小，但我不确定。

4

1 回答 1

1

我认为这与插值无关，并且与默认geom_raster情况下（的特殊情况geom_tile）使用x/y作为每个图块的中心这一事实有关。指定vjust = 1将瓷砖向下移动，但想法仍然相同。您看到的超出限制的geom_hline是第一行瓷砖的上半部分和最后一行瓷砖的下半部分。

如果您不想显示这些半图块，您可以将 y 轴限制设置为 incoord_cartesian而不是ylim（这会将限制应用于坐标系，而不是比例；?coord_cartesian有关更多详细信息，请参阅）以放大到特定数据范围：

p + geom_raster(aes(x=Var1, y=Var2, fill=value), interpolate = TRUE) +

  # limit coordinate system, without expansion
  coord_cartesian(ylim = c(1*10^-10, 1*10^-9), expand = FALSE)

协调限制法

或者，如果您有其他几何图层超出c(1*10^-10, 1*10^-9)且不想缩放，您可以添加一个遮罩层来覆盖半块：

p + geom_raster(aes(x=Var1, y=Var2, fill=value), interpolate = TRUE) +

  # add masking layers
  annotate(geom = "rect", xmin = -Inf, xmax = Inf, ymin = -Inf, ymax = 1*10^-10, fill = "white") +
  annotate(geom = "rect", xmin = -Inf, xmax = Inf, ymin = 1*10^-9, ymax = Inf, fill = "white")

掩蔽方法

于 2019-01-08T05:10:48.243 回答