幂律函数与R中数据的拟合

# 安装并加载powerlaw包
install.packages("powerlaw")
library(powerlaw)

# 生成幂律分布的随机数
shape <- 0.5  # 形状参数，对应幂律分布的指数部分
scale <- -1 / shape  # 计算对应的尺度参数
power_law_randoms <- rgamma(n = 1000, shape = shape, scale = scale)

# 拟合幂律分布数据
data <- read.csv("data.csv")  # 读取数据
fit <- fit_power_law(data$value)  # 拟合幂律分布

# 绘制拟合结果
hist(power_law_randoms, log="y", main = "Power Law Distribution Random Numbers", xlab = "Values", ylab = "Frequency (log)")