Matt Hoffman
unread,Mar 8, 2013, 5:55:44 PM3/8/13Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to stan-dev
Hi again,
Suppose we have as part of our model
β_n ~ Normal(μ, 1/sqrt(τ))
τ ~ Gamma(a, b).
We can integrate out τ and say
log p(β) = a*log(b) - log(Γ(a)) - 0.5*N*log(2*π) + log(Γ(a + N/2)) -
(a + N/2)*log(b + 0.5*Σ_n (β_n-μ)^2).
Is there a name for that distribution, akin to the beta-binomial or
Dirichlet-multinomial distribution? Inverse-gamma-normal would be the
analogous name, but there's already the normal-inverse-gamma, which is
something else.
I ask because I want to implement this distribution in Stan, and I
don't know what to call it.
Matt