I’m presenting some work based on this model to an audience with lots of Stan users next week, so it would be nice to know whether this model would be possible to run in Stan, where IIUC it’s impossible to sample from discrete variables? But it seems like it would be very difficult then.