幂律分布拟合

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)

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