covariance between residuals of a pair of models for a trait pair
calcCres12(imod1, imod2, T1, T2, beta, SSy, Sxy, xcovo, Mx, Nqq)
| 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 |