hierarchicalFDR: Hierachical testing with FDR control

Description Usage Arguments Details Value References See Also Examples

View source: R/hierarchicalFDR.R

Description

Apply hierarchical test for each hierarchy, and test external variables for FDR control at level alpha

Usage

1

Arguments

X

original data

y

associated response

group

vector with index of groups. group[i] contains the index of the group of the variable var[i].

var

vector whith the variables contained in each group. group[i] contains the index of the group of the variable var[i].

test

function for testing the nullity of a group of coefficients in linear regression. 3 parameters : X : design matrix, y response and varToTest : vector of variables to test; return a pvalue

Details

Version of the hierarchical testing procedure of Yekutieli for HCgglasso output. You can use th selFDR function to select groups at a desired level alpha.

Value

a list containing :

pvalues

pvalues of the different test (without correction)

adjPvalues

adjusted pvalues

groupId

Index of the group

hierMatrix

Matrix describing the hierarchical tree.

References

Yekutieli, Daniel. "Hierarchical False Discovery Rate-Controlling Methodology." Journal of the American Statistical Association 103.481 (2008): 309-16.

See Also

selFDR, hierarchicalFWER

Examples

1
2
3
4
5
set.seed(42)
X = simuBlockGaussian(50,12,5,0.7)
y = drop(X[,c(2,7,12)]%*%c(2,2,-2)+rnorm(50,0,0.5))
res = HCgglasso(X,y)
test = hierarchicalFDR(X, y, res$group[[20]], res$var[[20]])

HCgglasso documentation built on May 2, 2019, 4:54 p.m.