Nonstationary multivariate Gaussian processes for electronic health records

•Proposing a novel non-stationary multivariate Gaussian process model for EHRs.•Providing a computationally efficient and separable version of our model.•Illustrating promising fitting and prediction performance on synthetic and real data.•Discovering a statistically significant relationship between our model and a patient health metric.