Tutorials › Statistical Theory I - Probability and Distribution Theory › Conditional Laws Given an Order Statistic

Order Statistics and Sample Extremes · Tutorial 460 of 994

Conditional Laws Given an Order Statistic

Learn why, given the kth order statistic, the observations below and above it follow independent truncated versions of the parent distribution.

Intermediate 10 min read

What You'll Learn

  • State the conditions under which the conditional laws given an order statistic are defined.
  • Derive the conditional joint density of the ordered observations below and above the kth value.
  • Use the lower- and upper-truncated parent distributions to calculate conditional probabilities.
  • Find the conditional distribution of counts in intervals on either side of the order statistic.
  • Avoid confusing conditioning on the kth value with truncating all observations in the sample.

Why Conditioning Splits the Sample

In an iid continuous sample, knowing the value of the kth order statistic gives more information than simply knowing that one observation equals that value. It tells us that exactly \(k-1\) observations lie below it and \(n-k\) lie above it. The values on each side still have a tractable distribution: they behave like samples from the parent distribution restricted to the corresponding side of the observed order statistic.

The precise statement concerns the ordered observations below and above the kth value. Conditional on that value, the lower group has the distribution of \(k-1\) draws from the parent distribution truncated below the value, and the upper group has the distribution of \(n-k\) draws truncated above it. These two ordered groups are conditionally independent.

Truncated Parent Distributions

Let \(X_1,\ldots,X_n\) be iid with a continuous distribution function \(F\) and density \(f\). Write \(X_{(1)}\le\cdots\le X_{(n)}\) for their order statistics. Fix \(k\) with \(1\le k\le n\). At a value \(x\) with \(0<F(x)<1\), the parent distribution restricted to values below \(x\) has conditional CDF \(F(t)/F(x)\) for \(t<x\). Restricted to values above \(x\), it has conditional CDF \((F(t)-F(x))/(1-F(x))\) for \(t>x\).

Definition: The lower-truncated parent law at \(x\) is the conditional law of \(X_1\) given \(X_1<x\), with density \(f(t)/F(x)\) for \(t<x\) and zero otherwise. The upper-truncated parent law at \(x\) is the conditional law of \(X_1\) given \(X_1>x\), with density \(f(t)/(1-F(x))\) for \(t>x\) and zero otherwise. Continuity makes the choice of strict or non-strict inequality immaterial.

The conditioning event \(X_{(k)}=x\) has probability zero for a continuous sample. Conditional distributions are therefore understood through conditional densities, which specify a regular conditional law for almost every \(x\) at which the density of \(X_{(k)}\) is positive. The formulas below give a natural version at any such \(x\) with \(0<F(x)<1\) and \(f(x)>0\).

Conditional Joint Law

The central result identifies the full conditional law of both ordered groups, not just their individual marginal distributions.

Theorem: For \(x\) with \(f_{X_{(k)}}(x)>0\) and \(0<F(x)<1\), the conditional density of \((X_{(1)},\ldots,X_{(k-1)},X_{(k+1)},\ldots,X_{(n)})\), given \(X_{(k)}=x\), is
$$ \frac{(k-1)!\prod_{i=1}^{k-1} f(u_i)}{F(x)^{k-1}}\, \frac{(n-k)!\prod_{j=1}^{n-k} f(v_j)}{[1-F(x)]^{n-k}} $$

on the region \(u_1<\cdots<u_{k-1}<x<v_1<\cdots<v_{n-k}\), and is zero elsewhere. Empty products are one. Consequently, the ordered lower group and ordered upper group are conditionally independent. Each group has the order-statistic distribution from an iid sample of the relevant size from its corresponding truncated parent law.

Proof. The joint density of all \(n\) order statistics is \(n!\prod_{i=1}^n f(x_i)\) on \(x_1<\cdots<x_n\). Integrating over the \(k-1\) ordered coordinates below \(x\) and the \(n-k\) ordered coordinates above \(x\) gives the marginal density

$$ f_{X_{(k)}}(x) = \frac{n!}{(k-1)!(n-k)!} F(x)^{k-1}[1-F(x)]^{n-k}f(x). $$

Here the integral of \(\prod f(u_i)\) over \(u_1<\cdots<u_{k-1}<x\) is \(F(x)^{k-1}/(k-1)!\): the integral over the full \((k-1)\)-dimensional region with every coordinate below \(x\) is \(F(x)^{k-1}\), and its \((k-1)!\) orderings have equal integrals. The corresponding upper integral is \([1-F(x)]^{n-k}/(n-k)!\). Dividing the full joint density \(n!f(x)\prod f(u_i)\prod f(v_j)\) by this marginal density cancels \(n!f(x)\) and yields the displayed conditional density. Its factorization into a function of the \(u_i\)'s and a function of the \(v_j\)'s proves conditional independence. Each factor is the usual ordered-sample density formed from the appropriate truncated density. \(\square\)

In particular, for \(a<x\), one observation drawn from the lower-truncated law satisfies \(P(X_1\le a\mid X_1<x)=F(a)/F(x)\). For \(b>x\), one observation drawn from the upper-truncated law satisfies \(P(X_1>b\mid X_1>x)=[1-F(b)]/[1-F(x)]\). Conditional probabilities involving several observations use these probabilities together with the appropriate group size and order.

Conditional Interval Counts

An immediate consequence is a binomial law for counts within either truncated group. For example, let \(a<b<x\), and let \(N_-\) count how many of the \(k-1\) observations below \(X_{(k)}\) fall in \((a,b]\). Conditional on \(X_{(k)}=x\), each lower-group draw has probability \([F(b)-F(a)]/F(x)\) of falling in that interval.

Result: Under the conditions of the theorem,
$$ N_-\mid X_{(k)}=x \sim \operatorname{Bin}\!\left(k-1,\frac{F(b)-F(a)}{F(x)}\right), \qquad a<b<x. $$

For an interval \(x<a<b\) in the upper group, the corresponding count has law \(\operatorname{Bin}(n-k,[F(b)-F(a)]/[1-F(x)])\).

Derivation. The theorem says that the lower group has the distribution of the order statistics of \(k-1\) iid draws from the lower-truncated law. Counting observations in \((a,b]\) is unaffected by sorting. Each draw has probability \(p=[F(b)-F(a)]/F(x)\) of falling in the interval. Independence of the draws therefore gives \(P(N_-=r\mid X_{(k)}=x)=\binom{k-1}{r}p^r(1-p)^{k-1-r}\) for \(r=0,\ldots,k-1\), which is the stated binomial law. The upper-group result follows by the same argument using the upper-truncated law. Since the two ordered groups are conditionally independent, counts formed separately from the lower and upper groups are conditionally independent as well. \(\square\)

Worked Conditional Probabilities

Worked Example: Requiring the Lower Observations to Exceed a Threshold

Let \(X_1,\ldots,X_5\) be iid exponential random variables with rate \(2\), so \(F(t)=1-e^{-2t}\) for \(t\ge0\). Find \(P(X_{(1)}>0.4\mid X_{(3)}=1)\).

Given \(X_{(3)}=1\), the two observations below 1 have the lower-truncated parent law at 1. For one such draw, the probability of exceeding \(0.4\) is

$$ \frac{F(1)-F(0.4)}{F(1)} = \frac{(1-e^{-2})-(1-e^{-0.8})}{1-e^{-2}} = \frac{e^{-0.8}-e^{-2}}{1-e^{-2}} \approx 0.3632. $$

The event \(X_{(1)}>0.4\) requires both lower observations to exceed \(0.4\). They are conditionally independent draws before sorting, so

$$ P(X_{(1)}>0.4\mid X_{(3)}=1) = \left(\frac{e^{-0.8}-e^{-2}}{1-e^{-2}}\right)^2 \approx 0.1319. $$

The result is the probability that both of the two lower observations fall in \((0.4,1)\), given that they are below 1.

Worked Example: A Conditional Probability in a Normal Sample

Suppose \(X_1,\ldots,X_4\) are iid \(N(0,1)\). Given \(X_{(2)}=0\), find the probability that the largest observation exceeds 1.

There are \(n-k=2\) observations above \(X_{(2)}\). Each has the standard normal distribution conditional on being greater than 0. Let \(\Phi\) denote the standard normal CDF. For one upper-truncated draw,

$$ P(X_1\le 1\mid X_1>0) = \frac{\Phi(1)-\Phi(0)}{1-\Phi(0)} = \frac{0.84134-0.5}{0.5} \approx 0.68268. $$

The largest observation is at most 1 exactly when both upper observations are at most 1. Thus

$$ P(X_{(4)}>1\mid X_{(2)}=0) = 1-(0.68268)^2 \approx 0.53394. $$

The complement is essential here: the event that the largest exceeds 1 means at least one of the two upper observations exceeds 1.

Worked Example: Counting Lower Observations in an Interval

Let \(X_1,\ldots,X_4\) be iid with density \(f(t)=2t\) on \(0<t<1\), so \(F(t)=t^2\). Given \(X_{(2)}=0.5\), find the probability that the one observation below \(X_{(2)}\) is at most \(0.25\).

The lower group has size \(k-1=1\), and its parent law is truncated below \(0.5\). Therefore,

$$ P(X_{(1)}\le 0.25\mid X_{(2)}=0.5) = \frac{F(0.25)}{F(0.5)} = \frac{(0.25)^2}{(0.5)^2} = \frac{1/16}{1/4} = \frac14. $$

For a count version, the number of lower observations in \((0.25,0.5)\) is exactly one with probability \([F(0.5)-F(0.25)]/F(0.5)=3/4\). These probabilities sum to one, as they must, because the single lower observation lies either at most \(0.25\) or in \((0.25,0.5)\), up to a zero-probability boundary.

Common Mistake and Why the Result Matters

A common error is to treat all \(n\) observations as independent draws from the parent distribution truncated at \(x\). That cannot be right: given \(X_{(k)}=x\), exactly \(k-1\) observations are below \(x\), exactly \(n-k\) are above it, and the kth observation is fixed at \(x\). The lower and upper truncated laws apply to their respective groups, with these different group sizes.

A second pitfall is to use an ordinary conditional-probability ratio with \(P(X_{(k)}=x)\) in the denominator. That probability is zero in a continuous model. The valid calculation conditions through the density of \(X_{(k)}\), as in the derivation above. The resulting formulas make conditional questions about extremes, quantiles, and observations near a sample rank into ordinary calculations under truncated distributions.

Check Your Understanding

Use the conditional truncated laws to answer these questions.

  1. Given \(X_{(k)}=x\), how many observations are in the lower group, and what parent law governs each unsorted lower-group draw?
  2. For \(a<b<x\), state the conditional probability that one lower-group draw falls in \((a,b]\).
  3. If \(n=7\) and \(k=4\), what is the conditional distribution of the number of observations below \(X_{(4)}=x\) that fall in an interval \((a,b]\) with \(a<b<x\)?
  4. Why is it invalid to calculate \(P(A\mid X_{(k)}=x)\) by dividing \(P(A\cap\{X_{(k)}=x\})\) by \(P(X_{(k)}=x)\)?
  5. In the normal example, why was the probability that the maximum exceeds 1 calculated using a complement?