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

Question

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

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

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

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

提前谢谢。

Answer 1

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

Answer 2

调查一下scipy.optimize。如果您必须使用最大似然，请定义一个似然函数，并使用fmin()找到-likelihood的最小值(=似然的最大值)。如果其他方法是可接受的，请考虑curve_fit() (它使用最小二乘拟合)。