Multivariate Gaussian Process Model for Correlated Time Series in R