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"