variance of model residuals for trait T1 at all models that have joint effect estimates

Vres.all(Nqq, mbeta, SSy, Sxy)

Arguments

Nqq

matrix in which Nqq[i,j] = number of individuals measured in both trait i and trait j

mbeta

list of joint beta estimates for each trait

SSy

matrix of trait cross-products

Sxy

matrix with each column being the cross-product between SNPs and a trait