从R--ggplot2中混合模型的模型预测中添加置信区间？

我从我的数据中得到了平均ans置信区间的模型预测，我想将其添加到图表中。我知道如何绘制数据，但如何添加模型拟合平均值和置信区间？对于后一种情况，geom_丝带似乎不起作用

---

这里有一个方法。首先，我们使用expand.grid将要预测的每个值都列成行，并且仅列这些值。这样可以避免重复

plot_data <- expand.grid(level = c("low", "high"), fertilizer=c("N","P"))
plot_data <- cbind(plot_data, predict(mod1, plot_data, interval="confidence"))

相关的：

plot_data <- expand.grid(level = c("low", "high"), fertilizer=c("N","P"))
plot_data <- cbind(plot_data, predict(mod1, plot_data, interval="confidence"))

ggplot(plot_data, aes(fertilizer, color=fertilizer)) + 
  geom_point(aes(y=fit)) + 
  geom_errorbar(aes(ymin=lwr, ymax=upr)) + 
  scale_color_manual(values=c("#E69F00", "#1B9E77")) + 
  facet_grid(cols=vars(level))