The only difference is that our method uses marginal probabilities given by R instead of the p(x, y; θ) of standard HMMs.
only a forward-backward algorithm is required for the efficient calculation of our parameter estimation process.