The Potts-Ising model for discrete multivariate data