Dear Community, I'm having a particularly nice problem: I need to do a Bayesian regression based on a mixture of a Gaussian and two half-normal distributions. The scale ranges from 0 - 100 and the "modes" of the half-normal distributions need to be fixed on these values, hence only the variance needs to be estimated. Can R do this, or do I need to code it in Matlab? Best, Jochen