\begin{Schunk}
\begin{Sinput}
 plot(x1, y, main="Ugly Poisson Data")
 library(rockchalk)
 newdat <- expand.grid(x1=plotSeq(dat$x1, length.out=50), x2=mean(dat$x2))
 newdat$p1 <- predict(glm1, newdata=newdat, type="response")
 lines(newdat$x1, newdat$p1, lwd=3, col="red")
 newdat$lmp1 <- predict(lm1, newdata=newdat)
 lines(newdat$x1, newdat$lmp1, lwd=3, col="green")
 legend("topleft", legend=c("Exp. Poisson", "OLS"), lwd=c(3,3), col=c("red","green"))
 
\end{Sinput}
\end{Schunk}
