IDEAS home Printed from
MyIDEAS: Log in (now much improved!) to save this paper

Spatially Adaptive Bayesian P-Splines with Heteroscedastic Errors

Listed author(s):
  • Ciprian Crainiceanu

    (Johns Hokins Bloomberg School of Public Health, Department of Biostatistics)

  • David Ruppert

    (School of Operational Research & Industrial Engineering, Cornell University)

  • Raymond Carroll

    (Department of Statistics, Texas A&M University)

Registered author(s):

    An increasingly popular tool for nonparametric smoothing are penalized splines (P-splines) which use low-rank spline bases to make computations tractable while maintaining accuracy as good as smoothing splines. This paper extends penalized spline methodology by both modeling the variance function nonparametrically and using a spatially adaptive smoothing parameter. These extensions have been studied before, but never together and never in the multivariate case. This combination is needed for satisfactory inference and can be implemented effectively by Bayesian \mbox{MCMC}. The variance process controlling the spatially-adaptive shrinkage of the mean and the variance of the heteroscedastic error process are modeled as log-penalized splines. We discuss the choice of priors and extensions of the methodology,in particular, to multivariate smoothing using low-rank thin plate splines. A fully Bayesian approach provides the joint posterior distribution of all parameters, in particular, of the error standard deviation and penalty functions. In the multivariate case we produce maps of the standard deviation and penalty functions. Our methodology can be implemented using the Bayesian software WinBUGS.

    If you experience problems downloading a file, check if you have the proper application to view it first. In case of further problems read the IDEAS help page. Note that these files are not on the IDEAS site. Please be patient as the files may be large.

    File URL:
    Download Restriction: no

    Paper provided by Berkeley Electronic Press in its series Johns Hopkins University Dept. of Biostatistics Working Paper Series with number 1061.

    in new window

    Date of creation: 27 Sep 2004
    Handle: RePEc:bep:jhubio:1061
    Contact details of provider: Web page:

    References listed on IDEAS
    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.:

    in new window

    1. Andrew Gelman, 2004. "Prior distributions for variance parameters in hierarchical models," EERI Research Paper Series EERI_RP_2004_06, Economics and Econometrics Research Institute (EERI), Brussels.
    2. Andrew Gelman, 2004. "Prior distributions for variance parameters in hierarchical models," Econometrics 0404001, EconWPA.
    Full references (including those not matched with items on IDEAS)

    This item is not listed on Wikipedia, on a reading list or among the top items on IDEAS.

    When requesting a correction, please mention this item's handle: RePEc:bep:jhubio:1061. 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: (Christopher F. Baum)

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    If references are entirely missing, you can add them using this form.

    If the full references list an item that is present in RePEc, but the system did not link to it, you can help with this form.

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your profile, as there may be some citations waiting for confirmation.

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    This information is provided to you by IDEAS at the Research Division of the Federal Reserve Bank of St. Louis using RePEc data.