Gibbs update for dyadic variance
rs2_fc(Z, rho,offset=0,nu0=NULL,s20=NULL)
Z
rho
offset
Xbeta(X,beta+ U%*%t(V) + outer(a,b,"+")
nu0
s20
a new value of s2
Peter Hoff
Useful links