I am trying to do a MAP estimation of the model below. Because the model function has a normalcdf, the outcome of the model is in the range $0-1$. The observations are bound to the same range. However, using a Gaussian distribution for the likelihood means that the uncertainty bounds are sometimes outside the $0-1$ range. Statistically speaking, is this a problem, and if so, what would be a better solution?

$R_i \sim \mathcal{N}(\mu,\sigma)$

$\mu_i = normalcdf(\frac{1}{\alpha} * log(A/\beta) )$

$\alpha \sim \mathcal{N}(0.35, 0.1) $

$\beta\sim \mathcal{N}(70, 10) $

$\sigma \sim \mathcal{U}(0, 0.3) $