Empirical likelihood for a one sample mean vector hypothesis testing {mvhtests}

Empirical likelihood for a one sample mean vector hypothesis testing

Description

Empirical likelihood for a one sample mean vector hypothesis testing.

Usage

el.test1(x, mu, R = 1, ncores = 1, graph = FALSE)

Arguments

x

A matrix containing Euclidean data.

mu

The hypothesized mean vector.

R

If R is 1 no bootstrap calibration is performed and the classical p-value via the \chi^2 distribution is returned. If R is greater than 1, the bootstrap p-value is returned.

ncores

The number of cores to use, set to 1 by default.

graph

A boolean variable which is taken into consideration only when bootstrap calibration is performed. IF TRUE the histogram of the bootstrap test statistic values is plotted.

Details

The H_0 is that \pmb{\mu} = \pmb{\mu}_0 and the constraint imposed by EL is

\frac{1}{n}\sum_{i=1}^{n}\left\lbrace\left[1+\pmb{\lambda}^T\left({\bf x}_i-\pmb{\mu}_0 \right)\right]^{-1}\left({\bf x}_i-\pmb{\mu}_0\right)\right\rbrace={\bf 0},

where the \pmb{\lambda} is the Lagrangian parameter introduced to maximize the above expression. Note that the maximization of is with respect to the \pmb{\lambda}. The probabilities have the following form

p_i=\frac{1}{n}\left[1+\pmb{\lambda}^T \left({\bf x}_i-\pmb{\mu}_0 \right)\right]^{-1}.

The log-likelihood ratio test statistic can be written as

\Lambda=\sum_{i=1}^{n}\log{np_i}.

where d denotes the number of variables. Under H_0 \Lambda \sim \chi^2_d, asymptotically. Alternatively the bootstrap p-value may be computed.

Value

A list with the outcome of the function el.test which includes the -2 log-likelihood ratio, the observed P-value by chi-square approximation, the final value of Lagrange multiplier \lambda, the gradient at the maximum, the Hessian matrix, the weights on the observations (probabilities multiplied by the sample size) and the number of iteration performed. In addition the runtime of the procedure is reported. In the case of bootstrap, the bootstrap p-value is also returned.

Author(s)

Michail Tsagris.

R implementation and documentation: Michail Tsagris [e-mail].

References

Jing B.Y. and Wood A.T.A. (1996). Exponential empirical likelihood is not Bartlett correctable. Annals of Statistics, 24(1): 365–369.

Owen A. (1990). Empirical likelihood ratio confidence regions. Annals of Statistics, 18(1): 90–120.

Owen A.B. (2001). Empirical likelihood. Chapman and Hall/CRC Press.

See Also

eel.test1, hotel1T2, james, hotel2T2, maov, el.test2

Examples

x <- as.matrix(iris[, 1:4])
el.test1(x, mu = numeric(4) )
eel.test1(x, mu = numeric(4) )