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