HMT {MLGL} | R Documentation |
MLGL
objectHierarchical Multiple Testing procedure
Apply Hierarchical Multiple Testing procedure on a MLGL
object
HMT( res, X, y, control = c("FWER", "FDR"), alpha = 0.05, test = partialFtest, ... )
res |
|
X |
matrix of size n*p |
y |
vector of size n. |
control |
either "FDR" or "FWER" |
alpha |
control level for testing procedure |
test |
test used in the testing procedure. Default is partialFtest |
... |
extra parameters fpr selFDR |
a list containing :
lambda values maximizing the number of rejects
A vector containing the index of selected variables for the first lambdaOpt
value
A vector containing the values index of selected groups for the first lambdaOpt
value
Selected groups for the first lambdaOpt
value
Selected groups for all lambda values
Control level
Test used in the testing procedure
"FDR" or "FWER"
Elapsed time
list containing the output of the testing function for each lambda. Each element can be used with the selFWER or selFDR functions.
lambda path
Number of groups before testing
Numer of groups after testing
hierarchicalFWER hierarchicalFDR selFWER selFDR
set.seed(42) X <- simuBlockGaussian(50, 12, 5, 0.7) y <- X[,c(2,7,12)] %*% c(2,2,-2) + rnorm(50, 0, 0.5) res <- MLGL(X, y) # perform hierarchical testing with FWER control out <- HMT(res, X, y, alpha = 0.05) # test a new value of alpha for a specific lambda selFWER(out$hierTest[[60]], alpha = 0.1)