map2stan facilitates this form via the dmvnormNC density, which uses an internal Cholesky decomposition of the correlation matrix to build the varying effects. Here is the previous varying slopes model, now with the non-centered notation: f4nc <- alist( y ~ dnorm( mu , sigma ), mu <-...