Griddy Gibbs update for dyadic correlation

Source: R/rrho_fc.R

Simulation of dyadic correlation from its approximate full conditional distribution using griddy Gibbs sampling

rrho_fc(Z, Sab, s2 = 1, offset = 0, ngp = 100, asp = NULL)

Arguments

Z

n X n normal relational matrix
Sab

covariance of additive effects
s2

residual variance
offset

matrix of the same dimension as Z. It is assumed that Z-offset follows an SRM distribution, so the offset should contain any regression terms and multiplicative effects (such as Xbeta(X,beta+ U%*%t(V) )
ngp

the number of points for an unevenly-spaced grid on which to approximate the full conditional distribution
asp

use arc sine prior (TRUE) or uniform prior (FALSE)

Value

a value of rho

Author

Peter Hoff