Cumulative Link Models for Ordinal Regression with the R Package ordinal

This paper introduces the R-package ordinal for the analysis of ordinal data using cumulative link models. The model framework implemented in ordinal includes partial proportional odds, structured thresholds, scale effects and flexible link functions. The package also support cumulative link models with random effects which are covered in a future paper. A speedy and reliable regularized Newton estimation scheme using analytical derivatives provides maximum likelihood estimation of the model class. The paper describes the implementation in the package as well as how to use the functionality in the package for analysis of ordinal data including topics on model identifiability and customized modelling. The package implements methods for profile likelihood confidence intervals, predictions of various kinds as well as methods for checking the convergence of the fitted models.

[1]  F. Harrell,et al.  Partial Proportional Odds Models for Ordinal Response Variables , 1990 .

[2]  Hans Bruun Nielsen,et al.  Introduction to Numerical Computation--Analysis and MATLAB Illustrations , 2004 .

[3]  J. Pratt Concavity of the Log Likelihood , 1981 .

[4]  Per B. Brockhoff,et al.  lmerTest Package: Tests in Linear Mixed Effects Models , 2017 .

[5]  P. Brockhoff,et al.  Sensometrics: Thurstonian and Statistical Models , 2012 .

[6]  J. Randall,et al.  The Analysis of Sensory Data by Generalized Linear Model , 1989 .

[7]  P. McCullagh Regression Models for Ordinal Data , 1980 .

[8]  F. C. Genter,et al.  Goodness‐of‐link testing in ordinal regression models , 1985 .

[9]  Alan Agresti,et al.  Categorical Data Analysis , 1991, International Encyclopedia of Statistical Science.

[10]  J. Burridge,et al.  A Note on Maximum Likelihood Estimation for Regression Models using Grouped Data , 1981 .

[11]  Y. Pawitan In all likelihood : statistical modelling and inference using likelihood , 2002 .

[12]  C. Cox,et al.  Location-scale cumulative odds models for ordinal data: a generalized non-linear model approach. , 1995, Statistics in medicine.

[13]  Paul-Christian Bürkner,et al.  brms: An R Package for Bayesian Multilevel Models Using Stan , 2017 .

[14]  T. Yee The VGAM Package for Categorical Data Analysis , 2010 .

[15]  F. J. Aranda-Ordaz An extension of the proportional-hazards model for grouped data. , 1983, Biometrics.

[16]  Vern T. Farewell,et al.  A Study of Distributional Shape in Life Testing , 1977 .

[17]  W. Greene,et al.  Modeling Ordered Choices: A Primer , 2010 .

[18]  Brian D. Ripley,et al.  Modern applied statistics with S, 4th Edition , 2002, Statistics and computing.

[19]  S. Rabe-Hesketh,et al.  Generalized multilevel structural equation modeling , 2004 .