A Gibbs sampler for updating the multiplicative effect matrices U and V
rUV_fc(Z, U, V, Suv, rho, s2 = 1, offset = 0)
a new value of U
a new value of V
n X n normal relational matrix
current value of U
current value of V
covariance of (U V)
dyadic correlation
dyadic variance
a matrix of the same dimension as Z. It is assumed that
Z-offset is equal to the multiplicative effects plus dyadic noise, so the
offset should contain any additive effects (such as Xbeta(X,beta+
outer(a,b,"+")
)
Peter Hoff