powered by
dgamma for beta
logd_beta(b, beta_shape = 2, beta_scale = 2)
log dgamma
current beta
prior for the shape (gamma distibution) of beta
prior for the scale of beta