如何仅在预测线上绘制 CI答案

【问题标题】：How to plot a CI only in the predicted line如何仅在预测线上绘制 CI
【发布时间】：2021-01-11 21:51:25
【问题描述】：

我使用 GNM（广义非线性模型）进行预测并将所有数据存储在数据框中。我还存储了一个pred_up 和pred_down 变量，它们只是置信区间的上下部分，我想将其绘制为阴影区域。通过下面的代码，我得到以下图表

ggplot(plot_males, aes(x = year, y = real_data)) +
        geom_line(aes(color = "Obs"), na.rm=TRUE) +
        geom_line(aes(x = year, y = pred, color = "Prediction"), na.rm=TRUE) +
        geom_line(aes(x = year, y = pred_up, color = "CI"), na.rm=TRUE) +
        geom_line(aes(x = year, y = pred_down, color = "CI"), na.rm=TRUE)

我想要完成的是类似于下面的图表，置信区间区域为阴影。我知道 ggplot2 有一种直接绘制置信区间的方法，但它没有 GNM 方法，所以我将我的上下 CI 区间存储在数据集中。

有没有办法给区间加阴影？

【问题讨论】：

看看geom_ribbon。
您可以将任意模型函数和-参数传递给“ggplot2”的stat_summary 函数，因此您可以用它绘制一个GNM CI。 ……当然，如果你已经计算了 CI，那么就没有必要重新计算了。
请参阅this vignette，了解如何将 geom_ribbon() 与预先计算的置信范围一起使用

标签： r ggplot2

【解决方案1】：

尝试添加

geom_ribbon(aes(x = year, ymin = pred_down, ymax = pred_up), 
              fill = "red", alpha = 0.15, na.rm = TRUE)

【讨论】：