Установка обобщенного распределения Парето (GPD) ⇐ Python

[Anonymous]

Я пытаюсь установить GPD в Python. Стандартный GenPareto.fit от Scipy Works, однако я хотел бы получить ошибку параметров (чувствительность) и, к сожалению, она не предоставлена. Когда я стараюсь установить GPD самостоятельно с вероятностью журнала, он дает очень неудовлетворительные результаты (с Differial_evolution). Не могли бы вы помочь мне с написанием правильной функции подгонки для GPD или какой-то хитрости, чтобы получить чувствительность для параметров?from scipy.stats import genpareto
shape_fit, loc_fit, scale_fit = genpareto.fit(exceedances, floc=0)

def gpd_neg_log_likelihood(params, data):
xi, beta = params[0], params[1]
if beta

Подробнее здесь: https://stackoverflow.com/questions/794 ... bution-gpd