# GT.localfdr: Between and within group local fdr scores In GroupTest: Multiple Testing Procedure for Grouped Hypotheses

## Description

This function calculate the between-group and within-group local fdr scores for a given set of all the parameters.

## Usage

 1 GT.localfdr(TestStatistic, pi1, pi2.1, L, muL, sigmaL, cL)

## Arguments

 TestStatistic An array of list. Each element of the array corresponds to one group, containing the test statistic, stored as X, and the group size, stored as mg. L The number of Gaussian component under the alternative hypothesis. pi1 π_1, the probability that a group is significant. pi2.1 π_{2|1}, the probability that an individual null hypothesis is false given that the group is significant. muL a vector of means for all the components of the Gaussian mixture. sigmaL a vector of standard deviation of all the components of the Gaussian mixture. cL a vector of the probability for all the components of the Gaussian mixture.

## Value

This function returns an array of G lists where G is the number of groups.

 TSGroupTest[[g]] in each element, the individual conditional local fdr score ( P(θ_{gj}=0|x, θ_{g}=1)), the group-wise local fdr score (P(θ_g=0|x)), are stored.

## Examples

 1 2 3 data(GroupTest_simulate) GroupTest_simulate <- GT.localfdr( GroupTest_simulate, L=2, pi1=0.5, pi2.1=0.5, muL=c(-1, 1), sigmaL=c(1,2), cL=c(0.4,0.6) )

GroupTest documentation built on May 30, 2017, 12:10 a.m.