We discuss a Bayesian formalism which gives rise to a type of wavelet threshold estimation in non-parametric regression. A prior distribution is imposed on the wavelet coe#cients of the unknown response function, designed to capture the sparseness of wavelet expansion common to most applications. For the prior speci#ed, the posterior median yields a thresholding procedure. Our prior model for the underlying function can be adjusted to give functions falling in any speci#c Besov space. We ...