Thanks Tom. One of the limitations of these models is that they are mixture models which marginalize out the latent class (which, if there are K skills, there are 2^K latent classes) which becomes computationally difficult (for large K) for the existing
libraries to handle as they typically use EM algorithm.
I'm only just becoming familiar with EP. Do you know if this limitation would be a problem when fitting the model with EP as well?