\begin{Schunk}
\begin{Sinput}
 lmmod3 <- glm(yobs2 ~ I(1/xseq), family=gaussian(link="inverse"))
 summary(lmmod3)
\end{Sinput}
\begin{Soutput}
Call:
glm(formula = yobs2 ~ I(1/xseq), family = gaussian(link = "inverse"))

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.3914  -1.0866  -0.3662   0.5736  11.4282  

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  0.26355    0.03979   6.623 6.47e-11 ***
I(1/xseq)    2.85141    0.50428   5.654 2.18e-08 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for gaussian family taken to be 2.702469)

    Null deviance: 2320.6  on 799  degrees of freedom
Residual deviance: 2156.5  on 798  degrees of freedom
AIC: 3069.6

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