Variational Learning on Aggregate Outputs with Gaussian Processes