Multivariate Bayesian function estimation