digplanet beta 1: Athena
Share digplanet:

Agriculture

Applied sciences

Arts

Belief

Chronology

Culture

Education

Environment

Geography

Health

History

Humanities

Language

Law

Life

Mathematics

Nature

People

Politics

Science

Society

Technology

In statistics, a likelihood ratio test is a statistical test used to compare the fit of two models, one of which (the null model) is a special case of the other (the alternative model). The test is based on the likelihood ratio, which expresses how many times more likely the data are under one model than the other. This likelihood ratio, or equivalently its logarithm, can then be used to compute a p-value, or compared to a critical value to decide whether to reject the null model in favour of the alternative model. When the logarithm of the likelihood ratio is used, the statistic is known as a log-likelihood ratio statistic, and the probability distribution of this test statistic, assuming that the null model is true, can be approximated using Wilks' theorem.

In the case of distinguishing between two models, each of which has no unknown parameters, use of the likelihood ratio test can be justified by the Neyman–Pearson lemma, which demonstrates that such a test has the highest power among all competitors.[1]

Use

Each of the two competing models, the null model and the alternative model, is separately fitted to the data and the log-likelihood recorded. The test statistic (often denoted by D) is twice the difference in these log-likelihoods:

\begin{align} D & = -2\ln\left( \frac{\text{likelihood for null model}}{\text{likelihood for alternative model}} \right) \\ &= -2\ln(\text{likelihood for null model}) + 2\ln(\text{likelihood for alternative model}) \\ \end{align}

The model with more parameters will always fit at least as well (have a weakly greater log-likelihood). Whether it fits significantly better and should thus be preferred is determined by deriving the probability or p-value of the difference D. Where the null hypothesis represents a special case of the alternative hypothesis, the probability distribution of the test statistic is approximately a chi-squared distribution with degrees of freedom equal to df2 − df1 .[2] Symbols df1 and df2 represent the number of free parameters of models 1 and 2, the null model and the alternative model, respectively. The test requires nested models, that is: models in which the more complex one can be transformed into the simpler model by imposing a set of constraints on the parameters.[3]

For example: if the null model has 1 free parameter and a log-likelihood of −8024 and the alternative model has 3 degrees of freedom and a LL of −8012, then the probability of this difference is that of chi-squared value of +2·(8024 − 8012) = 24 with 3 − 1 = 2 degrees of freedom. Certain assumptions must be met for the statistic to follow a chi-squared distribution and often empirical p-values are computed.

Background

The likelihood ratio, often denoted by $\Lambda$ (the capital Greek letter lambda), is the ratio of the likelihood function varying the parameters over two different sets in the numerator and denominator. A likelihood ratio test is a statistical test for making a decision between two hypotheses based on the value of this ratio.

It is central to the NeymanPearson approach to statistical hypothesis testing, and, like statistical hypothesis testing in general, is both widely used and criticized.[citation needed]

Simple-vs-simple hypotheses

A statistical model is often a parametrized family of probability density functions or probability mass functions $f(x|\theta)$. A simple-vs-simple hypotheses test has completely specified models under both the null and alternative hypotheses, which for convenience are written in terms of fixed values of a notional parameter $\theta$:

\begin{align} H_0 &:& \theta=\theta_0 ,\\ H_1 &:& \theta=\theta_1 . \end{align}

Note that under either hypothesis, the distribution of the data is fully specified; there are no unknown parameters to estimate. The likelihood ratio test statistic can be written as:[4][5]

$\Lambda(x) = \frac{ L(\theta_0|x) }{ L(\theta_1|x) } = \frac{ f(x|\theta_0) }{ f(x|\theta_1) }$

or

$\Lambda(x)=\frac{L(\theta_0\mid x)}{\sup\{\,L(\theta\mid x):\theta\in\{\theta_0,\theta_1\}\}},$

where $L(\theta|x)$ is the likelihood function. Note that some references may use the reciprocal as the definition.[6] In the form stated here, the likelihood ratio is small if the alternative model is better than the null model and the likelihood ratio test provides the decision rule as:

If $\Lambda > c$, do not reject $H_0$;
If $\Lambda < c$, reject $H_0$;
Reject with probability $q$ if $\Lambda = c .$

The values $c, \; q$ are usually chosen to obtain a specified significance level $\alpha$, through the relation: $q\cdot P(\Lambda=c \;|\; H_0) + P(\Lambda < c \; | \; H_0) = \alpha$.[citation needed] The Neyman-Pearson lemma states that this likelihood ratio test is the most powerful among all level $\alpha$ tests for this problem.[1]

Definition (likelihood ratio test for composite hypotheses)

A null hypothesis is often stated by saying the parameter $\theta$ is in a specified subset $\Theta_0$ of the parameter space $\Theta$.

\begin{align} H_0 &:& \theta \in \Theta_0\\ H_1 &:& \theta \in \Theta_0^{\complement} \end{align}

The likelihood function is $L(\theta|x) = f(x|\theta)$ (with $f(x|\theta)$ being the pdf or pmf) is a function of the parameter $\theta$ with $x$ held fixed at the value that was actually observed, i.e., the data. The likelihood ratio test statistic is [7]

$\Lambda(x)=\frac{\sup\{\,L(\theta\mid x):\theta\in\Theta_0\,\}}{\sup\{\,L(\theta\mid x):\theta\in\Theta\,\}}.$

Here, the $\sup$ notation refers to the Supremum function.

A likelihood ratio test is any test with critical region (or rejection region) of the form $\{x|\Lambda \le c\}$ where $c$ is any number satisfying $0\le c\le 1$. Many common test statistics such as the Z-test, the F-test, Pearson's chi-squared test and the G-test are tests for nested models and can be phrased as log-likelihood ratios or approximations thereof.

Interpretation

Being a function of the data $x$, the LR is therefore a statistic. The likelihood ratio test rejects the null hypothesis if the value of this statistic is too small. How small is too small depends on the significance level of the test, i.e., on what probability of Type I error is considered tolerable ("Type I" errors consist of the rejection of a null hypothesis that is true).

The numerator corresponds to the maximum likelihood of an observed outcome under the null hypothesis. The denominator corresponds to the maximum likelihood of an observed outcome varying parameters over the whole parameter space. The numerator of this ratio is less than the denominator. The likelihood ratio hence is between 0 and 1. Low values of the likelihood ratio mean that the observed result was less likely to occur under the null hypothesis as compared to the alternative. High values of the statistic mean that the observed outcome was nearly as likely to occur under the null hypothesis as compared to the alternative, and the null hypothesis cannot be rejected.

Distribution: Wilks' theorem

If the distribution of the likelihood ratio corresponding to a particular null and alternative hypothesis can be explicitly determined then it can directly be used to form decision regions (to accept/reject the null hypothesis). In most cases, however, the exact distribution of the likelihood ratio corresponding to specific hypotheses is very difficult to determine. A convenient result, attributed to Samuel S. Wilks, says that as the sample size $n$ approaches $\infty$, the test statistic $-2 \log(\Lambda)$ for a nested model will be asymptotically $\chi^2$ distributed with degrees of freedom equal to the difference in dimensionality of $\Theta$ and $\Theta_0$.[8] This means that for a great variety of hypotheses, a practitioner can compute the likelihood ratio $\Lambda$ for the data and compare $-2\log(\Lambda)$ to the chi squared value corresponding to a desired statistical significance as an approximate statistical test.

Examples

Coin tossing

An example, in the case of Pearson's test, we might try to compare two coins to determine whether they have the same probability of coming up heads. Our observation can be put into a contingency table with rows corresponding to the coin and columns corresponding to heads or tails. The elements of the contingency table will be the number of times the coin for that row came up heads or tails. The contents of this table are our observation $X$.

 Heads Tails Coin 1 $k_{1H}$ $k_{1T}$ Coin 2 $k_{2H}$ $k_{2T}$

Here $\Theta$ consists of the parameters $p_{1H}$, $p_{1T}$, $p_{2H}$, and $p_{2T}$, which are the probability that coins 1 and 2 come up heads or tails. The hypothesis space $H$ is defined by the usual constraints on a distribution, $0 \le p_{ij} \le 1$, and $p_{iH} + p_{iT} = 1$. The null hypothesis $H_0$ is the subspace where $p_{1j} = p_{2j}$. In all of these constraints, $i = 1,2$ and $j = H,T$.

Writing $n_{ij}$ for the best values for $p_{ij}$ under the hypothesis $H$, maximum likelihood is achieved with

$n_{ij} = \frac{k_{ij}}{k_{iH}+k_{iT}}.$

Writing $m_{ij}$ for the best values for $p_{ij}$ under the null hypothesis $H_0$, maximum likelihood is achieved with

$m_{ij} = \frac{k_{1j}+k_{2j}}{k_{1H}+k_{2H}+k_{1T}+k_{2T}},$

which does not depend on the coin $i$.

The hypothesis and null hypothesis can be rewritten slightly so that they satisfy the constraints for the logarithm of the likelihood ratio to have the desired nice distribution. Since the constraint causes the two-dimensional $H$ to be reduced to the one-dimensional $H_0$, the asymptotic distribution for the test will be $\chi^2(1)$, the $\chi^2$ distribution with one degree of freedom.

For the general contingency table, we can write the log-likelihood ratio statistic as

$-2 \log \Lambda = 2\sum_{i, j} k_{ij} \log \frac{n_{ij}}{m_{ij}}.$

References

1. ^ a b Neyman, Jerzy; Pearson, Egon S. (1933). "On the Problem of the Most Efficient Tests of Statistical Hypotheses". Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 231 (694–706): 289–337. doi:10.1098/rsta.1933.0009. JSTOR 91247. edit
2. ^ Huelsenbeck, J. P.; Crandall, K. A. (1997). "Phylogeny Estimation and Hypothesis Testing Using Maximum Likelihood". Annual Review of Ecology and Systematics 28: 437–466. doi:10.1146/annurev.ecolsys.28.1.437. edit
3. ^ An example using phylogenetic analyses is described at Huelsenbeck, J. P.; Hillis, D. M.; Nielsen, R. (1996). "A Likelihood-Ratio Test of Monophyly". Systematic Biology 45 (4): 546. doi:10.1093/sysbio/45.4.546. edit
4. ^ Mood, A.M.; Graybill, F.A. (1963) Introduction to the Theory of Statistics, 2nd edition. McGraw-Hill ISBN 978-0070428638 (page 286)
5. ^ Kendall, M.G., Stuart, A. (1973) The Advanced Theory of Statistics, Volume 2, Griffin. ISBN 0852642156 (page 234)
6. ^ Cox, D. R. and Hinkley, D. V Theoretical Statistics, Chapman and Hall, 1974. (page 92)
7. ^ Casella, George; Berger, Roger L. (2001) Statistical Inference, Second edition. ISBN 978-0534243128 (page 375)
8. ^ Wilks, S. S. (1938). "The Large-Sample Distribution of the Likelihood Ratio for Testing Composite Hypotheses". The Annals of Mathematical Statistics 9: 60–62. doi:10.1214/aoms/1177732360. edit
 156 videos foundNext >
 Cycle 3 - Topic 6: Likelihood Ratio TestIntroduction to the Likelihood Ratio Test; Discussion over #37C from ASM manual, SOA #22 and #37.13 from ASM manual. What Are Likelihood Ratios and How Are They UsedVideo describing the role of likelihood ratios in diagnostic testing. Likelihood Ratio.mp41. What is likelihood ratio? 2. What is use of likelihood ratio? 3. What distribution does likelihood ratio follows? The Likelihood Ratio TestThe likelihood ratio test maximizes the probability of correctly deciding hypothesis H_1 is true for any given probability of deciding H_0 is false when the ... The Generalized Likelihood Ratio TestThere is no universally optimal test strategy for composite hypotheses (unknown parameters in the pdfs). The generalized likelihood ratio test (GLRT) is a li... Modern repeated measures analysis using mixed models in SPSS (1)Repeated measures analyse an introduction to the Mixed models (random effects) option in SPSS. Demonstrates different Covariance matrix types & how to use th... Statistics, Probability and Diagnostic MedicinePresented by: Jennifer Le-Rademacher, PhD, Division of Biostatistics, Medical College of Wisconsin Diagnostic devices are used to identify or to predict the ... Tips for Teachers of Evidence-based Medicine: Making Sense of Diagnostic Test Results Using LikelihFrom the Springer article: Tips for Teachers of Evidence-based Medicine: Making Sense of Diagnostic Test Results Using Likelihood Ratios http://www.springerl... Modern repeated measures analysis using mixed models in SPSS (2)This uses a Repeated measures analyse as an introduction to the Mixed models (random effects) option in SPSS. Demonstrates different Covariance matrix types ... restricted vs full model likelihood ratio test October 18
 156 videos foundNext >
 2 news items
 Hospital outcomes for paediatric pneumonia and diarrhoea patients admitted in ... BMC Pediatrics Fri, 10 May 2013 17:31:42 -0700 Variables were retained in the model if a likelihood ratio test (LRT, p value of <0.05) supported improved model fit. A goodness of fit test was carried out and an area under receiver operating curve (AROC) calculated to determine if the final model ... Medscape Risk of Ischaemic Heart Disease in Patients With Inflammatory Bowel Disease Medscape Mon, 22 Apr 2013 13:19:58 -0700 *p Value for homogeneity test performed by a likelihood ratio test. †Partly adjusted: excluding first-year cases and adjusting for age, gender, socioeconomic status and calendar year. ‡Partly adjusted plus additional adjustment for use (time-dependent ...
 Limit to books that you can completely read online Include partial books (book previews) .gsc-branding { display:block; }

Oops, we seem to be having trouble contacting Twitter