在准泊松回归图中插入置信区间

时间:2014-11-11 12:16:58

标签: r plot scatter-plot trendline

我使用以下代码进行了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

中的绘图功能来实现

由于

0 个答案:

没有答案