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.

[1]  G. C. Tiao,et al.  Bayesian inference in statistical analysis , 1973 .

[2]  Luke Tierney,et al.  A Space-Efficient Recursive Procedure for Estimating a Quantile of an Unknown Distribution , 1983 .

[3]  F. O’Sullivan A Statistical Perspective on Ill-posed Inverse Problems , 1986 .

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

[5]  Kerrie Mengersen,et al.  [Bayesian Computation and Stochastic Systems]: Rejoinder , 1995 .

[6]  J. Besag,et al.  Bayesian Computation and Stochastic Systems , 1995 .

[7]  H. Pruscha,et al.  Der Einfluß von Bestandeskenngrößen, Topographie, Standort und Witterung auf die Entwicklung des Kronenzustandes im Bereich des Forstamtes Rothenbuch , 1996, Forstwissenschaftliches Centralblatt vereinigt mit Tharandter forstliches Jahrbuch.

[8]  Paul H. C. Eilers,et al.  Flexible smoothing with B-splines and penalties , 1996 .

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

[10]  Paul H. C. Eilers,et al.  Direct generalized additive modeling with penalized likelihood , 1998 .

[11]  L. Fahrmeir,et al.  Bayesian Semiparametric Regression Analysis of Multicategorical Time-Space Data , 2001 .

[12]  Ludwig Fahrmeir,et al.  Semiparametric Analysis of the Socio-Demographic and Spatial Determinants of Undernutrition in Two African Countries , 2001 .

[13]  S. Chib,et al.  Marginal Likelihood From the Metropolis–Hastings Output , 2001 .

[14]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[15]  Eric R. Ziegel,et al.  The Elements of Statistical Learning , 2003, Technometrics.

[16]  S. Lang,et al.  Bayesian P-Splines , 2004 .

[17]  Leonhard Held,et al.  Simultaneous Posterior Probability Statements From Monte Carlo Output , 2004 .

[18]  T. Kneib,et al.  BayesX: Analyzing Bayesian Structural Additive Regression Models , 2005 .

[19]  S. Wood ON CONFIDENCE INTERVALS FOR GENERALIZED ADDITIVE MODELS BASED ON PENALIZED REGRESSION SPLINES , 2006 .

[20]  Andreas Brezger,et al.  Generalized structured additive regression based on Bayesian P-splines , 2006, Comput. Stat. Data Anal..

[21]  C. Holmes,et al.  Bayesian auxiliary variable models for binary and multinomial regression , 2006 .

[22]  G. D. Lin,et al.  THE HURWITZ ZETA DISTRIBUTION , 2006 .

[23]  S. Wood Generalized Additive Models: An Introduction with R , 2006 .

[24]  Alan Y. Chiang,et al.  Generalized Additive Models: An Introduction With R , 2007, Technometrics.

[25]  Paul A. Raschky,et al.  Working Papers in Economics and Statistics Charity Hazard -a Real Hazard to Natural Disaster Insurance University of Innsbruck Working Papers in Economics and Statistics Charity Hazard -a Real Hazard to Natural Disaster Insurance? , 2022 .

[26]  P. Raschky Estimating the effects of risk transfer mechanisms against floods in Europe and U.S.A.: A dynamic panel approach , 2007 .

[27]  P. Raschky The overprotective parent - Bureaucratic agencies and natural hazard management , 2007 .

[28]  M. Kocher,et al.  Working Papers in Economics and Statistics the Impact of Naïve Advice and Observational Learning in Beauty-contest Games University of Innsbruck Working Papers in Economics and Statistics the Impact of Naïve Advice and Observational Learning in Beauty-contest Games , 2022 .

[29]  Todd L. Cherry,et al.  Conditional cooperation on three continents , 2008 .

[30]  Ludwig Fahrmeir,et al.  Propriety of posteriors in structured additive regression models: Theory and empirical evidence , 2009 .

[31]  R. Schwarze,et al.  Working Papers in Economics and Statistics Can Minimum Prices Assure the Quality of Professional Services? University of Innsbruck Working Papers in Economics and Statistics Can Minimum Prices Assure the Quality of Professional Services? , 2022 .

[32]  Paul H. C. Eilers,et al.  Splines, knots, and penalties , 2010 .