r - 使用 poweRlaw 包的幂律拟合截止分布

Question

我目前正在尝试寻找一种方法来计算适合 MLE 截止分布的幂律。分布如下所示:

如您所见，我能够分别拟合整个分布(幂律拟合)和下限(exp-fit)。我没弄清楚的是如何拟合分布的上限(例如 8 < x < 100)。

有什么方法可以用 poweRlaw 包或任何其他 R 包来做到这一点？我所希望的是看起来像这样的东西(注意:它只是一个随机分布):

代码(如有必要):

#Power-Law
library("poweRlaw")
xmin1 <- 8
xmin2 <- 100
plf0 <- displ$new(deg)
plf0$setXmin(xmin1)
plf0_pars <- estimate_pars(plf0)
plf0$setPars(plf0_pars)

#Exponential
exp1 <- disexp$new(deg)
exp1$setXmin(xmin2)
exp1_pars <- estimate_pars(exp1)
exp1$setPars(exp1_pars)
plot(plf0)
lines(plf0, col="green")
lines(exp1, col="red")

Answer 1

你不能用 poweRlaw 包(我是包的作者)来适应这种类型的模型。

在near package中添加到package中的可能性不大。