There is the cdfquantreg R package that has likelihoods that put positive mass on 0 and / or 1, but it does not use Stan. I started to incorporate it into RStanArm but didn’t finish it. In the meantime you could transform your
(y * (n - 1) + .5) / n
to force it interior, where
n is the number of observations.