SEAGLE: Compute score-like test statistic and p-value for GxE test...

In SEAGLE: Scalable Exact Algorithm for Large-Scale Set-Based Gene-Environment Interaction Tests

Description

This function computes the score test statistic and corresponding p-value for the GxE test with the SEAGLE algorithm with input data that have been prepared with the prep.SEAGLE function

Usage

SEAGLE(obj.SEAGLE, init.tau = 0.5, init.sigma = 0.5, pv = "liu")

Arguments

obj.SEAGLE	Input data prepared with prep.SEAGLE function
init.tau	Initial estimate for tau (Default is 0.5)
init.sigma	Initial estimate for sigma (Default is 0.5)
pv	Method of obtaining p-value (Either "liu" or "davies", Default is liu)

Value

Score-like test statistic T for the GxE effect and corresponding p-value

Examples

dat <- makeSimData(H=cosihap, n=500, L=10, gammaG=1, gammaGE=0, causal=4, seed=1)
objSEAGLE <- prep.SEAGLE(y=dat$y, X=dat$X, intercept=1, E=dat$E, G=dat$G)
res <- SEAGLE(objSEAGLE, init.tau=0.5, init.sigma=0.5)

SEAGLE documentation built on Nov. 6, 2021, 1:06 a.m.