I’m certainly no expert on MCMC, but I suspect that HMC requires everything in your model to be differentiable, but the Integer components are not. What happens if you choose another sampler, like Metropolis Hastings?
2 Likes