I am currently working on a not-that-easy problem involving order statistics. As I am unsure as to how I could solve it, I thought it might already possess a solution. So here I am, my questions is: do you know of a pre-existing solution/approximation to (or do you see an easy way of solving) the probability below? $$\pi = \Pr\left[G > c\right],$$ where $$ G_1 = \max_{1 \leq j \leq K} \left(Z_{(j)} - Z_{(j - 1)}\right), Z_{(0)} = 0, Z_j \sim N(\mu, \sigma^2) \text{ and } Z_i \perp Z_j, \forall i \ne j. $$ Just to clarify things up, $Z_{(j)}$ is the $j$th order statistics of $Z_1, Z_2, \ldots, Z_K$. Also, if it helps, the problem can actually be reformulated this way (but I would prefer the one above):

$$G_2 = \max_{1 \leq j \leq K} \left(\tfrac{Z_{(j)}}{Z_{(j - 1)}}\right), Z_{(0)} = 1, Z_j \sim N(\mu, \sigma^2) \text{ and } Z_i \perp Z_j, \forall i \ne j.$$

I understand that mathematically speaking, those expressions are not equivalent, but for $\mu \in (0, 1)$ and $\sigma$ small enough, it should not matter.

Sorry if the exact same problem was posted before, sorry for my english.


ETA: OK, worked quite a bit on this damned problem, and found that for $\mu \in (0, 1)$, e.g. $\mu = 0.1$, and $\sigma$ small enough, e.g. $\sigma = 0.01$, $G$ actually simplifies to

$$G = \max_{1 \leq j \leq n} \left(Z_{(j)} - Z_{(j - 1)}\right) \approx Z_{(1)}.$$

Thus, the probability $\pi$ is easily computed

$$\pi = \Pr\left[G > c\right] \approx \Pr\left[Z_{(1)} > c\right] = {\left(1 - \Phi\left(\tfrac{c - \mu}{\sigma}\right)\right)}^K$$

The critical value $c$ is also easily extracted

$${\left(1 - \Phi\left(\tfrac{c - \mu}{\sigma}\right)\right)}^K = \alpha \quad \Leftrightarrow \quad c = \mu + \Phi^{-1} \left(1 - \alpha^{\tfrac{1}{K}}\right)\sigma.$$

All I got to do now (other than finding a more general solution to the above problem) is find values of $\mu$ and $\sigma$, which are context-dependent.

If anyone feels I should pile on and add some more details, please let me know. I am, after all, new to this forum.


closed as unclear what you're asking by kjetil b halvorsen, mdewey, whuber Dec 14 '18 at 16:29

Please clarify your specific problem or add additional details to highlight exactly what you need. As it's currently written, it’s hard to tell exactly what you're asking. See the How to Ask page for help clarifying this question. If this question can be reworded to fit the rules in the help center, please edit the question.

  • 3
    $\begingroup$ There is something strange about these formulations of your problem, because clearly they are not equivalent. For instance, about half the $Z_{(j)}$ will be negative and in that range the ratios of order statistics will be less than unity, so the second formulation effectively ignores the first half of the gaps on average! Are you perhaps assuming in the second formulation that the $Z_j$ are lognormal? Also, could you explain why you include a constant $Z_{(0)}$ in the list of order statistics? $\endgroup$ – whuber Jan 24 '14 at 20:11
  • 1
    $\begingroup$ Well, technically you are right, but here it is, the variables $Z_j, j = 1, 2, \ldots, K$ are in truth proportions (whose support $\Omega$ is $(0, 1)$, which renders both formulation equivalent) that are asymptotically distributed as normal distributions (as stated by the central limit theorem). Also, I included the variable $Z_{(0)}$ because it was easier this way. Otherwise, I would have had to state $G$ as: $$G = \max \left(Z_{(1)}, \max_{2 \leq j \leq K} \left(Z_{(j)} - Z_{(j - 1)} \right)\right).$$ Hope this clears things up. $\endgroup$ – RSMax Jan 25 '14 at 2:12
  • 3
    $\begingroup$ Largest/maximum spacing statistics are used in different statistical fields. An asymptotic result is e.g. jstor.org/discover/10.2307/… $\endgroup$ – Michael M Jan 25 '14 at 10:15
  • 2
    $\begingroup$ After your edit it's even less clear what you are doing. Your comments indicate that some kind of asymptotic result is needed, but the nature of the asymptotics is not in evidence. The question itself still is internally contradictory, because the distributions of $G_1$ and $G_2$ have no direct relationship to each other. $\endgroup$ – whuber Jan 25 '14 at 19:30
  • 2
    $\begingroup$ If they are reasonable proportion estimators, then they definitely are not asymptotically Normal: asymptotically, they converge to constants (the correct underlying proportion). You can obtain asymptotic Normality only by standardizing them. That process will push $Z_{(0)}$ out to $-\infty$, leading to a trivial (and useless) solution. $\endgroup$ – whuber Jun 24 '16 at 1:35