I modified the function VARBVSZBIN so that includes the additional
option of being able to compute the maximum likelihood estimate for
the prior variance of the additive effects ("sa").