相关疑难解决方法(0)

我想为一个配有这样的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模型生成的波段.

谢谢!

我手动创建了一个预期寿命数据集,伴随着95%的置信区间.我在时间尺度上绘制这些图,但更喜欢用阴影而不是虚线.代码显示:

p1 = ggplot()
p2 = p1 + geom_line(aes(x=pl$Time, y=pl$menle), colour="blue")
p3 = p2 + geom_line(aes(x=pl$Time, y=pl$menlelb), colour="blue", lty="dotted")
p4 = p3 + geom_line(aes(x=pl$Time, y=pl$menleub), colour="blue", lty="dotted")

是否有一种简单的方法来遮蔽间隔而不仅仅是线条？如果我遗漏了一些简单的事情,我会事先道歉,但我找不到任何表明这样做的简单方法.