我使用以下代码进行了quasipoisson回归:
plot(ori_brk ~ year_imd,pch=16,col="darkgreen",cex.lab=1.2,cex.axis=1.2,xlab="Year",ylab="No. of break days",ylim=c(0,30),cex=1.2)
并使用以下代码插入趋势线:
xv<-seq(1951,2007,1)
yv<-predict(ori_brk_qpmdl,list(Distance=xv))
lines(xv,exp(yv),col="darkgreen",lwd=2)
现在我需要在趋势线周围插入95%的置信区间。
我知道我可以使用以下代码使用ggplot2:
ggplot(ori,aes(x=year,y=brk))+ylim(0,35)+stat_smooth(method='glm',family='quasipoisson',lty="dotted") +geom_point()+theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),
panel.background = element_blank(), axis.line = element_line(colour = "black"))
但我特别希望使用R
中的绘图功能来实现由于