GBJ_objective: GBJ_objective.R

In GBJ: Generalized Berk-Jones Test for Set-Based Inference in Genetic Association Studies

Description

Calculates the GBJ objective function at given threshold points. Used both to calculate the observed GBJ statistic and also to find the boundary points for p-value calculation (through uniroot). Mostly for internal use.

Usage

GBJ_objective(t_vec, d, k_vec = NULL, pairwise_cors, offset = 0)

Arguments

t_vec	A scalar or vector of threshold points (magnitudes of test statistics) to calculate the objective at.
d	The total number of test statistics in the set.
k_vec	If t_vec is not a vector all of the test statistics, let us know which ordered objective functions we are calculating (can be a vector with the same length as t).
pairwise_cors	A vector of all d(d-1)/2 pairwise correlations between the test statistics, where d is total number of test statistics in the set.
offset	Used if we want to subtract the observed gbj value for uniroot.

Value

The GBJ objective at t (for given d, kkk, pairwise_cors) minus the offset

Examples

GBJ_objective(t_vec=1:5, d=5, k_vec=NULL, pairwise_cors=rep(0.2,10), offset=0)

GBJ documentation built on March 26, 2020, 6:05 p.m.