ФИНАНСЫ И БАНКОВСКОЕ ДЕЛО
DOI 10.26425/1816-4277-2018-1-141-148
ПОСТРОЕНИЕ КРИВОЙ ЛОРЕНЦА И ОЦЕНКА ИНДИКАТОРОВ ДИФФЕРЕНЦИАЦИИ ДЕНЕЖНЫХ ДОХОДОВ НАСЕЛЕНИЯ НА ОСНОВЕ ЭКСПОНЕНЦИАЛЬНОГО РАСПРЕДЕЛЕНИЯ
Аннотация. В статье обсуждаются вопросы применения двухпараметрической функции экспоненциального распределения для построения кривой Лоренца и вычисления таких показателей дифференциации населения по доходам, как коэффициент Джини, коэффициент Рейнбоу, коэффициент фондов и др. Показано, что все эти коэффициенты зависят от отношения двух параметров: минимального и среднедушевого денежных доходов населения. Для количественной оценки параметров функции распределения использованы численные методы решения задачи минимизации различных критериев близости теоретической кривой и статистических данных.
Ключевые слова: денежные доходы, население, функция распределения, кривая Лоренца, коэффициент Джини, коэффициент Рейнбоу, коэффициент фондов.
THE CONSTRUCTION OF THE LORENZ CURVE AND ESTIMATION OF INDICATORS OF DIFFERENTIATION OF MONETARY INCOMES OF POPULATION BASED ON EXPONENTIAL DISTRIBUTION
Annotation. The article discusses the use of the two-parameter exponential distribution function to construct the Lorenz curve and to calculate such indicators of differentiation ofpopulation by income, as the Gini coefficient, Rainbow coefficient, Funds coefficient etc. It is shown that all these coefficients depend on the ratio of two parameters: the minimum and the average per capita monetary incomes of the population. For a quantitative estimation of the parameters of the distribution function numerical methods of solving the problem of minimization of different criteria of the proximity of theoretical curve and statistical data were used. Keywords: incomes, population, distribution function, Lorenz curve, Gini coefficient, Rainbow coefficient, Funds coefficient.
Вопросам анализа уровня жизни населения посвящено большое количество работ [1; 2; 4; 5]. Для количественной оценки показателей уровня жизни населения используются различные подходы и источники информации, в том числе данные о распределении населения по величине среднедушевых денежных доходов, ежегодно публикуемых Росстатом [6]. Анализ этой информации осложняется следующими обстоятельствами. Во-первых, исходной относительной точностью статистических данных, которые получают на основе выборочного обследования бюджетов домашних хозяйств и макроэкономического показателя денежных доходов населения. Вторая причина неточности числовых показателей дифференциации связана с отсутствием в исходных данных Росстата информации о распределении по доходам как малообеспеченных, так и высокообеспеченных групп населения. Методик прямого измерения доходов в этих диапазонах не существует [4]. Поэтому значения индикаторов дифференциации населения России
© Лебедев В.В., Лебедев К.В., 2018
Работа выполнена при финансовой поддержке РФФИ (грант №16-06-00280)
УДК 519.254 JEL J39
Лебедев В.В. Лебедев К.В.
Lebedev Valery Lebedev Konstantin
по денежным доходам существенно зависят от выбора той или иной гипотезы о распределении как низкодоходных, так и высокодоходных групп населения. Сказанное означает, что при исследовании дифференциации населения по доходам используется концепция «мягких» измерений, которая предполагает применение различных подходов к определению числовых значений индикаторов дифференциации и их сравнительный анализ, позволяющий получить обоснованную оценку значений этих индикаторов с некоторой точностью [3].
Подготовлено авторами по материалам исследования
Рис. 1. Распределения населения России по денежным доходам в 2014 г. (слева) и в 2015 г. (справа)
100
У „Л- —
«О
70 >'
60 50
40 А
/
30 д
20 А
» /
0
10 20 30 40 50 60 70
Информация Росстата о распределении населения по величине среднедушевых денежных доходов в месяц может быть приведена в таблице вида (х;; у;.) [6]. Здесь хк - значение дохода (в тыс. руб.), ук - доля населения в %, имеющего денежный доход в месяц ниже уровня хк, к = 1, 2,..., N. На рис. 1 треугольниками отмечены точки (х,; у,), соответствующие статистическим данным о распределении населения России по среднедушевым денежным доходам в месяц в 2014 г. и 2015 г. Такое распределение населения по доходам характерно для всего периода 2009-2016 гг.
Обычно используемые функции распределения у = 1'(х. а.....а ) зависят от т параметров а 1.....а . Для
определения параметров теоретической функции распределения требуется решить задачу сглаживания, которая заключается в построении на плоскости Оху аппроксимирующей кривой, проходящей вблизи таблично заданных точек (хк; ук) , где к 1. 2,..., N. При этом задача аппроксимации решается в некотором классе функций. Так, наиболее часто используемое логнормальное распределение является двухпараметрическим: у = /'Ух, ц, <з) .
Формально задача аппроксимации сводится к минимизации конкретного критерия близости теоретической
линии у 1< (х. а 1.....а т) и статистических данных (х;; ук), например, суммы квадратов отклонений расчетных
значений от табличных: д,
(1)
А-=1
Распределение населения по величине среднедушевых денежных доходов в месяц х по методологии, применяемой Росстатом, описывается двухпараметрическим логарифмически нормальным распределением [1; 4]. Параметрами этого распределения служат математическое ожидание и среднее квадратическое отклонение нормального распределения величины 1пх.
Для анализа денежных доходов населения применяют и другие подходы, например, трехпараметрическое распределение населения по величине среднедушевых денежных доходов, частным случаем которой является логарифмически нормальное распределение, или логистическе распределение, которое также имеет три параметра [4; 5]. В этой работе предполагается, что функция распределения у = /<' (х, х0,"/.. а) при х > х0 является решением модифицированного дифференциального уравнения Ферхюльста с начальным условием у(х0) = 0 . Здесь х0,"/. и а- положительные параметры.
^ = Щ-У){1 + А. (2)
ах \ а ]
Использование модели (2) позволило получить формулы для вычисления различных индикаторов дифференциации денежных доходов населения от значений параметров х0, л и а . Количественные оценки показателей дифференциации населения по доходам существенно зависят как от выбора критерия близости теоретической кривой и статистических данных, так и от точности аппроксимации [5]. Согласно выполненным расчетам, несущественное снижение точности может привести к значительному изменению значений индекса Джини, коэффициента фондов и среднемесячного денежного дохода населения. Несмотря на относительную точность, использование модели (2) позволило, тем не менее, «уловить» тенденцию роста среднего месячного денежного дохода населения России в 2009-2016 гг. со среднегодовым темпом прироста около 9 % [5]. Таким образом, этот подход показал свою дееспособность: он может эффективно использоваться при аппроксимации таблично заданных функций, графики которых имеют ^-образную форму с горизонтальной асимптотой.
Отметим важное обстоятельство: таблично заданные функции распределения населения по доходам зачастую дают основание принять гипотезу о том, что кривая у = F (х, ах,..., a J не имеет точек перегиба. В этом случае в качестве функции распределения можно использовать двухпараметрическую функцию у = 1< (х. х0, И), являющуюся при х > х0 решением дифференциального уравнения
Ф = (3)
dx h
с начальным условием v(x0) = 0. Целесообразность перехода от модели логистического распределения (2) к модели
(3) определяется тем, что последняя является предельным случаем модели (2) при h = \! X и а —> +со . Выпишем функцию распределения у = F(x, х0, И):
ГО, х < х0;
Здесь параметрами служат х0 - минимальный уровень денежного дохода населения, и число h - параметр, определяющий скорость роста функции у = /< (х. х0, И). Для сокращения записи в дальнейшем вместо формулы
(4) будем писать /< (х) = 1 - е <Л л"1 Л, х > х0, помня, что при х < х0 значения этой функции равны нулю. Отметим, что это замечание относится ко всем встречающимся в тексте функциям, зависящим от переменной х.
Итак, будем считать, что задача аппроксимации (1) решена, и параметры х0 и h функции (4), определяющей долю населения, денежный доход которого меньше х, найдены. На рис. 1, кроме статистических данных (xk',vk) распределения населения России по величине среднедушевых денежных доходов в месяц в 2014 и 2015 гг. приведены графики соответствующих теоретических функций (4). Здесь параметры экспоненциальной функции распределения (4) получены на основе решения численными методами задачи минимизации среднего квадратического отклонения теоретической кривой от эмпирических точек (для 2014 г. х0=5,78 и /7=21,05; для 2015 г. х0=6,12 и /7=23,65).
При найденных значениях параметров х0 и h функции (4) можно получить несложные расчетные формулы для вычисления среднемесячного денежного дохода населения, коэффициента Джини, коэффициента Рейнбоу, коэффициента фондов и других показателей, характеризующих неравномерность распределения доходов, а также уравнение кривой Лоренца. Для этого построим сначала функцию плотности распределения /(х) = F'(x):
/(х) = е~(х_х°)//г / h, х > х0 • (5)
Отметим, что в силу уравнения (5) параметр h определяет наклон графика функции (4) при х = х0, т. к. F'(x 0) = 1/ И.. Ниже мы увидим, что параметр h имеет еще одну важную трактовку.
Для дальнейшего изложения удобно ввести еще один параметр: N() - общую численность населения. Тогда функция N(x) = N0F(x) определяет число жителей, имеющих денежный доход в месяц ниже уровня х. Из сказанного следует, что функция Nt)G(x), где G(x) равна:
= , (6)
■Ч)
определяет совокупный месячный доход трудящихся, заработная плата которых меньше х. А это значит, что совокупный месячный доход всего населения равен:
ыов{лсо) = ыо (7)
и поэтому значение среднедушевого месячного дохода населения составит:
хт = С(+со) • (8)
Для используемого здесь экспоненциального распределения населения по доходам (4) из формул (5)-(8) следует:
С(х) = х0+И-(х + И)е~{х~х°]/1\ х > х0; (9)
хт,=х0+И. (Ю)
Поэтому, согласно модели (3), среднедушевой доход населения России был равен в 2014 г. 26,83 тыс. руб., а в 2015 г. равен 29,77 тыс. руб.: среднедушевой доход за один год увеличился на 9,88 %.
Итак, в силу (10) параметр И равен разности между среднедушевым денежным доходом всего населения в месяц и его минимальным (расчетным) значением. Вычислим долю населения, которая имеет среднедушевой доход менее хаг. Имеем:
Р{хт,) = \-Ие- 0,632. (11)
Поэтому в случае экспоненциального распределения по доходам 63,2% населения имеет доход ниже среднедушевого месячного дохода. Обозначим ^(хш,) = ут,.
Доля совокупного месячного дохода населения, у которого доход меньше х, в совокупном денежном доходе всего населения составляет
Я(х) = вдх) = ед (12)
Для экспоненциального распределения эта функция определяется уравнением
Н(х) = -, Х>Х0. (13)
х0 +п
В силу уравнения (13) имеем:
(х0 + И + И)е~1
НЫ = 1--——, (14)
т. е. Н(хаг) = 1 - (1 + g) / е, где g = h / хаг. Поэтому доля трудящихся, денежный доход которых меньше хш, , превосходит долю совокупного месячного дохода этой части населения в совокупном денежном доходе всего населения на величину
Р(хт)-Щхт) = 8/е • (15)
Как видим, параметр g содержит важную информацию: из формулы (15) следует, что чем больше значение этого параметра, тем меньшую долю совокупного дохода получают 63,2 % населения, имеющего доход ниже среднедушевого месячного дохода хаг. Далее мы увидим, что все используемые показатели дифференциации распределения населения по доходам являются функциями этого параметра. Так как ^ = V / (1 + V), где V = И / х0, то параметр V естественно назвать базовым параметром дифференциации денежных доходов населения или
Подготовлено авторами по материалам исследования Рис. 2. Кривая Лоренца (линия 1) и линия «полного равенства» (линия 2)
базовым параметром неравенства. Важно, что этот параметр V равен отношению двух параметров функции (4), которые находятся из решения соответствующей задачи аппроксимации. При этом в силу неравенства g'(v) > О параметр g возрастает при увеличении параметра V. Кроме того, £ —» 0 при V —> 0. Поэтому чем больше расчетное значение базового параметра неравенства V, тем выше дифференциация населения по доходам. Для данных 2014 г. и 2015 г. имеем: параметр V равен 3,64 и 3, 86; параметр g - 0,785 и 0,794, соответственно. Как видим, приросты параметров V и g за год составили 5,7 % и 1,13 %, соответственно.
Построим теперь функцию Лоренца г = Ь(у), которая используется не только для визуализации неравенства в распределении доходов населения, но и для вычисления основных индикаторов дифференциации. Функция Лоренца выражает количественную зависимость накопленной (кумулятивной) доли дохода населения г от соответствующей кумулятивной его доли у .
С формальной точки зрения функция Лоренца - это функция, заданная параметрически следующим образом:
Ш <16>
где х - параметр. Из сказанного следует, что
2 = Н{Р-\У)). (17)
График функции (17) называют кривой Лоренца (рис. 2). Так как для экспоненциального распределения функция х = К_1 (у) принимает вид
х = х0 — И 1п(1 — у), (18)
то (после соответствующих преобразований) для функции Лоренца г = Ь(у) получаем
л ч [;; + &(1-;;)1п(1- У), если уе [0,1);
Ц}') = \Л . (19)
[1, если у = 1.
Из формулы (19) следует, что отклонение кривой Лоренца от прямой г = у стремится к нулю при g —> 0 . Поэтому линию г = у называют линией «полного равенства» (при g —> 0 все население получает доход на уровне х0, т. к. в этом случае /? —> О . а х,п, х0). Найдем максимум отклонения графика функции (17) от прямой г = у . Для этого надо решить уравнение Ь\у) = 1. Так как Ь'{у) = \- g\ri{\ - у) - g, то из условия Ь'{ у) = 1 следует 1п(1 —= -1, откуда получаем у = 1-6' 1. т.е. у = усп,- В силу формулы (15) максимум отклонения ут, - /-(>'„,,) = 1' {хаг) - Н(хау) = ^ / е. Отметим также, что //(0) = х0 / хаг. т. е. //(0) = 1 - # < 1. Как видим, «провисание» кривой Лоренца сокращается при уменьшении параметра g и, следовательно, при уменьшении базового
параметра неравенства распределения населения по доходам V . На рис. 2 построены линия «полного равенства» и кривая Лоренца (19) для 2015 г. (здесь значение параметра g =0,794). Согласно выполненным расчетам, кривые Лоренца для 2014 г. и 2015 г. практически совпадают, т. к. параметры g различались в эти годы всего лишь на 1,13 % (в 2014 т. g = 0,785).
Вычислим теперь некоторые индикаторы неравенства распределения населения по доходам. Коэффициент Джини КС по определению равен удвоенной площади области, лежащей на рис. 2 между линией «полного равенства» (2) и кривой Лоренца (1). Поэтому
1
КО = 1-2\Ь(уУ}у. (20)
о
Из определения следует, что значение коэффициента Джини изменяется в диапазоне от 0 до 1. При этом чем больше значение коэффициента Джини отличается от нуля, тем в большей степени совокупные доходы распределены неравномерно. Подставив функцию (18) в формулу (19) и вычислив по частям соответствующий несобственный интеграл, для коэффициента Джини получим:
КС = 8/ 2. (21)
Коэффициент Рейнбоу определяется следующим образом: КЯ = х„ / х». Здесь х„ - величина такого дохода, который имеют не менее 10 % самых высокообеспеченных членов общества; х» - величина дохода, менее которого имеют 10 % членов общества с самыми низкими доходами [2]. Поэтому КЯ = 1< 1 (0.9) / /• 1 (0.1), т.е. здесь в силу (17) имеем:
1 + е 1п 10
КЯ =---. (22)
1 + £1п(10/9)
Коэффициент фондов - это отношение среднедушевых доходов 10 % членов общества с самыми высокими доходами к среднедушевым доходам 10 % членов общества с самыми низкими доходами. Поэтому коэффициент фондов определяется по формуле КЬ' = (1 -¿(0,9)) /¿(0,1) . Для экспоненциального распределения (4) в силу формулы (18) следует:
ДГ= 1 + ^П1° . (23)
1-9£1п(10/9)
Коэффициент дифференциации удельных полярных доходов К1)а - это отношение среднего денежного дохода двух частей общества. Первая из них имеет денежные доходы выше среднедушевого месячного дохода всего населения, вторая - ниже этого уровня. Поэтому
ш _ \ -Н(хт) Н(хт) а \-FixJ Р{хт)'
В силу (11) и (15) получаем:
кв =(1 + я)/е.1-(1 + я)/е = (1 + я)(е-1) 1 /<? 1-1 /<? (<?-1)-£
Для данных Росстата 2015 г. рассчитанное выше значение параметра g = 0,794 . Поэтому коэффициент Джини равен 0,397, коэффициент фондов - 11,47, коэффициент Рейнбоу - 2,61, коэффициент дифференциации удельных полярных доходов К1)а -3,34. Выше было установлено, что среднедушевой доход в 2015 г. увеличился на 9,88 %. Несмотря на это все значения индикаторов дифференциации увеличились незначительно, т. к. прирост параметра g составил 1,13 %.
На рис. 3 приведены графики полученных выше зависимостей от величины базового параметра неравенства населения V следующих индикаторов дифференциации населения по доходам: коэффициента Джини (линия 1), коэффициента фондов К/' (линия 2), коэффициента дифференциации удельных полярных доходов К1)а (линия 3),
Подготовлено авторами по материалам исследования
Рис. 3. Зависимость индикаторов дифференциации населения от базового параметра V = И /х
коэффициента Рейнбоу КЯ (линия 4). Для наглядности коэффициент Джини увеличен здесь в 10 раз (при V = 4 коэффициент Джини равен 0,4). Выполненные расчеты этих и других индикаторов дифференциации населения по доходам дают основание сделать следующий вывод: наиболее чувствительными к изменению базового параметра V является коэффициент фондов, а наименее чувствительным - коэффициент Рейнбоу. Так, например, при увеличении параметра V от 3 до 4 (на 33,3 %) коэффициент Джини увеличивается от 0,375 до 0,4 (на 6,7 %), а коэффициент Рейнбоу - от 2,527 до 2,621 (на 3,7 %).
Использованный подход позволил получить не только уравнение функции Лоренца, но и несложные формулы для расчета индикаторов дифференциации населения по доходам. Более того, несмотря на относительную точность, он «уловил» тенденцию роста среднего месячного денежного дохода населения России в течение 20092015 гг. со среднегодовым темпом прироста на уровне 9-10%, что соответствует расчетам на основании использования других подходов. При этом, однако, согласно модели (3) росло и расчетное значение параметра хп, вследствие чего значение базового параметра дифференциации населения V в течение этого периода мало отличалось от 4, и, соответственно, значение коэффициента Джини мало отличалось от 0,4.
Библиографический список
1. Айвазян, С. А. Анализ качества и образа жизни населения (эконометрический подход) / С. А. Айвазян. - М. : Наука, 2012.-402 с.
2. Гречаный, С. А. Коэффициент Рейнбоу и возможности введения прогрессивного налога в России / С. А. Гречаный, В. А. Родин // Вестник ВГУ, Серия «Экономика и управление». - 2008. - № 2. - С. 44-47.
3. Клейнер, Г.Б. Системная экономика, экономическая кибернетика, мягкие измерения и проблемы развития общества / Г. Б. Клейнер // Экономика и управление: проблемы, решения. - 2017. - № 3. - Том 4 (63). С. 3-7.
4. Колмаков, И.Б. Методы и модели прогнозирования показателей дифференциации денежных доходов населения / И. Б. Колмаков. - М. : Институт микроэкономики, 2004. - 168 с.
5. Лебедев, В. В. О применении уравнения Ферхюльста для анализа дифференциации денежных доходов населения России / В. В. Лебедев, К. В. Лебедев // Экономика и управление: проблемы, решения. - 2017. - № 6. - Том 3 (66). - С. 218-224.
6. Российский статистический ежегодник. 2016: Стат. сб. / Росстат. - М. : 2016 - 725 с.
References
1. Ajvazyan S. A. Analiz kachestva i obraza zliizni naseleniya (ehkonometricheskij podhod) [Analyzing the quality and lifestyle of the population (econometric approachJ], Moscow, Nauka, 2012. 402 p.
2. Grechanyj S. A., Rodin V A. Koehliicient Rejnbou i vozmozhnosti wedeniya progressivnogo nalo-ga v Rossii [Rainbow coefficient and the possibility of introducing a progressive tax in Russia], //Vestnik VGU, Seriya: ehkonomika i upravlenie, 2008,1. 2. pp. AA-A1.
BecmuuK yuueepcumema № 1, 2018
3. Klejner G.B. Sistemnaya ehkonomika, ehkonomicheskaya kibernetika, myagkie izmereniya i problemy razvitiya obshchestva
[System Economics, economic Cybernetics, soft measurement and problems of development of society]. Ekonomika i upravlenie: problemy, resheniya, 2017, I. 3, Vol. 4 (63), pp. 3-7.
4. Kolmakov I.B. Metody i modeli prognozirovaniya pokazatelej differenciacii denezhnyh doho-dov naseleniya [Methods and models of forecasting of indicators of differentiation of monetary incomes of the population]. Moscow, Institut mikroehkonomiki, 2004, 168 p.
5. Lebedev V V, Lebedev K.V O primenenii uravneniya Ferhyul 'sta dlya analiza differenciacii de-nezhnyh dohodov naseleniya Rossii [On the application of the Verhulst equation for the analysis of differentiation of monetary incomes of the population of Russia] // Ekonomika i upravlenie: problemy, resheniya, 2017, I. 6, Vol. 3 (66), pp. 218-224.
6. Rossijskij statisticheskij ezhegodnik [Russian statistical Yearbook. 2016]. Moscow, 2016. 725 p.