In this paper, we consider a number of algorithms for estimating the parameters of ME models, including Generalized Iterative Scaling and Improved Iterative Scaling, as well as general purpose optimization techniques such as gradient ascent, conjugate gradient, and variable metric methods.
for the comparison, we implemented both Generalized and Improved Iterative Scaling in C++ using the primitives provided by PETSc.