Full-conditional draw of the dyadic variance for replicated relational data

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.

Usage

rs2_rep_fc(E.T, rho)

Arguments

E.T

numeric array of dimension n x n x N holding the residual matrices along the third margin (one square slice per replicate).

rho

current value of the within-dyad correlation.

Value

a single numeric draw of s2.

Author

lame authors