An latent Dirichlet allocation (with variational inference and Gibbs sampling) implementation in Julia