Variational Inference of Joint Models using Multivariate Gaussian Convolution Processes