Pythonを使用してべき乗則分布の指数を推定する方法は？

Question

私は現在、いくつかのネットワーク解析を行っています。平均クラ​​スタリング係数対ノード次数のべき乗則分布の指数を推定したいと思います。好ましい方法は、最尤推定（MLE）

Iは同様の符号がここに実装されている知っている：http://tuvalu.santafe.edu/~aaronc/powerlaws/

しかしながら、これらは実装しているコードは、x対曲線のPr（> = X X）をプロットする特異的ですxは与えられたシーケンス内のノードです。言い換えれば、Pr（X> = x）はxに等しいかそれより大きいXにノードが存在する確率である。

私は数日間グーグルで行きましたが、私はまだこれを行うツールを見つけることができませんでした。 2つの値のシーケンスを与えられたべき乗則分布の指数を推定する。

ありがとうございます。

Answer 1

scipy.optimizeを参照してください。 Maximum Likelihoodを使用する必要がある場合は、尤度関数を定義し、fmin（）を使用して最小値（=尤度の最大値）を求めます。他の方法が受け入れられる場合は、curve_fit（）（最小二乗適合を使用する）を検討してください。

Answer 2

これは遅れるかもしれないが、最近リリースされた 'powerlaw'パッケージを調べることができます。それは、実用的な例と数字でiPython notebookを持っています。