An Upper Bound for the Average of Top Order Statistics - MathOverflow most recent 30 from http://mathoverflow.net 2013-05-19T06:44:52Z http://mathoverflow.net/feeds/question/42922 http://www.creativecommons.org/licenses/by-nc/2.5/rdf http://mathoverflow.net/questions/42922/an-upper-bound-for-the-average-of-top-order-statistics An Upper Bound for the Average of Top Order Statistics Santiago 2010-10-20T18:40:54Z 2010-10-22T11:55:55Z <p>The following problem arises when we try to bound the expected offline optimal value of a simple online assignment problem with random values and unit weights, by its deterministic approximation.</p> <h2>The Problem</h2> <p>Consider a sequence <code>$\{X_i\}_{i=1}^n$</code> of non-negative integrable i.i.d. random variables with absolutely continuous c.d.f. <code>$F(x)$</code>. Let <code>$X_{(i)}$</code> be the <code>$i^{\rm th}$</code> order statistics, so that <code>$X_{(1)}$</code> is the minimum of the sequence and <code>$X_{(n)}$</code> is the maximum. Now, let <code>$T_k$</code> be the average of the top <code>$k^{\rm th}$</code> order statistics, that is, <code>$T_k = \frac 1 k \sum_{i=n-k+1}^n X_{(i)}$</code>. We would like to show that the expected value of the average of the top order statistics is upper bounded as follows</p> <p><code>$$\mathbb{E} T_{k} \le \mathbb{E} \left[ X | F(X) \ge 1 - k/n \right],$$</code></p> <p>where the right hand size is the conditional expectation of X given that it is larger than the <code>$(1-k/n)$</code>-percentile. In order to keep things simple, we may assume that the c.d.f. <code>$F(\cdot)$</code> is strictly increasing in its domain.</p> <p>Moreover, we may we fix <code>$\rho \in (0,1)$</code>, and set <code>$k=[\rho n]$</code>; so that we are interested in the average of the top <code>$\rho$</code> fraction of the sequence. If we scale both <code>$n$</code> and <code>$k$</code> to infinity while keeping the ratio <code>$\rho$</code> fixed, it seems to be the case that the bound is asymptotically tight:</p> <p><code>$$\lim_{n \rightarrow \infty} \mathbb{E} T_{[\rho n]} = \mathbb{E} \left[ X | F(X) \ge 1 - \rho \right].$$</code></p> <p>Can you show if these results hold? I have done some numerical experiments with a couple of distributions (uniform, truncated normal, exponential) that confirm these results. Any help or pointer would be appreciated.</p> <p>Thanks in advance.</p> <h2>More motivation</h2> <p>This is a stripped down version of a more complicated problem. Suppose that we have an incoming inventory of <code>$n$</code> different items with unknown value arriving in an online fashion. We can only keep <code>$k$</code> items of the <code>$n$</code> total. The decision to keep an item has to be made at the moment of arrival, and once we decide to keep an item we need to stick to this decision.</p> <p>The value of i-th item, denoted by <code>$X_i$</code>, is unknown, and revealed when it arrives (before a decision needs to be made). However, we do have a prior for the values; they are drawn independently from an identical distribution <code>$F(\cdot)$</code>. The objective of the problem is designing an online policy that maximizes the expected total value of the assignment.</p> <p>A useful benchmark, when comparing online policies, is the offline optimal solution. Given a realization of the values <code>$X=\{X_i\}_{i=1}^n$</code>, the optimal value of the offline problem, denoted by <code>$P(X)$</code>, is<br> <code>\begin{align} P(X) = \max_y &amp;\; \sum_{i=1}^n y_i X_i \\ \text{s.t.} &amp; \sum_{i=1}^n y_i = k, \\ &amp; y_i \in \{0,1\} \end{align}</code> In this simple case, the offline optimal solution is to keep the items with the k-th highest values. Equivalently, we have that <code>$P(X) = \sum_{i=n-k+1}^n X_{(i)}$</code>. We are interested in the expected optimal value of offline optimal solution, which is given by <code>$\mathbb{E} P(X) = k \mathbb{E} T_k$</code>.</p> <p>A simple online policy could instruct us to keep those items with value larger than the <code>$(1-k/n)$</code>-percentile. Such policy can be shown to attain an expected value of <code>$n \mathbb{E} \left[ X \mathbf{1} \{ F(X) \ge 1 - k/n\} \right]$</code> when <code>$n$</code> is large.</p> <p>The results that we want to prove would allows us to upper bound the expected optimal value of the offline problem by a bound that could be attained, asymptotically, by an online policy. This would confirm that our policy is good.</p> http://mathoverflow.net/questions/42922/an-upper-bound-for-the-average-of-top-order-statistics/43140#43140 Answer by jat for An Upper Bound for the Average of Top Order Statistics jat 2010-10-22T07:31:08Z 2010-10-22T07:31:08Z <p>There's a variational principle, attributed to Ky Fan, that provides an alternative description of the sum of the $k$ largest numbers in a nonnegative sequence: $$\frac{1}{k} \sum_{i=1}^k x^\circ_i = \inf_{x = y+ z} \ \Vert y \Vert_\infty + \frac{1}{k} \Vert z \Vert_1$$ The infimum is attained by a 1-parameter family (sometimes I get this backward, so check it!): $$y_i = \begin{cases} x_i, &amp; x_i \leq t \newline t, &amp; \text{else} \end{cases} \quad\text{and}\quad z_i = \begin{cases} x_i - t, &amp; x_i > t \newline 0, &amp; \text{else} \end{cases}$$ Then the infimum occurs over $t > 0$. The advantage of this formulation for your problem is that you can draw the expectation inside the infimum, compute the expectation (which is easy), and then take the infimum. You should get the value you noted above.</p> http://mathoverflow.net/questions/42922/an-upper-bound-for-the-average-of-top-order-statistics/43159#43159 Answer by zhoraster for An Upper Bound for the Average of Top Order Statistics zhoraster 2010-10-22T10:23:43Z 2010-10-22T11:55:55Z <p>Write $$E[T_k]=E[E[T_k|X_{(n-k)}]]$$ The distribution of $X_{(n-k+1)},\dots,X_{(n)}$ given $X_{(n-k)}=x$ is the same as the conditional distribution of a monotone arrangement of $n-k$ independent r.v.'s with cdf $F$ given they all are not less than $x$. (This can be proved easily e.g. by using the quantile transform.)</p> <p>Now the mean value of monotone arrangement of a sequence is equal to the mean value of the sequence itself. So by linearity of conditional expectation we can write $$E[E[T_k|X_{(n-k)}]]=E[E[X|X\ge X_{(n-k)}]],$$ where $X$ is an independent of $X_{(n-k)}$ r.v. with cdf $F$.</p> <p>Now it is again about the quantile transform. Let $U=F(X)$, $U_{(n-k)}=F(X_{(n-k)})$. We have (by the independence of $U$ and $U_{(n-k)}$) $$E[X|X\ge X_{(n-k)}]=E[F^{(-1)}(U)|U>U_{(n-k)}] = \frac{1}{1-U_{(n-k)}}\phi(U_{(n-k)}),$$ where $\phi(u) = \int_u^1 F^{(-1)}(x)dx$. But we do know the density of $U_{(n-k)}$, so the expectation of this expression is equal to $$\frac{n!}{(n-k-1)k!}\int_0^1 u^{n-k-1}(1-u)^{k-1}\phi(u)du = \frac{n}{k}E[\phi(B)],$$ where $B$ has $\mathrm{Beta}(n-k,k)$ distribution. But $\phi(x)$ is easily seen to be concave, so by Jensen's inequality $E[\phi(B)]\le \phi(E[B]) = \phi(1-k/n)$, as required.</p>