\begin{Schunk}
\begin{Sinput}
 msem.gotBully.03 <- '
  Level: 1
  gotBully.w =~ gotBu2_i + gotBu4_i + gotBu5_i + gotBu7_i
  gotBully.w ~ Age
  
  Level: 2
  gotBu2_i ~~ gotBu2_i
  gotBu4_i ~~ gotBu4_i
  gotBu5_i ~~ gotBu5_i
  gotBu7_i ~~ gotBu7_i '
 ## Fitting the model
 msem.gotBully.03.fit <- 
  sem(model = msem.gotBully.03, data = hbsc.2, mimic = "Mplus", std.lv = TRUE, cluster = "schl_id")
 ## Requesting an estimation summary 
 summary(msem.gotBully.03.fit, fit.measures = TRUE, standardized = TRUE)
\end{Sinput}
\begin{Soutput}
lavaan 0.6-3 did NOT end normally after 3 iterations
** WARNING ** Estimates below are most likely unreliable

  Optimization method                           NLMINB
  Number of free parameters                         17

                                                  Used       Total
  Number of observations                           242         333
  Number of clusters [schl_id]                       8

  Estimator                                         ML
  Model Fit Test Statistic                          NA
  Degrees of freedom                                NA
  P-value                                           NA

Parameter Estimates:

  Information                                 Observed
  Observed information based on                Hessian
  Standard Errors                             Standard


Level 1 [within]:

Latent Variables:
                   Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
  gotBully.w =~                                                         
    gotBu2_i          0.874       NA                      0.875    0.675
    gotBu4_i          1.008       NA                      1.009    0.738
    gotBu5_i          0.842       NA                      0.843    0.722
    gotBu7_i          0.978       NA                      0.979    0.737

Regressions:
                   Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
  gotBully.w ~                                                          
    Age               0.035       NA                      0.035    0.050

Intercepts:
                   Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
   .gotBu2_i          0.000                               0.000    0.000
   .gotBu4_i          0.000                               0.000    0.000
   .gotBu5_i          0.000                               0.000    0.000
   .gotBu7_i          0.000                               0.000    0.000
   .gotBully.w        0.000                               0.000    0.000

Variances:
                   Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
   .gotBu2_i          0.917       NA                      0.917    0.545
   .gotBu4_i          0.850       NA                      0.850    0.455
   .gotBu5_i          0.654       NA                      0.654    0.479
   .gotBu7_i          0.805       NA                      0.805    0.457
   .gotBully.w        1.000                               0.997    0.997


Level 2 [schl_id]:

Intercepts:
                   Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
   .gotBu2_i          1.915       NA                      1.915    2.045
   .gotBu4_i          2.027       NA                      2.027    2.086
   .gotBu5_i          1.622       NA                      1.622    1.951
   .gotBu7_i          1.900       NA                      1.900    2.028

Variances:
                   Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
   .gotBu2_i          0.877       NA                      0.877    1.000
   .gotBu4_i          0.944       NA                      0.944    1.000
   .gotBu5_i          0.691       NA                      0.691    1.000
   .gotBu7_i          0.878       NA                      0.878    1.000
\end{Soutput}
\end{Schunk}
