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

时间:2012-01-28 18:24:26

标签: python social-networking

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

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

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

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

提前致谢。

2 个答案:

答案 0 :(得分:1)

我知道这可能会迟到,但你可以查看最近发布的`powerlaw'包。它有iPython notebook个工作示例和数据。

答案 1 :(得分:0)

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