A Bayesian mixed-effects statistical model where the random effect is parameterized by a singular value decomposition (SVD). The matrices of the SVD representation are estimated using a Bayesian framework, enabling uncertainty quantification of the singular values and vectors in the decomposition.