alpha     ~ Gam(10,10)
mu_{k}    ~ NorMV(vec(0.0,1,dim), 1)     , k \in [1,K]
si2       ~ IG(10,10)
pi        ~ DirSym(alpha, K)
z_{n}     ~ Mult(pi)                     , n \in [1,N]
x_{n}     ~ NorMV(mu_{z_{n}}, si2)       , n \in [1,N]

mu_{k,d} \in R , k \in [1,K] , d \in [1,dim]
x_{n,d}  \in R , n \in [1,N] , d \in [1,dim]