Create change_e.r

Function to change how M and h are estimated in the EM.