\begin{Schunk}
\begin{Sinput}
 agam <- glm(yobs ~ I(1/xseq), family = Gamma(link = "identity"), control = glm.control(maxit=100),start=c(1,65))
 library(MASS)
 myshape <- gamma.shape(agam)
 gampred <- predict(agam,type="response", se=T, dispersion=1/myshape$alpha)
 
\end{Sinput}
\end{Schunk}
