\begin{Schunk}
\begin{Sinput}
 par(mfcol=c(2,1))
 plot(xseq,gampred$fit,type="l",xlab="x",ylab="y", main="Gamma, link=identity")
 lines(xseq, gampred$fit+2*gampred$se,lty=2,col="green")
 lines(xseq, gampred$fit-2*gampred$se,lty=2,col="green")
 plot(xseq,lmpred$fit,type="l",xlab="x",ylab="y", main="Gaussian, link=identity")
 lines(xseq, lmpred$fit+2*lmpred$se,lty=2,col="green")
 lines(xseq, lmpred$fit-2*lmpred$se,lty=2,col="green")
 
\end{Sinput}
\end{Schunk}
