gcontrol {gap} | R Documentation |
the genomic control statistics of Devlin and Roeder, the list of parameters after data are optional (with default value)
gcontrol(data,zeta,kappa,tau2,epsilon,ngib,burn,idum)
data |
the data matrix |
zeta |
with default value 1000 |
kappa |
with default value 4 |
tau2 |
with default value 1 |
epsilon |
with default value 0.01 |
ngib |
number of Gibbs steps, with default value 500 |
burn |
number of burn-ins with default value 50 |
idum |
seed for pseudorandom number sequence |
The returned value is a list containing:
deltot |
the probability |
x2 |
the statistic |
A |
the A vector |
Devlin B, Roeder K (1999) Genomic control for association studies. Biometrics 55:997-1004
Adapted from gcontrol.c Bobby Jones and Kathryn Roeder, use -Dexecutable for standalone program, function getnum in the original code needs %*s to skip id string
Bobby Jones, Jing hua Zhao
## Not run: test<-c(1,2,3,4,5,6, 1,2,1,23,1,2, 100,1,2,12,1,1, 1,2,3,4,5,61, 1,2,11,23,1,2, 10,11,2,12,1,11) test<-t(matrix(test,nrow=6)) gcontrol(test) ## End(Not run)