我想为装有这样的 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 模型生成的波段。
谢谢!