Although genome-wide expression data sets from multiple species are now more commonly generated, there have been few studies on how to best integrate this type of correlated data into models. Starting with a single-species, linear regression model that predicts transcription factor binding sites as a case study, we investigated how best to take into account the correlated expression data when extending this model to multiple species. Using a multivariate regression model, the phylogenetic relationships among the species were accounted for in two ways: 1) a repeated-measures model, where the error term is constrained, and 2) a Bayesian hierarchical model, where the prior distributions of the regression coefficients are constrained. We show that both multiple-species models improve predictive performance over the single-species model. When compared with each other, the repeated-measures model out-performed the Bayesian model. We suggest a possible explanation for the better performance of the model with the constrained error term.