getlogitp: Get the p-values of the multiple SNPs based on the Wald test...

Description Usage Arguments Value See Also Examples

Description

Get p-value using logistic regresion for each of the multiple SNPs

Usage

1
getlogitp(Y, X)

Arguments

Y

Response or phenotype data. It can be a disease indicator; =0 for controls, =1 for cases.

X

Genotype or other data; each row for a subject, and each column for an SNP (or a predictor). The value of each SNP is the # of the copies for an allele. A matrix with dimension n by p (n : number of observation, p : number of SNPs (or predictors) ).

Value

p-values for each SNPs.

See Also

GATES2 GatesSimes Hyst

Examples

1
2
3
4
simula <- simPathAR1Snp(nGenes=20, nGenes1=1, nSNPlim=c(1, 20), nSNP0=1:3,
                           LOR=.2, rholim=c(0,0),
                           n=100, MAFlim=c(0.05, 0.4), p0=0.05)
logitp <- getlogitp(simula$Y, simula$X)

aSPU documentation built on May 1, 2019, 7:04 p.m.