Last data update: 2014.03.03

R: Calculation of the quantiles for the mht procedure
quantilemhtR Documentation

Calculation of the quantiles for the mht procedure

Description

Calculation of the quantiles for the mht procedure

Usage

quantilemht(data,k,alpha,IT,maxq,sigma)

Arguments

data

Input matrix of dimension n*p; each row is an observation vector.

k

Positive integer.

alpha

A user supplied type I error sequence. Default is (0.1,0.05).

IT

Number of simulations in the calculation of the quantile. Default is 1000.

maxq

Number of maximum multiple hypotheses testing to do. Default is log(min(n,p)-1,2).

sigma

Value of the variance if it is known; 0 otherwise. Default is 0.

Details

Calculation of the quantile of the statistic that tests the null hypothesis Hk:"E(Y)=X_1β_1+X_2β_2+…+X_kβ_k", i.e there is no more signal remaining in the variables (k+1,...,p).

Value

quantile

The estimated quantiles.

nbrprob

Number of times there was not enough simulation to estimate the alpha-quantile.

Results