Multivariate time series prediction