On estimation of multivariate prediction regions in partial least squares regression