internal function for calcAdjPP for that gives list of covariance matrix of residuals for all trait pairs

allC12(M, nummods, beta, SSy, Sxy, xcovo, Mx, Nqq)

Arguments

M

number of traits

nummods

list where component i is the number of models for trait i

beta

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

xcovo

covariance matrix of c("one"=1,X)

Mx

vector of SNP means

Nqq

has N[i,j] = no. with both trait i and trait j measured