
It has been shown in [1]} that the optimal solution of (REF ) is given by
\(Q^\dagger ({\Theta }_j)=\frac{\exp (\mathbb {E}_{\prod _{i\ne j}Q({\Theta }_i)}\ln p(\Theta ,\lbrace \mathcal {Y}_l\rbrace _{l=1}^L))}{\int \mathbb {E}_{\prod _{i\ne j}Q({\Theta }_i)}\ln p(\Theta ,\lbrace \mathcal {Y}_l\rbrace _{l=1}^L)d\Theta _j}, \forall j.\) 
