Bivariate random effect model using skew‐normal distribution with application to HIV‐RNA

Correlated data arise in a longitudinal studies from epidemiological and clinical research. Random effects models are commonly used to model correlated data. Mostly in the longitudinal data setting we assume that the random effects and within subject errors are normally distributed. However, the normality assumption may not always give robust results, particularly if the data exhibit skewness. In this paper, we develop a Bayesian approach to bivariate mixed model and relax the normality assumption by using a multivariate skew-normal distribution. Specifically, we compare various potential models and illustrate the procedure using a real data set from HIV study.