Choosing the largest of several means can be a demanding problem, especially in the presence of a covariate. A hierarchical Bayesian approach to ranking and selection, as well as estimation of related means in the presence of a covariate, is considered. For the multiple slopes model we compute, in addition to the posterior means and standard deviations of the parameters, the posterior probabilities that each mean, at a given value of the covariate, is the largest. The vector of posterior probabilities thus obtained provides an easily understandable answer to the selection problem. Although calculation of the posterior probabilities may involve four-dimensional numerical integration in the difficult unbalanced design and unknown variance case, an efficient Monte Carlo method of evaluation has been developed and is given in the article. By reanalyzing a well-known data set on the breaking strength and thickness of starch films, we illustrate how our Bayesian approach produces meaningful conclusions, some of which would perhaps be difficult to obtain otherwise. For the starch film example, we found that it took only 1.4 seconds to compute the quantities of interest using an IBM 3090–600S machine. Because the computation time is quite small, it is apparent that the Bayesian procedure can be implemented for everyday use.