We consider the nonparametric regression problem with multiple predictors and an additive error, where function is assumed to be coordinatewise nondecreasing. propose a Bayesian approach make inferences on multivariate monotone function, obtain posterior contraction rate, construct universally consistent testing procedure for monotonicity. To facilitate analysis, we temporarily set aside shape ...