covariance between residuals of a pair of models for a trait pair

calcCres12(imod1, imod2, T1, T2, beta, SSy, Sxy, xcovo, Mx, Nqq)

Arguments

imod1

index of model from STR for trait T1, to be called from beta

imod2

index of model from STR for trait T2, to be called from beta

T1

index of first trait

T2

index of second trait

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