Description Usage Arguments Value Author(s) References See Also Examples
Return exact minP test p-value for multiple traits - single SNP association.
1 | minP(Zi, r)
|
Zi |
a vector of summary Z-scores for single SNP |
r |
estimated correlation matrix based on the summary Z-scores (output of estcov) |
return exact minP test
Junghi Kim, Yun Bai and Wei Pan
Junghi Kim, Yun Bai and Wei Pan (2015) An Adaptive Association Test for Multiple Phenotypes with GWAS Summary Statistics, Genetic Epidemiology, 8:651-663
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | # -- n.snp: number of SNPs
# -- n.trait: number of traits
# -- n.subject: number of subjects
n.snp <- 100
n.traits <- 10
n.subjects <- 1000
traits <- matrix(rnorm(n.subjects*n.traits), n.subjects, n.traits)
v <- cov(traits)
allZ <- rmvnorm(n.snp, Sigma=v)
colnames(allZ) <- paste("trait", 1:n.traits, sep="")
rownames(allZ) <- paste("snp", 1:n.snp, sep="")
r <- estcov(allZ)
MTaSPUs(Z = allZ, v = r, B = 100, pow = c(1:4, Inf), transform = FALSE)
MTaSPUs(Z = allZ[1,], v = r, B = 100, pow = c(1:4, Inf), transform = FALSE)
minP(Zi= allZ[1,], r = r)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.