The MATLAB function "varbvs" now includes an option to update the
variance of the residual (sigma).