Вычисление процентиля на основе данных Survey of Consumer Finance

Я знаю, что есть пакет под названием survey, делающий подобные вещи. Я не знаю, есть ли какой-либо код или способ вычисления процентиля в R.

Что я хотел бы получить, так это относительный собственный капитал, которым владеют верхние x% домохозяйств. Таким образом, это будет соотношение:

  • (Общее богатство верхнего x процентиля домохозяйства (от 0,99 квантиля до максимума) / (общее богатство в экономике) *

Я нашел svrqsr в пакете survey, но это дает мне относительную долю по сравнению с домохозяйствами с более низким 20 процентилем.

Любое предложение или ссылка на код, которую я могу использовать?


r percentile inequality
person user9620331    schedule 01.05.2018    source источник
comment
См. здесь советы по созданию воспроизводимого сообщения, которое люди могут легко ответить. Я предполагаю, что не могу видеть ваши данные или код, но quantile в базе R не работает для ваших целей?   -  person camille    schedule 01.05.2018
comment
Я заметил, что могу использовать svylorenz для вычисления чисел. Я пробую это.   -  person user9620331    schedule 01.05.2018

Ответы (2)


arrow_upward
2
arrow_downward

следуйте инструкциям на http://asdfree.com/survey-of-consumer-finances-scf.html, пока не дойдете до

library(convey)
scf_design$designs <- lapply( scf_design$designs , convey_prep )

затем используйте

scf_MIcombine( with( scf_design , svyqsr( ~ networth , alpha1 = 0.01 ) ) )
person Anthony Damico    schedule 15.05.2018

arrow_upward
0
arrow_downward

Как сказал @camille, quantile может дать вам ваш порог (в вашем случае собственный капитал домохозяйства 99-го процентиля). После этого остается только суммировать значения, превышающие пороговое значение.

amt<-10^rnorm(1e5,4,1) # Assuming a log normal distribution
amt_99<-amt[amt>quantile(amt,0.99)]  #Top 1%
sum(amt_99)/sum(amt)  # Relative net worth
#Output : [1] 0.4521971
person Rohit    schedule 02.05.2018