Skip Navigation

Biometrika 2001 88(4):1055-1071; doi:10.1093/biomet/88.4.1055
© 2001 by Biometrika Trust
This Article
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Dimatteo, I.
Right arrow Articles by Kass, R. E.
Right arrow Search for Related Content
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Bayesian curve-fitting with free-knot splines

Ilaria Dimatteo1, Christopher R.Genovese1 and Robert E.Kass1

1 Department of Statistics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, U.S.Adimatteo{at}stat.cmu.edugenovese{at}stat.cmu.edukass{at}stat.cmu.edu

We describe a Bayesian method, for fitting curves to data drawn from an exponential family, that uses splines for which the number and locations of knots are free parameters.The method uses reversible-jump Markov chain Monte Carlo to change the knot configurations and a locality heuristic to speed up mixing. For nonnormal models, we approximate the integrated likelihood ratios needed to compute acceptance probabilities by using the Bayesian information criterion, BIC, under priors that make this approximation accurate. Our technique is based on a marginalised chain on the knot number and locations, but we provide methods for inference about the regression coefficients, and functions of them, in both normal and nonnormal models. Simulation results suggest that the method performs well, and we illustrate the method in two neuroscience applications.

Key Words: BIC; Generalised linear model; Nonparametric regression; Reversible-jump Markov chain Monte Carlo; Smoothing; Unit-information prior

Received February 2000. Revised May 2001


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
 BiostatisticsHome page
K. Berhane and N.-T. Molitor
A Bayesian approach to functional-based multilevel modeling of longitudinal data: applications to environmental epidemiology
Biostat., March 18, 2008; (2008) kxm059v1.
[Abstract] [Full Text] [PDF]

Home page
 Neural Comput.Home page
H. Shimazaki and S. Shinomoto
A Method for Selecting the Bin Size of a Time Histogram
Neural Comput., June 1, 2007; 19(6): 1503 - 1527.
[Abstract] [Full Text] [PDF]

Home page
 Neural Comput.Home page
W. Truccolo and J. P. Donoghue
Nonparametric Modeling of Neural Point Processes via Stochastic Gradient Boosting Regression.
Neural Comput., March 1, 2007; 19(3): 672 - 705.
[Abstract] [Full Text] [PDF]

Home page
 J. Neurophysiol.Home page
V. Ventura, C. Cai, and R. E. Kass
Trial-to-Trial Variability and Its Effect on Time-Varying Dependency Between Two Neurons
J Neurophysiol, October 1, 2005; 94(4): 2928 - 2939.
[Abstract] [Full Text] [PDF]

Home page
 J. Neurophysiol.Home page
V. Ventura, C. Cai, and R. E. Kass
Statistical Assessment of Time-Varying Dependency Between Two Neurons
J Neurophysiol, October 1, 2005; 94(4): 2940 - 2947.
[Abstract] [Full Text] [PDF]

Home page
 J. Neurophysiol.Home page
R. E. Kass, V. Ventura, and E. N. Brown
Statistical Issues in the Analysis of Neuronal Data
J Neurophysiol, July 1, 2005; 94(1): 8 - 25.
[Abstract] [Full Text] [PDF]

Home page
 BioinformaticsHome page
V. N. Minin, K. S. Dorman, F. Fang, and M. A. Suchard
Dual multiple change-point model leads to more accurate recombination detection
Bioinformatics, July 1, 2005; 21(13): 3034 - 3042.
[Abstract] [Full Text] [PDF]


Disclaimer:
Please note that abstracts for content published before 1996 were created through digital scanning and may therefore not exactly replicate the text of the original print issues. All efforts have been made to ensure accuracy, but the Publisher will not be held responsible for any remaining inaccuracies. If you require any further clarification, please contact our Customer Services Department.