3
$\begingroup$

Consider a random vector $X^{(m)} = (X^{(m)}_1,\dots,X^{(m)}_m)$ where, for fixed $m$, the elements of $X^{(m)}$ are i.i.d. $\mathcal{N}(0,\sigma^2 / m)$. Define $$Z_m =\max_{k=1,\dots,m}X^{(m)}_k.$$ What can we say about the distribution of $Z_m$ as $m$ gets large? Ideally I'd like a result such of the form $Z_m=\mathcal{O}_p(f(Z))$. Note that this problem is different than the standard problem where the variance is fixed, as is discussed here and here.

$\endgroup$
10
  • $\begingroup$ Your description of the problem is ambiguous. Are you trying to say that the variance of $X_i$ is $\sigma^2/i$ for all $i$? In that case "iid" makes no sense. Or are you asking about the maxima of $m$ iid vectors, all of which have variance $\sigma^2/m$? $\endgroup$
    – whuber
    Commented Feb 14, 2020 at 20:06
  • $\begingroup$ @Aksakal It says "i.i.d." $\endgroup$
    – Henry
    Commented Feb 14, 2020 at 20:42
  • $\begingroup$ @whuber I'm asking about the later. For a given value of $m$, the elements of the random $m$-vector are i.i.d. with a Gaussian distribution with variance $\sigma^2 /m$ $\endgroup$
    – nothing
    Commented Feb 14, 2020 at 20:48
  • $\begingroup$ Hopefully the new notation clarifies this... $\endgroup$
    – nothing
    Commented Feb 14, 2020 at 20:54
  • 1
    $\begingroup$ Thank you; the new notation helps. It still leaves the interpretation open, though, because $Z_m$ converges to zero in probability and if you wish to standardize it to get some nontrivial answer, you're right back in the classical case of the extreme-value distribution associated with the standard Normal. The meaning of "fundamentally different" therefore remains obscure. $\endgroup$
    – whuber
    Commented Feb 14, 2020 at 20:57

1 Answer 1

4
$\begingroup$

The lesson of extreme value theory is that when you rescale the maximum of $m$ iid standard Normal variables by an amount proportional to $\sqrt{2\log m}$ and shift that to a location near $\sqrt{2\log m},$ it converges to the standard Gumbel distribution as $m$ increases.

If first you rescale your $m$ iid Normal variables by $\sqrt{m},$ they become standard Normal, and the preceding applies.

Consequently, for suitable constants $\alpha$ and $\beta,$ $$Z_m \beta \sqrt{2m\log m} - \alpha \sqrt{2\log m}$$ converges to a Gumbel distribution.

In particular, this means we may approximate $Z_m$ closely (in distribution) by

$$Z_m \approx \frac{Y }{\beta \sqrt{2m\log m}} + \frac{\alpha}{\beta\sqrt{m}}$$

for a Gumbel variate $Y.$ Since the denominators diverge, the (unnormalized) limiting distribution of $Z_m$ is $0$ (in probability). Notice that it does so by being squeezed narrowly around a positive value $\alpha / (\beta\sqrt{m})$ while this central value creeps slowly down towards zero. Here's a histogram with $m=10000$ based on $4000$ such samples:

Figure

The red curve plots the limiting Gumbel density, $f(y) = \exp(-e^{-y} - y).$

$\endgroup$
2
  • $\begingroup$ The link "lesson ..." seems wrong. I guess it should define $a$ and $b$. $\endgroup$
    – Yves
    Commented Feb 16, 2020 at 16:45
  • $\begingroup$ @Yves Oops, you're right: the last digit "5" of the thread id was cut off during the pasting. I'll fix that right away. $\endgroup$
    – whuber
    Commented Feb 16, 2020 at 18:37

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.