GHC_objective: GHC_objective.R

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

GHC_objective.R

Description

Internal function to calculate the value of the kth GHC objective (possibly minus an offset) when the kth p-value order statistic is x, for a set of size d.

Usage

GHC_objective(x, k, d, offset, pairwise_cors)

Arguments

x	The p-value of the kth p-value order statistic.
k	Which objective to use.
d	The size of the set.
offset	Used to zero the correct value when we put this into uniroot
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.

Value

The value of the kth HC objective

Examples

GHC_objective(x=0.1, k=2, d=5, offset=0, pairwise_cors=rep(0.2,10))

ryanrsun/GBJ documentation built on Feb. 6, 2024, 1:21 a.m.