Update hamimc.R