Подобрать степенной закон к эмпирическим данным в Python ⇐ Python

[Anonymous]

Я экспериментирую с сопоставлением степенного закона с эмпирическими данными с помощью модуля степенного закона. Я создал следующие данные, которые соответствуют степенному закону распределения показателя 2:

Код: Выделить всё

x = range(1,1000)
y = []

for i in x:
y.append(i**(-2))

Я ожидаю, что показатель подобранной степени будет равен 2. Однако полученный показатель степени сильно отклоняется от теоретического значения:

Код: Выделить всё

    fitted_pl = powerlaw.Fit(y)

fitted_pl.alpha
Out[115]: 1.4017584065981563

Не могли бы вы объяснить, почему это происходит, или указать, что я здесь сделал не так?

Спасибо за добрый ответ!

Подробнее здесь: https://stackoverflow.com/questions/181 ... -in-python