Posterior Contraction and Credible Sets for Filaments of Regression Functions

by   Wei Li, et al.

A filament consists of local maximizers of a smooth function f when moving in a certain direction. Filamentary structures are important features of the shape of objects and are also considered as important lower dimensional characterization of multivariate data. There have been some recent theoretical studies of filaments in the nonparametric kernel density estimation context. This paper supplements the current literature in two ways. First, we provide a Bayesian approach to the filament estimation in regression context and study the posterior contraction rates using a finite random series of B-splines basis. Compared with the kernel-estimation method, this has theoretical advantage as the bias can be better controlled when the function is smoother, which allows obtaining better rates. Assuming that f: R^2 R belongs to an isotropic Hölder class of order α≥ 4, with the optimal choice of smoothing parameters, the posterior contraction rates for the filament points on some appropriately defined integral curves and for the Hausdorff distance of the filament are both (n/ n)^(2-α)/(2(1+α)). Secondly, we provide a way to construct a credible set with sufficient frequentist coverage for the filaments. Our valid credible region consists of posterior filaments that have frequentist interpretation. We demonstrate the success of our proposed method in simulations and application to earthquake data.


Adaptive Bayesian semiparametric density estimation in sup-norm

We investigate the problem of deriving adaptive posterior rates of contr...

A Theoretical Framework for Bayesian Nonparametric Regression: Orthonormal Random Series and Rates of Contraction

We develop a unifying framework for Bayesian nonparametric regression to...

Posterior Contraction Rates for Graph-Based Semi-Supervised Classification

This paper studies Bayesian nonparametric estimation of a binary regress...

Contraction rates for sparse variational approximations in Gaussian process regression

We study the theoretical properties of a variational Bayes method in the...

Posterior Consistency in the Binomial (n,p) Model with Unknown n and p: A Numerical Study

Estimating the parameters from k independent Bin(n,p) random variables, ...

A Divide and Conquer Algorithm of Bayesian Density Estimation

Data sets for statistical analysis become extremely large even with some...

Frank-Wolfe Bayesian Quadrature: Probabilistic Integration with Theoretical Guarantees

There is renewed interest in formulating integration as an inference pro...

Please sign up or login with your details

Forgot password? Click here to reset