5
$\begingroup$

The horseshoe prior for, say, an unknown mean $\mu$ is expressed $$ \begin{align*} \mu_i | \lambda_i, \tau &\sim N(0, \lambda_i^2 \tau^2) \\ \lambda_i &\sim C^+(0,1) \end{align*} $$ where the prior on $\lambda_i$ can be equivalently stated as $\kappa_i := \frac{1}{1+\lambda_i^2} \sim$ Beta(.5,.5), and it is this distribution on $\kappa_i$ that has the horseshoe shape for which the scheme is named. In the normal mean model $y | \mu \sim N_d (\mu, I)$, the $\kappa_i$ say how much the MLE $y$ should be shrunk, and the philosophy of the horseshoe is that you want $\kappa_i$ to be near 0 to produce no shrinkage in signal components, or near 1 to produce strong shrinkage to 0 of noise components. Hence you'd like a prior on $\kappa_i$ that has its mass toward 0 and 1, like Beta(.5,.5).

However, Beta(b,b) for b in (0,1) will have the same horseshoe shape. $b$ approaching 0 will shift mass out of the middle toward 0 and 1, with limiting distribution Bernoulli(.5). $b$ approaching 1 will shift mass from the ends into the middle, approaching U(0,1).

Is there any work on using values of $b$ other than .5 in horseshoe priors? If so, what guides such choices?

$\endgroup$
1
  • 2
    $\begingroup$ The Bayesian way would probably be to put another prior on your $b$ and see what posterior Horseshoe you'd end up with given your data... :-) $\endgroup$ Commented Nov 24, 2021 at 20:55

1 Answer 1

3
$\begingroup$

In the subjective Bayesian paradigm, you always have a choice of priors, so obviously you can use that alternative prior if you want. The only real issue here is whether or not you call that other prior a "horseshoe", which is really an aesthetic/terminology judgment, not a substantive statistical issue.$^\dagger$

As to what would guide the use of other types of beta priors with $0<b<1$ (compared to using $b = \tfrac{1}{2}$), presumably we would use $\tfrac{1}{2}<b<1$ when we want to be closer to a uniform distribution, and we would use a value $0<b<\tfrac{1}{2}$ when we want greater prior mass at the extremities. In the context of the model you mention for signals, this would be a contextual judgment of how much we wish to favour values near zero or one a priori.


$^\dagger$ Incidentally, a rather silly aspect of the whole "horseshoe" terminology is the implicit assumption that the horseshoe can only be looked at in one direction (i.e., held so that it appears convex). If you rotate the horseshoe 180-degrees it will then appear concave, so we could just as easily call concave priors "horseshoe" priors. In fact, none of these priors actually look like horseshoes, since actual horseshoes rotate around to an extent that they do not even look like a function.

$\endgroup$
1
  • $\begingroup$ The limit as b goes to 0 is Bernoulli(.5) which produces a spike and slab prior, albeit a more strict version as the membership probability is strictly 50-50 rather than specified by its own prior. The spike and slab scheme suffers from exponential explosion of the search space (or so I understand) so I wonder if using small $b$ starts getting into these kinds of computational difficulties. $\endgroup$ Commented Nov 25, 2021 at 3:59

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.