### R code from vignette source 'HEM.Rnw'