\begin{Schunk}
\begin{Sinput}
 sem.01.ord <- '
  ## the measurement model
  gotBully =~ NA*gotBu1_o + gotBu2_o + gotBu3_o
               + gotBu4_o + gotBu5_o + gotBu6_o
               + gotBu7_o + gotBu8_o + gotBu9_o
  gotBully ~~ 1*gotBully
  alcohol =~ NA*alc1_o + alc2_o + alc3_o
              + alc4_o + alc5_o
  alcohol ~~ 1*alcohol
  # regress 
  alcohol ~ gotBully
  '
 ## Run on each imputed data set
 sem.01.ord.fits <- list()
 for(i in seq_along(imps)){
      sem.01.ord.fits[[i]] <- sem(model = sem.01.ord, mimic = "Mplus",
                        estimator = "DWLS", data = imps[[i]])
  }
 
\end{Sinput}
\end{Schunk}
