r - R给定分布中幂律的最大似然估计（而不是样本）

Question

我有一个数据框，其中 xy 值表示值及其计数，例如 (1, 1000)、(2, 100)、(3, 10) 等。我想使用 MLE 对这个分布拟合幂律。

我可以使用 power.law.fit 或 power.law 库，但这些库似乎采用了特定的数据样本，而不是表示值及其计数的 xy 值。

有没有其他图书馆可以做这项工作？谢谢！

Answer 1

你可以使用这个poweRlaw包——它只是有点笨重。只需将您的值和计数扩展为单个向量，例如

dd = data.frame(x=1:3, counts = 3:1)
x = rep(dd$x, dd$counts)
library(poweRlaw)
m = displ$new(x)