r - 準ポアソン回帰プロットに信頼区間を挿入

Question

次のコードを使用して準ポアソン回帰を実行しました。

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のプロット関数を使用してそれを行うことを特に検討しています

ありがとう