答案 0 :(得分:2)
你的原始代码中几乎已经有了它。使用interval = "confidence"
重复预测,并将这些行添加到您的绘图中。
# Right after you get the prediction intervals:
conf.int <- predict(Mod, interval = "confidence")
# Change the cbind line to this:
mydata <- data.frame(D3S, pred.int, conf.int) # checks names and adds .1 to dupes
# Right after you added prediction interval, do this:
myplot <- myplot +
geom_line(aes(y = lwr.1), color="red") +
geom_line(aes(y = upr.1), color="red")
现在像以前一样制作你的最终情节。