Added EM algorithm for estimating average yield of sub regions