I try to verify the following result of extreme value statistics that the distribution of the block maxima of Gaussian random variables converges to a Gumbel distribution. To this end, I first generate Gaussian random varaibles, divide them in blocks of length 100, compute the maximum value of each block and fit a Gumbel distribution
data=RandomVariate[NormalDistribution[0, 1], 100000];
blocks=Table[data[[i + 1 ;; i + 100]], {i, 0, 99900}];
max=Table[blocks[[k]]//Max, {k, 1, 99901}];
est = EstimatedDistribution[max,GumbelDistribution[a, b]];
Show[Histogram[max,{.1},"PDF",PlotRange->{{0, 8},All},Frame->True, Axes->False],Plot[PDF[est, x], {x, 0, 6}]]
As you see there is a poor agreement. Am I doing something wrong?
ExtremeValueDistribution
for maxima,GumbelDistribution
is for minima. This is explained in the documentation forGumbelDistribution
$\endgroup$