### R code from vignette source 'vignettes/HEM/inst/doc/HEM.Rnw'