如何使用 Python 估计幂律分布的指数？

Question

我目前正在做一些网络分析，我想估计平均聚类系数与节点度的幂律分布的指数。首选方法是最大似然估计（MLE）

我知道这里已经实现了类似的代码：http : //tuvalu.santafe.edu/~aaronc/powerlaws/

但是，他们实现的代码专门用于绘制曲线 Pr(X >= x) vs x，其中 x 是给定序列中的节点。换句话说，Pr(X >= x) 是 X 中存在等于或大于 x 的节点的概率。

我已经用谷歌搜索了好几天，但我仍然没有找到任何可以做到这一点的工具，例如在给定两个值序列的情况下估计幂律分布的指数。

提前致谢。

Answer 1

查看scipy.optimize。如果必须使用最大似然，请定义似然函数并使用 fmin() 查找 - 似然的最小值（=似然的最大值）。如果可以接受其他方法，请考虑 curve_fit() （使用最小二乘拟合）。