Good Evening,

I am new to SEM in stata. I have the following path diagram, in my empirical dataset, I have count data, binary and continuous endogenous variables.
in my exogenous variables, there are group identity I want to control in the main equation. (Y = endo1 + endo2 + endo3 + endo4 + i.exo1 + exo2 + exo3 + exo4)
what is the best command that can estimate this model? Do gsem and gllamm both work for my dataset?
I would really appreciate if anyone can help point me to a stata code example which is similar to the model diagram.

Thank you in advance.