However, that is computationally impractical because a complicated integral computation is needed.
we estimate an approximate distribution of P (π|w, y, α, θ) using Variational Bayes Method(H. Attias, 1999).