Расчет логарифмической вероятности распределений в Python

Question

Каков простой способ для вычисления логарифмической вероятности любого распределения для данных?

Answer 1

Решение по OP.

Python имеет 82 стандартных дистрибутива, которые можно найти здесь и в scipy.stats.distributions

Предположим, вы найдетепараметры, такие что функция плотности вероятности (pdf) соответствует данным следующим образом:

dist = getattr(stats.stats, 'distribution name')
params = dist.fit(data)

Тогда, поскольку это стандартное распределение, включенное в библиотеку SciPy, pdf и logpdf могут бытьнайти и использовать очень легко следующим образом:

LLH = dist.logpdf(data,*params).sum()

Обратите внимание, что это соответствует функции правдоподобия, определенной здесь .