We consider model-based clustering methods for continuous, correlated data that account for external information available in the presence of mixed-type fixed covariates by proposing the MoEClust suite of models. These allow covariates influence the component weights and/or component densities by modelling the parameters of the mixture as functions of the covariates. A familiar range of constrained eigen-decomposition parameterisations of the component covariance matrices are also accommodated. This paper thus addresses the equivalent aims of including covariates in Gaussian Parsimonious Clustering Models and incorporating parsimonious covariance structures into the Gaussian mixture of experts framework. The MoEClust models demonstrate significant improvement from both perspectives in applications to univariate and multivariate data sets.Keywords: Model-based clustering, mixtures of experts, EM algorithm, parsimony, multivariate response, covariates arXiv:1711.05632v2 [stat.ME] 10 Dec 2018 cluster membership indicator vector, where z ig = 1 if observation i belongs to cluster g and z ig = 0 otherwise, the first approach assumes that z i affects the distribution of x i . In probabilistic terms, this means to replace the actual group-specific conditional dis-The name 'cluster-weighted model' (CWM) is frequently given to this approach, e.g. Dang et al. (2017) and Ingrassia et al. (2015); the latter provides a recent extension allowing for mixed-type covariates. The use of the alternative term 'mixtures of regressions with random covariates' to describe CWMs (e.g. Hennig (2000)), provides opportunity to clarify that the remainder of this paper focuses on the second approach, with fixed covariates affecting cluster membership via f (yThis is achieved using the mixture of experts (MoE) paradigm (Jacobs et al., 1991), in which the parameters of the mixture are modelled as functions of fixed, potentially mixed-type covariates. However, for multivariate, continuous, correlated responses, a unifying framework combining the special cases of the Gaussian MoE model with the flexibility afforded by the various parsimonious covariance parameterisations within the Gaussian Parsimonious Clustering Models (GPCM) family of finite mixture models (Banfield & Raftery, 1993;Celeux & Govaert, 1995) has to date been lacking. Indeed, the main contribution of this paper is in addressing the aim of incorporating covariates into the GPCM family and the equivalent aim of incorporating GPCM covariance constraints into the Gaussian MoE framework, by proposing the MoEClust family of models: the name comes from the interest in employing MoE models chiefly for clustering purposes. From both perspectives, MoEClust models demonstrate significant improvement in applications to univariate and multivariate response data.A software implementation for the full suite of MoEClust models is provided by the associated R package MoEClust (Murphy & Murphy, 2018), which is available from www.r-project.org (R Core Team, 2018), with which all results were obtaine...