Description Usage Arguments Value Author(s) References See Also Examples
Get the p-value of GATES. Usually it is used to get genomewise p-values. This function is taken from postgwas package. There is a little modification of the code GATES in postgwas package. 1) The approximated matrix may have negative eigen value, we modified it not to have negative values; 2) we added one more return (the key gene location) for Hyst method.
1 | GATES2(ldmatrix, p)
|
ldmatrix |
numeric. A correlation matrix of SNPs, dimensions matching the p and snps arguments. |
p |
p-value for each SNPs. |
A p-value of GATES and the key gene location (to be used by Hyst).
Milan Hiersche(taken from pastgwas package), Il-Youp Kwak(modified a little)
Miao-Xin Li, Hong-Sheng Gui, Johnny S.H. Kwan and Pak C. Sham (2011) GATES: A Rapid and Powerful Gene-Based Association Test Using Extended Simes Procedure The American Journal of Human Genetics 88, 283-293
1 2 3 4 5 6 7 8 9 10 11 | #simula <- simPathAR1Snp(nGenes=20, nGenes1=1, nSNPlim=c(1, 20), nSNP0=1:3,
# LOR=.2, rholim=c(0,0),
# n=30, MAFlim=c(0.05, 0.4), p0=0.05)
#Ps <- getlogitp(simula$Y, simula$X)
## get correlation of SNPs using controls
#ldmat <- cor(simula$X[ simula$Y == 0, ])
#o.pvec = order(Ps)
# ldmat <- ldmat[o.pvec, o.pvec]
#(gatesp <- GATES2(ldmat, sort(Ps))[1])
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.