Normality of random effects and error terms is a routine assumption for linear effects mixed models. However, such assumption may be unrealistic, obscuring important features of withinand among-unit variation. A simple and robust Bayesian parametric approach that relaxes this assumption by using a multivariate Skew elliptical distributions, which includes the Skew-t, Skew-normal, t-Student, and...