Simulate a and Sab from full conditional distributions under bin likelihood
raSab_bin_fc(Z, Y, a, b, Sab, Sab0=NULL, eta0=NULL, SS = round(sqrt(nrow(Z))))
new value of Z
new value of Sab
new value of a
a square matrix, the current value of Z
square binary relational matrix
current value of row effects
current value of column effects
current value of Cov(a,b)
prior (inverse) scale matrix for the prior distribution
prior degrees of freedom for the prior distribution
number of iterations
Peter Hoff