2016-07-22 12 views
1

私はデータのサンプルを持っています。サンプルサイズは約500,000です。私がしようとすると、しかし、問題が来る法則分布フィッティング

pl_rg <- conpl$new(a) 
estimate_xmin(pl_rg, xmax = 100) 

:私は現在 だから、これはそのために私のコードです

R.

poweRlawパッケージを使用して、べき乗分布を試料に合うようにしようとしています x<-estimate_xmin(pl_rg,xmax = 100)のコードを実行します。私はこのコードが非常に非常に時間がかかることがわかった。私は今までにこのコードを最大5時間実行しており、まだ実行中です。

このように、法則やその他の方法を加速する方法はありますか?

+0

再現可能な例http://stackoverflow.com/help/mcveを提供してください。スピードアップしたい作業コードがある場合は、質問をhttp://codereview.stackexchange.comに移動することを検討してください –

答えて

1

estimate_xminが遅い理由は、潜在的な可能性のあるデータの可能なすべての値を試しているということです。xminの値です。値を指定することで、機能を大幅に高速化することができます。

estimate_xmin(pl_rg, xmins = 1:10) 

詳細については、ヘルプページを参照してください。

+0

ありがとうございます。 xminの初期値はどのように指定できますか?データの確率密度関数に従っていくつかの値を選択するだけで、プロットがログログ軸で線形になるように開始できますか? –

+0

上記のコードでは、 'xmins'引数を使ってxminの値を選択しています。 – csgillespie