Empirical likelihood

In probability theory and statistics, empirical likelihood (EL) is a nonparametric method for estimating the parameters of statistical models. It requires fewer assumptions about the error distribution while retaining some of the merits in likelihood-based inference. The estimation method requires that the data are independent and identically distributed (iid). It performs well even when the distribution is asymmetric or censored. ^[1] EL methods can also handle constraints and prior information on parameters. Art Owen pioneered work in this area with his 1988 paper. ^[2]

Definition

Given a set of ${\displaystyle n}$ i.i.d. realizations ${\displaystyle y_{i}}$ of random variables ${\displaystyle Y_{i}}$, then the empirical distribution function is ${\displaystyle {\hat {F}}(y):=\sum _{i=1}^{n}\pi _{i}I(Y_{i}<y)}$, with the indicator function ${\displaystyle I}$ and the (normalized) weights ${\displaystyle \pi _{i}}$. Then, the empirical likelihood is: ^[3]

${\displaystyle L:=\prod _{i=1}^{n}{\frac {{\hat {F}}(y_{i})-{\hat {F}}(y_{i}-\delta y)}{\delta y}},}$

where ${\displaystyle \delta y}$ is a small number (potentially the difference to the next smaller sample).

Empirical likelihood estimation can be augmented with side information by using further constraints (similar to the generalized estimating equations approach) for the empirical distribution function. E.g. a constraint like the following can be incorporated using a Lagrange multiplier ${\displaystyle E[h(Y;\theta )]=\int _{-\infty }^{\infty }h(y;\theta )dF=0}$ which implies ${\displaystyle {\hat {E}}[h(y;\theta )]=\sum _{i=1}^{n}h(y_{i};\theta )\pi _{i}=0}$.

With similar constraints, we could also model correlation.

Discrete random variables

The empirical-likelihood method can also be also employed for discrete distributions. ^[4] Given ${\displaystyle \ p_{i}:={\hat {F}}(y_{i})-{\hat {F}}(y_{i}-\delta y),\ i=1,...,n}$ such that ${\displaystyle p_{i}\geq 0{\text{ and }}\sum _{i=1}^{n}\ p_{i}=1.}$

Then the empirical likelihood is again ${\displaystyle L(p_{1},...,p_{n})=\prod _{i=1}^{n}\ p_{i}}$.

Using the Lagrangian multiplier method to maximize the logarithm of the empirical likelihood subject to the trivial normalization constraint, we find ${\displaystyle p_{i}=1/n}$ as a maximum. Therefore, ${\displaystyle {\hat {F}}}$ is the empirical distribution function.

Estimation Procedure

EL estimates are calculated by maximizing the empirical likelihood function (see above) subject to constraints based on the estimating function and the trivial assumption that the probability weights of the likelihood function sum to 1. ^[5] This procedure is represented as:

${\displaystyle \max _{\pi _{i},\theta }\ln(L)=\max _{\pi _{i},\theta }\sum _{i=1}^{n}\ln \pi _{i}}$

subject to the constraints

${\displaystyle s.t.\sum _{i=1}^{n}\pi _{i}=1,\sum _{i=1}^{n}\pi _{i}h(y_{i};\theta )=0,\forall i\in [1..n]\quad 0\leq \pi _{i}.}$ ^{[6] : Equation (73)}

The value of the theta parameter can be found by solving the Lagrangian function

${\displaystyle {\mathcal {L}}=\sum _{i=1}^{n}\ln \pi _{i}+\mu (1-\sum _{i=1}^{n}\pi _{i})-n\tau '\sum _{i=1}^{n}\pi _{i}h(y_{i};\theta ).}$ ^{[6] : Equation (74)}

There is a clear analogy between this maximization problem and the one solved for maximum entropy.

The parameters ${\displaystyle \pi _{i}}$ are nuisance parameters.

Empirical Likelihood Ratio (ELR)

An empirical likelihood ratio function is defined and used to obtain confidence intervals parameter of interest θ similar to parametric likelihood ratio confidence intervals. ^{[7] [8]} Let L(F) be the empirical likelihood of function ${\displaystyle F}$, then the ELR would be:

${\displaystyle R(F)=L(F)/L(F_{n})}$.

Consider sets of the form

${\displaystyle C=\{T(F)|R(F)\geq r\}}$.

Under such conditions a test of ${\displaystyle T(F)=t}$ rejects when t does not belong to ${\displaystyle C}$, that is, when no distribution F with ${\displaystyle T(F)=t}$ has likelihood ${\displaystyle L(F)\geq rL(F_{n})}$.

The central result is for the mean of X. Clearly, some restrictions on ${\displaystyle F}$ are needed, or else ${\displaystyle C=\mathbb {R} ^{p}}$ whenever ${\displaystyle r<1}$. To see this, let:

${\displaystyle F=\epsilon \delta _{x}+(1-\epsilon )F_{n}}$

If ${\displaystyle \epsilon }$ is small enough and ${\displaystyle \epsilon >0}$, then ${\displaystyle R(F)\geq r}$.

But then, as ${\displaystyle x}$ ranges through ${\displaystyle \mathbb {R} ^{p}}$, so does the mean of ${\displaystyle F}$, tracing out ${\displaystyle C=\mathbb {R} ^{p}}$. The problem can be solved by restricting to distributions F that are supported in a bounded set. It turns out to be possible to restrict attention t distributions with support in the sample, in other words, to distribution ${\displaystyle F\ll F_{n}}$. Such method is convenient since the statistician might not be willing to specify a bounded support for ${\displaystyle F}$, and since ${\displaystyle t}$ converts the construction of ${\displaystyle C}$ into a finite dimensional problem.

Other Applications

The use of empirical likelihood is not limited to confidence intervals. In efficient quantile regression, an EL-based categorization ^[9] procedure helps determine the shape of the true discrete distribution at level p, and also provides a way of formulating a consistent estimator. In addition, EL can be used in place of parametric likelihood to form model selection criteria. ^[10] Empirical likelihood can naturally be applied in survival analysis ^[11] or regression problems ^[12]

See also

• Bootstrapping (statistics)
• Jackknife (statistics)
• Nelson–Aalen estimator

Literature

• Nordman, Daniel J., and Soumendra N. Lahiri. "A review of empirical likelihood methods for time series." Journal of Statistical Planning and Inference 155 (2014): 1-18. https://doi.org/10.1016/j.jspi.2013.10.001

References

1. Owen, Art B. (2001). Empirical likelihood. Boca Raton, Fla. ISBN   978-1-4200-3615-2. OCLC   71012491.
2. Owen, Art B. (1988). "Empirical likelihood ratio confidence intervals for a single functional". Biometrika. 75 (2): 237–249. doi:10.1093/biomet/75.2.237. ISSN   0006-3444.
3. ${\displaystyle {\frac {{\hat {F}}(y_{i})-{\hat {F}}(y_{i}-\delta y)}{\delta y}}}$ is an estimate of the probability density, compare histogram
4. Wang, Dong; Chen, Song Xi (2009-02-01). "Empirical likelihood for estimating equations with missing values". The Annals of Statistics. 37 (1). arXiv: 0903.0726 . doi:10.1214/07-aos585. ISSN   0090-5364. S2CID   5427751.
5. Mittelhammer, Judge, and Miller (2000), 292.
6. Bera, Anil K.; Bilias, Yannis (2002). "The MM, ME, ML, EL, EF and GMM approaches to estimation: a synthesis". Journal of Econometrics. 107 (1–2): 51–86. doi:10.1016/S0304-4076(01)00113-0.
7. Owen, Art (1990-03-01). "Empirical Likelihood Ratio Confidence Regions". The Annals of Statistics. 18 (1). doi: 10.1214/aos/1176347494 . ISSN   0090-5364.
8. Dong, Lauren Bin; Giles, David E. A. (2007-01-30). "An Empirical Likelihood Ratio Test for Normality". Communications in Statistics - Simulation and Computation. 36 (1): 197–215. doi:10.1080/03610910601096544. ISSN   0361-0918. S2CID   16866055.
9. Chen, Jien; Lazar, Nicole A. (2010-01-27). "Quantile estimation for discrete data via empirical likelihood". Journal of Nonparametric Statistics. 22 (2): 237–255. doi:10.1080/10485250903301525. ISSN   1048-5252. S2CID   119684596.
10. Chen, Chixiang; Wang, Ming; Wu, Rongling; Li, Runze (2022). "A Robust Consistent Information Criterion for Model Selection Based on Empirical Likelihood". Statistica Sinica. arXiv: 2006.13281 . doi:10.5705/ss.202020.0254. ISSN   1017-0405. S2CID   220042083.
11. Zhou, M. (2015). Empirical Likelihood Method in Survival Analysis (1st ed.). Chapman and Hall/CRC. https://doi.org/10.1201/b18598
12. Chen, Song Xi, and Ingrid Van Keilegom. "A review on empirical likelihood methods for regression." TEST volume 18, pages 415–447, (2009) https://doi.org/10.1007/s11749-009-0159-5