3

我对 geom_ribbon 有疑问，我认为这是一个错误，而不是一个功能。我想放大我的情节的“有趣”部分，但我不希望 ggplot 仅仅因为整个事情不适合情节而排除任何东西。为此，我使用 coord_cartesian 进行限制。它适用于线和点，可能还有许多其他东西（条），但不适用于 geom_ribbon。所以这里有一个例子：

# Load libraries
library(ggplot2)

# Create data:
set.seed(1234)
LineA=c(seq(1,20,0.1))
LineB=c(seq(1,25,0.1))
LineC=c(seq(1,19,0.1))
LineD=c(seq(1,60,0.1))

my_df=data.frame(Mean = c(sort(sample(LineA,40)),sort(sample(LineB,40)),sort(sample(LineC,40)), 
                           sort(sample(LineD,40))))

my_df$Names=c(rep("Line-A",40),rep("Line-B",40),rep("Line-C",40),rep("Line-D",40))
my_df$SD=c(runif(n = 120, min = 1, max = 5),runif(n = 40, min = 1, max = 20))
my_df$Time=c(1:40,1:40,1:40,1:40)
my_df$Mean_low=my_df$Mean-my_df$SD
my_df$Mean_low[my_df$Mean_low<0]=0
my_df$Mean_hi=my_df$Mean+my_df$SD
head(my_df)

# Plot
# Ribbon visible:
ggplot(my_df, aes(x=Time, y=Mean)) + geom_line(aes(colour = Names), size = 1) +
    geom_point(size = 2, aes(shape = Names, color = Names))+
    geom_ribbon(aes(x = Time, y=NULL, ymin = Mean_low, ymax = Mean_hi, fill = Names),
                show.legend = F, linetype = 0, alpha = 0.1, na.rm = T) +
    geom_hline(yintercept = 20, linetype = "dotdash", color = "red", size = 1)+
    theme_classic()+
    scale_y_continuous("Mean value", breaks = seq(0, 100, 2), expand = expansion(mult = c(0, 0.01))) +
    scale_x_continuous("Days", breaks = seq(0, max(my_df$Time),2),
                       expand = expansion(mult = c(0.01, 0.005))) +
    coord_cartesian(ylim = c(0, 100), xlim = c(0, 50))

丝带可见 如果允许所有功能都适合绘图，则功能区可见，但当我限制 y 轴时，功能区完全缺失 Line-D，如下所示：

ggplot(my_df, aes(x=Time, y=Mean)) + geom_line(aes(colour = Names), size = 1) +
    geom_point(size = 2, aes(shape = Names, color = Names))+
    geom_ribbon(aes(x = Time, y=NULL, ymin = Mean_low, ymax = Mean_hi, fill = Names),
                show.legend = F, linetype = 0, alpha = 0.1, na.rm = T) +
    geom_hline(yintercept = 20, linetype = "dotdash", color = "red", size = 1)+
    theme_classic()+
    scale_y_continuous("Mean value", breaks = seq(0, 100, 2), expand = expansion(mult = c(0, 0.01))) +
    scale_x_continuous("Days", breaks = seq(0, max(my_df$Time),2),
                       expand = expansion(mult = c(0.01, 0.005))) +
    coord_cartesian(ylim = c(0, 30), xlim = c(0, 50))

没有丝带

我发现只有一种解决方法，如这里所述：通过手动删除数据（值的 NA）来扩展 geom_ribbon 中的范围以获取超出限制的值，但这是一种解决方法，而不是解决方案。限制/缩放适用于大多数其他 geom 选项，那么为什么不适用 geom_ribbon 呢？有谁知道更优雅的解决方案？它是一个错误吗？我应该尝试让 ggplot 的人知道吗？

谢谢！！

4

1 回答 1

0

只需安装 ragg 库 [library(ragg)] 即可在导出/保存绘图时显示功能区。但是，在 R-studio 绘图中放大时，截止带仍然不可见。这可能是 ggplot 中的错误。

于 2021-12-13T15:41:49.330 回答