Two-stage Bayesian regression enforcing a fixed distribution (not Just Hierarchical regression)

[junpenglao]

This will become a stochastic likelihood and HMC does not differentiate over these stochastic node. You don’t need a custom sampler for this.