estimateEM function

Perform One EM