numpy.random.gumbel

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

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

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

Примечание

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

Параметры

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,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0nVVRGLTgnPz4KPCEtLSBUaGlzIGZpbGUgd2FzIGdlbmVyYXRlZCBieSBkdmlzdmdtIDIuOS4xIC0tPgo8c3ZnIHZlcnNpb249JzEuMScgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJyB4bWxuczp4bGluaz0naHR0cDovL3d3dy53My5vcmcvMTk5OS94bGluaycgd2lkdGg9JzEzOS4xMTk0MDRwdCcgaGVpZ2h0PScyOC45MjkwNDlwdCcgdmlld0JveD0nMTYzLjU2NjA3MSA3OC43MDQ4NTcgMTM5LjExOTQwNCAyOC45MjkwNDknPgo8ZGVmcz4KPHBhdGggaWQ9J2cyLTEyJyBkPSdNNC4wMDQ5ODEtMy40MDEyNDVDNC4wMDQ5ODEtMy44NDM1ODcgMy41OTI1MjgtNC4yMTQxOTcgMi45MjkwMTYtNC4yMTQxOTdDMi43NDk2ODktNC4yMTQxOTcgMi41NTg0MDYtNC4xNjAzOTkgMi4zOTcwMTEtNC4wNzY3MTJDMS44NDcwNzMtMy43OTU3NjYgMS40NTI1NTMtMy4yMDM5ODUgMS4zMjEwNDYtMi42NzE5OEwuMzgyNTY1IDEuMDgxOTQzQy4zNjQ2MzMgMS4xNDE3MTkgLjQ0MjM0MSAxLjE1OTY1MSAuNTAyMTE3IDEuMTU5NjUxQy41NTU5MTUgMS4xNTk2NTEgLjU5MTc4MSAxLjE1MzY3NCAuNjAzNzM2IDEuMDk5ODc1Qy42MDk3MTQgMS4wOTM4OTggLjkzMjUwMy0uMTk3MjYgLjk5ODI1Ny0uNDY2MjUyQzEuMTc3NTg0LS4xMDE2MTkgMS41NTQxNzIgLjA1OTc3NiAxLjk0MjcxNSAuMDU5Nzc2QzIuNzY3NjIxIC4wNTk3NzYgMy43OTU3NjYtLjUwODA5NSAzLjc5NTc2Ni0xLjQzNDYyQzMuNzk1NzY2LTEuNzU3NDEgMy42NzYyMTQtMi4xMDQxMSAzLjMxMTU4Mi0yLjM2MTE0NkMzLjgxMzY5OS0yLjY3Nzk1OCA0LjAwNDk4MS0zLjAyNDY1OCA0LjAwNDk4MS0zLjQwMTI0NVpNMi43MDc4NDYtMi4zNjcxMjNDMi42NDgwNy0yLjM1NTE2OCAyLjQ5ODYzLTIuMzQzMjEzIDIuNDM4ODU0LTIuMzQzMjEzUzIuMjIzNjYxLTIuMzQzMjEzIDIuMTgxODE4LTIuMzYxMTQ2QzIuMjc3NDYtMi4zNzMxMDEgMi4zNzMxMDEtMi4zNzkwNzggMi40Njg3NDItMi4zNzkwNzhDMi41NDA0NzMtMi4zNzkwNzggMi42MTgxODItMi4zNzMxMDEgMi43MDc4NDYtMi4zNjcxMjNaTTMuNTgwNTczLTMuNDYxMDIxQzMuNTgwNTczLTMuMDM2NjEzIDMuMzIzNTM3LTIuNjYwMDI1IDMuMDQyNTktMi40OTg2M0MyLjg2OTI0LTIuNTU4NDA2IDIuNjg5OTEzLTIuNTc2MzM5IDIuNTEwNTg1LTIuNTc2MzM5QzIuMzYxMTQ2LTIuNTc2MzM5IDEuOTA2ODQ5LTIuNTg4Mjk0IDEuOTA2ODQ5LTIuMzQ5MTkxQzEuOTA2ODQ5LTIuMTI4MDIgMi4yODk0MTUtMi4xNDU5NTMgMi40Mzg4NTQtMi4xNDU5NTNDMi41ODIzMTYtMi4xNDU5NTMgMi44OTkxMjgtMi4xNjk4NjMgMy4wMjQ2NTgtMi4yNDE1OTRDMy4xOTIwMy0yLjEyMjA0MiAzLjMyOTUxNC0xLjk0ODY5MiAzLjMyOTUxNC0xLjU5MDAzN0MzLjMyOTUxNC0uNzc3MDg2IDIuNzc5NTc3LS4xMzc0ODQgMS45MzY3MzctLjEzNzQ4NEMxLjUwMDM3NC0uMTM3NDg0IDEuMTI5NzYzLS40MTg0MzEgMS4xMjk3NjMtLjg3ODcwNUMxLjEyOTc2My0uOTMyNTAzIDEuMTI5NzYzLS45ODAzMjQgMS4xNDE3MTktMS4wMjgxNDRMMS41NDgxOTQtMi42NDgwN0MxLjY5MTY1Ni0zLjIyNzg5NSAyLjI2NTUwNC00LjAxNjkzNiAyLjkzNDk5NC00LjAxNjkzNkMzLjI5MzY0OS00LjAxNjkzNiAzLjU4MDU3My0zLjg1NTU0MiAzLjU4MDU3My0zLjQ2MTAyMVonLz4KPHBhdGggaWQ9J2cyLTIyJyBkPSdNMi43OTE1MzItLjg0MjgzOUMyLjc2MTY0NC0uNzI5MjY1IDIuNzYxNjQ0LS43MTczMSAyLjY1NDA0Ny0uNTc5ODI2QzIuNTIyNTQtLjQxODQzMSAyLjI0MTU5NC0uMTM3NDg0IDEuODcwOTg0LS4xMzc0ODRDMS4zNDQ5NTYtLjEzNzQ4NCAxLjM0NDk1Ni0uNTQ5OTM4IDEuMzQ0OTU2LS42Njk0ODlDMS4zNDQ5NTYtLjgyNDkwNyAxLjM2Mjg4OS0uOTAyNjE1IDEuNDEwNzEtMS4wOTk4NzVDMS40NjQ1MDgtMS4yOTcxMzYgMS40NzY0NjMtMS4zNDQ5NTYgMS41MzAyNjItMS41NjYxMjdTMS41OTYwMTUtMS44MzUxMTggMS42NDk4MTMtMi4wMzIzNzlDMS42Nzk3MDEtMi4xNjk4NjMgMS43Mzk0NzctMi40MDI5ODkgMS43Mzk0NzctMi40MzI4NzdDMS43Mzk0NzctMi41NzYzMzkgMS42MjU5MDMtMi42MzYxMTUgMS41MzAyNjItMi42MzYxMTVDMS4zMDkwOTEtMi42MzYxMTUgMS4yNTUyOTMtMi40MjY4OTkgMS4yMzczNi0yLjM0OTE5MUwuNDAwNDk4IDEuMDA0MjM0Qy4zODg1NDMgMS4wNTIwNTUgLjM4MjU2NSAxLjA2OTk4OCAuMzgyNTY1IDEuMDg3OTJDLjM4MjU2NSAxLjE1OTY1MSAuNDM2MzY0IDEuMjg1MTgxIC41OTc3NTggMS4yODUxODFDLjgwNjk3NCAxLjI4NTE4MSAuODYwNzcyIDEuMTA1ODUzIC44OTA2NiAuOTgwMzI0Qy45MjA1NDggLjg4NDY4MiAxLjAwNDIzNCAuNTE0MDcyIDEuMTcxNjA2LS4xMTk1NTJDMS40MTY2ODcgLjA1OTc3NiAxLjc4MTMyIC4wNTk3NzYgMS44NDcwNzMgLjA1OTc3NkMyLjMwMTM3IC4wNTk3NzYgMi41NzYzMzktLjIwOTIxNSAyLjc0MzcxMS0uMzc2NTg4QzIuODU3Mjg1IC4wMTE5NTUgMy4yNjM3NjEgLjA1OTc3NiAzLjM5NTI2OCAuMDU5Nzc2QzMuNjEwNDYxIC4wNTk3NzYgMy43NjU4NzgtLjA1OTc3NiAzLjg3OTQ1Mi0uMjQ1MDgxQzQuMDI4ODkyLS40ODQxODQgNC4xMTI1NzgtLjgzMDg4NCA0LjExMjU3OC0uODYwNzcyQzQuMTEyNTc4LS44NzI3MjcgNC4xMDY2LS45NDQ0NTggMy45OTkwMDQtLjk0NDQ1OEMzLjkwOTM0LS45NDQ0NTggMy44OTczODUtLjkwMjYxNSAzLjg3MzQ3NC0uODA2OTc0QzMuNzc3ODMzLS40NDIzNDEgMy42NTIzMDQtLjEzNzQ4NCAzLjQxOTE3OC0uMTM3NDg0QzMuMjE1OTQtLjEzNzQ4NCAzLjE5ODAwNy0uMzUyNjc3IDMuMTk4MDA3LS40NDIzNDFDMy4xOTgwMDctLjUyMDA1IDMuMjE1OTQtLjU4NTgwMyAzLjI2Mzc2MS0uNzgzMDY0QzMuMzIzNTM3LTEuMDEwMjEyIDMuMzIzNTM3LTEuMDIyMTY3IDMuMzc3MzM1LTEuMjQzMzM3TDMuNjEwNDYxLTIuMTU3OTA4QzMuNjI4Mzk0LTIuMjI5NjM5IDMuNjU4MjgxLTIuMzQ5MTkxIDMuNjU4MjgxLTIuMzc5MDc4QzMuNjU4MjgxLTIuNDUwODA5IDMuNjA0NDgzLTIuNTc2MzM5IDMuNDQzMDg4LTIuNTc2MzM5QzMuMzQxNDY5LTIuNTc2MzM5IDMuMjM5ODUxLTIuNTEwNTg1IDMuMTk4MDA3LTIuNDI2ODk5QzMuMTc0MDk3LTIuMzg1MDU2IDMuMTMyMjU0LTIuMTk5NzUxIDMuMTAyMzY2LTIuMDg2MTc3TDIuNzkxNTMyLS44NDI4MzlaJy8+CjxwYXRoIGlkPSdnMi02MScgZD0nTTMuMjE1OTQtNC4xOTYyNjRDMy4yNTc3ODMtNC4yOTE5MDUgMy4yNTc3ODMtNC4zMDM4NjEgMy4yNTc3ODMtNC4zMjc3NzFDMy4yNTc3ODMtNC40MTE0NTcgMy4xODYwNTItNC40ODMxODggMy4xMDIzNjYtNC40ODMxODhDMy4wMDA3NDctNC40ODMxODggMi45NzA4NTktNC40MTE0NTcgMi45MzQ5OTQtNC4zMjc3NzFMLjU5Nzc1OCAxLjIwNzQ3MkMuNTU1OTE1IDEuMzAzMTEzIC41NTU5MTUgMS4zMTUwNjggLjU1NTkxNSAxLjMzODk3OUMuNTU1OTE1IDEuNDIyNjY1IC42Mjc2NDYgMS40OTQzOTYgLjcxMTMzMyAxLjQ5NDM5NkMuODEyOTUxIDEuNDk0Mzk2IC44NDI4MzkgMS40MjI2NjUgLjg3ODcwNSAxLjMzODk3OUwzLjIxNTk0LTQuMTk2MjY0WicvPgo8cGF0aCBpZD0nZzItMTIwJyBkPSdNMy4zNjUzOC0yLjM0OTE5MUMzLjE1NjE2NC0yLjI4OTQxNSAzLjEwODM0NC0yLjExNjA2NSAzLjEwODM0NC0yLjAyNjQwMUMzLjEwODM0NC0xLjgzNTExOCAzLjI2Mzc2MS0xLjc5MzI3NSAzLjM0NzQ0Ny0xLjc5MzI3NUMzLjUyMDc5Ny0xLjc5MzI3NSAzLjY5NDE0Ny0xLjkzNjczNyAzLjY5NDE0Ny0yLjE2OTg2M0MzLjY5NDE0Ny0yLjQ5MjY1MyAzLjM0MTQ2OS0yLjYzNjExNSAzLjAzNjYxMy0yLjYzNjExNUMyLjY0MjA5Mi0yLjYzNjExNSAyLjQwMjk4OS0yLjMzMTI1OCAyLjMzNzIzNS0yLjIxNzY4NEMyLjI1OTUyNy0yLjM2NzEyMyAyLjAzMjM3OS0yLjYzNjExNSAxLjU4NDA2LTIuNjM2MTE1Qy44OTY2MzgtMi42MzYxMTUgLjQ5MDE2Mi0xLjkyNDc4MiAuNDkwMTYyLTEuNzE1NTY3Qy40OTAxNjItMS42ODU2NzkgLjUxNDA3Mi0xLjYzMTg4IC41OTc3NTgtMS42MzE4OFMuNjk5Mzc3LTEuNjY3NzQ2IC43MTczMS0xLjcyMTU0NEMuODY2NzUtMi4yMDU3MjkgMS4yNzMyMjUtMi40Mzg4NTQgMS41NjYxMjctMi40Mzg4NTRTMS45NDg2OTItMi4yNDc1NzIgMS45NDg2OTItMi4wNTAzMTFDMS45NDg2OTItMS45Nzg1OCAxLjk0ODY5Mi0xLjkyNDc4MiAxLjkxNDcyLTEuNzM5NDc3QzEuNzYzMzg3LTEuMTgzNTYyIDEuNjMxMTg4LS42Mzk2MDEgMS42MDE5OTMtLjU2Nzg3QzEuNTEyMzI5LS4zNDA3MjIgMS4yOTcxMzYtLjEzNzQ4NCAxLjA0NjA3Ny0uMTM3NDg0QzEuMDEwMjEyLS4xMzc0ODQgLjg0MjgzOS0uMTM3NDg0IC43MDUzNTUtLjIyNzE0OEMuOTM4NDgxLS4zMDQ4NTcgLjk2MjM5MS0uNTAyMTE3IC45NjIzOTEtLjU0OTkzOEMuOTYyMzkxLS43MDUzNTUgLjg0MjgzOS0uNzgzMDY0IC43MjMyODgtLjc4MzA2NEMuNTU1OTE1LS43ODMwNjQgLjM3NjU4OC0uNjUxNTU3IC4zNzY1ODgtLjQwNjQ3NkMuMzc2NTg4LS4wNjU3NTMgLjc1MzE3NiAuMDU5Nzc2IDEuMDM0MTIyIC4wNTk3NzZDMS4zNzQ4NDQgLjA1OTc3NiAxLjYxOTkyNS0uMTczMzUgMS43MzM0OTktLjM1ODY1NUMxLjg1MzA1MS0uMTA3NTk3IDIuMTM5OTc1IC4wNTk3NzYgMi40ODA2OTcgLjA1OTc3NkMzLjE4NjA1MiAuMDU5Nzc2IDMuNTgwNTczLS42NjM1MTIgMy41ODA1NzMtLjg2MDc3MkMzLjU4MDU3My0uODcyNzI3IDMuNTc0NTk1LS45NDQ0NTggMy40NjY5OTktLjk0NDQ1OEMzLjM4MzMxMy0uOTQ0NDU4IDMuMzcxMzU3LS45MDI2MTUgMy4zNTM0MjUtLjg0ODgxN0MzLjE4MDA3NS0uMzI4NzY3IDIuNzU1NjY2LS4xMzc0ODQgMi41MDQ2MDgtLjEzNzQ4NEMyLjI3NzQ2LS4xMzc0ODQgMi4xMTYwNjUtLjI2ODk5MSAyLjExNjA2NS0uNTIwMDVDMi4xMTYwNjUtLjYzMzYyNCAyLjE0NTk1My0uNzY1MTMxIDIuMTk5NzUxLS45NzQzNDZMMi4zOTEwMzQtMS43NTE0MzJDMi40NTA4MDktMS45ODQ1NTggMi40ODA2OTctMi4wOTIxNTQgMi42MDYyMjctMi4yMzU2MTZDMi42ODk5MTMtMi4zMjUyOCAyLjgzMzM3NS0yLjQzODg1NCAzLjAyNDY1OC0yLjQzODg1NEMzLjA1NDU0NS0yLjQzODg1NCAzLjIzMzg3My0yLjQzODg1NCAzLjM2NTM4LTIuMzQ5MTkxWicvPgo8cGF0aCBpZD0nZzUtNDAnIGQ9J00yLjM3OTA3OCAxLjM1MDkzNEMxLjcyMTU0NCAuODQ4ODE3IDEuMTgzNTYyLS4wOTU2NDEgMS4xODM1NjItMS40OTQzOTZDMS4xODM1NjItMi44NzUyMTggMS43MDk1ODktMy44MzE2MzEgMi4zNzkwNzgtNC4zMzk3MjZDMi4zNzkwNzgtNC4zNDU3MDQgMi4zOTcwMTEtNC4zNjk2MTQgMi4zOTcwMTEtNC4zOTk1MDJDMi4zOTcwMTEtNC40MzUzNjcgMi4zNzMxMDEtNC40ODMxODggMi4zMDEzNy00LjQ4MzE4OEMyLjE5Mzc3My00LjQ4MzE4OCAuNzUzMTc2LTMuNTQ0NzA3IC43NTMxNzYtMS40OTQzOTZDLjc1MzE3NiAuNTQ5OTM4IDIuMTg3Nzk2IDEuNDk0Mzk2IDIuMzAxMzcgMS40OTQzOTZDMi4zNzMxMDEgMS40OTQzOTYgMi4zOTcwMTEgMS40NDY1NzUgMi4zOTcwMTEgMS40MTA3MUMyLjM5NzAxMSAxLjM4MDgyMiAyLjM3OTA3OCAxLjM1NjkxMiAyLjM3OTA3OCAxLjM1MDkzNFonLz4KPHBhdGggaWQ9J2c1LTQxJyBkPSdNMi4xMTYwNjUtMS40OTQzOTZDMi4xMTYwNjUtMy41Mzg3MyAuNjgxNDQ1LTQuNDgzMTg4IC41Njc4Ny00LjQ4MzE4OEMuNDk2MTM5LTQuNDgzMTg4IC40NzgyMDctNC40MjkzOSAuNDc4MjA3LTQuMzk5NTAyQy40NzgyMDctNC4zNjM2MzYgLjQ4NDE4NC00LjM1NzY1OSAuNTkxNzgxLTQuMjYyMDE3QzEuMTI5NzYzLTMuODAxNzQzIDEuNjg1Njc5LTIuOTE3MDYxIDEuNjg1Njc5LTEuNDk0Mzk2QzEuNjg1Njc5LS4yMjcxNDggMS4yNDMzMzcgLjcyMzI4OCAuNTY3ODcgMS4yOTExNThDLjQ4NDE4NCAxLjM2ODg2NyAuNDc4MjA3IDEuMzc0ODQ0IC40NzgyMDcgMS40MTA3MUMuNDc4MjA3IDEuNDQwNTk4IC40OTYxMzkgMS40OTQzOTYgLjU2Nzg3IDEuNDk0Mzk2Qy42NzU0NjcgMS40OTQzOTYgMi4xMTYwNjUgLjU1NTkxNSAyLjExNjA2NS0xLjQ5NDM5NlonLz4KPHBhdGggaWQ9J2czLTEyJyBkPSdNNC43ODIwNjctNC41MzQ5OTRDNC43ODIwNjctNS4yNjAyNzQgNC4xOTIyNzktNS42MTg5MjkgMy41MzA3Ni01LjYxODkyOUMyLjUwMjYxNS01LjYxODkyOSAxLjczNzQ4NC00LjQ1NTI5MyAxLjUxNDMyMS0zLjU3ODU4TC4yNTUwNDQgMS40NTA1NkMuMjM5MTAzIDEuNTA2MzUxIC4zMDI4NjQgMS41NDYyMDIgLjM1MDY4NSAxLjU0NjIwMkMuNDE0NDQ2IDEuNTQ2MjAyIC40NzgyMDcgMS41NDYyMDIgLjQ5NDE0NyAxLjQ5MDQxMUwxLjAzNjExNS0uNjQ1NTc5QzEuMjE5NDI3LS4yMTUxOTMgMS41OTQwMjIgLjA3MTczMSAyLjE2Nzg3IC4wNzE3MzFDMy4yNjc3NDYgLjA3MTczMSA0LjQ3MTIzMy0uNzA5MzQgNC40NzEyMzMtMS45NzY1ODhDNC40NzEyMzMtMi41MDI2MTUgNC4yNzE5OC0yLjg1MzMgMy44NzM0NzQtMy4xODAwNzVDMy45MjEyOTUtMy4xODgwNDUgNC4wNjQ3NTctMy4yOTE2NTYgNC4xMDQ2MDgtMy4zMjM1MzdDNC41MDMxMTMtMy41OTQ1MjEgNC43ODIwNjctNC4wNDg4MTcgNC43ODIwNjctNC41MzQ5OTRaTTMuMjExOTU1LTMuMTgwMDc1QzMuMDY4NDkzLTMuMTU2MTY0IDIuOTI1MDMxLTMuMTQwMjI0IDIuNzczNTk5LTMuMTQwMjI0QzIuNjg1OTI4LTMuMTQwMjI0IDIuNTk4MjU3LTMuMTQwMjI0IDIuNTEwNTg1LTMuMTcyMTA1QzIuNjYyMDE3LTMuMjExOTU1IDIuODIxNDItMy4yMTE5NTUgMi45NzI4NTItMy4yMTE5NTVDMy4wNTI1NTMtMy4yMTE5NTUgMy4xMzIyNTQtMy4yMDM5ODUgMy4yMTE5NTUtMy4xODAwNzVaTTQuMjc5OTUtNC42ODY0MjZDNC4yNzk5NS00LjI4NzkyIDQuMTA0NjA4LTMuODE3Njg0IDMuODI1NjU0LTMuNTM4NzNDMy43ODU4MDMtMy40OTg4NzkgMy42MTg0MzEtMy4zNDc0NDcgMy41NzA2MS0zLjMzOTQ3N0MzLjM0NzQ0Ny0zLjM3OTMyOCAzLjE3MjEwNS0zLjQzNTExOCAyLjk0MDk3MS0zLjQzNTExOEMyLjczMzc0OC0zLjQzNTExOCAyLjIzMTYzMS0zLjQ1MTA1OSAyLjIzMTYzMS0zLjE0MDIyNEMyLjIzMTYzMS0yLjkwMTEyMSAyLjY4NTkyOC0yLjkxNzA2MSAyLjg0NTMzLTIuOTE3MDYxQzMuMDkyNDAzLTIuOTE3MDYxIDMuMzIzNTM3LTIuOTQwOTcxIDMuNTU0NjctMy4wMjg2NDNDMy43OTM3NzMtMi44MzczNiAzLjkxMzMyNS0yLjU4MjMxNiAzLjkxMzMyNS0yLjE4MzgxMUMzLjkxMzMyNS0xLjcxMzU3NCAzLjc1MzkyMy0xLjE4NzU0NyAzLjUyMjc5LS44Njg3NDJDMy4yMjc4OTUtLjQ0NjMyNiAyLjczMzc0OC0uMTUxNDMyIDIuMTM1OTktLjE1MTQzMkMxLjYwOTk2My0uMTUxNDMyIDEuMTk1NTE3LS41NTc5MDggMS4xOTU1MTctMS4xMTU4MTZDMS4xOTU1MTctMS4yMTk0MjcgMS4yMDM0ODctMS4zMjMwMzkgMS4yMjczOTctMS40MDI3NEwxLjc1MzQyNS0zLjUyMjc5QzEuOTI4NzY3LTQuMjI0MTU5IDIuNjA2MjI3LTUuMzk1NzY2IDMuNDk4ODc5LTUuMzk1NzY2QzMuOTA1MzU1LTUuMzk1NzY2IDQuMjc5OTUtNS4xOTY1MTMgNC4yNzk5NS00LjY4NjQyNlonLz4KPHBhdGggaWQ9J2czLTIyJyBkPSdNMS45Mjg3NjctMi44MTM0NUMxLjk2ODYxOC0yLjk2NDg4MiAyLjAzMjM3OS0zLjIyNzg5NSAyLjAzMjM3OS0zLjI2Nzc0NkMyLjAzMjM3OS0zLjQzNTExOCAxLjkwNDg1Ny0zLjUxNDgxOSAxLjc2OTM2NS0zLjUxNDgxOUMxLjQ5ODM4MS0zLjUxNDgxOSAxLjQzNDYyLTMuMjUxODA2IDEuNDAyNzQtMy4xNDAyMjRMLjI5NDg5NCAxLjI4MzE4OEMuMjYzMDE0IDEuNDEwNzEgLjI2MzAxNCAxLjQ1MDU2IC4yNjMwMTQgMS40NjY1MDFDLjI2MzAxNCAxLjY2NTc1MyAuNDIyNDE2IDEuNzEzNTc0IC41MTgwNTcgMS43MTM1NzRDLjU1NzkwOCAxLjcxMzU3NCAuNzQxMjIgMS43MDU2MDQgLjg0NDgzMiAxLjQ5ODM4MUMuODY4NzQyIDEuNDM0NjIgMS4wOTk4NzUgLjQ4NjE3NyAxLjI1OTI3OC0uMTUxNDMyQzEuMzk0NzctLjA0NzgyMSAxLjY2NTc1MyAuMDc5NzAxIDIuMDcyMjI5IC4wNzk3MDFDMi43MjU3NzggLjA3OTcwMSAzLjE1NjE2NC0uNDU0Mjk2IDMuMTgwMDc1LS40ODYxNzdDMy4zMjM1MzcgLjA2Mzc2MSAzLjg2NTUwNCAuMDc5NzAxIDMuOTYxMTQ2IC4wNzk3MDFDNC4zMjc3NzEgLjA3OTcwMSA0LjUxMTA4My0uMjIzMTYzIDQuNTc0ODQ0LS4zNTg2NTVDNC43MzQyNDctLjY0NTU3OSA0Ljg0NTgyOC0xLjEwNzg0NiA0Ljg0NTgyOC0xLjEzOTcyNkM0Ljg0NTgyOC0xLjE4NzU0NyA0LjgxMzk0OC0xLjI0MzMzNyA0LjcxODMwNi0xLjI0MzMzN1M0LjYwNjcyNS0xLjE5NTUxNyA0LjU1ODkwNC0uOTk2MjY0QzQuNDQ3MzIzLS41NTc5MDggNC4yOTU4OS0uMTQzNDYyIDMuOTg1MDU2LS4xNDM0NjJDMy44MDE3NDMtLjE0MzQ2MiAzLjczMDAxMi0uMjk0ODk0IDMuNzMwMDEyLS41MTgwNTdDMy43MzAwMTItLjY1MzU0OSAzLjgxNzY4NC0uOTk2MjY0IDMuODczNDc0LTEuMjI3Mzk3TDQuMDgwNjk3LTIuMDQwMzQ5QzQuMTI4NTE4LTIuMjQ3NTcyIDQuMTY4MzY5LTIuNDE0OTQ0IDQuMjMyMTMtMi42NTQwNDdDNC4yNzE5OC0yLjgyOTM5IDQuMzUxNjgxLTMuMTQwMjI0IDQuMzUxNjgxLTMuMTg4MDQ1QzQuMzUxNjgxLTMuMzg3Mjk4IDQuMTkyMjc5LTMuNDM1MTE4IDQuMDk2NjM4LTMuNDM1MTE4QzMuODE3Njg0LTMuNDM1MTE4IDMuNzY5ODYzLTMuMjM1ODY2IDMuNjgyMTkyLTIuODc3MjFMMy41MTQ4MTktMi4yMTU2OTFMMy4yNjc3NDYtMS4yMTk0MjdMMy4xODgwNDUtLjkwMDYyM0MzLjE3MjEwNS0uODUyODAyIDIuOTU2OTEyLS41NDk5MzggMi43NzM1OTktLjM5ODUwNkMyLjYzODEwNy0uMjk0ODk0IDIuNDA2OTc0LS4xNDM0NjIgMi4xMTIwOC0uMTQzNDYyQzEuNzM3NDg0LS4xNDM0NjIgMS40ODI0NDEtLjM0MjcxNSAxLjQ4MjQ0MS0uODM2ODYyQzEuNDgyNDQxLTEuMDQ0MDg1IDEuNTQ2MjAyLTEuMjgzMTg4IDEuNTk0MDIyLTEuNDc0NDcxTDEuOTI4NzY3LTIuODEzNDVaJy8+CjxwYXRoIGlkPSdnMy02MScgZD0nTTMuNzA2MTAyLTUuNjQyODM5QzMuNzUzOTIzLTUuNzU0NDIxIDMuNzUzOTIzLTUuNzcwMzYxIDMuNzUzOTIzLTUuNzk0MjcxQzMuNzUzOTIzLTUuODk3ODgzIDMuNjc0MjIyLTUuOTc3NTg0IDMuNTcwNjEtNS45Nzc1ODRDMy40NDMwODgtNS45Nzc1ODQgMy40MTEyMDgtNS44ODE5NDMgMy4zNzkzMjgtNS44MDIyNDJMLjUxODA1NyAxLjY1Nzc4M0MuNDcwMjM3IDEuNzY5MzY1IC40NzAyMzcgMS43ODUzMDUgLjQ3MDIzNyAxLjgwOTIxNUMuNDcwMjM3IDEuOTEyODI3IC41NDk5MzggMS45OTI1MjggLjY1MzU0OSAxLjk5MjUyOEMuNzgxMDcxIDEuOTkyNTI4IC44MTI5NTEgMS44OTY4ODcgLjg0NDgzMiAxLjgxNzE4NkwzLjcwNjEwMi01LjY0MjgzOVonLz4KPHBhdGggaWQ9J2czLTEwMScgZD0nTTEuNjAxOTkzLTEuODE3MTg2QzEuNzc3MzM1LTEuODE3MTg2

где $\mu$ является параметром положения, а $\beta$ является параметром масштаба.

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

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

Функция имеет среднее значение $\mu + 0.57721\beta$ и дисперсией

$\frac{\pi^2}{6}\beta^2$ .

Ссылки

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()

../../../_images/numpy-random-gumbel-1_00_00.png

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

>>> 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()

<img alt="../../../_images/numpy-random-gumbel-1_01_00.png" src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAeUAAAEsCAMAAADKCVssAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAACH1BMVEX///8AeAD//v40NDT/AAAAgAAdcKofd7Tx8fEAAAD4+PgedrMBfwXU1NTk5OQeHh4qcKovLy8jdLD7/PsqKiokJCTwBggWixb+/v3f398XFxfc7tw4ODiAgIDq6uozmTNzSW8zbKP6AgN/Q2YMhgzJ5Mn/CQn9AAESemwCfw8Qe18QEBCbMgAbd6Sgz6Du9u7/8fH/1tYed7Bra2uv16/l8uWMxYwqlSr/EhISeWbqCg8GgwYOe1XMzMz/6OgEfhr/trWrq6vCwsKhMUn2BAcaeJ7/+fn0+PTjDhU7Z5zcEhv/Hh6+IjMcd6q83bz4+/gGBgb/SEj/X18JfTgNfE6ysrJCY5WgoKBRqVFarVqEwYTT6dMgkCCaNFCQOli/v78GfST/n54UeXk8njwWeYP/KChKX5FlUXvVFSBGokUHfS+1Jjv/qKgZd5bOGSP/cnL/VFT/goL09PRNUEwLfERttm3/PT0+Pj6IPl/IHSuuKj+Wypb/aWkXeI3/4eFUVFSnp6dISEhlsmXc3Nz/wMBVWYf/NjaKiooMe1rY2NhCQkL/z89jY2N3d3dzuXP/yMheVYH/e3vu7u7/j4//lJRQXIv/iopVWASVN1R8vXz/Li6Tk5MMegGcnJzOzs62trY0ZgbR0dGmXzJdXV2jLQAqbVZhUCxhk0SzonuxQxuUVR+bQxGuKABUaRPwOjLkgHMccKiaMgCUNQCbNAKINlOnLJDXAAAWD0lEQVR42uyd61cURxrG64xvnRZxuV9aiAgIQ9SViwqyEA3xMqAoCPGG4AWVeIiK4AreWWOMJ6uJl9XV6LobTbKb7J69nN2/cKu6uocBEWemq2emm+f5Mj345Tn+pp6uy1tVjEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQ5FpmzlJIt3LMDKOckwXpV06GUV6ahcasvSlnLc04ykvx2gr8/ykogzIEyhAog/JictTbC8pBd2T21LSvun29theUA+yotoZbal81WT9YCMpea2Nu1q/k56GsAylzFCov5byUOyqtBWWPtS0nV1LO+eCD1FFuHuL8VG9nbcvAZoW5AZS9lkX5iy2XU0a5cxXnm22wDT0t4tsAKKeC8p++ZVHKyzyec+0d5XyoOaZli3f0JlD2nvKK/IszlFda6yeeOTIHREQPxv6lfnFkdtop3yrLzd1X9ueUtOUHomc9u7tlHl4UmZ0B72XGUvRebuFv5/PiyOw0j6TKilNHWaQzr5/vr6XNoBwURz3tnB+Z5+8isw+boBwMR4Mimm/PR7OhdL4mDsp+dGSKofHheeczQ5s4r2kG5SA46uH8xF7D1px/E+OrUROUA+BolPOTNuO8YxeG70fmZPZ1UPa/o0HOhzYoyAV7SGpk/PG5pki0odeUg7LvHYlQvm435aMUo6km+c4O3eb8lAnKPnfU3M5LOxXk1UTVT5/13w0rzN3T8t87h4Kd2YuD8gPOXzAL8qM6okqraz3dMdHfTfSy0c7s0l5Q9rWjhhre3mxR/qyE6DvVxw6FGNsqMPdF7EivB2VfO2rh/DaTlKu6iM5mx4ykmlqJbshXc614M4Oynx31ipFSuaScN0Y0ljdrvHxHvJufSMybOR8EZR87quf8sPgwCs4SdVXNmRUZFpg7xOd1q72Dsl8dyWYqV5XlGKrk4FtzX4+JwvdFN7uG13SCsl8dmT32YLhSjKEevT3DGRpX46nbAR5MLYKvK/K6gQ6qrcUYav+e67MjF1rPDXI+WYTlH3qSPSehwpFbzpMdHT+1QprPGWyUyrYQdmPjgZUEk+JMVTB/JTt8VR9cEvA0r+34pvi/Mt/8c5RuZrVaiLanfeOlUcWkuOpp2LA1d4Ayrql9lYcKGQHir1zNCknNxm7QLTaeBdlazwVnn7AeQso65ddw1mxz/TKkZzcFA000krhqgUosxtE/aLZD5mg7BXllb+1v+uvxz6iJjvOEY0bC1FuFK/mh6I33gPKHlE+VFxhf9e+t8Ke3GR9RB0LUg49Jbq7SU2SgbIHlG8Vr2BetWV7cnOaqDu0IGUWGSEaFr+JZlD2gnIMZO2OnMnNCaIJtjBl1iF+CS9UVw2U9Y6kyorZvtz8/PwKbxzZ64mhbqLp91GWqf4fMewqBGWfOZq0agPkcLiPvZeyGFKHR4O5ayrQlAtLebtcZxonOvd+ynJM/QPno6DsK0dyNUr2vbaGqTXyfsqmGE3VDak+OSj7x9GAyt8nRBdYHG1Zjqb+wfkkKPvJkSwM6BXwXhI1zaL8Lt5iNFXSHsRigiBT3qTmvR4S3WVxUZajqedBLOYMMGV5moScsLxBNPwuynN599EY56tA2T+OGtqtwW9jmMKN8VJuorrXASzmDDDl65w/YNaiYj+Ll7IYTf0vgPNfAaZsV/hMET2Mg7Kz7bWkLYDrj8Gl3Kyq9ZqIRkLxUzZW0+vg1X8Fl3JLtEjkCUuAcsHu/6qkB2U/OFplzWJFWom2JkLZOFOyPHCRHVjKg2pEdIborJGQsseeBy6yA0v5hSrV20P0KDHKxqOfAjfLGVTK5pBV9jFdR10FCVLO/veroK0yB5VyrVpCnCA6ZiSqjudBq/ILKmW7fqCb6nYlTDn0A+c/g7Imqb0VFV8UF23T7ahQ7ZO4T3Q1YciG2fGKv+qd3S0H5WSl9lZ8s5JtyS3U7MiuH3gWu6MibrHQvzj/BZR1SVLeV8HYkuOaHdn1A1NEu5KgzH7h/LkJyvooV+wTn787pL7rqsfuVfUDjURdRjKUO0VkD4OyR5Q17a0w7fqBDqLtSVFmP3P+9wgoe5PYutryTP3AGSMpnRSRXblAUQkoJ0j5R/29L7t+gHUT5SVHecNy/mp3ASjrGUmVFbOKe8VFx/U6qldTlE1EU0aSusT5PytBOZMdjarlhqdEz5KlfFpEduzcKChnmqNmu9qjj+h+spQ3rOGvdlSCcuY6susHImFqDSVL2fia859iGjMoZ5qjU6oK844s6zPcRDZVgnKGOjJlwZd8uCD3wBluIrtkpjGDcoY5arFvB3spa4GSpmzsF5E9MwsOyhnmyA7srdbGmaQhZ9+Ukd1VAMoZ6ajcDuxhoscuKBvZa/jykmhjBuXMcuQE9jjRHTeUZWT/mtoKQDkTHakKXWa2ztmanrBEZL+Ork+DckY5Krf3LD4kGmeuKGd/zHkbteWBcuY5cgJ7guipO8rGR5z/jegKKGeeIzuwZZnItEvKIrJP1FFJHihnmiMnsBuJRphLyjKyt9uNGZQzyZET2OfUiTGuKMvI/mu1asygnEmOnMC+oe4Pckf5vIjsa6puH5QzyJET2KHuuI6ZeJ9EZB9UjRmUM8iRE9hN1pmMrimv5XztdnkzJCi71vdL8osO6XHkBPYTomcaKMvI3lVNO6pA2a3MdW/Yin1LdThyAluWiTRpoGyc4HzndnlFESi7pryNvSlbpsPRERXYoUiYukM6KMvIrpKNGZTd6vi63HWqitNtPbYT2B3O6U9uKe8UkS3vhzwKyi5VeG8b21L2B/nocm9FNLCtMhEdlGVkn6/aQdXToOxOW4qYs7vCZVs+4lwINULUqIeyjGzjO19I5+j/5uH198sB6fS+3gI6U7a9d0H3q12x+e58p+p033l7sF51D1c+8vZ5o9yE6rN0/P3w8839k5V85q/L/2x47zD9C33oP3q12x+e5/l3sH9y0d593c7Y4zB1b28x/3b2x78t82N8t2/dD23sB17q/P+N1y+e/42+2N9z/T63N51rR2uB1z2P9J3Y479i521y+o9J3nV/e3+j/h3T6wP93R7c8sJ96eO8/V16O47J34i+c51669J7/9m9q2O88w5u7b479rX9v6V6P3p1b74j575z8f37n+4k/cYvG/8719Y1e47l/jK86k20/Y3v5x9T1K9b6p+2F/N6/9e7l/a+9jD1F4p15n/m/3W3/e4x9b5+9H4x7L1P33+p9a6y7Y3m7+d116f39X83+uK6/o9k93fWq37z/H27/x6Z/h4/c2865N9F+9J+x/H1oT286r2r/v1/o1sF3j6r5fT2l/b4h7s5P33747H1i9d9+8R+v+e/f6+82/6h/X/q1t/sX9yL/u7/q1Z2W8uR3vQ1p95n/a/w+uS/3/b/o5X4c4v7R5t7u7m+3yN96q7+x/t1t/aT/97/e+v5k7/Q/p5/2a68p75x+b5f/4n6d3X5/7995T379H1v6W/7/e6X4e2+R8e9O38927b3g9H5/c07d/tK1t5H618/T1p1jX6J1775P5q+3t955Z9d/W5t/L7d1X/7nW4b/sL0t9+R+8/u6/zH/3P/9b96f74t/v/t5+x95f6r/o7a8o6J/V/7664f7/i5/e/iO3/Z75L/994t0+k54/7n/q77/e4/27947f3+5O051X7o7a/X/4K7jY5N/9p/+N/15l/T1/7L+p9/82O+X1X8J51/vN6/Yf88r/w/6L9/L/99v4f/l+P/9/b39/yD24y/wK/v4d/j/v3/uY7e87/vJ2s/+a/4L/7/+80/99X/9/o1a706H0H3v8t3/41/4H/b7e/9v9+/uX7a39P+z9vH69/t7/9T2L3r7/o+/d73737zZ3/3i/P6f/v7/eW8/9r5/V/x43+7+847+7/vZ9P9/Y/9e/c/vH4f+u/zO5/25/t/7+6X1P6/vW+99w3u++z+/v//iW/Y7/L6b9/6+v++9+/v//+3r/f8+/vv/e77/9777b+/e17b62e66+/uXv8/u+/zD3/+/c85t/z/+66eP3//7+7b3a499f31r1w29/57+z+n9/+n2a2/v/4z/f78m/1/8t/+H/9H9z2z88x/95367z651/9l9X79T3f/491j/y/z2/+N/76/v56L/45H/++v5L//b0/L+7//7d/m5t/j//++v+v/e//8u+/3j//8196//f7f/++/wA540n7d+4f+/f89D/76f+6+v/3n27v+/6X//7X/j74D/p2l+0eN/b7a++/vZ+585T83u/9z73o9799f+/v/v7z7jO/g85/9j+v+7/30f/78v/f99/wAAAAgAAdcKofd7Tx8fEAAAD4+PgedrMBfwXU1NTk5OQeHh4qcKovLy8jdLD7/PsqKiokJCTwBggWixb+/v3f398XFxfc7tw4ODiAgIDq6uozmTNzSW8zbKP6AgN/Q2YMhgzJ5Mn/CQn9AAESemwCfw8Qe18QEBCbMgAbd6Sgz6Du9u7/8fH/1tYed7Bra2uv16/l8uWMxYwqlSr/EhISeWbqCg8GgwYOe1XMzMz/6OgEfhr/trWrq6vCwsKhMUn2BAcaeJ7/+fn0+PTjDhU7Z5zcEhv/Hh6+IjMcd6q83bz4+/gGBgb/SEj/X18JfTgNfE6ysrJCY5WgoKBRqVFarVqEwYTT6dMgkCCaNFCQOli/v78GfST/n54UeXk8njwWeYP/KChKX5FlUXvVFSBGokUHfS+1Jjv/qKgZd5bOGSP/cnL/VFT/goL09PRNUEwLfERttm3/PT0+Pj6IPl/IHSuuKj+Wypb/aWkXeI3/4eFUVFSnp6dISEhlsmXc3Nz/wMBVWYf/NjaKiooMe1rY2NhCQkL/z89jY2N3d3dzuXP/yMheVYH/e3vu7u7/j4//lJRQXIv/iopVWASVN1R8vXz/Li6Tk5MMegGcnJzOzs62trY0ZgbR0dGmXzJdXV2jLQAqbVZhUCxhk0SzonuxQxuUVR+bQxGuKABUaRPwOjLkgHMccKiaMgCUNQCbNAKINlOnLJDXAAAWD0lEQVR42uyd61cURxrG64xvnRZxuV9aiAgIQ9SViwqyEA3xMqAoCPGG4AWVeIiK4AreWWOMJ6uJl9XV6LobTbKb7J69nN2/cKu6uocBEWemq2emm+f5Mj345Tn+pp6uy1tVjEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQ5FpmzlJIt3LMDKOckwXpV06GUV6ahcasvSlnLc04ykvx2gr8/ykogzIEyhAog/JictTbC8pBd2T21LSvun29theUA+yotoZbal81WT9YCMpea2Nu1q/k56GsAylzFCov5byUOyqtBWWPtS0nV1LO+eCD1FFuHuL8VG9nbcvAZoW5AZS9lkX5iy2XU0a5cxXnm22wDT0t4tsAKKeC8p++ZVHKyzyec+0d5XyoOaZli3f0JlD2nvKK/IszlFda6yeeOTIHREQPxv6lfnFkdtop3yrLzd1X9ueUtOUHomc9u7tlHl4UmZ0B72XGUvRebuFv5/PiyOw0j6TKilNHWaQzr5/vr6XNoBwURz3tnB+Z5+8isw+boBwMR4Mimm/PR7OhdL4mDsp+dGSKofHheeczQ5s4r2kG5SA46uH8xF7D1px/E+OrUROUA+BolPOTNuO8YxeG70fmZPZ1UPa/o0HOhzYoyAV7SGpk/PG5pki0odeUg7LvHYlQvm435aMUo6km+c4O3eb8lAnKPnfU3M5LOxXk1UTVT5/13w0rzN3T8t87h4Kd2YuD8gPOXzAL8qM6okqraz3dMdHfTfSy0c7s0l5Q9rWjhhre3mxR/qyE6DvVxw6FGNsqMPdF7EivB2VfO2rh/DaTlKu6iM5mx4ykmlqJbshXc614M4Oynx31ipFSuaScN0Y0ljdrvHxHvJufSMybOR8EZR87quf8sPgwCs4SdVXNmRUZFpg7xOd1q72Dsl8dyWYqV5XlGKrk4FtzX4+JwvdFN7uG13SCsl8dmT32YLhSjKEevT3DGRpX46nbAR5MLYKPK6gQ6qrcUYav+e67MjF1rPDXI+WYTlH3qSPSehwpFbzpMdHT+1QprPGWyUyrYQdmPjgZUEk+JMVTB/JTt8VR9cEvA0r+34pvi/Mt/8c5RuZrVaiLanfeOlUcWkuOpp2LA1d4Ayrql9lYcKGQHir1zNCknNxm7QLTaeBdlazwVnn7AeQso65ddw1mxz/TKkZzcFA000krhqgUosxtE/aLZD5mg7BXllb+1v+uvxz6iJjvOEY0bC1FuFK/mh6I33gPKHlE+VFxhf9e+t8Ke3GR9RB0LUg49Jbq7SU2SgbIHlG8Vr2BetWV7cnOaqDu0IGUWGSEaFr+JZlD2gnIMZO2OnMnNCaIJtjBl1iF+CS9UVw2U9Y6kyorZvtz8/PwKbxzZ64mhbqLp91GWqf4fMewqBGWfOZq0agPkcLiPvZeyGFKHR4O5ayrQlAtLebtcZxonOvd+ynJM/QPno6DsK0dyNUr2vbaGqTXyfsqmGE3VDak+OSj7x9GAyt8nRBdYHG1Zjqb+wfkkKPvJkSwM6BXwXhI1zaL8Lt5iNFXSHsRigiBT3qTmvR4S3WVxUZajqedBLOYMMGV5moScsLxBNPwuynN599EY56tA2T+OGtqtwW9jmMKN8VJuorrXASzmDDDl65w/YNaiYj+Ll7IYTf0vgPNfAaZsV/hMET2Mg7Kz7bWkLYDrj8Gl3Kyq9ZqIRkLxUzZW0+vg1X8Fl3JLtEjkCUuAcsHu/6qkB2U/OFplzWJFWom2JkLZOFOyPHCRHVjKg2pEdIborJGQsseeBy6yA0v5hSrV20P0KDHKxqOfAjfLGVTK5pBV9jFdR10FCVLO/veroK0yB5VyrVpCnCA6ZiSqjudBq/ILKmW7fqCb6nYlTDn0A+c/g7Imqb0VFV8UF23T7ahQ7ZO4T3Q1YciG2fGKv+qd3S0H5WSl9lZ8s5JtyS3U7MiuH3gWu6MibrHQvzj/BZR1SVLeV8HYkuOaHdn1A1NEu5KgzH7h/LkJyvooV+wTn787pL7rqsfuVfUDjURdRjKUO0VkD4OyR5Q17a0w7fqBDqLtSVFmP3P+9wgoe5PYutryTP3AGSMpnRSRXblAUQkoJ0j5R/29L7t+gHUT5SVHecNy/mp3ASjrGUmVFbOKe8VFx/U6qldTlE1EU0aSusT5PytBOZMdjarlhqdEz5KlfFpEduzcKChnmqNmu9qjj+h+spQ3rOGvdlSCcuY6susHImFqDSVL2fia859iGjMoZ5qjU6oK844s6zPcRDZVgnKGOjJlwZd8uCD3wBluIrtkpjGDcoY5arFvB3spa4GSpmzsF5E9MwsOyhnmyA7srdbGmaQhZ9+Ukd1VAMoZ6ajcDuxhoscuKBvZa/jykmhjBuXMcuQE9jjRHTeUZWT/mtoKQDkTHakKXWa2ztmanrBEZL+Ork+DckY5Krf3LD4kGmeuKGd/zHkbteWBcuY5cgJ7guipO8rGR5z/jegKKGeeIzuwZZnItEvKIrJP1FFJHihnmiMnsBuJRphLyjKyt9uNGZQzyZET2OfUiTGuKMvI/mu1asygnEmOnMC+oe4Pckf5vIjsa6puH5QzyJET2KHuuI6ZeJ9EZB9UjRmUM8iRE9hN1pmMrimv5XztdnkzJCi71vdL8osO6XHkBPYTomcaKMvI3lVNO6pA2a3MdW/Yin1LdThyAluWiTRpoGyc4HzndnlFESi7pryNvSlbpsPRERXYoUiYukM6KMvIrpKNGZTd6vi63HWqitNtPbYT2B3O6U9uKe8UkS3vhzwKyi5VeG8b21L2B/nocm9FNLCtMhEdlGVkn6/aQdXToOxOW4qYs7vCZVs+4lwINULUqIeyjGzjO3UVNygnr1DFut+zL9flaHDkBPa0fb6qBsoisj82RGOmaVB2p1tF+fm3NDiKBvawvIxED2UrsmVj7gflzHAUDexx51BGDZR/w/lHRl6JMzID5XQ7cgI70kqtpi7K60VkZxtXnLk0UE6vIzMa2HaZiB7Kxucysgu67dUPUE6zo2hgT9jHnuuh/KGMbOMc0UsTlNPvaLNzfcxUTI/Y0BTZobsxPx1QTpujaGA7ZSKaKMvIvmnIXVetjaCcZkchJ7BDHbFzGIaWyN5vWB33CVBOsyN5cae6b/VCbEdJA+W9y0VkqwsCt4Jyeh3ZF3ey2OlNPZStyFa/nn5QTq+jSeceqJjpTU2Urcg2jKoSqvvMnyVgQaHca1/c6eyC00lZRPaaDeLzCtGebFBOo6NNnA+op35rF5xOysannJ+WV1p0WUeqg3LaHMlz3FSpfBvtKNBLOfsk519bV1oQ7S4A5bQ5ss5xs5AcFLFqaG7LTmRnjxFVgnLaHFn3wNn3QdFq3ZSdyDb+SNSWB8rpcrRKnjJhAdlDdFA7ZRHZl6yHs0THQDlNjtS1QtalXzuoLVs7ZRHZy/fKB1mcPQ3K6XGkrhWSGL4iumbo1yXOP7QejvpyaiQTKC/bWFz0oxtH9rVCksIxokoPKN/k/HP7ts9ZAzVQjlsbN5osx42jHnWtkKRwlajKA8pyK/N5p3c3EgHlhHVx3UWXjuxbKmRTq6bdhhdaa9USyBf/2Ky5NVCOU2+Kv11y2c3eik51SwVT8xbbPaG83h4yWwNyug/KiWpL1iH2pqxCPia3t6Je3VIhKW9P+g6DOIbMJ9XTY6K7JignqIqsQnd7K+xDzyXl3VSd5w3l05x/qp4iI9ErnUE5ft37nq0oS35vhXPouaBcRTTmDWR5xi5frw6FukMUngblBPXl5aL8T5J3ZB96LilXJnN5Z/xHyay1j/66QdQXAuVUOjKjtZuGcY3oK68o77RqOS7Kjd3RPZWgnBpHM6VAatXRK8qyMOi0/Si68m1VPlpp9j/lSc6v25StVUfPFF2yUKsW10A5dY7sO3gtyquTuqI1gf6XWrJQNWD0CJRT5mimFIipVUfvtN9ZsjCsfl5XHiinytFhNbsppVYdvZM8Fyo6sX3VT4eN+J2yc6WjlEerjnN2rNvatYPqHoJyahxZpUC2+PFp1fKswW+kK0csIKKfEkX3hzP/bO/OnNq47gL/C2y5GBJAAWVw2kggYUsxlD4cxwkZykEDhDK1bm2DAsiwD5jJF1PiIwRDbpL7IeOpJ00knmWnzS2faP7Dv7a6kldC5qz0Uvp8fPJaQHl/48P3qnbscy8qsOh7f5ccvTvnp1d5y4+YlOW6ZDJa7QisHbUqtOkZtGZkLP/DZIksjYFnJiK6Lur2KrTpGLVkEIo/IwM3sBsuKR1TNiEqocquOUVtGXJFH2xj7TWBZ6YgehFYQWEVXHeNsGeE3gXlC5QMsKxdRaz1T3xIe2mAcVFoyvfxEpHiw7LRT6NaDZeUi+iPDPGBFs1EzilsOn7KIfFOnDywrGdHoaYapZkULCEfKW54LbxkJdwY8pWBZuYgs5xnmD+HDE6VObDYpb5kuWbhEj+mo+Z4BLKdD7alPJFjuovu9VBxHhbaMHIgfu83Rx+/AciIuPy2RYPkCt31AVLBfq2G5pSg6memfl+0ILKeksKZSiuV5bjUqUrB7TawWyUyPTg31g+VUHH5EguVMduoKk5sqF+w4yWwK6n6znw4sX31hCVnOZNf9df6iQJGCvaeOZZrMU1FPWHv1fqhGB5b3q0pKTlX9kGEuVzPM6dGwZfUKNp/MLVHznnuNodtXgeVkZP65/EC4C1x4dkKlgs3vDIpOZrpuoetrseaqZTq52RqxrGLBjpPMLLuKdb1DO1dnRejkJgpbLrWpV7DjJnPefazny6TnqGVucjNiWdWCHS+ZUd4Exg27YDmrEXGTmxHLI2oWbD6Z+6ItI8sKxt0vwXIWI+InN8OWacE2qGnZFZvMJJBr40RzO1jOXkQXwmejkAYFm3AQm8z0by2I8ZCbBcvZimg+stVeKNhH6lp2He9mE81+jD1WsJyliL4Q7dzkC/aQgdU+mVmrB2N/G1jOSkSjXaErnvOW60Qbo9VMZu+xZ91DGD/vB8vZiIiMlectIsvqF+xEycxO92I83A6WZUdkqb4kOk9BLyaiQcFOlMysj2ju3hUN8nSxxTMHLd8VHY0SLiaiesHmz1nESWbWSrpgttjjWmA504jOM8yNz8WWScH2aWDZxTBNcZKZLSXx4BkDWJYTUevpyCiKs6xNweaTeSPe86ZBonnbBJZlRPR1eG5TsEwK9gctJNNPZtHZOPF6s70R46AVLEuOiAyV6y9GWQ5qU7D548xNLXG/cseJsWcaLEuMiA6V/xL1jLURewzaWDZMMkwg/vemXW3zHliWFhEZKp8ti3rGrlXBpiuQTaKDtTHzI7SrPWMCyxIiih4qc2hXsPnzzEVL8b/EdbWXp8GyQOF6hfFpZToRxQyVKf3aFWxhBmy2I/6XTDM2ks5bJrDMW/7Bgj6+SCeimKEy5b0qJx0