alpha     ~ Gam(1,1)
eta       ~ Gam(1,1)

beta_{k}    ~ DirSym(eta, V)              , k \in [1,K]
theta_{d,k} ~ DirSym(alpha, K)            , d \in [1,D] , k \in [1,K]

z_{d,n}     ~ Mult(theta_{d,z_{d,(n-1)}}) , d \in [1,D] , n \in [1,N_{d}]
w_{d,n}     ~ Mult(beta_{z_{d,n}})        , d \in [1,D] , n \in [1,N_{d}]