For an ordinal response scale, the multinomial variables are simulated either under a marginal cumulative link model (rmult.clm
) or under a marginal continuation ratio model (rmult.crm
).
For a nominal response scale, the multinomial responses are simulated under a marginal baseline category logit model (rmult.bcl
).
Correlated binary responses are simulated using the function rbin
.
See vignette for more details about the threshold approaches and the implemented marginal models.
Li, S. T. and Hammond, J. L. (1975) Generation of pseudorandom numbers with specified univariate distributions and correlation coefficients. IEEE Transacions on Systems, Man and Cybernetics, 5, 557--561.
McCullagh, P. (1980) Regression models for ordinal data. Journal of the Royal Statistical Society B, 42, 109--142.
McFadden, D. (1974) Conditional logit analysis of qualitative choice behavior. New York: Academic Press, 105--142.
Touloumis, A., Agresti, A. and Kateri, M. (2013) GEE for multinomial responses using a local odds ratios parameterization. Biometrics, 69, 633--640.
Tutz, G. (1991) Sequential models in categorical regression. Computational Statistics & Data Analysis, 11, 275--295.