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