inits_MCMC_hierarchical <- function(){list("mu.alpha.prob" = 1, "mu.gamma.prob" = 2, "mu.delta.prob" = 2, "mu.rho.log" = -.0.5)}