Tolerance regions for random vectors and best linear predictors

In the case that vectors X and Y have a joint multivariate normal distribution, tolerance regions are found for the best linear predictor of Y using X if samples are used to estimate the regression coeffierante. Tolerance regions are also found for Y. In addition, simultaneous tolerance intervals for all linear functions of Y or of the best linear predictor of Y using X are found.