Dear everyone,


I am looking for something similar to Heckman selection model/svysemlog with a modification.

I have a selectiion variable with two values (0 and 1) in the first step, and a mulitnomial non-ordinal categorical variable (with six categories) in the second step.
I am interested only in positive (1) values in the first step (around 30% of the total sample).

What I did at the first place was a. logit analysis for the first step b. multinomial logit for the second step. However, I was advised to use the Heckman selection model for multiple reasons.
However, if I am not mistaken, Heckman (and svysemlog) cannot be used if the outcome variable is a non-ordinal variable.


I have two questions:


a. Is there any Stata package that adresses my problem?
b. Do you have any advice how to proceed, in case there is no ready-made solution in Stata?


Thanks in advance!