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