\begin{Schunk}
\begin{Sinput}
 agam2<-glm(yobs2~ I(1/xseq), family=Gamma(link="inverse"), control=glm.control(maxit=100),start=c(2,4))
 library(MASS)
 myshape2 <- gamma.shape(agam2)
 gampred2 <- predict(agam2,type="response", se=T, dispersion=1/myshape2$alpha)
 
\end{Sinput}
\end{Schunk}
