The variational inference algorithm (varbvszbin.m) appears to be
working---at least it recovers the same solution as varbvsbin.m in the
special case when we have a single covariate, the intercept.