A Gibbs sampler for updating the multiplicative effect matrices U and V, assuming they are the same across replicates.
rUV_rep_fc(E.T,U,V,rho,s2=1,shrink=TRUE)
a new value of U
a new value of V
Array of square residual relational matrix series with additive effects and covariates subtracted out. The third dimension of the array is for different replicates. Each slice of the array according to the third dimension is a square residual relational matrix.
current value of U
current value of V
dyadic correlation
dyadic variance
adaptively shrink the factors with a hierarchical prior
Peter Hoff, Yanjun He