Bayesian Inference on Order‐Constrained Parameters in Generalized Linear Models

In biomedical studies, there is often interest in assessing the association between one or more ordered categorical predictors and an outcome variable, adjusting for covariates. For a k-level predictor, one typically uses either a k-1 degree of freedom (df) test or a single df trend test, which requires scores for the different levels of the predictor. In the absence of knowledge of a parametric form for the response function, one can incorporate monotonicity constraints to improve the efficiency of tests of association. This article proposes a general Bayesian approach for inference on order-constrained parameters in generalized linear models. Instead of choosing a prior distribution with support on the constrained space, which can result in major computational difficulties, we propose to map draws from an unconstrained posterior density using an isotonic regression transformation. This approach allows flat regions over which increases in the level of a predictor have no effect. Bayes factors for assessing ordered trends can be computed based on the output from a Gibbs sampling algorithm. Results from a simulation study are presented and the approach is applied to data from a time-to-pregnancy study.

[1]  L. Wasserman,et al.  Computing Bayes Factors by Combining Simulation and Asymptotic Approximations , 1997 .

[2]  Govind S. Mudholkar,et al.  A Simple Approach to Testing Homogeneity of Order-Constrained Means , 1993 .

[3]  C. Tietze Fertility after discontinuation of intrauterine and oral contraception. , 1968, International journal of fertility.

[4]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[5]  Joseph G. Ibrahim,et al.  Monte Carlo Methods in Bayesian Computation , 2000 .

[6]  J. Albert A Bayesian Approach to Estimation of GPAs of University of Iowa Freshmen Under Order Restrictions , 1994 .

[7]  F. T. Wright,et al.  Order restricted statistical inference , 1988 .

[8]  P. Diggle,et al.  Additive isotonic regression models in epidemiology. , 2000, Statistics in medicine.

[9]  A. Gelfand,et al.  Nonparametric Bayesian bioassay including ordered polytomous response , 1991 .

[10]  D. J. Bartholomew,et al.  A Test of Homogeneity of Means Under Restricted Alternatives , 1961 .

[11]  Peter Bacchetti,et al.  Additive Isotonic Models , 1989 .

[12]  Reduced fertility among women employed as dental assistants exposed to high levels of nitrous oxide. , 1992 .

[13]  Shyamal D. Peddada,et al.  Confidence Interval Estimation Subject to Order Restrictions , 1994 .

[14]  Christian P. Robert,et al.  Maximum likelihood estimation under order restrictions by the prior feedback method , 1996 .

[15]  Audris Mockus,et al.  A nonparametric Bayes method for isotonic regression , 1995 .

[16]  Adrian F. M. Smith,et al.  Bayesian Analysis of Constrained Parameter and Truncated Data Problems , 1991 .

[17]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[18]  Hari Mukerjee,et al.  Order-restricted inferences in linear regression , 1995 .

[19]  C. Weinberg,et al.  Reduced fertility among women employed as dental assistants exposed to high levels of nitrous oxide. , 1992, The New England journal of medicine.

[20]  S. Chib,et al.  Bayesian analysis of binary and polychotomous response data , 1993 .

[21]  A. Agresti,et al.  Order-restricted tests for stratified comparisons of binomial proportions. , 1996, Biometrics.

[22]  Adrian F. M. Smith,et al.  Bayesian Inference for Generalized Linear and Proportional Hazards Models Via Gibbs Sampling , 1993 .

[23]  D D Baird,et al.  Incidence of early loss of pregnancy. , 1988, The New England journal of medicine.

[24]  Purushottam W. Laud,et al.  Nonparametric Bayesian bioassay with prior constraints on the shape of the potency curve , 1993 .

[25]  Ming-Hui Chen,et al.  Monte Carlo methods for Bayesian analysis of constrained parameter problems , 1998 .

[26]  Andrew Thomas,et al.  WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility , 2000, Stat. Comput..