variance of model residuals for trait T1 at model index imod
Vres.hat(imod, T1, SSy, Sxy, mbeta, Nqq)
| imod | model index |
|---|---|
| T1 | index of trait |
| SSy | matrix of trait cross-products |
| Sxy | matrix with each column being the cross-product between SNPs and a trait |
| mbeta | list of joint beta estimates for each trait |
| Nqq | matrix in which Nqq[i,j] = number of individuals measured in both trait i and trait j |