Latent Dirichlet Allocation (MCMC with model selection)

Does latent dirichlet allocation

Do LDA by Gibbs sampling
USAGE: [Phi,Theta,L,LogLAll]=LDA_Gibbs(WS,DS, K, Alpha, Beta, Iter, BURNIN, LAG)
This function was inspired by the following code...
A Matlab demo of collapsed Gibbs sampling for learning LDA by Xiaodong Yu
http://xiaodong-yu.blogspot.jp/2009/08/matlab-demo-of-gibbs-sampling-for.html
I'd be happy to have any suggestions or comments!
本プログラムは、日本語用のテキストマイニングを行うツールボックスの中の一つとして作られました。
https://bitbucket.org/IkkoKimura/tmt/downloads
もしよろしければ、こちらの方もご参照ください。