我有一个数据框,其x-y值代表值及其数量,例如: (1,1000),(2,100),(3,10)等我想用MLE拟合这个分布的幂律。
我可以使用power.law.fit或poweRlaw库,但看起来这些库接收数据的特定样本而不是表示值及其计数的x-y值。
是否有其他图书馆可以完成这项工作?谢谢!
答案 0 :(得分:1)
您可以使用poweRlaw
套餐 - 它只是有点笨拙。只需将您的值扩展并计入单个向量,例如
dd = data.frame(x=1:3, counts = 3:1)
x = rep(dd$x, dd$counts)
library(poweRlaw)
m = displ$new(x)