我正在处理不同国家的家庭总净收入的微观数据,在计算不平等指数之前,我想调整右尾以考虑这些指数对最高收入和最高收入的敏感性。因此,我正在尝试使用 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)
知道我做错了什么吗?谢谢!