gap (version 0.4)

gcontrol: 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

source

http://www.stat.cmu.edu

References

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

Examples

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)