Как оценить показатель степени степенного распределения с использованием Python?

Question

В настоящее время я делаю некоторый анализ сети и хочу оценить показатель степени степенного распределения среднего коэффициента кластеризации в зависимости от степени узла.Предпочтительным методом является оценка максимального правдоподобия (MLE)

Я знаю, что подобный код был реализован здесь: http://tuvalu.santafe.edu/~aaronc/powerlaws/

Однако код, который они реализовали, специально предназначен для построения кривой Pr (X).> = x) против x, где x - узел в данной последовательности.Другими словами, Pr (X> = x) - это вероятность того, что в X есть узел, который равен или больше, чем x.

У меня есть Google в течение нескольких дней, и мне до сих пор не удалось найти какой-либоинструменты, которые могли бы сделать это, например, оценка степени степенного распределения с учетом двух последовательностей значений.

Заранее спасибо.

Answer 1

Я понимаю, что это может быть поздно, но вы можете взглянуть на пакет ` powerlaw ', который был выпущен недавно.Он имеет ноутбук iPython с рабочими примерами и рисунками.

Answer 2

Посмотрите на scipy.optimize . Если вы должны использовать максимальное правдоподобие, определите функцию правдоподобия и используйте fmin (), чтобы найти минимум-правдоподобия (= максимум правдоподобия). Если другие методы приемлемы, рассмотрите функцию curve_fit () (которая использует подгонку методом наименьших квадратов).