\begin{Schunk}
\begin{Sinput}
 summary(agam2, dispersion=1/myshape2$alpha)
\end{Sinput}
\begin{Soutput}
Call:
glm(formula = yobs2 ~ I(1/xseq), family = Gamma(link = "inverse"), 
    start = c(2, 4), control = glm.control(maxit = 100))

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.5708  -0.8150  -0.2532   0.3325   2.5156  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)  0.26736    0.03352   7.976 1.51e-15 ***
I(1/xseq)    2.80621    0.34735   8.079 6.54e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for Gamma family taken to be 0.6756965)

    Null deviance: 654.43  on 799  degrees of freedom
Residual deviance: 599.09  on 798  degrees of freedom
AIC: 2479.4

Number of Fisher Scoring iterations: 7
\end{Soutput}
\begin{Sinput}
 
\end{Sinput}
\end{Schunk}
