[R] Bayesian Regression with half-normal distributions

[Jochen Mattes]

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