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.
Arguments
- Z
n x n normal relational matrix.
- rho
current value of rho.
- s2
current value of the dyadic variance.
- offset
matrix matching
Z;Z - offsetis treated as dyadic noise (so it should absorb any additive/multiplicative effects).- asp
logical; use an arc-sine prior (
TRUE, the default whenNULL) or a uniform prior (FALSE).