Bayesian modeling of the dependence in longitudinal data via partial autocorrelations and marginal variances

Many parameters and positive-definiteness are two major obstacles in estimating and modeling a correlation matrix for longitudinal data. In addition, when longitudinal data is incomplete, incorrectly modeling the correlation matrix often results in bias in estimating mean regression parameters. In this paper, we introduce a flexible and parsimonious class of regression models for a covariance matrix parameterized using marginal variances and partial autocorrelations. The partial autocorrelations can freely vary in the interval ( − 1 , 1 ) while maintaining positive definiteness of the correlation matrix so the regression parameters in these models will have no constraints. We propose a class of priors for the regression coefficients and examine the importance of correctly modeling the correlation structure on estimation of longitudinal (mean) trajectories and the performance of the DIC in choosing the correct correlation model via simulations. The regression approach is illustrated on data from a longitudinal clinical trial.