Sample additive row effects and their covariance for the binary family

Joint Gibbs update of the social-relations additive row effects a and the additive-effect covariance Sab under a probit (binary) likelihood. Given the latent normal scores Z and the observed 0/1 matrix Y, each row effect is redrawn from its Gaussian full conditional (conditioned on the column effects b through Sab) truncated to the interval that keeps every latent score sign-compatible with the observed responses. The covariance Sab is then refreshed from its inverse-Wishart full conditional.

Usage

raSab_bin_fc(
  Z,
  Y,
  a,
  b,
  Sab,
  Sab0 = NULL,
  eta0 = NULL,
  SS = round(sqrt(nrow(Z)))
)

Arguments

Z

square matrix of current latent normal scores

Y

square binary relational matrix (NA entries allowed)

a

current value of the row effects

b

current value of the column effects

Sab

current 2x2 covariance of the additive effects

Sab0

prior inverse-scale matrix (defaults to the 2x2 identity)

eta0

prior degrees of freedom (defaults to 4)

SS

number of inner Gibbs sweeps

Value

list with the updated Z, a and Sab

Author

lame authors