reconstruct joint model ioutput for two snp models
qq.twosnp(b1, b2, g1, g2, vb1, vb2, vg1, vg2, fx, fz, rho, v1, v2, eta, N)
| b1 | coef for snp X, trait 1 |
|---|---|
| b2 | coef for snp X, trait 2 |
| g1 | coef for snp Z, trait 1 |
| g2 | coef for snp Z, trait 2 |
| vb1 | var(coef) for snp X, trait 1 |
| vb2 | var(coef) for snp X, trait 2 |
| vg1 | var(coef) for snp Z, trait 1 |
| vg2 | var( coef ) for snp Z, trait 2 |
| fx | MAF of X |
| fz | MAF of Z |
| rho | cor(X,Z) |
| v1 | variance of trait 1 |
| v2 | variance of trait 2 |
| eta | cor(Y1,Y2) |
| N | number of samples |
list of coefficients and their (co-)variances that would be expected from a joint model
The SNPs are X and Z, the traits are 1 and 2