Removed state variable x for Gibbs sampler
Added guards in nbinom in case there are no DE genes
Different initial conditions for eta and zeta
Update zeta using method by West