data("down")
fit.glm=glm(cases/births~age,weight=births,family=binomial,data=down)
fit.seg=segmented(fit.glm,seg.Z=~age,psi=25)
dev.new()
plot(down$age,log(down$cases/down$births),pch=20,cex=1.5,ylab="logit(cases/births)",xlab="Mother Age")
par(new=TRUE)
plot(fit.seg,col="red",lwd=2,ylab="",xlab="",ylim=c(min(log(down$cases/down$births)),max(log(down$cases/down$births))))
第一张图是自己画的,第二张图是示例。两张图的横坐标对比后发现,我绘制的图中横坐标上多出一组黑色的线段,每个黑色线段对应一个点,如果用:plot(fit.seg,add=TRUE,col="red",lwd=2,ylab="",xlab="") 就不会出现这些黑色线段了。所以想问一下在只绘制拟合结果(红色线)的时候怎样删掉它们?