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)
a value of rho
n X n normal relational matrix
covariance of additive effects
residual variance
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)
)
the number of points for an unevenly-spaced grid on which to approximate the full conditional distribution
use arc sine prior (TRUE) or uniform prior (FALSE)
Peter Hoff