\begin{Schunk}
\begin{Sinput}
 glm1 <- glm(y ~ x1 + x2, data=dat, family=poisson(link=log))
 summary(glm1)
\end{Sinput}
\begin{Soutput}
Call:
glm(formula = y ~ x1 + x2, family = poisson(link = log), data = dat)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.1831  -0.6880   0.1172   0.7351   2.5707  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)  9.978e-02  2.032e-04     491   <2e-16 ***
x1           7.500e-01  7.970e-06   94100   <2e-16 ***
x2          -1.500e-01  5.053e-06  -29686   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for poisson family taken to be 1)

    Null deviance: 2.0717e+10  on 98  degrees of freedom
Residual deviance: 1.0285e+02  on 96  degrees of freedom
  (1 observation deleted due to missingness)
AIC: 1664.3

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