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

Question

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

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

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

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

提前谢谢。在

Answer 1

查看scipy.optimize。如果必须使用最大似然函数，请定义一个似然函数并使用fmin（）来查找最小似然值（=似然最大值）。如果可以接受其他方法，请考虑曲线拟合（）（它使用最小二乘拟合）。在

Answer 2

我知道这可能晚了，但是你可以看看最近发布的`powerlaw'包。它有一个iPython notebook，带有工作示例和图形。在