Empirical and exponential empirical likelihood tests for one sample
Empirical and exponential empirical likelihood tests for one sample.
eel.test1(x, mu, tol = 1e-09, logged = FALSE) el.test1(x, mu, tol = 1e-07, logged = FALSE)
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. |
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.
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 λ, the likelihood ratio test statistic and the relevant p-value. If no convergence occured, the value of the λ 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. |
Michail Tsagris
R implementation and documentation: Michail Tsagris <mtsagris@yahoo.gr> and Manos Papadakis <papadakm95@gmail.com>.
Owen A. B. (2001). Empirical likelihood. Chapman and Hall/CRC Press.
x <- rnorm(500) system.time(a1 <- eel.test1(x, 0) ) system.time(a2 <- el.test1(x, 0) )
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.