Monotone likelihood ratio

A monotonic likelihood ratio in distributions

f(x)

and

g(x)

The ratio of the density functions above is increasing in the parameter $x$ , so $f(x)/g(x)$ satisfies the monotone likelihood ratio property.

In statistics, the monotone likelihood ratio property is a property of the ratio of two probability density functions (PDFs). Formally, distributions ƒ(x) and g(x) bear the property if

{\text{for every }}x_{1}>x_{0},\quad {\frac {f(x_{1})}{g(x_{1})}}\geq {\frac {f(x_{0})}{g(x_{0})}}

that is, if the ratio is nondecreasing in the argument $x$ .

If the functions are first-differentiable, the property may sometimes be stated

{\frac {\partial }{\partial x}}\left({\frac {f(x)}{g(x)}}\right)\geq 0

For two distributions that satisfy the definition with respect to some argument x, we say they "have the MLRP in x." For a family of distributions that all satisfy the definition with respect to some statistic T(X), we say they "have the MLR in T(X)."

Intuition

The MLRP is used to represent a data-generating process that enjoys a straightforward relationship between the magnitude of some observed variable and the distribution it draws from. If $f(x)$ satisfies the MLRP with respect to $g(x)$ , the higher the observed value $x$ , the more likely it was drawn from distribution $f$ rather than $g$ . As usual for monotonic relationships, the likelihood ratio's monotonicity comes in handy in statistics, particularly when using maximum-likelihood estimation. Also, distribution families with MLR have a number of well-behaved stochastic properties, such as first-order stochastic dominance and increasing hazard ratios. Unfortunately, as is also usual, the strength of this assumption comes at the price of realism. Many processes in the world do not exhibit a monotonic correspondence between input and output.

Example: Working hard or slacking off

Suppose you are working on a project, and you can either work hard or slack off. Call your choice of effort $e$ and the quality of the resulting project $q$ . If the MLRP holds for the distribution of q conditional on your effort $e$ , the higher the quality the more likely you worked hard. Conversely, the lower the quality the more likely you slacked off.

Choose effort $e\in \{H,L\}$ where H means high, L means low
Observe $q$ drawn from $f(q\mid e)$ . By Bayes' law with a uniform prior,
$Pr[e=H\mid q]={\frac {f(q\mid H)}{f(q\mid H)+f(q\mid L)}}$
Suppose $f(q\mid e)$ satisfies the MLRP. Rearranging, the probability the worker worked hard is

{\frac {1}{1+f(q\mid L)/f(q\mid H)}}

which, thanks to the MLRP, is monotonically increasing in

q

. Hence if some employer is doing a "performance review" he can infer his employee's behavior from the merits of his work.

Families of distributions satisfying MLR

Statistical models often assume that data are generated by a distribution from some family of distributions and seek to determine that distribution. This task is simplified if the family has the Monotone Likelihood Ratio Property (MLRP).

A family of density functions $\{f_{\theta }(x)\}_{{\theta \in \Theta }}$ indexed by a parameter $\theta$ taking values in an ordered set $\Theta$ is said to have a monotone likelihood ratio (MLR) in the statistic $T(X)$ if for any $\theta _{1}<\theta _{2}$ ,

{\frac {f_{{\theta _{2}}}(X=x_{1},x_{2},x_{3},\dots )}{f_{{\theta _{1}}}(X=x_{1},x_{2},x_{3},\dots )}}

is a non-decreasing function of

T(X)

Then we say the family of distributions "has MLR in $T(X)$ ".

List of families

Family	$T(X)$ in which $f_{\theta }(X)$ has the MLR
Exponential $[\lambda ]$	$\sum x_{i}$ observations
Binomial $[n,p]$	$\sum x_{i}$ observations
Poisson $[\lambda ]$	$\sum x_{i}$ observations
Normal $[\mu ,\sigma ]$	if $\sigma$ known, $\sum x_{i}$ observations

Hypothesis testing

If the family of random variables has the MLRP in $T(X)$ , a uniformly most powerful test can easily be determined for the hypotheses $H_{0}:\theta \leq \theta _{0}$ versus $H_{1}:\theta >\theta _{0}$ .

Example: Effort and output

Example: Let $e$ be an input into a stochastic technology --- worker's effort, for instance --- and $y$ its output, the likelihood of which is described by a probability density function $f(y;e).$ Then the monotone likelihood ratio property (MLRP) of the family $f$ is expressed as follows: for any $e_1,e_2$ , the fact that $e_{2}>e_{1}$ implies that the ratio $f(y;e_{2})/f(y;e_{1})$ is increasing in $y$ .

Relation to other statistical properties

Monotone likelihoods are used in several areas of statistical theory, including point estimation and hypothesis testing, as well as in probability models.

Exponential families

One-parameter exponential-families have monotone likelihood-functions. In particular, the one-dimensional exponential-family of probability density functions or probability mass functions with

f_\theta(x) = c(\theta)h(x)\exp(\pi(\theta)T(x))

has a monotone non-decreasing likelihood ratio in the sufficient statistic T(x), provided that $\pi(\theta)$ is non-decreasing.

Most powerful tests: The Karlin–Rubin theorem

Monotone likelihood-functions are used to construct uniformly most powerful tests, according to the Karlin–Rubin theorem.^[1] Consider a scalar measurement having a probability density function parameterized by a scalar parameter θ, and define the likelihood ratio $l(x) = f_{\theta_1}(x) / f_{\theta_0}(x)$ . If $l(x)$ is monotone non-decreasing, in $x$ , for any pair $\theta_1 \geq \theta_0$ (meaning that the greater $x$ is, the more likely $H_{1}$ is), then the threshold test:

\phi(x) = \begin{cases} 1 & \text{if } x > x_0 \\ 0 & \text{if } x < x_0 \end{cases}

where

x_{0}

is chosen such that

\operatorname{E}_{\theta_0}\phi(X)=\alpha

is the UMP test of size α for testing $H_0: \theta \leq \theta_0 \text{ vs. } H_1: \theta > \theta_0 .$

Note that exactly the same test is also UMP for testing $H_0: \theta = \theta_0 \text{ vs. } H_1: \theta > \theta_0 .$

Median unbiased estimation

Monotone likelihood-functions are used to construct median-unbiased estimators, using methods specified by Johann Pfanzagl and others.^[2]^[3] One such procedure is an analogue of the Rao--Blackwell procedure for mean-unbiased estimators: The procedure holds for a smaller class of probability distributions than does the Rao--Blackwell procedure for mean-unbiased estimation but for a larger class of loss functions.^[4]

Lifetime analysis: Survival analysis and reliability

If a family of distributions $f_{\theta }(x)$ has the monotone likelihood ratio property in $T(X)$ ,

the family has monotone decreasing hazard rates in $\theta$ (but not necessarily in $T(X)$ )
the family exhibits the first-order (and hence second-order) stochastic dominance in $x$ , and the best Bayesian update of $\theta$ is increasing in $T(X)$ .

But not conversely: neither monotone hazard rates nor stochastic dominance imply the MLRP.

Proofs

Let distribution family $f_{\theta }$ satisfy MLR in x, so that for $\theta _{1}>\theta _{0}$ and $x_{1}>x_{0}$ :

{\frac {f_{{\theta _{1}}}(x_{1})}{f_{{\theta _{0}}}(x_{1})}}\geq {\frac {f_{{\theta _{1}}}(x_{0})}{f_{{\theta _{0}}}(x_{0})}},

or equivalently:

f_{{\theta _{1}}}(x_{1})f_{{\theta _{0}}}(x_{0})\geq f_{{\theta _{1}}}(x_{0})f_{{\theta _{0}}}(x_{1}).\,

Integrating this expression twice, we obtain:

1. To $x_{1}$ with respect to $x_{0}$ $\int _{{\min _{x}\in X}}^{{x_{1}}}f_{{\theta _{1}}}(x_{1})f_{{\theta _{0}}}(x_{0})\,dx_{0}$ $\geq \int _{{\min _{x}\in X}}^{{x_{1}}}f_{{\theta _{1}}}(x_{0})f_{{\theta _{0}}}(x_{1})\,dx_{0}$ integrate and rearrange to obtain ${\frac {f_{{\theta _{1}}}}{f_{{\theta _{0}}}}}(x)\geq {\frac {F_{{\theta _{1}}}}{F_{{\theta _{0}}}}}(x)$		2. From $x_{0}$ with respect to $x_{1}$ $\int _{{x_{0}}}^{{\max _{x}\in X}}f_{{\theta _{1}}}(x_{1})f_{{\theta _{0}}}(x_{0})\,dx_{1}$ $\geq \int _{{x_{0}}}^{{\max _{x}\in X}}f_{{\theta _{1}}}(x_{0})f_{{\theta _{0}}}(x_{1})\,dx_{1}$ integrate and rearrange to obtain ${\frac {1-F_{{\theta _{1}}}(x)}{1-F_{{\theta _{0}}}(x)}}\geq {\frac {f_{{\theta _{1}}}}{f_{{\theta _{0}}}}}(x)$

1. To

x_{1}

with respect to

x_{0}

\int _{{\min _{x}\in X}}^{{x_{1}}}f_{{\theta _{1}}}(x_{1})f_{{\theta _{0}}}(x_{0})\,dx_{0}

\geq \int _{{\min _{x}\in X}}^{{x_{1}}}f_{{\theta _{1}}}(x_{0})f_{{\theta _{0}}}(x_{1})\,dx_{0}

integrate and rearrange to obtain

{\frac {f_{{\theta _{1}}}}{f_{{\theta _{0}}}}}(x)\geq {\frac {F_{{\theta _{1}}}}{F_{{\theta _{0}}}}}(x)

2. From

x_{0}

with respect to

x_{1}

\int _{{x_{0}}}^{{\max _{x}\in X}}f_{{\theta _{1}}}(x_{1})f_{{\theta _{0}}}(x_{0})\,dx_{1}

\geq \int _{{x_{0}}}^{{\max _{x}\in X}}f_{{\theta _{1}}}(x_{0})f_{{\theta _{0}}}(x_{1})\,dx_{1}

integrate and rearrange to obtain

{\frac {1-F_{{\theta _{1}}}(x)}{1-F_{{\theta _{0}}}(x)}}\geq {\frac {f_{{\theta _{1}}}}{f_{{\theta _{0}}}}}(x)

First-order stochastic dominance

Combine the two inequalities above to get first-order dominance:

F_{{\theta _{1}}}(x)\leq F_{{\theta _{0}}}(x)\ \forall x

Monotone hazard rate

Use only the second inequality above to get a monotone hazard rate:

{\frac {f_{{\theta _{1}}}(x)}{1-F_{{\theta _{1}}}(x)}}\leq {\frac {f_{{\theta _{0}}}(x)}{1-F_{{\theta _{0}}}(x)}}\ \forall x

Uses

Economics

The MLR is an important condition on the type distribution of agents in mechanism design. Most solutions to mechanism design models assume a type distribution to satisfy the MLR to take advantage of a common solution method.

References

↑ Casella, G.; Berger, R.L. (2008), Statistical Inference, Brooks/Cole. ISBN 0-495-39187-5 (Theorem 8.3.17)
↑ Pfanzagl, Johann. "On optimal median unbiased estimators in the presence of nuisance parameters." The Annals of Statistics (1979): 187-193.
↑ Brown, L. D.; Cohen, Arthur; Strawderman, W. E. A Complete Class Theorem for Strict Monotone Likelihood Ratio With Applications. Ann. Statist. 4 (1976), no. 4, 712--722. doi:10.1214/aos/1176343543. http://projecteuclid.org/euclid.aos/1176343543.
↑ Page 713: Brown, L. D.; Cohen, Arthur; Strawderman, W. E. A Complete Class Theorem for Strict Monotone Likelihood Ratio With Applications. Ann. Statist. 4 (1976), no. 4, 712--722. doi:10.1214/aos/1176343543. http://projecteuclid.org/euclid.aos/1176343543.

Statistics

Descriptive statistics

Continuous data

Center	Mean arithmetic geometric harmonic Median Mode

Dispersion	Variance Standard deviation Coefficient of variation Percentile Range Interquartile range

Shape	Moments Skewness Kurtosis L-moments

Count data

Index of dispersion

Summary tables

Dependence

Graphics

Data collection

Study design	Population Statistic Effect size Statistical power Sample size determination Missing data

Survey methodology	Sampling Standard error stratified cluster Opinion poll Questionnaire

Controlled experiments	Design control optimal Controlled trial Randomized Random assignment Replication Blocking Interaction Factorial experiment

Uncontrolled studies	Observational study Natural experiment Quasi-experiment

Statistical inference

Statistical theory

Frequentist inference

Point estimation	Estimating equations Maximum likelihood Method of moments M-estimator Minimum distance Unbiased estimators Mean-unbiased minimum-variance Rao–Blackwellization Lehmann–Scheffé theorem Median unbiased Plug-in

Interval estimation	Confidence interval Pivot Likelihood interval Prediction interval Tolerance interval Resampling Bootstrap Jackknife

Testing hypotheses	1- & 2-tails Power Uniformly most powerful test Permutation test Randomization test Multiple comparisons

Parametric tests	Likelihood-ratio Wald Score

Specific tests

Z (normal) Student's t-test F

Goodness of fit	Chi-squared Kolmogorov–Smirnov Anderson–Darling Normality (Shapiro–Wilk) Likelihood-ratio test Model selection Cross validation AIC BIC

Rank statistics	Sign Sample median Signed rank (Wilcoxon) Hodges–Lehmann estimator Rank sum (Mann–Whitney) Nonparametric anova 1-way (Kruskal–Wallis) 2-way (Friedman) Ordered alternative (Jonckheere–Terpstra)

Bayesian inference

Correlation	Pearson product–moment Partial correlation Confounding variable Coefficient of determination

Regression analysis	Errors and residuals Regression model validation Mixed effects models Simultaneous equations models Multivariate adaptive regression splines (MARS)

Linear regression	Simple linear regression Ordinary least squares General linear model Bayesian regression

Non-standard predictors	Nonlinear regression Nonparametric Semiparametric Isotonic Robust Heteroscedasticity Homoscedasticity

Generalized linear model	Exponential families Logistic (Bernoulli) / Binomial / Poisson regressions

Partition of variance	Analysis of variance (ANOVA, anova) Analysis of covariance Multivariate ANOVA Degrees of freedom

Categorical / Multivariate / Time-series / Survival analysis

Categorical

Multivariate

Time-series

General	Decomposition Trend Stationarity Seasonal adjustment Exponential smoothing Cointegration Structural break Granger causality

Specific tests	Dickey–Fuller Johansen Q-statistic (Ljung–Box) Durbin–Watson Breusch–Godfrey

Time domain	Autocorrelation (ACF) partial (PACF) Cross-correlation (XCF) ARMA model ARIMA model (Box–Jenkins) Autoregressive conditional heteroskedasticity (ARCH) Vector autoregression (VAR)

Frequency domain	Spectral density estimation Fourier analysis Wavelet

Survival

Survival function	Kaplan–Meier estimator (product limit) Proportional hazards models Accelerated failure time (AFT) model First hitting time

Hazard function	Nelson–Aalen estimator

Test	Log-rank test

Applications

Biostatistics	Bioinformatics Clinical trials / studies Epidemiology Medical statistics

Engineering statistics	Chemometrics Methods engineering Probabilistic design Process / quality control Reliability System identification

Social statistics	Actuarial science Census Crime statistics Demography Econometrics National accounts Official statistics Population statistics Psychometrics

Spatial statistics	Cartography Environmental statistics Geographic information system Geostatistics Kriging

Category
Portal
Commons
WikiProject

Theory of probability distributions

probability mass function (pmf) probability density function (pdf) cumulative distribution function (cdf) quantile function

raw moment central moment mean variance standard deviation skewness kurtosis L-moment

moment-generating function (mgf) characteristic function probability-generating function (pgf) cumulant combinant

This article is issued from Wikipedia - version of the 11/22/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.