I have the "outer loop" implemented for the variational inference
method allowing for arbitary covariates (Z). Next, I need to implement
the "inner loop" (varbvszbin.m).