Multivariate normal integrals and contingency tables with ordered categories

Ak-dimensional multivariate normal distribution is made discrete by partitioning thek-dimensional Euclidean space with rectangular grids. The collection of probability integrals over the partitioned cubes is ak-dimensional contingency table with ordered categories. It is shown that loglinear model with main effects plus two-way interactions provides an accurate approximation for thek-dimensional table. The complete multivariate normal integral table is computed via the iterative proportional fitting algorithm from bivariate normal integral tables. This approach imposes no restriction on the correlation matrix. Comparisons with other numerical integration algorithms are reported. The approximation suggests association models for discretized multivariate normal distributions and contingency tables with ordered categories.The contingency-table approach occurred to me while I was collaborating with Paul Holland of the Educational Testing Service in 1985 on bivariate dependence functions. Holland maintains a belief that the continuous can learn from the discrete. This work is a reassertion of his claim.This research was sponsored by the National Science Council, Republic of China.