Hi,

I was wondering how to fit a two part multilevel model using gsem. The first part is a binary multilevel model (fit for the probability of observing a positive-versus-zero outcome) and the second is a multinomial multilevel model (with gsem).
How can I now combine these two parts and calculate marginal effects? I guess the twopm command is not appropriate here.

Thank you already in advance.