Multivariate Gaussian process regression for nonlinear modelling with colored noise