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

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