We present an efficient and accurate method to evaluate multivariate normal probabilities with arbitrary singular correlation matrices. The new method is applied to the construction of simultaneous confidence intervals and simultaneous all pairwise confidence intervals for multinomial proportions when the sample size is sufficiently large.
Key Words: multivariate normal, singular distribution, numerical integration, statistical computation.