GT.localfdr {GroupTest}

Between and within group local fdr scores

Description

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

Usage

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

\pi_1, the probability that a group is significant.

pi2.1

\pi_{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(\theta_{gj}=0|x, \theta_{g}=1)), the group-wise local fdr score (P(\theta_g=0|x)), are stored.

Examples

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) )