shad/ML1: add em.R