powered by
dgamma for gamma
logd_gamma(g, gamma_shape = 2, gamma_scale = 2)
log dgamma
current gamma
prior for the shape (gamma distibution) of gamma
prior for the scale of gamma