\begin{Schunk}
\begin{Sinput}
 MLRModel <- "
  wjspl ~ edlevel + newschl + suspend + expelled + haveld + female + age
  wjspl ~ 1
  "
 output <- sem(model = MLRModel, data = dat)
 summary(output)
\end{Sinput}
\begin{Soutput}
lavaan (0.5-23.1097) converged normally after  58 iterations

                                                  Used       Total
  Number of observations                           313         322

  Estimator                                         ML
  Minimum Function Test Statistic                0.000
  Degrees of freedom                                 0

Parameter Estimates:

  Information                                 Expected
  Standard Errors                             Standard

Regressions:
                   Estimate  Std.Err  z-value  P(>|z|)
  wjspl ~                                             
    edlevel           1.162    0.324    3.584    0.000
    newschl           0.063    0.747    0.085    0.932
    suspend          -0.052    0.773   -0.067    0.947
    expelled         -2.758    1.098   -2.510    0.012
    haveld           -6.974    0.987   -7.063    0.000
    female            0.720    0.792    0.909    0.363
    age               0.412    0.207    1.994    0.046

Intercepts:
                   Estimate  Std.Err  z-value  P(>|z|)
   .wjspl            21.953    4.414    4.974    0.000

Variances:
                   Estimate  Std.Err  z-value  P(>|z|)
   .wjspl            39.904    3.190   12.510    0.000
\end{Soutput}
\end{Schunk}
