Я хотел бы создать доверительный диапазон для модели с такими gls, как это:
require(ggplot2)
require(nlme)
mp <-data.frame(year=c(1990:2010))
mp$wav <- rnorm(nrow(mp))*cos(2*pi*mp$year)+2*sin(rnorm(nrow(mp)*pi*mp$wav))+5
mp$wow <- rnorm(nrow(mp))*mp$wav+rnorm(nrow(mp))*mp$wav^3
m01 <- gls(wow~poly(wav,3), data=mp, correlation = corARMA(p=1))
mp$fit <- as.numeric(fitted(m01))
p <- ggplot(mp, aes(year, wow))+ geom_point()+ geom_line(aes(year,fit))
p
Это отображает только установленные значения и данные, и мне хотелось бы что-то в стиле
p <- ggplot(mp, aes(year, wow))+ geom_point()+ geom_smooth()
p
но с полосами, генерируемыми моделью gls.
Спасибо!