In the setting of inference with two-step monotone incomplete data drawn from N d (µ, ∑), a multivariate normal population with mean µ and covariance matrix ∑, we derive a stochastic representation for the exact distribution of a generalization of Hotelling's T 2 -statistic, thereby enabling the construction of exact level ellipsoidal confidence regions for µ. By applying the equivariance of μ andˆ, Σ the maximum likelihood estimators of µ and ∑, respectively, we show that the T 2 -statistic is invariant under affine transformations. Further, as a consequence of the exact stochastic representation, we derive upper and lower bounds for the cumulative distribution function of the T 2 -statistic. We apply these results to construct simultaneous confidence regions for linear combinations of µ, and we apply these results to analyze a dataset consisting of cholesterol measurements on a group of Pennsylvania heart disease patients.