Learn R Programming

aSPU (version 1.38)

GEEaSPU: GEEaSPU test for association with multiple traits in generalized estimation equations.

Description

It gives p-values of the GEESPU tests and GEEaSPU test.

Usage

GEEaSPU(traits, geno, Z = NULL, model = c("binomial", "gaussian"),
  gamma = c(1:8, Inf), n.sim = 1000)

Arguments

traits
trait matrix. The row for individuals and the column for traits.
geno
A matrix of genetic information.
Z
covariates.
model
Use "gaussian" for a quantitative trait, and use "binomial" for a binary trait.
gamma
power used in GEEaSPU test. A vector of the powers.
n.sim
number of simulations

Value

  • p-values for the GEE-SPU and GEE-aSPU test.

References

Yiwei Zhang, Zhiyuan Xu, Xiaotong Shen, Wei Pan (2014) Testing for association with multiple traits in generalized estimation equations, with application to neuroimaging data. Neuroimage. 1;96:309-25

Examples

Run this code
traits <- matrix(rnorm(100*5, 0,1), ncol=5)
Z <- rnorm(100, 2, 0.5)
geno <- rbinom(100, 2, 0.5)
out <- GEEaSPU(traits, geno, Z = NULL, model = "gaussian",
		  gamma = c(1:8,Inf), n.sim = 100)

Run the code above in your browser using DataLab