numpy.random.RandomState.gumbel

метод

RandomState.gumbel(loc=0.0, scale=1.0, size=None)

Генерирует выборку из распределения Гумбеля.

Генерирует выборку из распределения Гумбеля с заданными параметрами сдвига и масштаба. Более подробную информацию о распределении Гумбеля см. в примечаниях и ссылках ниже.

Примечание

Новый код должен использовать метод gumbel экземпляра default_rng(); см. random-quick-start.

Параметры

locfloat or array_like of floats, optional

Сдвиг, определяющий положение моды распределения. По умолчанию 0.

scalefloat or array_like of floats, optional

Масштабный параметр распределения. По умолчанию 1. Должен быть неотрицательным.

sizeint or tuple of ints, optional

Форма выходных данных. Если заданная форма, например, (m, n, k), то генерируется m * n * k выборок. Если размерность равна None (по умолчанию), возвращается одно значение, если loc и scale оба являются скалярными величинами. В противном случае генерируется np.broadcast(loc, scale).size выборок.

Возвращаемое значение

outndarray or scalar

Сгенерированная выборка из распределения Гумбеля с заданными параметрами.

См. также

scipy.stats.gumbel_l, scipy.stats.gumbel_r, scipy.stats.genextreme, weibull

Generator.gumbel

который следует использовать в новом коде.

Примечания

Распределение Гумбеля (или распределение наименьших экстремальных значений (SEV) или распределение наименьших экстремальных значений типа I) является одним из класса распределений обобщенных экстремальных значений (GEV), используемых в задачах моделирования экстремальных значений. Распределение Гумбеля является частным случаем распределения экстремальных значений типа I для максимальных значений из распределений с "экспоненциально подобными" хвостами.

Плотность вероятности для распределения Гумбеля:

<img alt="p(x) = \frac{e^{-(x - \mu)/ \beta}}{\beta} e^{ -e^{-(x - \mu)/ \beta}}," src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0nVVRGLTgnPz4KPCEtLSBUaGlzIGZpbGUgd2FzIGdlbmVyYXRlZCBieSBkdmlzdmdtIDIuNi4zIC0tPgo8c3ZnIHZlcnNpb249JzEuMScgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJyB4bWxuczp4bGluaz0naHR0cDovL3d3dy53My5vcmcvMTk5OS94bGluaycgd2lkdGg9JzEzOS4xMTk0MDRwdCcgaGVpZ2h0PScyOC45MjkwNDlwdCcgdmlld0JveD0nMTYzLjU2NjA3MSA3OC43MDQ4NTcgMTM5LjExOTQwNCAyOC45MjkwNDknPgo8ZGVmcz4KPHBhdGggaWQ9J2c1LTQwJyBkPSdNMi4zNzkwNzggMS4zNTA5MzRDMS43MjE1NDQgMC44NDg4MTcgMS4xODM1NjIgLTAuMDk1NjQxIDEuMTgzNTYyIC0xLjQ5NDM5NkMxLjE4MzU2MiAtMi44NzUyMTggMS43MDk1ODkgLTMuODMxNjMxIDIuMzc5MDc4IC00LjMzOTcyNkMyLjM3OTA3OCAtNC4zNDU3MDQgMi4zOTcwMTEgLTQuMzY5NjE0IDIuMzk3MDExIC00LjM5OTUwMkMyLjM5NzAxMSAtNC40MzUzNjcgMi4zNzMxMDEgLTQuNDgzMTg4IDIuMzAxMzcgLTQuNDgzMTg4QzIuMTkzNzczIC00LjQ4MzE4OCAwLjc1MzE3NiAtMy41NDQ3MDcgMC43NTMxNzYgLTEuNDk0Mzk2QzAuNzUzMTc2IDAuNTQ5OTM4IDIuMTg3Nzk2IDEuNDk0Mzk2IDIuMzAxMzcgMS40OTQzOTZDMi4zNzMxMDEgMS40OTQzOTYgMi4zOTcwMTEgMS40NDY1NzUgMi4zOTcwMTEgMS40MTA3MUMyLjM5NzAxMSAxLjM4MDgyMiAyLjM3OTA3OCAxLjM1NjkxMiAyLjM3OTA3OCAxLjM1MDkzNFonLz4KPHBhdGggaWQ9J2c1LTQxJyBkPSdNMi4xMTYwNjUgLTEuNDk0Mzk2QzIuMTE2MDY1IC0zLjUzODczIDAuNjgxNDQ1IC00LjQ4MzE4OCAwLjU2Nzg3IC00LjQ4MzE4OEMwLjQ5NjEzOSAtNC40ODMxODggMC40NzgyMDcgLTQuNDI5MzkgMC40NzgyMDcgLTQuMzk5NTAyQzAuNDc4MjA3IC00LjM2MzYzNiAwLjQ4NDE4NCAtNC4zNTc2NTkgMC41OTE3ODEgLTQuMjYyMDE3QzEuMTI5NzYzIC0zLjg5MTc0MyAxLjY4NTY3OSAtMi45MTcwNjEgMS42ODU2NzkgLTEuNDk0Mzk2QzEuNjg1Njc5IC0wLjIyNzE0OCAxLjI0MzMzNyAwLjc0MzI4OCAwLjc1MzE3NiAxLjI5MTExOFMwLjc1MzE3NiAwLjU0OTkzOCAyLjE4NzU5MiAxLjQ5NDM5NiAyLjMzMTEgMS40OTQzOTZaJy8+CjxwYXRoIGlkPSdnMC0wJyBkPSdNNC43NTgxNTcgLTEuMzM4OTc5QzQuODUzNzk4IC0xLjMzODk3OSA1LjAwMzIzOCAtMS4zMzg5NzkgNS4wMDMyMzggLTEuNDk0Mzk2UzQuODUzNzk4IC0xLjY0OTgxMyA0Ljc1ODE1NyAtMS42NDk4MTNIMC45OTIyNzlDMC44OTY2MzggLTEuNjQ5ODEzIDAuNzQ3MTk4IC0xLjY0OTgxMyAwLjc0NzE5OCAtMS40OTQzOTZTMC44OTY2MzggLTEuMzM4OTc5IDAuOTkyMjc5IC0xLjMzODk3OUg0Ljc1ODE1N1onLz4KPHBhdGggaWQ9J2cyLTEyJyBkPSdNNC4wMDQ5ODEgLTMuNDAxMjQ1QzQuMDA0OTgxIC0zLjg0MzU4NyAzLjU5MjUyOCAtNC4yMTQxOTcgMi45MjkwMTYgLTQuMjE0MTk3QzIuNzQ5Njg5IC00LjIxNDE5NyAyLjU1ODQwNiAtNC4xNjAzOTkgMi4zOTcwMTEgLTQuMDc2NzEyQzEuODQ3MDczIC0zLjc5NTc2NiAxLjQ1MjU1MyAtMy4yMDM5ODUgMS4zMjEwNDYgLTIuNjcxOThMMC4zODI1NjUgMS4wODE5NDNDMC4zNjQ2MzMgMS4xNDE3MTkgMC40NDIzNDEgMS4xNTk2NTEgMC41MDIxMTcgMS4xNTk2NTFDMC41NTU5MTUgMS4xNTk2NTEgMC41OTE3ODEgMS4xNTM2NzQgMC42MDM3MzYgMS4wOTk4NzVDMC42MDk3MTQgMS4wOTM4OTggMC45MzI1MDMgLTAuMTk3MjYgMC45OTgyNTcgLTAuNDY2MjUyQzEuMTc3NTg0IC0wLjEwMTYxOSAxLjU1NDE3MiAwLjA1OTc3NiAxLjk0MjcxNSAwLjA1OTc3NkMyLjc2NzYyMSAwLjA1OTc3NiAzLjc5NTc2NiAtMC41MDgwOTUgMy43OTU3NjYgLTEuNDM0NjJDMy43OTU3NjYgLTEuNzU3NDEgMy42NzYyMTQgLTIuMTA0MTEgMy4zMTE1ODIgLTIuMzYxMTQ2QzMuODEzNjk5IC0yLjY3Nzk1OCA0LjAwNDk4MSAtMy4wMjQ2NTggNC4wMDQ5ODEgLTMuNDAxMjQ1Wk0yLjcwNzg0NiAtMi4zNjcxMjNDMi42NDgwNyAtMi4zNTUxNjggMi40OTg2MyAtMi4zNDMyMTMgMi40Mzg4NTQgLTIuMzQzMjEzUzIuMjIzNjYxIC0yLjM0MzIxMyAyLjE4MTgxOCAtMi4zNjExNDZDMi4yNzc0NiAtMi4zNzMxMDEgMi4zNzMxMDEgLTIuMzc5MDc4IDIuNDY4NzQyIC0yLjM3OTA3OEMyLjU0MDQ3MyAtMi4zNzkwNzggMi42MTgxODIgLTIuMzczMTAxIDIuNzA3ODQ2IC0yLjM2NzEyM1pNMy41ODA1NzMgLTMuNDYxMDIxQzMuNTgwNTczIC0zLjAzNjYxMyAzLjMyMzUzNyAtMi42NjAwMjUgMy4wNDI1OSAtMi40OTg2M0MyLjg2OTI0IC0yLjU1ODQwNiAyLjY4OTkxMyAtMi41NzYzMzkgMi41MTA1ODUgLTIuNTc2MzM5QzIuMzYxMTQ2IC0yLjU3NjMzOSAxLjkwNjg0OSAtMi41ODgyOTQgMS45MDY4NDkgLTIuMzQ5MTkxQzEuOTA2ODQ5IC0yLjEyODAyIDIuMjg5NDE1IC0yLjE0NTk1MyAyLjQzODg1NCAtMi4xNDU5NTNDMi41ODIzMTYgLTIuMTQ1OTUzIDIuODk5MTI4IC0yLjE2OTg2MyAzLjAyNDY1OCAtMi4yNDE1OTRDMy4xOTIwMyAtMi4xMjIwNDIgMy4zMjk1MTQgLTEuOTQ4NjkyIDMuMzI5NTE0IC0xLjU5MDAzN0MzLjMyOTUxNCAtMC43NzcwODYgMi43Nzk1NzcgLTAuMTM3NDg0IDEuOTM2NzM3IC0wLjEzNzQ4NEMxLjUwMDM3NCAtMC4xMzc0ODQgMS4xMjk3NjMgLTAuNDE4NDMxIDEuMTI5NzYzIC0wLjg3ODcwNUMxLjEyOTc2MyAtMC45MzI1MDMgMS4xMjk3NjMgLTAuOTgwMzI0IDEuMTQxNzE5IC0xLjAyODE0NEwxLjU0ODE5NCAtMi42NDgwN0MxLjY5MTY1NiAtMy4yMjc4OTUgMi4yNjU1MDQgLTQuMDE2OTM2IDIuOTM0OTk0IC00LjAxNjkzNkMzLjI5MzY0OSAtNC4wMTY5MzYgMy41ODA1NzMgLTMuODU1NTQyIDMuNTgwNTczIC0zLjQ2MTAyMVonLz4KPHBhdGggaWQ9J2cyLTIyJyBkPSdNMi43OTE1MzIgLTAuODQyODM5QzIuNzYxNjQ0IC0wLjcyOTI2NSAyLjc2MTY0NCAtMC43MTczMSAyLjY1NDA0NyAtMC41Nzk4MjZDMi41MjI1NCAtMC40MTg0MzEgMi4yNDE1OTQgLTAuMTM3NDg0IDEuODcwOTg0IC0wLjEzNzQ4NEMxLjM0NDk1NiAtMC4xMzc0ODQgMS4zNDQ5NTYgLTAuNTQ5OTM4IDEuMzQ0OTU2IC0wLjY2OTQ4OUMxLjM0NDk1NiAtMC44MjQ5MDcgMS4zNjI4ODkgLTAuOTAyNjE1IDEuNDEwNzEgLTEuMDk5ODc1QzEuNDY0NTA4IC0xLjI5NzEzNiAxLjQ3NjQ2MyAtMS4zNDQ5NTYgMS41MzAyNjIgLTEuNTY2MTI3UzEuNTk2MDE1IC0xLjgzNTExOCAxLjY0OTgxMyAtMi4wMzIzNzlDMS42Nzk3MDEgLTIuMTY5ODYzIDEuNzM5NDc3IC0yLjQwMjk4OSAxLjczOTQ3NyAtMi40MzI4NzdDMS43Mzk0NzcgLTIuNTc2MzM5IDEuNjI1OTAzIC0yLjYzNjExNSAxLjUzMDI2MiAtMi42MzYxMTVDMS4zMDkwOTEgLTIuNjM2MTE1IDEuMjU1MjkzIC0yLjQyNjg5OSAxLjIzNzM2IC0yLjM0OTE5MUMwLjQwMDQ5OCAxLjAwNDIzNEMwLjM4ODU0MyAxLjA1MjA1NSAwLjM4MjU2NSAxLjA2OTk4OCAwLjM4MjU2NSAxLjA4NzkyQzAuMzgyNTY1IDEuMTU5NjUxIDAuNDM2MzY0IDEuMjg1MTgxIDAuNTk3NzU4IDEuMjg1MTgxQzAuODA2OTc0IDEuMjg1MTgxIDAuODYwNzcyIDEuMTA1ODUzIDAuODkwNjYgMC45ODAzMjRDMC45MjA1NDggMC44ODQ2ODIgMS4wMDQyMzQgMC41MTQwNzIgMS4xNzE2MDYgLTAuMTE5NTUyQzEuNDE2Njg3IDAuMDU5Nzc2IDEuNzgxMzIgMC4wNTk3NzYgMS44NDcwNzMgMC4wNTk3NzZDMi4zMDEzNyAwLjA1OTc3NiAyLjU3NjMzOSAtMC4yMDkyMTUgMi43NDM3MTEgLTAuMzc2NTg4QzIuODU3Mjg1IDAuMDExOTU1IDMuMjYzNzYxIDAuMDU5Nzc2IDMuMzk1MjY4IDAuMDU5Nzc2QzMuNjEwNDYxIDAuMDU5Nzc2IDMuNzY1ODc4IC0wLjA1OTc3NiAzLjg3OTQ1MiAtMC4yNDUwODFDNC4wMjg4OTIgLTAuNDg0MTg0IDQuMTEyNTc4IC0wLjgzMDg4NCA0LjExMjU3OCAtMC44NjA3NzJDNC4xMTI1NzggLTAuODcyNzI3IDQuMTA2NiAtMC45NDQ0NTggMy45OTkwMDQgLTAuOTQ0NDU4QzMuOTA5MzQgLTAuOTQ0NDU4IDMuODk3Mzg1IC0wLjkwMjYxNSAzLjg3MzQ3NCAtMC44MDY5NzRDMy43Nzc4MzMgLTAuNDQyMzQxIDMuNjUyMzA0IC0wLjEzNzQ4NCAzLjQxOTE3OCAtMC4xMzc0ODRDMy4yMTU5NCAtMC4xMzc0ODQgMy4xOTgwMDcgLTAuMzUyNjc3IDMuMTk4MDA3IC0wLjQ0MjM0MUMzLjE5ODAwNyAtMC41MjAwNSAzLjIxNTk0IC0wLjU4NTgwMyAzLjI2Mzc2MSAtMC43ODMwNjRDMy4zMjM1MzcgLTEuMDEwMjEyIDMuMzIzNTM3IC0xLjAyMjE2NyAzLjM3NzMzNSAtMS4yNDMzMzdMMy42MTA0NjEgLTIuMTU3OTA4QzMuNjI4Mzk0IC0yLjIyOTYzOSAzLjY1ODI4MSAtMi4zNDkxOTEgMy42NTgyODEgLTIuMzc5MDc4QzMuNjU4MjgxIC0yLjQ1MDgwOSAzLjYwNDQ4MyAtMi41NzYzMzkgMy40NDMwODggLTIuNTc2MzM5QzMuMzQxNDY5IC0yLjU3NjMzOSAzLjIzOTg1MSAtMi41MTA1ODUgMy4xOTgwMDcgLTIuNDI2ODk5QzMuMTc0MDk3IC0yLjM4NTA1NiAzLjEzMjI1NCAtMi4xOTk3NTEgMy4xMDIzNjYgLTIuMDg2MTc3TDIuNzkxNTMyIC0wLjg0MjgzOVonLz4KPHBhdGggaWQ9J2cyLTYxJyBkPSdNMy4yMTU5NCAtNC4xOTYyNjRDMy4yNTc3ODMgLTQuMjkxOTA1IDMuMjU3NzgzIC00LjMwMzg2MSAzLjI1Nzc4MyAtNC4zMjc3NzFMMC41OTc3NTggMS4yMDc0NzJDMC41NTU5MTUgMS4zMDMxMTMgMC41NTU5MTUgMS4zMTUwNjggMC41NTU5MTUgMS4zMzg5NzlDMC41NTU5MTUgMS40MjI2NjUgMC42Mjc2NDYgMS40OTQzOTYgMC43MTEzMzMgMS40OTQzOTZDMC44MTI5NTEgMS40OTQzOTYgMC44NDI4MzkgMS40MjI2NjUgMC44Nzg3MDUgMS4zMzg5NzlMMy4yMTU5NCAtNC4xOTYyNjRaJy8+CjxwYXRoIGlkPSdnMi0xMjAnIGQ9J00zLjM2NTM4IC0yLjM0OTE5MUMzLjE1NjE2NCAtMi4yODk0MTUgMy4xMDgzNDQgLTIuMTE2MDY1IDMuMTA4MzQ0IC0yLjAyNjQwMUMzLjEwODM0NCAtMS44MzUxMTggMy4yNjM3NjEgLTEuNzkzMjc1IDMuMzQ3NDQ3IC0xLjc5MzI3NUMzLjUyMDc5NyAtMS43OTMyNzUgMy42OTQxNDcgLTEuOTM2NzM3IDMuNjk0MTQ3IC0yLjE2OTg2M0MzLjY5NDE0NyAtMi40OTI2NTMgMy4zNDE0NjkgLTIuNjM2MTE1IDMuMDM2NjEzIC0yLjYzNjExNUMyLjY0MjA5MiAtMi42MzYxMTUgMi40MDI5ODkgLTIuMzMxMjU4IDIuMzM3MjM1IC0yLjIxNzY4NEMyLjI1OTUyNyAtMi4zNjcxMjMgMi4wMzIzNzkgLTIuNjM2MTE1IDEuNTg0MDYgLTIuNjM2MTE1QzAuODk2NjM4IC0yLjYzNjExNSAwLjQ5MDE2MiAtMS45MjQ3ODIgMC40OTAxNjIgLTEuNzE1NTY3QzAuNDkwMTYyIC0xLjY4NTY3OSAwLjUxNDA3MiAtMS42MzE4OCAwLjcwOTM0IC0xLjYzMTY4OEMwLjc2NTEzMSAtMS43MjE1NDNDMC44NjY3NSAtMi4yMDU3MjkgMS4yNzMyMjUgLTIuNDM4ODU0IDEuNTY2MTI3IC0yLjQzODg1NFMxLjk0ODY5MiAtMi4yNDc1NzIgMS45NDg2OTIgLTIuMDUwMzExQzEuOTQ4NjkyIC0xLjk3ODU4IDEuOTQ4NjkyIC0xLjkyNDc4MiAxLjkwMDg3MiAtMS43Mzk0NzdDMS43NjMzODcgLTEuMTgzNTYyIDEuNjMxODggLTAuNjM5NjAxIDEuNjAxOTkzIC0wLjU2Nzg3QzEuNTEyMzI5IC0wLjM0MDcyMiAxLjI5NzEzNiAtMC4xMzc0ODQgMS4wNDYwNzcgLTAuMTM3NDg0QzEuMDEwMjEyIC0wLjEzNzQ4NCAwLjg0MjgzOSAtMC4xMzc0ODQgMC43MDUzNTUgLTAuMjI3MTQ4QzAuOTM4NDgxIC0wLjMwNDg1NyAwLjk2MjM5MSAtMC41MDIxMTcgMC45NjIzOTEgLTAuNTQ5OTM4QzAuOTYyMzkxIC0wLjcwNTM1NSAwLjg0MjgzOSAtMC43ODMwNjQgMC43MjMyODggLTAuNzgzMDY0QzAuNTU1OTE1IC0wLjc4MzA2NCAwLjM3NjU4OCAtMC42NTE1NTcgMC4zNzY1ODggLTAuNDA2NDc2QzAuMzc2NTg4IC0wLjA2NTc1MyAwLjc1MzE3NiAwLjA1OTc3NiAxLjAzNDEyMiAwLjA1OTc3NkMxLjM3NDg0NCAwLjA1OTc3NiAxLjYxOTkyNSAtMC4xNzMzNSAxLjczMzQ5OSAtMC4zNTg2NTVDMS44NTMwNTEgLTAuMTA3NTk3IDIuMTM5OTc1IDAuMDU5Nzc2IDIuNDgwNjk3IDAuMDU5Nzc2QzMuMTg2MDUyIDAuMDU5Nzc2IDMuNTgwNTczIC0wLjY2MzUxMiAzLjU4MDU3MyAtMC44NjA3NzJDMy41ODA1NzMgLTAuODcyNzI3IDMuNTc0NTk1IC0wLjk0NDQ1OCAzLjQ2Njk5OSAtMC45NDQ0NThDMy4zODMzMTMgLTAuOTQ0NDU4IDMuMzcxMzU3IC0wLjkwMjYxNSAzLjM1MzQyNSAtMC44NDg4MTdDMy4xODAwNzUgLTAuMzI4NzY3IDIuNzU1NjY2IC0wLjEzNzQ4NCAyLjUwNDYwOCAtMC4xMzc0ODRDMi4yNzc0NiAtMC4xMzc0ODQgMi4xMTYwNjUgLTAuMjY4OTkxIDIuMTE2MDY1IC0wLjUyMDA1QzIuMTE2MDY1IC0wLjYzMzYyNCAyLjE0NTk1MyAtMC43NjUxMzEgMi4xOTk3NTEgLTAuOTc0MzQ2TDIuMzkxMDM0IC0xLjc1MTQzMkMyLjQ1MDgwOSAtMS45ODQ1NTggMi40ODA2OTcgLTIuMDkyMTU0IDIuNjA2MjI3IC0yLjIzNTYxNkMyLjY4OTkxMyAtMi4zMjUyOCAyLjgzMzM3NSAtMi40Mzg4NTQgMy4wMjQ2NTggLTIuNDM4ODU0QzMuMDU0NTQ1IC0yLjQzODg1NCAzLjIzMzg3MyAtMi40Mzg4NTQgMy4zNjUzOCAtMi4zNDkxOTFaJy8+CjxwYXRoIGlkPSdnNy00MCcgZD0nTTMuODg1NDMgMi45MDUxMDZDMy44ODU0MyAyLjg2OTI0IDMuODg1NDMgMi44NDUzMyAzLjY4MjE5MiAyLjY0MjA5MkMyLjQ4NjY3NSAxLjQzNDYyIDEuODE3MTg2IC0wLjUzNzk4MyAxLjgxNzE4NiAtMi45NzY4MzdDMS44MTcxODYgLTUuMjk2MTM5IDIuMzc5MDc4IC03LjI5MjY1MyAzLjc2NTg3OCAtOC43MDMzNjJDMy44ODU0MyAtOC44MTA5NTkgMy44ODU0MyAtOC44MzQ4NjkgMy44ODU0MyAtOC44NzA3MzVDMy44ODU0MyAtOC45NDI0NjYgMy44MjU2NTQgLTguOTY2Mzc2IDMuNzc3ODMzIC04Ljk2NjM3NkMzLjYyMjQxNiAtOC45NjYzNzYgMi42NDIwOTIgLTguMTA1NjA0IDIuMDU2Mjg5IC02LjkzMzk5OEMxLjQ0NjU3NSAtNS43MjY1MjYgMS4xNzE2MDYgLTQuNDQ3MzIzIDEuMTcxNjA2IC0yLjk3NjgzN0MxLjE3MTYwNiAtMS45MTI4MjcgMS4zMzg5NzkgLTAuNDkwMTYyIDEuOTYwNjQ4IDAuNzg5MDQxQzIuNjY2MDAyIDIuMjIzNjYxIDMuNjQ2MzI2IDMuMDAwNzQ3IDMuNzc3ODMzIDMuMDAwNzQ3QzMuODI1NjU0IDMuMDAwNzQ3IDMuODg1NDMgMi45NzY4MzcgMy44ODU0MyAyLjkwNTEwNlonLz4KPHBhdGggaWQ9J2c3LTQxJyBkPSdNMy4zNzEzNTcgLTIuOTc2ODM3QzMuMzcxMzU3IC0zLjg4NTQzIDMuMjUxODA2IC01LjM2Nzg3IDIuNTgyMzE2IC02Ljc1NDY3QzEuODc2OTYxIC04LjE4OTI5IDAuODk2NjM4IC04Ljk2NjM3NiAwLjc2NTEzMSAtOC45NjYzNzZDMC43MTczMSAtOC45NjYzNzYgMC42NTc1MzQgLTguOTQyNDY2IDAuNjU3NTM0IC04Ljg3MDczNUMwLjY1NzUzNCAtOC44MzQ4NjkgMC42NTc1MzQgLTguODEwOTU5IDAuODYwNzcyIC04LjYwNzcyMUMyLjA1NjI4OSAtNy40MDAyNDkgMi43MjU3NzggLTUuNDI3NjQ2IDIuNzI1Nzc4IC0yLjk4ODc5MkMyLjcyNTc3OCAtMC42Njk0ODkgMi4xNjM4ODUgMS4zMjcwMjQgMC43NzcwODYgMi43Mzc3MzNDMC42NTc1MzQgMi44NDUzMyAwLjY1NzUzNCAyLjg2OTI0IDAuNjU3NTM0IDIuOTA1MTA2QzAuNjU3NTM0IDIuOTc2ODM3IDAuNzE3MzEgMy4wMDA3NDcgMC43NjUxMzEgMy4wMDA3NDdDMC45MjA1NDggMy4wMDA3NDcgMS45MDA4NzIgMi4xMzk5NzUgMi40ODY2NzUgMC45NjgzNjlDMy4wOTYzODkgLTAuMjUxMDU5IDMuMzcxMzU3IC0xLjU0MjIxNyAzLjM3MTM1NyAtMi45NzY4MzdaJy8+CjxwYXRoIGlkPSdnNy02MScgZD0nTTguMDY5NzM4IC0zLjg3MzQ3NEM4LjIzNzExMSAtMy44NzM0NzQgOC40NTIzMDQgLTMuODczNDc0IDguNDUyMzA0IC00LjA4ODY2N0M4LjQ1MjMwNCAtNC4zMTU4MTYgOC4yNDkwNjYgLTQuMzE1ODE2IDguMDY5NzM4IC00LjMxNTgxNkgxLjAyODE0NEMwLjg2MDc3MiAtNC4zMTU4MTYgMC42NDU1NzkgLTQuMzE1ODE2IDAuNjQ1NTc5IC00LjEwMDYyM0MwLjY0NTU3OSAtMy44NzM0NzQgMC44NDg4MTcgLTMuODczNDc0IDEuMDI4MTQ0IC0zLjg3MzQ3NEg4LjA2OTczOFpNOC4wNjk3MzggLTEuNjQ5ODEzQzguMjM3MTExIC0xLjY0OTgxMyA4LjQ1MjMwNCAtMS42NDk4MTMgOC40NTIzMDQgLTEuODY1MDA2QzguNDUyMzA0IC0yLjA5MjE1NCA4LjI0OTA2NiAtMi4wOTIxNTQgOC4wNjk3MzggLTIuMDkyMTU0SDEuMDI4MTQ0QzAuODYwNzcyIC0yLjA5MjE1NCAwLjY0NTU3OSAtMi4wOTIxNTQgMC42NDU1NzkgLTEuODc2OTYxQzAuNjQ1NTc5IC0xLjY0OTgxMyAwLjg0ODgxNyAtMS42NDk4MTMgMS4wMjgxNDQgLTEuNjQ5ODEzSDguMDY5NzM4WicvPgo8cGF0aCBpZD0nZzMtMTInIGQ9J000Ljc4MjA2NyAtNC41MzQ5OTRDNC43ODIwNjcgLTUuMjYwMjc0IDQuMTkyMjc5IC01LjYxODkyOSAzLjUzMDc2IC01LjYxODkyOUMyLjUwMjYxNSAtNS42MTg5MjkgMS43Mzc0ODQgLTQuNDU1MjkzIDEuNTE0MzIxIC0zLjc3ODU4TDAuMjU1MDQ0IDEuNDUwNTZDMC4yMzkxMDMgMS41MDYzNTEgMC4zMDI4NjQgMS41NDYyMDIgMC4zNTA2ODUgMS41NDYyMDJDMC40MTQ0NDYgMS41NDYyMDIgMC40NzgyMDcgMS41NDYyMDIgMC40OTQxNDcgMS40OTA0MTFMMS4wMzYxMTUgLTAuNjQ1NTc5QzEuMjE5NDI3IC0wLjIxNTE5MyAxLjU5NDAyMiAwLjA3MTczMSAyLjE2Nzg3IDAuMDcxNzMxQzMuMjY3NzQ2IDAuMDcxNzMxIDQuNDcxMjMzIC0wLjcwOTM0IDQuNDcxMjMzIC0xLjk3NjU4OEM0LjQ3MTIzMyAtMi41MDI2MTUgNC4yNzE5OCAtMi44NTMzIDMuODczNDc0IC0zLjE4MDA3NUMzLjkyMTI5NSAtMy4xODgwNDUgNC4wNjQ3NTcgLTMuMjkxNjU2IDQuMTA0NjA4IC0zLjMyMzUzN0M0LjUwMzExMyAtMy41OTQ1MjEgNC43ODIwNjcgLTQuMDQ4ODE3IDQuNzgyMDY3IC00LjUzNDk5NFpNMy4yMTE5NTUgLTMuMTgwMDc1QzMuMDY4NDkzIC0zLjE1NjE2NCAyLjkyNTAzMSAtMy4xNDAyMjQgMi43NzM1OTkgLTMuMTQwMjI0QzIuNjg1OTI4IC0zLjE0MDIyNCAyLjU5ODI1NyAtMy4xNDAyMjQgMi41MTA1ODUgLTMuMTcyMTA1QzIuNjYyMDE3IC0zLjIxMTk1NSAyLjgyMTQyIC0zLjIxMTk1NSAyLjk3Mjg1MiAtMy4yMTE5NTVDMy4wNTI1NTMgLTMuMjExOTU1IDMuMTMyMjU0IC0zLjIwMzk4NSAzLjIxMTk1NSAtMy4xODAwNzVaTTE0LjI3OTk1IC00LjY4NjQyNkMzLjgyNTY1NCAtMy4zMjUyOCAyLjgyOTM5IC0zLjQzNTExOCAyLjk0MDk3MSAtMy40MzUxMThDMi43MzM3NDggLTMuNDM1MTE4IDIuMjMxNjMxIC0zLjQ1MTA1OSAyLjIzMTYzMSAtMy4xNDAyMjRDMi4yMzE2MzEgLTIuOTAxMTIxIDIuNjg1OTI4IC0yLjkxNzA2MSAyLjg0NTMzIC0yLjkxNzA2MUMzLjA5MjQwMyAtMi45MTcwNjEgMy4zMjM1MzcgLTIuOTQwOTcxIDMuNTU0NjcgLTMuMDI4NjQzQzMuNzkzNzczIC0yLjgzNzM2IDMuOTEzMzI1IC0yLjU4MjMxNiAzLjkxMzMyNSAtMi4xODM4MTFDMy45MTMzMjUgLTEuNzEzNTc0IDMuNzUzOTIzIC0xLjE4NzU0NyAzLjUyMjc5IC0wLjg2ODc0MkMzLjIyNzg5NSAtMC40NDYzMjYgMi43MzM3NDggLTAuMTUxNDMyIDIuMTM1OTkgLTAuMTUxNDMyQzEuNjA5OTYzIC0wLjE1MTQzMiAxLjE5

где это параметр сдвига, а это параметр масштаба.

Распределение Гамбеля (названное в честь немецкого математика Эмиля Юлиуса Гамбеля) было использовано очень рано в литературе по гидрологии для моделирования наводнений. Оно также используется для моделирования максимальной скорости ветра и интенсивности осадков. Это распределение с «толстыми хвостами» — вероятность события в хвосте распределения больше, чем если бы использовалась гауссова функция, поэтому 100-летние наводнения происходят удивительно часто. Наводнения первоначально моделировались как гауссовский процесс, что недооценивало частоту экстремальных событий.

Это одно из семейства распределений экстремальных значений, распределения обобщённых экстремальных значений (GEV), которое также включает в себя распределения Вейбулла и Фреше.

Функция имеет среднее значение

и дисперсия .

Литература

1

Гамбел, Е. Дж., «Статистика экстремумов», Нью-Йорк: Издательство Колумбийского университета, 1958.

2

Райсс, Р.-Д. и Томас, М., «Статистический анализ экстремальных значений из страхования, финансов, гидрологии и других областей», Базель: Birkhauser Verlag, 2001.

Примеры

Генерация выборок из распределения:

>>> mu, beta = 0, 0.1 # location and scale
>>> s = np.random.gumbel(mu, beta, 1000)

Отображение гистограммы выборок вместе с функцией плотности вероятности:

>>> import matplotlib.pyplot as plt
>>> count, bins, ignored = plt.hist(s, 30, density=True)
>>> plt.plot(bins, (1/beta)*np.exp(-(bins - mu)/beta)
...          * np.exp( -np.exp( -(bins - mu) /beta) ),
...          linewidth=2, color='r')
>>> plt.show()

Покажите, как распределение экстремальных значений может возникнуть из гауссовского процесса и сравните его с гауссовым:

>>> means = []
>>> maxima = []
>>> for i in range(0,1000) :
...    a = np.random.normal(mu, beta, 1000)
...    means.append(a.mean())
...    maxima.append(a.max())
>>> count, bins, ignored = plt.hist(maxima, 30, density=True)
>>> beta = np.std(maxima) * np.sqrt(6) / np.pi
>>> mu = np.mean(maxima) - 0.57721*beta
>>> plt.plot(bins, (1/beta)*np.exp(-(bins - mu)/beta)
...          * np.exp(-np.exp(-(bins - mu)/beta)),
...          linewidth=2, color='r')
>>> plt.plot(bins, 1/(beta * np.sqrt(2 * np.pi))
...          * np.exp(-(bins - mu)**2 / (2 * beta**2)),
...          linewidth=2, color='g')
>>> plt.show()