R中拟合单参数帕累托分布出现错误码100的解决办法咨询

我来帮你拆解这个问题——你遇到的报错核心原因是初始参数设置违反了单参数帕累托分布的基本约束，咱们一步步来搞定：

错误原因分析 #

单参数帕累托分布（actuar包中的pareto1）有一个硬性要求：所有观测数据必须大于等于阈值参数min（这个参数是分布的下限阈值，只有当x ≥ min时，概率密度才有非零定义）。

你设置的初始min=12e6（1200万），但看你的数据里有大量数值远小于这个值（比如459745、803526、200多万的样本比比皆是），这就导致这些数据点在该初始参数下的对数似然值为负无穷，optim函数根本无法计算目标函数，直接抛出"无法在初始参数处评估函数"的错误。

解决步骤 #

1. 确定合理的min初始值 #

首先计算你的数据最小值，min的初始值必须小于等于这个最小值：

min_val <- min(lrg_clms)
# 输出结果：459745

所以我们可以直接用数据最小值作为min的初始值，或者稍小一点的数值（比如400000），确保所有数据都满足x ≥ min的要求。

2. 调整初始参数重新拟合 #

修改start参数，用合理的min值重新运行拟合代码：

require(actuar)
require(fitdistrplus)

# 定义你的数据
lrg_clms <- c(20000000.4, 2943143.86, 2223020, 3156203.17, 2581621.69, 45842102.64, 2303396.3, 74924503.71, 2107205.28, 803526, 115984539.2, 2904623, 6e+06, 4782296.13, 2354548.46, 2389383.84, 2288983.99, 4436729, 7500000, 1045758.15, 2495385.48, 126794168, 5546141.72, 2238181.12, 2471143.02, 4839064.67, 15293527.76, 24888213.71, 459745, 1790472.95, 4730418.92, 7241093, 14574723.75)

# 使用数据最小值作为min初始值，shape初始值设为1
fit2 <- fitdist(lrg_clms, "pareto1", start = list(shape = 1, min = min(lrg_clms)))

# 查看拟合结果
summary(fit2)

3. 可选：用矩估计优化shape初始值（提升拟合稳定性） #

如果想让初始参数更贴合数据，可以用矩估计先计算shape的初始值。单参数帕累托的矩估计公式为：当shape > 1时，均值μ = min * shape / (shape - 1)，反推可得：

mu <- mean(lrg_clms)
init_shape <- mu / (mu - min_val)

# 使用矩估计的shape初始值拟合
fit2 <- fitdist(lrg_clms, "pareto1", start = list(shape = init_shape, min = min_val))

这样调整后，拟合应该就能正常运行了。

内容的提问来源于stack exchange，提问作者user33484