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