A different, also very explicit way to go about this type of problem, that also generalizes fully, is to use a Bayesian hierarchical model of a dirichlet proces

A different, also very explicit way to go about this type of problem, that also generalizes fully, is to use a Bayesian hierarchical model of a dirichlet process and sub-isotonic regressions. Gelman et al have written a lot about this, and they have a proposed general workflow [1] [1] Bayesian Workflow https://arxiv.org/abs/2011.01808 https://arxiv.org/abs/2011.01808