Draws a new value of the dyadic correlation rho from its full conditional using a random-walk-style Metropolis-Hastings step. Each dyad's pair of directed residuals is modelled as bivariate normal with unit variances (after scaling by the dyadic variance) and correlation rho; the proposal is a normal truncated to (-1, 1) and the acceptance ratio combines the bivariate-normal log-likelihood, an optional arc-sine prior, and the truncated-proposal normaliser.
rrho_mh(Z, rho, s2 = 1, offset = 0, asp = NULL)a new (scalar) value of rho.
n x n normal relational matrix.
current value of rho.
current value of the dyadic variance.
matrix matching Z; Z - offset is treated as
dyadic noise (so it should absorb any additive/multiplicative effects).
logical; use an arc-sine prior (TRUE, the default when
NULL) or a uniform prior (FALSE).
lame authors