Additive proportional odds model for ordinal data using Laplace P-splines. The combination of Laplace approximations and P-splines enable fast and flexible inference in a Bayesian framework. Specific approximations are proposed to account for the asymmetry in the marginal posterior distributions of non-penalized parameters. For more details, see Lambert and Gressani (2023) <doi:10.1177/1471082X231181173> ; Preprint: <arXiv:2210.01668>).

Calculates ordinated diet breadth with some plotting functions.

Implementation of cumulative link (mixed) models also known
    as ordered regression models, proportional odds models, proportional
    hazards models for grouped survival times and ordered logit/probit/...
    models. Estimation is via maximum likelihood and mixed models are fitted
    with the Laplace approximation and adaptive Gauss-Hermite quadrature.
    Multiple random effect terms are allowed and they may be nested, crossed or
    partially nested/crossed. Restrictions of symmetry and equidistance can be
    imposed on the thresholds (cut-points/intercepts). Standard model
    methods are available (summary, anova, drop-methods, step,
    confint, predict etc.) in addition to profile methods and slice
    methods for visualizing the likelihood function and checking
    convergence.

Provides a function for fitting various penalized Bayesian
    cumulative link ordinal response models when the number of parameters
    exceeds the sample size. These models have been described in 
    Zhang and Archer (2021) <doi:10.1186/s12859-021-04432-w>.

A regression framework for response variables which are continuous
    self-rating scales such as the Visual Analog Scale (VAS) used in pain
    assessment, or the Linear Analog Self-Assessment (LASA) scales in quality
    of life studies. These scales measure subjects' perception of an intangible
    quantity, and cannot be handled as ratio variables because of their inherent
    non-linearity. We treat them as ordinal variables, measured on a continuous
    scale. A function (the g function) connects the scale with an underlying
    continuous latent variable. The link function is the inverse of the CDF of the
    assumed underlying distribution of the latent variable. A variety of
    link functions are currently implemented. Such models are described in Manuguerra et al (2020) <doi:10.18637/jss.v096.i08>.