Samples the dyadic variance s2 from its Gibbs full conditional when
the residual cube stacks several independent replicate networks that share a
common within-dyad correlation rho. Each ordered dyad pair is whitened
by the symmetric inverse square root of the 2x2 correlation matrix and the
pooled whitened sum of squares drives an inverse-gamma draw.
rs2_rep_fc(E.T, rho)a single numeric draw of s2.
numeric array of dimension n x n x N holding the residual matrices along the third margin (one square slice per replicate).
current value of the within-dyad correlation.
lame authors