This is the method proposed by Van_MGAS
MGAS(z_vector, est_genetic_cor, est_pheno_cor)
Arguments
- z_vector
Column vectorized data matrix with rows represent
phenotype and columns represent genotype.
- est_genetic_cor
Estimated phenotype correlation matrix.
- est_pheno_cor
Estimated genotype correlation matrix.
Value
A numeric value represents the p value of MGAS.
Examples
z_vector<-MASS::mvrnorm(1,mu=rep(0,9),Sigma = diag(nrow = 9, ncol = 9))
genotype_covariance<-diag(nrow = 3,ncol = 3)
phenotype_covariance<-diag(nrow = 3,ncol = 3)
ThreeWayTest::MGAS(z_vector=z_vector,
est_genetic_cor = genotype_covariance,
est_pheno_cor =phenotype_covariance)
#> [1] 0.5115277