r - 幂律分布拟合

Question

我有一个数据样本，样本量约为 500,000。我目前正在尝试使用幂律分布拟合样本poweRlawR 中的包来拟合具有幂律分布的样本。

所以这是我的代码：

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

但是，当我尝试运行代码时出现问题x<-estimate_xmin(pl_rg,xmax = 100). 我发现这段代码非常非常耗时。到目前为止，我已经运行了这段代码长达 5 个小时，它仍在运行。

那么有什么方法可以加速这个过程或其他符合幂律的方法吗？

Answer 1

缓慢的原因estimate_xmin是它正在尝试数据的每个可能值以获得潜在xmin值。您可以通过指定值来显着加快功能，例如

estimate_xmin(pl_rg, xmins = 1:10)

有关更多详细信息，请参阅帮助页面。