I am trying to build a simple bernoulli_logit model for estimating model parameters with missing data, including missing categorical outcome data (choice). Since it is impossible to define missing data as int and since bernoulli_logit MUST have int as an outcome, my missing data formulation fails.
Is there a way around it?