r - 适合收入分配的幂律
问题描述
我正在处理不同国家的家庭总净收入的微观数据，在计算不平等指数之前，我想调整右尾以考虑这些指数对最高收入和最高收入的敏感性。因此，我正在尝试使用 powerlaw 包来拟合右尾的幂律分布。但是，我对 Xmin 有一些奇怪的估计。我使用了以下代码：
m_c1<-conpl$new(country2_1$hi100_eur)
est=estimate_xmin(m_c1)
m_c1$setXmin(est)
结果 Xmin 为：99863，但对于相关国家/地区而言，这意味着大约是第 25 个百分位。（在这种情况下，平均收入为 223290，丹麦电晕）。实际上，绘制 qqnorm 可以观察到右尾的偏差开始远高于估计的 Xmin，但在第 95 个百分位附近（红线）：
最后，还用获得的最小值绘制估计的幂律的 cdf，这显然不是最佳拟合：
鉴于数据集很大，我报告了我用来拟合幂律的列向量（家庭收入）：
dput(head(country2_1$hi100_eur))
country2_1$hi100_eur<-c(97000, 126550, 90000, 249500, 249500, 156166)
知道我做错了什么吗？谢谢！
解决方案
感谢 Gillespie 教授的澄清，我意识到问题出在 xmax 的定义中。默认情况下它设置为 1e5，但我的数据大于 xmax=1e5。修改 xmax 范围（或设置 xmin）拟合幂律没有任何奇怪的结果。