amen (version 1.0)

raSab_bin_fc: Simulate a and Sab from full conditional distributions under bin likelihood

Description

Simulate a and Sab from full conditional distributions under bin likelihood

Usage

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

Arguments

Z

a square matrix, the current value of Z

Y

square binary relational matrix

a

current value of row effects

b

current value of column effects

Sab

current value of Cov(a,b)

SS

number of iterations

Value

• Z
  new value of Z
• Sab
  new value of Sab
• a
  new value of a

Examples

## The function is currently defined as
raSab_bin_fc<-function(Z, Y, a, b, Sab, SS = round(sqrt(nrow(Z))))
{
    E <- Z - a %*% t(rep(1, nrow(Z)))
    MEL <- MEU <- -E
    MEL[Y == 0] <- -Inf
    MEU[Y == 1] <- Inf
    diag(MEU) <- Inf
    diag(MEL) <- -Inf
    lba <- apply(MEL, 1, max)
    lba[is.na(lba)] <- -Inf
    uba <- apply(MEU, 1, min)
    uba[is.na(uba)] <- Inf
    for (ss in 1:SS)
    {
        ea <- b * Sab[1, 2]/Sab[2, 2]
        sa <- sqrt(Sab[1, 1] - Sab[1, 2]^2/Sab[2, 2])
        a <- ea + sa * qnorm(runif(nrow(Z), pnorm((lba - ea)/sa),
            pnorm((uba - ea)/sa)))
        Sab <- solve(rwish(solve(diag(2) + crossprod(cbind(a, b))), 3 + nrow(Z)))
    }
    list(Z = E + a %*% t(rep(1, nrow(Z))), a = a, Sab = Sab)
  }