幂律分布拟合

Question

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

这是我的代码：

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)

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