Simultaneous probability statements for Bayesian P-splines
P-splines are a popular approach for fitting nonlinear effects of continuous covariates in semiparametric regression models. Recently, a Bayesian version for P-splines has been developed on the basis of Markov chain Monte Carlo simulation techniques for inference. In this work we adopt and generalize the concept of Bayesian contour probabilities to additive models with Gaussian or multicategorical responses. More specifically, we aim at computing the maximum credible level (sometimes called Bayesian p-value) for which a particular parameter vector of interest lies within the corresponding highest posterior density (HPD) region. We are particularly interested in parameter vectors that correspond to a constant, linear or more generally a polynomial fit. As an alternative to HPD regions simultaneous credible intervals could be used to define pseudo contour probabilities. Efficient algorithms for computing contour and pseudo contour probabilities are developed. The performance of the approach is assessed through simulation studies. Two applications on the determinants of undernutrition in developing countries and the health status of trees show how contour probabilities may be used in practice to assist the analyst in the model building process.
|Date of creation:||May 2007|
|Date of revision:|
|Contact details of provider:|| Postal: |
Web page: http://www.uibk.ac.at/fakultaeten/volkswirtschaft_und_statistik/index.html.en
More information through EDIRC
Please report citation or reference errors to , or , if you are the registered author of the cited work, log in to your RePEc Author Service profile, click on "citations" and make appropriate adjustments.:
- Andreas Brezger & Thomas Kneib & Stefan Lang, . "BayesX: Analyzing Bayesian Structural Additive Regression Models," Journal of Statistical Software, American Statistical Association, vol. 14(i11).
- David J. Spiegelhalter & Nicola G. Best & Bradley P. Carlin & Angelika van der Linde, 2002. "Bayesian measures of model complexity and fit," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 64(4), pages 583-639.
- Marx, Brian D. & Eilers, Paul H. C., 1998. "Direct generalized additive modeling with penalized likelihood," Computational Statistics & Data Analysis, Elsevier, vol. 28(2), pages 193-209, August.
- John Maindonald, . "Generalized Additive Models: An Introduction with R," Journal of Statistical Software, American Statistical Association, vol. 16(b03).
- Ludwig Fahrmeir & Stefan Lang, 2001. "Bayesian Semiparametric Regression Analysis of Multicategorical Time-Space Data," Annals of the Institute of Statistical Mathematics, Springer, vol. 53(1), pages 11-30, March.
- Chib S. & Jeliazkov I., 2001. "Marginal Likelihood From the Metropolis-Hastings Output," Journal of the American Statistical Association, American Statistical Association, vol. 96, pages 270-281, March.
When requesting a correction, please mention this item's handle: RePEc:inn:wpaper:2007-08. See general information about how to correct material in RePEc.
For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Janette Walde)
If references are entirely missing, you can add them using this form.