\begin{Schunk}
\begin{Sinput}
 summary(agam,dispersion=1/myshape$alpha)
\end{Sinput}
\begin{Soutput}
Call:
glm(formula = yobs ~ I(1/xseq), family = Gamma(link = "identity"), 
    start = c(1, 65), control = glm.control(maxit = 100))

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.5702  -0.8164  -0.2555   0.3364   2.5154  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)   0.7738     0.4524    1.71   0.0872 .  
I(1/xseq)    67.7209     5.4534   12.42   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for Gamma family taken to be 0.6757266)

    Null deviance: 786.66  on 799  degrees of freedom
Residual deviance: 599.11  on 798  degrees of freedom
AIC: 4756

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