internal function for calcAdjPP for the 2-trait case
calcD12( mod1, mod2, T1, T2, beta, SSy, Sxy, xcovo, Mx, Nqq, Vres, covY, nsnpspermodel )
| mod1 | vector of model indices for trait 1 |
|---|---|
| mod2 | vector of model indices for trait 2 |
| T1 | index of trait 1 |
| T2 | index of trait 2 |
| 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 |
| Vres | list of variance residuals |
| covY | covariance matrix of traits |
| nsnpspermodel | list of number of SNPs per model for each model in mod1, mod2 |
lbf12-lbf1-lbf2