added option to estimate sigma in EM (need testing)