Gibbs sampling for Bayesian nonconjugate models using auxilliary variables