您将如何使用ANCOVA模型绘制CI

Question

我正在使用一个数据集，该数据集涉及特定鸟类的响应变量（质量），协变量（长度）和分类独立变量（位置）。此数据集如下所示：

在这种情况下，有三个地点（奥兰多，坦帕和另一个佛罗里达城），每个地点有10个观测值。我决定使用lm函数运行ANCOVA。

m1<- lm(mass.d ~ length.d + location)

我能够用ggplot绘制这个模型：

predm1<- predict(m1)
ggplot(YOUR_DATA, aes(length.d, mass.d, color = pop)) + geom_point() + 
geom_line(aes(y= predm1))

现在我想增加95％的置信区间。

我有置信区间，但我对如何将它们输入到我的ggplot图表中的想法很少，因为我基本上有3个不同位置的三个不同的斜率。我一直在玩geom_ribbon，但一直没能找到答案。