Renext (version 3.0-0)

LRExp.test: Likelihood Ratio test of exponentiality vs. GPD

Description

Likelihood Ratio test of exponentiality vs. GPD.

Usage

LRExp.test(x,
              alternative = c("lomax", "GPD", "gpd", "maxlo"),
              method = c("num", "sim", "asymp"),
              nSamp = 15000,
              simW = FALSE)

Arguments

x

Numeric vector of sample values.

alternative

Character string describing the alternative distribution.

method

Method used to compute the $p$-value.

nSamp

Number of samples for a simulation, if method is "sim".

simW

Logical. If this is set to TRUE and method is "sim", the simulated values are returned as an element W in the list.

Value

• A list of results with elements statistic, p.value and method. Other elements are
• alternative
  Character describing the alternative hypothesis.
• W
  If simW is TRUE and method is "sim" only. A vector of nSamp simulated values of the statistic $W := -2 \log \textrm{LR}$.

Details

The Lomax and maxlo alternatives correspond to a GPD alternative with positive shape parameter $\xi > 0$ (Lomax) and GPD with $\xi < 0$ (maxlo). The asymptotic distribution of the Likelihood-ratio statistic is known. For the GPD alternative, this is a chi-square distribution with one df. For the Lomax alternative, this is the distribution of a product $BC$ where $B$ and $C$ are two independent random variables following a Bernoulli distribution with probability parameter $p = 0.5$ and a chi-square distribution with one df.

• Whenmethodis"num", a numerical approximation of the distribution is used. This method is not unlike that used by Kozubowski et al., but a different approximation is used. However, ifxhas a length$n > 500$, the method is turned to"asymp".
• Whenmethodis"sim",nSampsamples of the exponential distribution with the same size asxare drawn and the LR statistic is computed for each sample. The$p$-value is simply the estimated probability that a simulated LR is greater than the observed LR.
• Finally whenmethodis"asymp", the asymptotic distribution is used.

References

T.J. Kozubowski, A. K. Panorska, F. Qeadan, A. Gershunov and D. Rominger (2009) "Testing Exponentiality Versus Pareto Distribution via Likelihood Ratio" Comm. Statist. Simulation Comput. 38(1), pp. 118-139.

The approximation method used is described in the Renext Computing Details report.

See Also

Lomax, Maxlo, GPD for the alternatives used here.

Examples

set.seed(1234)
x <- rGPD(n = 50, loc = 0, scale = 1, shape = 0.1)
LRExp.test(x, method = "num")$p.value
LRExp.test(x, method = "asymp")$p.value
## requires much time
LRExp.test(x, method = "asymp")$p.value