Abstract:We consider a high-dimensional sparse normal means model where the goal is to estimate the mean vector assuming the proportion of non-zero means is unknown. We model the mean vector by a one-group global-local shrinkage prior belonging to a broad class of such priors that includes the horseshoe prior. We address some questions related to asymptotic properties of the resulting posterior distribution of the mean vector for the said class priors. We consider two ways to model the global parameter in this paper. Firstly by considering this as an unknown fixed parameter and then by an empirical Bayes estimate of it. In the second approach, we do a hierarchical Bayes treatment by assigning a suitable non-degenerate prior distribution to it. We first show that for the class of priors under study, the posterior distribution of the mean vector contracts around the true parameter at a near minimax rate when the empirical Bayes approach is used. Next, we prove that in the hierarchical Bayes approach, the corresponding Bayes estimate attains the minimax risk asymptotically under the squared error loss function. We also show that the posterior contracts around the true parameter at a near minimax rate. These results generalize those of van der Pas et al. (2014) \cite{van2014horseshoe}, (2017) \cite{van2017adaptive}, proved for the horseshoe prior. We have also studied in this work the asymptotic Bayes optimality of global-local shrinkage priors where the number of non-null hypotheses is unknown. Here our target is to propose some conditions on the prior density of the global parameter such that the Bayes risk induced by the decision rule attains Optimal Bayes risk, up to some multiplicative constant. Using our proposed condition, under the asymptotic framework of Bogdan et al. (2011) \cite{bogdan2011asymptotic}, we are able to provide an affirmative answer to satisfy our hunch.

On posterior distribution of Bayesian wavelet thresholding

Bayesian wavelet shrinkage with logistic prior

Kurtosis control in wavelet shrinkage with generalized secant hyperbolic prior

Deviation Bounds for Wavelet Shrinkage

Posterior Contraction rate and Asymptotic Bayes Optimality for one-group shrinkage priors in sparse normal means problem

Uniform almost sure convergence and asymptotic distribution of the wavelet-based estimators of partial derivatives of multivariate density function under weak dependence

Adapting to Unknown Smoothness by Aggregation of Thresholded Wavelet Estimators

Low and high dimensional wavelet thresholds for matrix-variate normal distribution

Wavelet methods in statistics: Some recent developments and their applications

Bayesian nonparametric regression using complex wavelets

Strong uniform convergence rates of the linear wavelet estimator of a multivariate copula density

Horseshoe priors for edge-preserving linear Bayesian inversion

Bayesian estimation for mean vector of multivariate normal distribution on the linear and nonlinear exponential balanced loss based on wavelet decomposition

Bernstein -- von Mises theorems for statistical inverse problems II: Compound Poisson processes

Oracle-type Posterior Contraction Rates in Bayesian Inverse Problems

Strong posterior contraction rates via Wasserstein dynamics

Tail-adaptive Bayesian shrinkage

A posterior contraction for Bayesian inverse problems in Banach spaces

Cross-validated wavelet block thresholding for non-Gaussian errors

A new approach to posterior contraction rates via Wasserstein dynamics