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