Python - Подгонка распределения GEV из этих значений
Я очень новичок в Python, и я огляделся по интернету, но не смог найти ничего логичного, что могло бы помочь мне с моей проблемой.
У меня есть значения осадков на графике, и теперь мне нужно установить распределение GEV из этих значений на графике. Каждое значение равно максимальному значению года, начиная с 1974 по 2017 год (так что всего 43 значения).
Это значения:
max_precip = [9.4, 38.0, 12.5, 35.3, 17.6, 12.9, 12.4, 19.6, 15.0, 13.2, 12.3, 16.9, 16.9, 29.4, 13.6, 11.1, 8.0, 16.6, 12.0, 13.1, 9.1, 9.7, 21.0, 11.2, 14.4, 18.8, 14.0, 19.9, 12.4, 10.8, 21.6, 15.4, 17.4, 14.8, 22.7, 11.5, 10.5, 11.8, 12.4, 16.6, 11.7, 12.9, 17.8]
Я обнаружил, что мне нужно использовать gev.fit, поэтому я подумал использовать следующее:
t = np.linspace(1,43,43)
fit = gev.fit(max_precip,loc=3)
pdf = gev.pdf(t, *fit)
plt.plot(t,pdf)
plt.plot(t, max_precip, "o")
Но это только печатает точки max_precip в графе, а не в распределении GEV.
Кто-нибудь может мне помочь? Извините, если этот вопрос уже задан, я не смог найти ничего подобного.
Я использовал эти импортные товары:
import csv
import matplotlib.pyplot as plt
import numpy as np
from dateutil.rrule import rrule, YEARLY
import datetime
from matplotlib.dates import DateFormatter
from scipy.stats import genextreme as gev
from scipy.stats import genpareto as gpd
from scipy.optimize import minimize
1 ответ
Я попытался подогнать ваши данные
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import genextreme as gev
def main(rvs):
shape, loc, scale = gev.fit(rvs)
return shape, loc, scale
if __name__ == '__main__':
rvs = [9.4, 38.0, 12.5, 35.3, 17.6, 12.9, 12.4, 19.6, 15.0, 13.2, 12.3, 16.9, 16.9, 29.4, 13.6, 11.1, 8.0, 16.6, 12.0, 13.1, 9.1, 9.7, 21.0, 11.2, 14.4, 18.8, 14.0, 19.9, 12.4, 10.8, 21.6, 15.4, 17.4, 14.8, 22.7, 11.5, 10.5, 11.8, 12.4, 16.6, 11.7, 12.9, 17.8]
shape, loc, scale = main(rvs)
print(shape)
print(loc)
print(scale)
l = loc + scale / shape
xx = np.linspace(l+0.00001, l+0.00001+35, num=71)
yy = gev.pdf(xx, shape, loc, scale)
hist, bins = np.histogram(rvs, bins=12, range=(-0.5, 23.5), density=True)
plt.bar(bins[:-1], hist, width = 2, align='edge')
plt.plot(xx, yy, 'ro')
plt.show()
но я вернулся
-0.21989526255575445
12.749780017954315
3.449061347316184
для shape, loc и scale. Если вы посмотрите на распределение GEV, как определено в scipy, когда форма отрицательная, допустимый интервал [loc + scale/shape... +infinity]. Я вычислил последнее значение и равен
-2.935417290135696
должно сработать...
Python3, Anaconda, scipy 1.1, Windows 10 64bit
ОБНОВИТЬ
Хорошо, я обновил код и добавил график, выглядит несколько разумным. Это то, что вы ищете? В принципе, трюк заключается в гистограмме и компоновке ячеек плотности, перекрывающихся с PDF
Ещё вопросы
- 0Добавить фиктивные данные в mySql Result
- 1Нет Access-Control-Allow-Origin с использованием настройки OneSignal
- 1Лучший способ определить, была ли вызвана функция с правильным количеством аргументов в JavaScript
- 0javascript / jquery получить значение переменной [дубликата]
- 0Set Session в PHP электронной почте не работает
- 0eJabberd с PHP Extauth
- 0какой тип содержимого установить, если содержимое страницы содержит смесь php, html, css и javascript / jquery… и расширение файла было php?
- 0Moodle 2.7 Нужно уметь создавать пользовательское меню для разных групп пользователей
- 1Manifest Merger Ошибка в проекте
- 1React Native - отключить воспроизведение звука в фоновом режиме через определенное время
- 0Jquery / CSS подменю
- 1Переходит на URL, но не получает ответа
- 0Перенести значение из выпадающего в текстовое поле
- 0Как снять угловой наблюдатель с направляющего элемента?
- 0Ошибка java-скрипта при перемещении ссылок jquery на главную страницу asp.net
- 0C ++ ошибка Valgrind
- 0Сортировка списка с помощью jquery UI
- 0Кнопка браузера не запускает маршруты в Backbone
- 0Параметр выпадающего не выбран - angularjs
- 0AngularJS: почему поиск не работает должным образом?
- 0Результат запроса MySQL: что означает показанное время?
- 1Как изменить значение логического значения при нажатии кнопки JButton
- 1Частичное переопределение метода: как?
- 0new-выражение и delete-выражение для константной ссылки и константного указателя
- 0Как выполнить функцию PHP crypt в Nodejs
- 0Программа перестает работать при вызове функции.
- 1Сглаживание неравномерного массива - ошибка исключения вне границ
- 0обеспечение единого #defines среди всех объектных файлов
- 0Google Map API не загружается, когда скрипт на той же HTML-странице
- 1Не удается создать NotificationCompat.Builder в Android O
- 0Что означает выбор * во вкладке «Администрирование-клиент» рабочей среды
- 1Блокнот Jupyter не будет импортировать matplotlib, хотя он, кажется, установлен
- 1в Цезии динамический вращающийся компас сбрасывается назад до 0, когда он проходит 360
- 1Использование нескольких регулярных выражений в C # web crawler [duplicate]
- 03 столбца делятся не по горизонтали
- 1Получить все возможные строки в печатном ASCII
- 0Как обновить количество базы данных из действия php
- 1Котлин как проблема мультиплатформенности
- 0Проблема отправки писем в php [дубликаты]
- 1Некоторые пользовательские маркеры не отображаются на Android-native-картах реакции
- 1Фоновые геозоны действительно работают на Android 8+?
- 0Создать стену для поста из блоков
- 1Javascript, имя функции с интерполяцией строк [дубликаты]
- 1Можно ли передать HTML в таблицу прокрутки в Gojs?
- 0ошибка: у класса 'Software' нет поля с именем 'ptr'
- 0назначить адрес узлу * в списке ссылок?
- 0Ошибка типа Javascript: oauth не является конструктором
- 1double.ToString («G»), но всегда включать точку?
- 1pandas.read_csv () интерпретирует TRUE как логическое значение, мне нужна строка
- 0почему при использовании сервиса не отображается предупреждение?
genextreme, вы увидите, что параметры формы имеют противоположные знаки. Таким образом, поддержка подгонки, которую вы вычислили: (loc + scale / shape, + ∞) = (-2.9354, + ∞).