Summary This article proposes a two‐stage simultaneous confidence procedure for the comparisons of k pairs of population means, without using multiplicity adjustment of more than two populations. The proposed procedure can be broadly applied to parametric or nonparametric models. It is robust and versatile because its derivation only utilizes a partitioning approach in conjunction with a bivariate adjustment, without any assumption on the underlying distribution. To elucidate the application, the proposed procedure is intertwined with the estimation of the therapeutic window of a drug. It provides confidence limits for the efficacy and the toxicity of the effective doses, highest ineffective dose, safe doses, and lowest unsafe dose, simultaneously. Such estimation information facilitates follow‐up studies in clinical trials. As an illustrative example, the new procedure is applied to analyze a data set on molecular cancer therapeutics regarding the apoptotic killing effects of different chemical compounds on two leukemia cell lines.