Two Stage Estimation for Generalized Structural Equation Models R package TSQLEM version 0.1.0 Provides a framework to estimate high dimensional generalized structural equation models using two stage quasi-likelihood expectation-maximization. The structural model supports binomial (logit and probit), Poisson, negative binomial, and gamma distributions for the outcome variable. Hattab (2026) "A Two Stage Quasi-Likelihood Estimation Method for High Dimensional Generalized Structural Equation Models" <doi:10.48550/arXiv.2608.16017>. Software Depends: R (>= 4.6.0) Imports: lavaan, Matrix, MASS, mvtnorm Suggests: psych Mohammad Hattab <mwhattab@unm.edu> Comprehensive R Archive Network (CRAN) Mohammad Hattab [aut, cre, cph] (ORCID: <https://orcid.org/0000-0003-0100-740X>) GPL (>= 3) 2026-09-28 application/tgz https://CRAN.R-project.org/package=TSQLEM doi:10.32614/CRAN.package.TSQLEM