“…This renders FIML computationally infeasible when the number of ordinal variables is large. As a result, two-and three-stage limited-information least squares (3S-LS) estimation and testing theory have been proposed in the literature (Jöreskog, 1990(Jöreskog, , 1994Jöreskog & Sörbom, 1996;Lee et al, 1990bLee et al, , 1992Muthén, 1984;Satorra, 2000;Satorra & Bentler, 2010, 1988Asparouhov & Muthén, 2006, 2010 and implemented in software such as LISREL (Jöreskog & Sörbom, 1996), Mplus , EQS (Bentler, 2006), and the R package lavaan (Rosseel, 2012;Rosseel et al, 2012). Bayesian estimation methods of estimation, testing and model selection have also been developed (see e.g.…”