gcontrol {gap}R Documentation

genomic control

Description

the genomic control statistics of Devlin and Roeder, the list of parameters after data are optional (with default value)

Usage

gcontrol(data,zeta,kappa,tau2,epsilon,ngib,burn,idum)

Arguments

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

Value

The returned value is a list containing:

deltot the probability
x2 the statistic
A the A vector

References

Devlin B, Roeder K (1999) Genomic control for association studies. Biometrics 55:997-1004

Note

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

Author(s)

Bobby Jones, Jing hua Zhao

Source

http://www.stat.cmu.edu

Examples

## 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)

[Package Contents]