powered by
Metropolis update for dyadic correlation
rrho_mh(E, rho, s2 = 1)
square residual relational matrix
current value of rho
current value of s2
a new value of rho