Nonparametric forecasting of multivariate probability density functions