(2015) provide a blueprint for Gibbs sampling with potential functions without additional parameters, our model has additional parameters of ρ S→T and ρ T →S so we need to optimize them.
we use stochastic EM (Celeux, 1985).