Gibbs update for dyadic variance with independent replicate relational data
rs2_rep_fc(E.T, rho)Array of square residual relational matrix series. 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 rho
a new value of s2