Simulates from the joint full conditional distribution of (a,b,beta), assuming same additive row and column effects and regression coefficient across replicates.
rbeta_ab_rep_fc(Z.T,Sab,rho,X.T,s2=1)
regression coefficients
additive row effects
additive column effects
n x n x T array, with the third dimension for replicates. Each slice of the array is a (latent) normal relational matrix, with multiplicative effects subtracted out
row and column covariance
dyadic correlation
n x n x p x T covariate array
dyadic variance
Peter Hoff, Yanjun He