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.
-
$\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$– HenryCommented 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$– nothingCommented Feb 14, 2020 at 20:48
-
$\begingroup$ Hopefully the new notation clarifies this... $\endgroup$– nothingCommented 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
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:
The red curve plots the limiting Gumbel density, $f(y) = \exp(-e^{-y} - y).$
-
$\begingroup$ The link "lesson ..." seems wrong. I guess it should define $a$ and $b$. $\endgroup$– YvesCommented 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