Statistics and Its Interface

Volume 2 (2009)

Number 1

Bayesian adaptive nonparametric M-regression

Pages: 71 – 81

DOI: http://dx.doi.org/10.4310/SII.2009.v2.n1.a7

Author

Colin Chen (SAS Institute Inc., Cary, North Carolina, U.S.A.)

Abstract

Nonparametric regression has been popularly used in curve fitting, signal denosing, and image processing. In such applications, the underlying functions (or signals) may vary irregularly, and it is very common that data are contaminated with outliers. Adaptive and robust techniques are needed to extract clean and accurate information. In this paper, we develop adaptive nonparametric M-regression with a Bayesian approach. This general approach fits M-regression using piecewise polynomial functions with an unknown number of knots at unknown locations, all treated as parameters to be inferred through Reversible Jump Markov Chain Monte Carlo (RJMCMC) of Green (1995). The Bayesian solution presented in this paper with computational details can be considered as an approximation to the general optimal solution for M-regression with free knots as described in Stone (2005). Numerical results show that the Bayesian approach performs well in various cases, especially with discontinuous underlying functions.

Keywords

least informative distribution, monotone score function, m-regression, noninformative prior, reversible jump Markov chain Monte Carlo, splines

2010 Mathematics Subject Classification

60K35

Full Text (PDF format)