Fixed bug in Ksi estimation in FactorEM, and also added iterative
FactorEM2 (will rename later)