R/flashfm.R
allC12.Rdinternal 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)
| 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 |