APaSPU: The function for adaptive sum of powered score (aSPU) test

In jasonzyx/APaSPU: Adaptive Sum of Powered Test (aSPU) Based on Approximated Score Vector

Description

It gives the simulation-based p-values of the sum of powered score (SPU) with different power and aSPU test.

Usage

APaSPU(betas, vcov_mat, pow = c(1:8, Inf), B)

Arguments

betas	the vector (length = p) of maximum likelihood estimates; p is the number of paramters to be tested.
vcov_mat	the p by p variance-covariance matrix of MLE.
pow	the power vector for aSPU tests.
B	number of simulations

Details

the suggested power term is c(1:8, Inf)

Value

simulation-based p-values of the sum of powered score (SPU) with different power and aSPU test

Author(s)

Zhiyuan (Jason) Xu, Yiwei Zhang and Wei Pan

References

Zhiyuan Xu, Wei Pan (2014) Approximate score-based testing with application to multivariate trait association analysis. Genetic Epidemiology. 39(6): 469-479.

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. 96:309-25

Examples

library(lme4)
data("exdat_GLMM")
fit1<-glmer(Y~.-ID + (1|ID),family="binomial",
data=exdat_GLMM,control=glmerControl(optimizer="bobyqa"))
betas<-fixef(fit1)[-1]
Cov<-vcov(fit1)[-1,-1]
APaSPU(betas,Cov,pow=c(1:8,Inf),B=1000)

jasonzyx/APaSPU documentation built on May 18, 2019, 5:55 p.m.