added EM cluster analysis