\begin{Schunk}
\begin{Sinput}
 mlm.01 <- lmer(gotBu_ave ~ 1 + Age + (1|schl_id), data = hbsc.2, REML = FALSE)
 summary(mlm.01)
\end{Sinput}
\begin{Soutput}
Linear mixed model fit by maximum likelihood  ['lmerMod']
Formula: gotBu_ave ~ 1 + Age + (1 | schl_id)
   Data: hbsc.2

     AIC      BIC   logLik deviance df.resid 
   609.8    623.7   -300.9    601.8      234 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-1.5547 -0.7147 -0.2137  0.3208  4.0431 

Random effects:
 Groups   Name        Variance Std.Dev.
 schl_id  (Intercept) 0.2124   0.4608  
 Residual             0.6851   0.8277  
Number of obs: 238, groups:  schl_id, 8

Fixed effects:
            Estimate Std. Error t value
(Intercept)  1.12275    0.77532   1.448
Age          0.05864    0.05962   0.984

Correlation of Fixed Effects:
    (Intr)
Age -0.973
\end{Soutput}
\end{Schunk}
