Make gibbs.jl a dependency for posterior sampling
Put hyperprior on the gammas
Modify for sampling gammas