LDA modeling algorithm