rrho_fc function

Griddy Gibbs update for dyadic correlation

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)

Returns

a value of rho

Author(s)

Peter Hoff