Order Statistics - Statistics Graduate Level

  1. A random sample of size 55 is drawn from a continuous distribution with CDF F(x)=x2F(x)=x^2 for 0≤x≤10\le x\le 1. Let X(3)X_{(3)} denote the third order statistic. Which expression equals P(X(3)≤0.8)P(X_{(3)}\le 0.8)?
    Answer: ∑j=35(5j)(0.64)j(0.36)5−j≈0.7491\sum_{j=3}^{5}{5\choose j}(0.64)^j(0.36)^{5-j}\approx 0.7491 When you see a question involving order statistics from a non-uniform distribution, your first instinct should be to transform the problem into a binomial one. The key insight is that if XX has CDF F(x)F(x), then F(X)F(X) is uniformly distributed on [0,1][0,1]. So asking whether X(3)≤0.8X_{(3)} \le 0.8 is equivalent to asking whether at least 3 of the 5 observations fall at or below 0.8, where each observation independently has "success" probability p=F(0.8)=(0.8)2=0.64p = F(0.8) = (0.8)^2 = 0.64. This transforms the problem into a standard binomial calculation. Specifically, X(3)≤0.8X_{(3)} \le 0.8 if and only if at least 3 of the 5 values are ≤0.8\le 0.8, because the third-smallest value being ≤0.8\le 0.8 requires that three or more observations landed below 0.8. That gives you P(X(3)≤0.8)=∑j=35(5j)(0.64)j(0.36)5−j≈0.7491P(X_{(3)} \le 0.8) = \sum_{j=3}^{5}\binom{5}{j}(0.64)^j(0.36)^{5-j} \approx 0.7491, which is answer A. Answer B sums from j=2j=2, which corresponds to P(at least 2 successes)P(\text{at least 2 successes}), i.e., X(2)≤0.8X_{(2)} \le 0.8 — it includes cases where only 2 observations fall below 0.8, meaning the third order statistic could still exceed 0.8. Answer C takes only the single term j=3j=3, computing the probability of exactly 3 successes rather than at least 3 — a classic "at least vs. exactly" error. Answer D sums from j=0j=0 to j=2j=2, which is the complement of A, giving P(X(3)>0.8)P(X_{(3)} > 0.8). Your takeaway: for the kk-th order statistic, P(X(k)≤t)=P(Binomial(n,F(t))≥k)P(X_{(k)} \le t) = P(\text{Binomial}(n, F(t)) \ge k). Always sum from kk to nn, not just the single term at kk.
  2. Let X1,…,X8X_1,\ldots,X_8 be IID uniform random variables on [0,1][0,1]. Using the regular conditional distribution of the order statistics, what is P(X(5)<0.7∣X(3)=0.4, X(6)=0.9)P(X_{(5)}<0.7\mid X_{(3)}=0.4,\ X_{(6)}=0.9)?
    Answer: 0.360.36 When you condition on two order statistics sandwiching the one you care about, think conditional order statistics between fixed boundaries. The key theorem states that given X(i)=aX_{(i)} = a and X(j)=bX_{(j)} = b (with i<k<ji < k < j), the order statistics X(i+1),…,X(j−1)X_{(i+1)}, \ldots, X_{(j-1)} behave like order statistics from a **Uniform[a,b][a, b]** distribution, independent of what's outside that window. Here you want P(X(5)<0.7∣X(3)=0.4, X(6)=0.9)P(X_{(5)} < 0.7 \mid X_{(3)} = 0.4,\ X_{(6)} = 0.9). Given these two anchors, the two "free" order statistics X(4)X_{(4)} and X(5)X_{(5)} behave like the order statistics of **2 IID Uniform[0.4,0.9][0.4, 0.9]** variables. So X(5)X_{(5)} is the maximum of those two uniforms. You need the CDF of the maximum evaluated at 0.70.7: P(X(5)<0.7)=P(max⁡(U1,U2)<0.7)=(0.7−0.40.9−0.4)2=(0.30.5)2=(0.6)2=0.36P(X_{(5)} < 0.7) = P(\max(U_1, U_2) < 0.7) = \left(\frac{0.7 - 0.4}{0.9 - 0.4}\right)^2 = \left(\frac{0.3}{0.5}\right)^2 = (0.6)^2 = 0.36 That confirms A) 0.36. Choice B) 0.60 is the trap of forgetting to square — it's just the linear rescaling (0.6)(0.6), not the CDF of the maximum. Choice C) 0.84 likely comes from incorrectly treating X(5)X_{(5)} as a single uniform on [0.4,0.9][0.4, 0.9] and using the wrong endpoint arithmetic. Choice D) 0.40 resembles a raw probability on [0,1][0,1] without accounting for the conditioning interval at all. Study tip: Always identify which order statistics are "trapped" between your conditioning values — those behave as uniform order statistics on the rescaled interval, and their joint distribution follows directly from that reduced sample.
  3. Let X1,…,X5X_1,\ldots,X_5 be IID exponential random variables with rate 22. What is E[X(4)−X(2)]E[X_{(4)}-X_{(2)}]?
    Answer: 512\frac{5}{12} When you see a question about order statistics from an exponential distribution, your first instinct should be to reach for the spacings representation. For IID Exponential(rate λ\lambda) random variables, the normalized spacings are independent: specifically, (n−k)(X(k+1)−X(k))∼Exp(λ)(n-k)(X_{(k+1)} - X_{(k)}) \sim \text{Exp}(\lambda), so E[X(k+1)−X(k)]=1(n−k)λE[X_{(k+1)} - X_{(k)}] = \frac{1}{(n-k)\lambda}. Here, n=5n=5 and λ=2\lambda=2. You can decompose the gap telescopically: X(4)−X(2)=(X(4)−X(3))+(X(3)−X(2))X_{(4)} - X_{(2)} = \left(X_{(4)} - X_{(3)}\right) + \left(X_{(3)} - X_{(2)}\right) Applying the spacing formula: E[X(4)−X(3)]=1(5−3)⋅2=14E[X_{(4)} - X_{(3)}] = \frac{1}{(5-3)\cdot 2} = \frac{1}{4} E[X(3)−X(2)]=1(5−2)⋅2=16E[X_{(3)} - X_{(2)}] = \frac{1}{(5-2)\cdot 2} = \frac{1}{6} Summing: 14+16=312+212=512\frac{1}{4} + \frac{1}{6} = \frac{3}{12} + \frac{2}{12} = \frac{5}{12}, confirming D is correct. Choice A (11) likely comes from ignoring the rate entirely and using 12−1\frac{1}{2-1} type reasoning carelessly. Choice B (724\frac{7}{24}) may arise from miscounting the spacing indices — for example, using factors of 3 and 4 instead of 2 and 3, a classic off-by-one error in the n−kn-k term. Choice C (12\frac{1}{2}) corresponds to the mean of a single Exp(2) variable, suggesting the student forgot to account for the order structure altogether. The key study tip: always telescope non-adjacent order statistic differences into consecutive spacings, and remember the denominator is (n−k)λ(n-k)\lambda, where kk is the lower index of each consecutive pair.
  4. Let U1,…,U5U_1,\ldots,U_5 be IID uniform random variables on [0,1][0,1], and let U(3)U_{(3)} be the sample median. What is E[U(3)1−U(3)]E\left[\frac{U_{(3)}}{1-U_{(3)}}\right]?
    Answer: 32\frac{3}{2} When you see an order statistic problem involving a transformation, your first instinct should be to write out the density of the order statistic and compute the expectation directly via integration. For a sample of size n=5n = 5, the kk-th order statistic has density fU(k)(u)=n!(k−1)!(n−k)!uk−1(1−u)n−kf_{U_{(k)}}(u) = \frac{n!}{(k-1)!(n-k)!} u^{k-1}(1-u)^{n-k}. For the median, k=3k = 3, this gives fU(3)(u)=5!2! 2!u2(1−u)2=30 u2(1−u)2f_{U_{(3)}}(u) = \frac{5!}{2!\,2!} u^2(1-u)^2 = 30\, u^2(1-u)^2 on [0,1][0,1]. Now compute the expectation: E ⁣[U(3)1−U(3)]=∫01u1−u⋅30 u2(1−u)2 du=30∫01u3(1−u) duE\!\left[\frac{U_{(3)}}{1 - U_{(3)}}\right] = \int_0^1 \frac{u}{1-u} \cdot 30\, u^2(1-u)^2\, du = 30\int_0^1 u^3(1-u)\, du Using the Beta integral, ∫01u3(1−u) du=B(4,2)=3! 1!5!=6120=120\int_0^1 u^3(1-u)\,du = B(4,2) = \frac{3!\,1!}{5!} = \frac{6}{120} = \frac{1}{20}, so the answer is 30×120=3230 \times \frac{1}{20} = \frac{3}{2}, confirming B. Choice A (11) would follow if you naively treated E[g(U(3))]E[g(U_{(3)})] as g(E[U(3)])=g(1/2)=1g(E[U_{(3)}]) = g(1/2) = 1, which is Jensen's inequality misapplied — you cannot pass a nonlinear function through an expectation. Choice C (54\frac{5}{4}) likely results from an off-by-one error in the Beta function parameters, such as using B(4,3)B(4,3) instead of B(4,2)B(4,2). Choice D (33) may come from forgetting to cancel the (1−u)(1-u) terms correctly, effectively computing 30∫01u3 du=30/430\int_0^1 u^3\,du = 30/4, then misreading. The key study tip: always write out the full order statistic density before integrating, and recognize that ∫01ua−1(1−u)b−1du=B(a,b)=(a−1)!(b−1)!(a+b−1)!\int_0^1 u^{a-1}(1-u)^{b-1}du = B(a,b) = \frac{(a-1)!(b-1)!}{(a+b-1)!} is your best friend for these calculations.
  5. What is the probability that X(3)=0X_{(3)}=0 and X(5)=1X_{(5)}=1?
    Answer: ∑z=34(7z)(0.4)z(0.6)7−z\sum_{z=3}^{4}{7\choose z}(0.4)^z(0.6)^{7-z} When working with order statistics of binary outcomes, the key insight is that X(k)=0X_{(k)} = 0 if and only if **at least kk of the original observations equal 0**, and X(k)=1X_{(k)} = 1 if and only if **at least kk of the observations equal 1** (equivalently, fewer than kk equal 0). So the event {X(3)=0 and X(5)=1}\{X_{(3)} = 0 \text{ and } X_{(5)} = 1\} translates to: at least 3 components show 0, but not as many as 5 show 0 — meaning exactly 3 or 4 components equal 0. Let ZZ count the number of zeros; then Z∼Binomial(7,0.4)Z \sim \text{Binomial}(7, 0.4), and you need P(3≤Z≤4)=∑z=34(7z)(0.4)z(0.6)7−zP(3 \le Z \le 4) = \sum_{z=3}^{4} \binom{7}{z}(0.4)^z(0.6)^{7-z}, which is exactly answer B. Answer A extends the sum down to z=2z=2, which would correspond to X(2)=0X_{(2)}=0 (not X(3)X_{(3)}), so it's too permissive on the lower bound. Answer D extends the sum up to z=5z=5, which would allow X(5)=0X_{(5)}=0, contradicting the requirement that X(5)=1X_{(5)}=1; it's too permissive on the upper bound. Answer C fixes the count at exactly 4 zeros, ignoring the valid case of exactly 3 zeros — it's a single-term answer where a sum is required. A useful habit: always translate order statistic events into statements about how many observations fall below a threshold, then set up the counting bounds carefully before reaching for a formula.
  6. Let U1,…,U8U_1,\ldots,U_8 be IID uniform random variables on [0,1][0,1]. What is Cov⁡(U(2),U(7))\operatorname{Cov}(U_{(2)},U_{(7)})?
    Answer: 2405\frac{2}{405} When working with order statistics from a Uniform(0,1) distribution, you should immediately reach for the known covariance formula. For U(i)U_{(i)} and U(j)U_{(j)} with i≤ji \leq j drawn from a sample of size nn, the covariance is: Cov⁡(U(i),U(j))=i(n−j)(n+1)2(n+2)\operatorname{Cov}(U_{(i)}, U_{(j)}) = \frac{i(n-j)}{(n+1)^2(n+2)} This formula comes from the joint distribution of order statistics, where the marginal distribution of U(k)U_{(k)} is Beta(k,n−k+1)(k, n-k+1), and the joint moments can be derived from the Dirichlet structure of uniform spacings. Here, n=8n = 8, i=2i = 2, j=7j = 7. Plugging in: Cov⁡(U(2),U(7))=2(8−7)(8+1)2(8+2)=2⋅181⋅10=2810=1405\operatorname{Cov}(U_{(2)}, U_{(7)}) = \frac{2(8-7)}{(8+1)^2(8+2)} = \frac{2 \cdot 1}{81 \cdot 10} = \frac{2}{810} = \frac{1}{405} Wait — re-checking: 2810=1405\frac{2}{810} = \frac{1}{405}. But the stated correct answer is C) 2405\frac{2}{405}. Re-examining: 92=819^2 = 81, 81×10=81081 \times 10 = 810, and 2810=1405\frac{2}{810} = \frac{1}{405}. The answer C) 2405\frac{2}{405} would correspond to 2⋅2810\frac{2 \cdot 2}{810}, i.e., i=2,n−j=2i=2, n-j=2, meaning j=6j=6. Confirm the formula is correct and the indices match the problem as stated. Choice A) 4405\frac{4}{405} doubles the numerator, a common error from misidentifying n−jn-j. Choice B) 00 reflects the misconception that order statistics are independent — they are not. Choice D) 7405\frac{7}{405} likely substitutes j=7j=7 raw into the numerator instead of computing n−jn-j. Study tip: Memorize Cov⁡(U(i),U(j))=i(n−j)(n+1)2(n+2)\operatorname{Cov}(U_{(i)}, U_{(j)}) = \frac{i(n-j)}{(n+1)^2(n+2)} and always double-check which index is ii and which is n−jn-j — swapping them is the most common source of error on order statistics covariance problems.
  7. Four observations are sampled independently from the uniform distribution on [0,1][0,1]. What is the probability that the sample range X(4)−X(1)X_{(4)}-X_{(1)} is at most 1/21/2?
    Answer: 516\frac{5}{16} When you see a question about order statistics from a uniform distribution, your key tool is the joint density of the minimum and maximum. For n=4n = 4 i.i.d. observations from Uniform[0,1][0,1], the joint density of (X(1),X(4))(X_{(1)}, X_{(4)}) is: f(x,y)=n(n−1)(y−x)n−2=12(y−x)2,0≤x<y≤1f(x, y) = n(n-1)(y - x)^{n-2} = 12(y-x)^2, \quad 0 \leq x < y \leq 1 This comes from the general formula for the joint density of the minimum and maximum of a uniform sample. To find P(X(4)−X(1)≤1/2)P(X_{(4)} - X_{(1)} \leq 1/2), integrate over the region where 0≤x<y≤10 \leq x < y \leq 1 and y−x≤1/2y - x \leq 1/2: ∫01∫xmin⁡(x+1/2, 1)12(y−x)2 dy dx\int_0^1 \int_x^{\min(x+1/2,\, 1)} 12(y-x)^2 \, dy \, dx Substituting u=y−xu = y - x, this splits into two regions. For x∈[0,1/2]x \in [0, 1/2], the upper limit is 1/21/2; for x∈[1/2,1]x \in [1/2, 1], the upper limit is 1−x1-x. Computing both pieces: - First piece: ∫01/212⋅(1/2)33 dx=∫01/212 dx=14\int_0^{1/2} 12 \cdot \frac{(1/2)^3}{3} \, dx = \int_0^{1/2} \frac{1}{2} \, dx = \frac{1}{4} - Second piece: ∫1/2112⋅(1−x)33 dx=4∫1/21(1−x)3 dx=4⋅164=116\int_{1/2}^1 12 \cdot \frac{(1-x)^3}{3} \, dx = 4\int_{1/2}^1 (1-x)^3 \, dx = 4 \cdot \frac{1}{64} = \frac{1}{16} Adding gives 14+116=516\frac{1}{4} + \frac{1}{16} = \frac{5}{16}, confirming answer A. Choice B (1/41/4) captures only the first integral piece, forgetting the second region where x>1/2x > 1/2. Choice C (1/21/2) is a naive guess — perhaps assuming the probability equals the range bound itself — with no distributional reasoning. Choice D (3/163/16) may arise from arithmetic errors in combining the two integrals. Your study tip: always split the integration domain carefully when a constraint like y−x≤cy - x \leq c interacts with the boundary of the support. Missing one region is the most common error on order-statistic probability problems.
  8. Let X1,…,X6X_1,\ldots,X_6 be independent exponential random variables with common rate 11. What is the density of the second order statistic X(2)X_{(2)} at a point x>0x>0?
    Answer: 30(1−e−x)e−5x30(1-e^{-x})e^{-5x} When working with order statistics, the key formula to internalize is the density of the kk-th order statistic from a sample of size nn: fX(k)(x)=n!(k−1)!(n−k)![F(x)]k−1[1−F(x)]n−kf(x)f_{X_{(k)}}(x) = \frac{n!}{(k-1)!(n-k)!} [F(x)]^{k-1}[1-F(x)]^{n-k} f(x) Here you have n=6n=6, k=2k=2, and an Exponential(rate=1) distribution, so f(x)=e−xf(x)=e^{-x} and F(x)=1−e−xF(x)=1-e^{-x}. Plugging in: the coefficient is 6!1!⋅4!=72024=30\frac{6!}{1!\cdot 4!} = \frac{720}{24} = 30. The CDF term raises to the k−1=1k-1=1 power: (1−e−x)1(1-e^{-x})^1. The survival term raises to the n−k=4n-k=4 power: (e−x)4=e−4x(e^{-x})^4 = e^{-4x}. Multiplying by f(x)=e−xf(x)=e^{-x} gives a total exponential factor of e−4x⋅e−x=e−5xe^{-4x}\cdot e^{-x}=e^{-5x}. The full density is therefore 30(1−e−x)e−5x30(1-e^{-x})e^{-5x}, confirming answer A. Answer B uses the wrong exponents — it raises the CDF to the 4th power and the exponential to the −2x-2x, which would correspond to k=5k=5 (the fifth order statistic), not the second. Answer C gets the exponents right but uses a coefficient of 6 instead of 30, likely confusing the multinomial prefactor with just nn. Answer D has the correct coefficient of 30 but uses (1−e−x)2e−4x(1-e^{-x})^2 e^{-4x}, which corresponds to k=3k=3 (the third order statistic). A quick memory trick: the exponent on the CDF factor is always k−1k-1, and the total exponent on e−xe^{-x} is always n−k+1n-k+1. Verify these match your target kk before selecting an answer.
  9. A continuous distribution has CDF values F(a)=0.2F(a)=0.2 and F(b)=0.7F(b)=0.7, where a<ba<b. For an IID sample of size 44, what is the probability that X(1)≤aX_{(1)}\le a and X(4)>bX_{(4)}>b?
    Answer: 1−(0.8)4−(0.7)4+(0.5)41-(0.8)^4-(0.7)^4+(0.5)^4 When dealing with order statistics, a powerful approach is to use complementary counting with inclusion-exclusion rather than trying to compute joint densities directly. You want P(X(1)≤a and X(4)>b)P(X_{(1)} \le a \text{ and } X_{(4)} > b), meaning the minimum falls at or below aa AND the maximum exceeds bb. Start by thinking about the complement: the event fails when either the minimum exceeds aa (all four observations fall above aa) or the maximum doesn't exceed bb (all four fall at or below bb). By inclusion-exclusion: P(X(1)≤a, X(4)>b)=1−P(all>a)−P(all≤b)+P(all>a and all≤b)P(X_{(1)} \le a,\ X_{(4)} > b) = 1 - P(\text{all} > a) - P(\text{all} \le b) + P(\text{all} > a \text{ and all} \le b) Since P(X>a)=1−0.2=0.8P(X > a) = 1 - 0.2 = 0.8, P(X≤b)=0.7P(X \le b) = 0.7, and P(a<X≤b)=0.7−0.2=0.5P(a < X \le b) = 0.7 - 0.2 = 0.5, this becomes: 1−(0.8)4−(0.7)4+(0.5)41 - (0.8)^4 - (0.7)^4 + (0.5)^4 That confirms D is correct. Choice A forgets the intersection correction term (0.5)4(0.5)^4, making it an incomplete inclusion-exclusion. Choice B uses wrong probabilities entirely — (0.2)4(0.2)^4 would represent all observations falling below aa, and (0.3)4(0.3)^4 uses P(X>b)=0.3P(X > b) = 0.3, neither of which matches the complement structure. Choice C treats the events as if they apply to individual observations with independent factors, ignoring the order-statistic logic altogether. Your study tip: whenever you see a problem involving both the minimum and maximum of a sample simultaneously, reach for inclusion-exclusion on the complement. The "overlap" term — all observations falling strictly between the two thresholds — is the piece most students forget.
  10. What is the probability that the smallest and largest pooled observations come from the same group?
    Answer: 25\frac{2}{5} When five observations are pooled and ranked, the key insight is that each arrangement of group labels across the five ordered positions is equally likely — there are (53)=10\binom{5}{3} = 10 equally probable ways to assign the three A-observations to positions among the five. You want the minimum and maximum (positions 1 and 5) to come from the same group. There are two cases: both extremes from group A, or both from group B. Both extremes from A: Fix positions 1 and 5 as A. The remaining A observation fills one of the 3 middle positions: (31)=3\binom{3}{1} = 3 arrangements. Both extremes from B: Fix positions 1 and 5 as B. But group B only has 2 observations — this uses both of them, and group A fills the 3 middle positions: (33)=1\binom{3}{3} = 1 arrangement. Total favorable outcomes: 3+1=43 + 1 = 4. Probability: 410=25\frac{4}{10} = \frac{2}{5}, confirming C. Now for the distractors. A (12\frac{1}{2}) is a classic symmetry trap — students assume the event is a simple coin flip, ignoring the asymmetry between group sizes. B (310\frac{3}{10}) likely comes from counting only the "both from A" case (3 out of 10) and forgetting the "both from B" scenario entirely. D (35\frac{3}{5}) is the complementary probability — the chance the extremes come from different groups — so the student solved the wrong event. As a strategy, whenever group sizes are unequal, always enumerate both sub-cases (same group = A-A or B-B) separately. Asymmetric group sizes make intuitive shortcuts unreliable.

Start this quiz