eel.test1: Empirical and exponential empirical likelihood tests for one...

View source: R/empirical_likelihood_tests.R

Empirical and exponential empirical likelihood tests for one sampleR Documentation

Empirical and exponential empirical likelihood tests for one sample

Description

Empirical and exponential empirical likelihood tests for one sample.

Usage

eel.test1(x, mu, tol = 1e-09, logged = FALSE)
el.test1(x, mu, tol = 1e-07, logged = FALSE) 

Arguments

x

A numerical vector.

mu

The hypothesised mean value.

tol

The tolerance value to stop the iterations of the Newton-Raphson.

logged

Should the logarithm of the p-value be returned? TRUE or FALSE.

Details

Exponential empirical likelihood is a non parametric method. In this case we use it as the non parametric alternative to the t-test. Newton-Raphson is used to maximise the log-likelihood ratio test statistic. In the case of no solution, NULL is returned. Despite the function having beeen written in R, it is pretty fast. As for the empirical likelihood ratio test, there is a condition for the range of possible values of mu. If mu is outside this range it is rejected immediately.

Value

iters

The number of iterations required by the Newton-Raphson algorithm. If no covnergence occured this is NULL. This is not returned for the empircial likelihood ratio test.

info

A vector with three elements, the value of the \lambda, the likelihood ratio test statistic and the relevant p-value. If no convergence occured, the value of the \lambda before is becomes NA, the value of test statistic is 10^5 and the p-value is 0. No convergence can be interpreted as rejection of the hypothesis test.

p

The estimated probabilities, one for each observation. If no covnergence occured this is NULL.

Author(s)

Michail Tsagris

R implementation and documentation: Michail Tsagris <mtsagris@uoc.gr> and Manos Papadakis <papadakm95@gmail.com>.

References

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

See Also

ftest, ttest1

Examples

x <- rnorm(500)
a1 <- eel.test1(x, 0)
a2 <- el.test1(x, 0)

Rfast documentation built on Nov. 9, 2023, 5:06 p.m.