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

Question

6

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

我知道类似的代码已经在这里实现了：http://tuvalu.santafe.edu/~aaronc/powerlaws/ 然而，他们实现的代码专门用于绘制曲线Pr(X >= x) vs x，其中x是给定序列中的一个节点。换句话说，Pr(X >= x)是X中有一个等于或大于x的节点的概率。

我已经谷歌了几天，仍然没有找到任何可以完成此操作的工具，例如给定两个值序列估计幂律分布的指数。

提前感谢。

- Mike

2个回答

0

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

- Jan

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- Dexter · Accepted Answer

我知道这可能有点晚，但你可以查看最近发布的“powerlaw”包。它带有一个iPython笔记本，其中包含工作示例和图表。请参考powerlaw和iPython notebook。