Я работаю с микроданными об общем чистом доходе домохозяйства для разных стран, и прежде чем вычислять индексы неравенства, я хочу скорректировать правильный хвост, чтобы учесть чувствительность этих индексов к верхним доходам и пропущенным.Поэтому я пытаюсь согласовать распределение степенного закона для правого хвоста с пакетом 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)
Есть идеи о том, что я делаю неправильно?Спасибо!