Abstract:A vast literature in statistics, biometrics, and econometrics is concerned with the analysis of binary and polychotomous response data. The classical approach fits a categorical response regression model using maximum likelihood, and inferences about the model are based on the associated asymptotic theory. The accuracy of classical confidence statements is questionable for small sample sizes. In this article, exact Bayesian methods for modeling categorical response data are developed using the idea of data augmentation. The general approach can be summarized as follows. The probit regression model for binary outcomes is seen to have an underlying normal regression structure on latent continuous data. Values of the latent data can be simulated from suitable truncated normal distributions. If the latent data are known, then the posterior distribution of the parameters can be computed using standard results for normal linear models. Draws from this posterior are used to sample new latent data, and the process is iterated with Gibbs sampling. This data augmentation approach provides a general framework for analyzing binary regression models. It leads to the same simplification achieved earlier for censored regression models. Under the proposed framework, the class of probit regression models can be enlarged by using mixtures of normal distributions to model the latent data. In this normal mixture class, one can investigate the sensitivity of the parameter estimates to the choice of “link function,” which relates the linear regression estimate to the fitted probabilities. In addition, this approach allows one to easily fit Bayesian hierarchical models. One specific model considered here reflects the belief that the vector of regression coefficients lies on a smaller dimension linear subspace. The methods can also be generalized to multinomial response models with J > 2 categories. In the ordered multinomial model, the J categories are ordered and a model is written linking the cumulative response probabilities with the linear regression structure. In the unordered multinomial model, the latent variables have a multivariate normal distribution with unknown variance-covariance matrix. For both multinomial models, the data augmentation method combined with Gibbs sampling is outlined. This approach is especially attractive for the multivariate probit model, where calculating the likelihood can be difficult.

Bayesian analysis of linear regression models with autoregressive symmetrical errors and incomplete data

Bayesian variable selection in linear regression models with non-normal errors

Bayesian diagnostics in a partially linear model with first-order autoregressive skew-normal errors

Bayesian Regression Analysis of Data with Random Effects Covariates from Nonlinear Longitudinal Measurements

A Bayesian hybrid method for the analysis of generalized linear models with missing-not-at-random covariates

Robust Bayesian model selection for heavy-tailed linear regression using finite mixtures

Bayesian estimation for longitudinal data in a joint model with HPCs

Bayesian analysis for nonlinear regression model under skewed errors, with application in growth curves

Robust Bayesian Inference for Berkson and Classical Measurement Error Models

Bayesian Analysis of Longitudinal Dyadic Data with Informative Missing Data Using a Dyadic Shared-Parameter Model

Bayesian Variance Changepoint Detection in Linear Models with Symmetric Heavy-Tailed Errors

Bayesian Restricted Likelihood Methods: Conditioning on Insufficient Statistics in Bayesian Regression

Censored autoregressive regression models with Student‐t innovations

Semiparametric Bayesian Analysis for Longitudinal Mixed Effects Models with Non-Normal AR(1) Errors

Bayesian Latent-Class Mixed-Effect Hybrid Models for Dyadic Longitudinal Data with Non-Ignorable Dropouts.

Bayesian Nonparametric Regression Analysis of Data with Random Effects Covariates from Longitudinal Measurements

Bayesian estimation of simultaneous equation model with lagged endogenous variables and first order serially correlated errors

Robust Bayesian Model Averaging for Linear Regression Models With Heavy-Tailed Errors

Bayesian Analysis of Binary and Polychotomous Response Data

Diagnostic analystics in the Bayesian vector autoregressive model

Full Bayesian analysis of double seasonal autoregressive models with real applications