Latent Dirichlet Allocation with Gibbs sampling