Joint segmentation of multivariate Gaussian processes using mixed linear models