Bayesian isotonic regression and trend analysis.
In many applications, the mean of a response variable can be assumed to be a nondecreasing function of a continuous predictor, controlling for covariates. In such cases, interest often focuses on estimating the regression function, while also assessing evidence of an association. This article proposes a new framework for Bayesian isotonic regression and order-restricted inference. Approximating the regression function with a high-dimensional piecewise linear model, the nondecreasing constraint is incorporated through a prior distribution for the slopes consisting of a product mixture of point masses (accounting for flat regions) and truncated normal densities. To borrow information across the intervals and smooth the curve, the prior is formulated as a latent autoregressive normal process. This structure facilitates efficient posterior computation, since the full conditional distributions of the parameters have simple conjugate forms. Point and interval estimates of the regression function and posterior probabilities of an association for different regions of the predictor can be estimated from a single MCMC run. Generalizations to categorical outcomes and multiple predictors are described, and the approach is applied to an epidemiology application.
Neelon
Brian
B
Department of Biostatistics, CB no. 7420, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA.
Dunson
David B
DB
D001499:Bayes Theorem; D001699:Biometry; D003633:Dichlorodiphenyl Dichloroethylene / Q000633:toxicity; D004812:Epidemiologic Methods; D005260:Female; D006801:Humans; D007231:Infant, Newborn; D007306:Insecticides / Q000633:toxicity; D016014:Linear Models; D008390:Markov Chains; D009010:Monte Carlo Method; D011247:Pregnancy; D047928:Premature Birth / Q000139:chemically induced; D012044:Regression Analysis
D007306:Insecticides; D003633:Dichlorodiphenyl Dichloroethylene