Model-based time-varying clustering of multivariate longitudinal data with covariates and outliers