Added an option to the R function varbvsoptimize ('update.sigma') that
computes the maximum likelihood estimate of the residual variance
rather than setting this value by hand.