Вставьте доверительный интервал в график квази-пуассоновской регрессии

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"))

Но я специально хочу сделать это с помощью функции plot в R

Спасибо