LDA function

Create a Latent Dirichlet Allocation model