Simulates an ordinal relational matrix whose marginal category frequencies match those of an observed ordinal matrix.
Details
A latent Gaussian relational array is drawn from the social relations
model with mean EZ and within-dyad correlation rho via
simZ. The latent values are then discretized into ordinal
categories by cutting them at empirical quantiles chosen so that the
proportion of entries falling in each category reproduces the marginal
distribution of the observed categories of Y.