Leistungsgesetz Verteilung passend

Ich habe eine Stichprobe von Daten, mit der Stichprobengröße etwa 500.000. Ich bin derzeit versucht, die Probe mit Potenzgesetz Verteilung zu passen, in R. das poweRlaw Paket mitLeistungsgesetz Verteilung passend

Das ist also mein Code zu diesem Zweck

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

jedoch, da kommt das Problem, wenn ich versuche, Führen Sie den Code x<-estimate_xmin(pl_rg,xmax = 100) aus. Ich fand diesen Code sehr zeitaufwendig. Ich habe diesen Code jetzt für bis zu 5 Stunden ausgeführt und es läuft noch.

Also jede Möglichkeit, den Prozess oder andere Möglichkeiten zur Anpassung der Potenzgesetz zu beschleunigen?

Quelle

2016-07-22 Burry Xie

Bitte stellen Sie ein reproduzierbares Beispiel http://StackOverflow.com/Help/Mcve bereit und für Fragen, bei denen Sie Arbeitscode haben, den Sie beschleunigen möchten, ziehen Sie in Betracht, Ihre Frage auf http://codereview.stackexchange.com zu verschieben –

Der Grund, dass estimate_xmin langsam ist, dass es jeden möglichen Wert Ihrer Daten für einen potentiellen xmin Wert versucht. Sie können die Funktion erheblich beschleunigen, indem Sie Werte angeben, z.

estimate_xmin(pl_rg, xmins = 1:10)

Weitere Informationen finden Sie auf der Hilfeseite.

Quelle

2016-07-22 10:35:45 csgillespie

Danke. aber wie kann ich die initialen Werte für xmin angeben? Kann ich einfach mehrere Werte entsprechend der Wahrscheinlichkeitsdichtefunktion der Daten wählen, wenn der Plot in der Log-Log-Achse linear zu sein beginnt? –

In dem obigen Code wählen Sie xmin Werte mit dem 'xmins' Argument. – csgillespie

Antwort

Verwandte Themen