Locally adaptive Bayesian P-splines with a Normal-Exponential-Gamma prior

  • Scheipl F
  • Kneib T
  • 11


    Mendeley users who have this article in their library.
  • 12


    Citations of this article.


An implementation of locally adaptive penalized spline smoothing using a class of heavy-tailed shrinkage priors for the estimation of functional forms with highly varying curvature or discontinuities is presented. These priors utilize scale mixtures of normals with locally varying exponential-gamma distributed variances for the differences of the P-spline coefficients. A fully Bayesian hierarchical structure is derived with inference about the posterior being based on Markov Chain Monte Carlo techniques. Three increasingly flexible and automatic approaches are introduced to estimate the spatially varying structure of the variances. An extensive simulation study for Gaussian, Poisson, and Binomial responses shows that the performance of this approach on a number of benchmark functions is competitive to that of previous approaches. Results from applications with Gaussian and Poisson responses support the simulation results. © 2009 Elsevier B.V. All rights reserved.

Get free article suggestions today

Mendeley saves you time finding and organizing research

Sign up here
Already have an account ?Sign in

Find this document


  • Fabian Scheipl

  • Thomas Kneib

Cite this document

Choose a citation style from the tabs below

Save time finding and organizing research with Mendeley

Sign up for free