我目前正在做一些网络分析,我想估计平均聚类系数与节点度的幂律分布的指数。首选的方法是最大似然估计(MLE)。
我知道这里已经实现了类似的代码:http://tuvalu.santafe.edu/~aaronc/powerlaws/
然而,他们实现的代码专门用于绘制曲线Pr(X >= x) vs x,其中x是给定序列中的节点。换句话说,Pr(X >= x)是x中存在等于或大于x的节点的概率。
我已经用谷歌搜索了几天,但我仍然没有找到任何可以做到这一点的工具,例如,在给定两个值序列的情况下,估计幂律分布的指数。
提前谢谢。
发布于 2012-01-29 06:24:17
调查一下scipy.optimize。如果您必须使用最大似然,请定义一个似然函数,并使用fmin()找到-likelihood的最小值(=似然的最大值)。如果其他方法是可接受的,请考虑curve_fit() (它使用最小二乘拟合)。
https://stackoverflow.com/questions/9047667
复制相似问题