mvord: An R Package for Fitting Multivariate Ordinal Regression Models

Abstract

The R package mvord implements composite likelihood estimation in the class of multivariate ordinal regression models with a multivariate probit and a multivariate logit link. A flexible modeling framework for multiple ordinal measurements on the same subject is set up, which takes into consideration the dependence among the multiple observations by employing different error structures. Heterogeneity in the error structure across the subjects can be accounted for by the package, which allows for covariate dependent error structures. In addition, different regression coefficients and threshold parameters for each response are supported. If a reduction of the parameter space is desired, constraints on the threshold as well as on the regression coefficients can be specified by the user. The proposed multivariate framework is illustrated by means of a credit risk application.

The link between attention deficit hyperactivity disorder (ADHD) symptoms and obesity-related traits: genetic and prenatal explanations

“…This method has less restrictive assumptions than the ordinary least squares regression and is invariant to the choice of transformation of the outcome [ 48 ]. Using a model for multiple outcomes allows a direct comparison of the effect size estimates of BMI PRS on inattention and hyperactivity symptoms while taking their mutual correlation into account [ 49 ]. We report the effect estimates of BMI PRS on the multivariate outcomes adjusted for sex and the first ten genetic PCs.…”

Multivariate ordinal probit analysis in the skeletal assessment of sex

“…We then fit a bivariate ordinal probit model to the ventral arc and subpubic concavity traits. We fit the model using the package “mvord” (Hirk, Hornik, Vana, & Genz, ) which uses composite likelihood estimation (Varin, Reid, & Firth, ). The estimated threshold parameters were 1.025, 1.807, 1.930, and 2.071 for VA and 1.606, 1.976, 2.000, and 2.185 for SPC.…”

New Approaches to Age Estimation Using Palatal Suture Fusion

“…The multivariate models were generated in order to provide the most mathematically correct approach to combining indicators in practice (39) as the palate aging indicators are biologically correlated to one another. The multivariate models use composite likelihood estimation (40,41), which differentiates them from the univariate models in this paper that employ maximum likelihood estimation. Consistent with Konigsberg et al (36,42,43), a Lagrange multiplier test for more than two categories (44–46) was applied to test the assumption of ordinal probit models that the unobserved error term is log‐normally distributed.…”