Single C++ function replacing the bipartite beta/a/b update block in lame.R
rbeta_ab_bip_gibbs_cpp(
Z,
Xlist,
UV_eff,
a_current,
b_current,
s2,
g_prior,
va,
vb,
rvar,
cvar
)List with beta, a, b
3D array (nA x nB x T)
List of T arrays (nA x nB x p)
nA x nB matrix (UGV' or U*V')
Current row effects (length nA)
Current column effects (length nB)
Dyadic variance
G-prior parameter
Row effect variance (diagonal element of Sab)
Column effect variance (diagonal element of Sab)
Whether to update row effects
Whether to update column effects