Internal function used compute the products in cross-sectional VB algorithm and Gibbs sampler
Xt_siginv_X(tx, siginv, y = NULL)
transpose of the X design matrix
inverse variance matrix
outcome matrix. if NULL
, function computes
first product; if not, function computes second product.
Jeff Goldsmith ajg2202@cumc.columbia.edu