Корреляционно-регрессионный анализ взаимосвязей социально-экономических явлений

КОРРЕЛЯЦИОННО-РЕГРЕССИОННЫЙ АНАЛИЗ ВЗАИМОСВЯЗЕЙ СОЦИАЛЬНО-ЭКОНОМИЧЕСКИХ ЯВЛЕНИЙ

План

1. Понятие, виды и формы факторных связей.

1.1. Понятие о функциональных и стохастических связях. Корреляционная

связь.

1.2.Графическое представлениекорреляционной связи

1.3.Виды и формы выражения корреляционной связи

1.4. Свойства корреляционной связи

2. Этапы корреляционно-регрессионного анализа, требования к исходному

статистическому материалу

3. Корреляционный анализ взаимосвязи признаков

3.1. Основные статистические методы выявления наличия

и направления корреляционной связи

3.2. Измерение тесноты корреляционной связи

4. Регрессионный анализ взаимосвязи признаков.

4.1. Выбор аналитической формы связи.

4. 2. Расчет числовые значения параметров математической функции связи.

4. 3. Анализ качества построенной модели.

5. Непараметические методы обнаружения взаимосвязей социально-

экономических явлений.

Пример 1.

Пусть первичные данные статистического наблюдения представлены в табл. 1, где каждому значению x_i фактора X соответствуют несколько значений результативного признака Y.

Таблица 1

Первичные данные статистического наблюдения

xi

yi

 

 

Зависимость Y от X можно представить в виде соответствий:

для x₁=8: ;

для x₂=9:

и т.д.

 

Средние арифметические значения распределений признака Y, вычисленные при условии, что фактор X принимает некоторое определённое значение x_i, называют условными средними. С изменением значения x_i меняется и средняя .

Для рассмотренного примера зависимость между значениями x_i и условными средними определяет табл. 2.

Таблица 2

Корреляционный анализ взаимосвязи признаков

Основные статистические методы выявления наличия

Пример.

 

x
y	-7	-3
z									-3	-11
w				-5						-3

 

Обращаясь к исходным данным примера, можно говорить о наличии прямой корреляционной связи между признаками Х,Y, обратной связи - между признаками Х, Z и отсутствии связи между Х,W.

Метод параллельных рядов обычно используется для установления характера связи при относительно небольшом объёме исходных данных.

2. Метод аналитической группировки

 

В статистической практике изучение взаимосвязей явлений проводится, как правило, по достаточно большому числу наблюдений, а значения наблюдаемых признаков X и Y представляются в сгруппированном виде

Наиболее удобной формой представления корреляционных зависимостей при большом числе наблюдений являются групповые аналитические таблицы, отражающие результаты аналитической группировки статистической совокупности по факторному признаку.

 

При построении аналитической таблицы для каждой выделенной j -ой группы подсчитывается численность составляющих её факторных значений x, а также суммарное и среднее групповые значения результативного признака

 

Пример.

Зависимость суммы прибыли банков от объема кредитных вложений.

Номер группы	Группы банков по объему кредитных вложений, млн. руб., х	Число банков, fj	Сумма прибыль, млн. руб.
всего	в среднем на один банк
	40 – 90		50,4	16,800
	90 – 140		241,3	40,217
	140 – 190		711,4	59,283
	190 – 240		704,0	78,222
	Итого		1707,1

 

.Если при переходе от одной группы к другой средние значения будут изменяться с определенной закономерностью - систематически возрастать или убывать, то между признаками X и Y существует корреляционная связь ( как в приведенной таблице ).

 

3. Метод корреляционной таблицы

 

При построении корреляционной таблицы значения признаков X и Y ранжируются в порядке возрастания, факторные значения x_i располагаются, как правило, в строках таблицы, результативные – в столбцах (графах), а на пересечении строк и столбцов проставляются числа, указывающие частоту появления различных результативных значений y_j при фиксированном факторном значении x_i. При таком построении таблицы каждая i -ая строка представляет распределение признака Y при условии X= x_i. В итоговой строке таблицы проставляются частоты повторения в эмпирических данных результативного значения y_j, а в итоговом столбце – частоты повторения значений x_i.

Пример.

Зависимость суммы прибыли банков от объема кредитных вложений.

Группы банков по объему кредитных вложений, млн. руб.	Группы банков по сумме прибыли, млн. руб.
6,2 – 27,2	27,2 – 48,2	48,2 – 69,2	69,2 – 90,2	Итого
40 – 90
90 – 140
140 – 190
190 – 240
Итого

 

По корреляционной таблице прослеживают визуально, как именно расположена в таблице основная масса частот повторения в эмпирических данных сочетаний (x_i, y_j):

- концентрация частот вдоль диагонали от левого верхнего угла таблицы к правому нижнему (т.е. большему значению X соответствует большее значение Y) означает наличие прямой корреляционной связи между признаками (как в таблице данного примера).

- концентрация частот вдоль диагонали от левого нижнего угла к правому верхнему (когда большему значению Х соответствует меньшее значение Y) означает наличие обратной корреляционной связи между признаками

Интенсивная концентрация частот около диагонали таблицы указывает на факт тесной корреляционной связи. Так, в корреляционной таблице данного примера наблюдается тесная связь между рассматриваемыми признаками.

Корреляционная таблица дает более правильную характеристику тесноты связи в случае, если число выделенных групп одинаково для обоих изучаемым признаков X и Y (как в приведенной таблице).

ПРИМЕР ПОСТРОЕНИЯ ЛИНЕЙНОГО УРАВНЕНИЯ РЕГРЕСИИ

 

Имеются следующие условные данные о размере страховой суммы и страховых выплат на автотранспортные средства одной из страховых компаний

Таблица 7

Расчётная таблица для определения параметровуравнения регрессии

№ регистрации авто	Объём страховых выплат, тыс. долл., у	Стоимость застрахованных авто, тыс. долл., x
1.	0,15	17,6	309,76	2,64	0,077
2.	1,95	18,8	353,44	36,66	0,542
3.	0,15	20,0	400,00	3,00	1,007
4.	3,90	21,2	449,44	82,68	1,472
5.	0,15	22,0	484,00	3,30	1,782
6.	0,45	23,8	566,44	10,71	2,480
7.	6,90	25,4	645,16	175,26	3,100
8.	0,45	27,0	729,00	12,15	3,720
9.	0,60	31,0	961,00	18,60	5,270
10.	10,95	33,4	1115,56	365,73	6,200
Итого	25,65	240,2	6013,80	710,73	25,650

Решение. Система нормальных уравнений (22)

 

для данных таблицы принимает вид:

 

Решение системы методом исключений или по формулам (23), (24) дает:

Тогда модель связи в виде линейного уравнения регрессий имеет вид:

 

 

Расчётные значения страховых выплат , полученных по построенной модели, представлены в последнем столбце табл.7.

Наблюдаются отклонения между фактическими и теоретическими значениями y_i. Однако суммарные значения совпадают (26,65=25,650). Это свидетельствует, что из всех возможных линейных моделей построенная модель является наилучшей.

Это объясняется тем, что она получена методом наименьших квадратов (МНК) путём решения задачи на экстремум.

 

Коэффициент a₁ при факторе Х называется коэффициентом регрессии. Величина коэффициента регрессии a₁ показывает, на сколько в среднем изменяется значения результативного признака Y при изменении фактора Х на единицу его измерения. Знак при a₁ показывает направление этого изменения.

В примере единицей измерения X и Y является 1000 $. Следовательно, коэффициент регрессии в построенной модели показывает, что при увеличении стоимости автомобиля на 1000 долларов объём страховых выплат в среднем вырастает на 0,3875 тыс. долларов.

 

При изучении многофакторных корреляционных связей методология их моделирования уравнениями регрессии аналогична рассмотренной.

 

Уравнения многофакторной регрессии имеют вид:

= f(x_1, x_{2, …,} x_m)

и позволяют приближенно оценить меру влияния на результативный признак Y каждого из включенных в модель факторов X.

Наиболее простыми для построения и анализа являются многофакторные линейные модели, которые содержат независимые переменные только первой степени:

.

 

Методом наименьших квадратов для нее получна систему нормальных уравнений:

 

.

 

Параметры уравнения множественной регрессии показывают степень влияния каждого фактора на анализируемый показатель при фиксированном (среднем) значении всех других факторов.

С изменением факторного признака a_k, на единицу его измерения (при неизменном значении других факторов) результативный признак Y изменяется в среднем на величину параметра a_k.

 

Необходимо отметить, что параметры связи a_k, в силу их расчета по МНК, являются усредненными по всей совокупности наблюдений (x_i, y_i). Они отражают взаимосвязь признаков X и Y только в общем итоге, по всей совокупности в целом (для каждой индивидуальной пары (x_i, y_i) значения a_k остаются неизвестными.

 

3. Анализ качества построенной модели

Этот этап заключается в анализе адекватности и практической пригодности построенной регрессионной модели.

Анализ качества регрессионной модели ставит задачу – определить, насколько построенная теоретическая модель взаимосвязи признаков отражает фактическую зависимость между этими признаками. Анализ выполняется в 3 этапа:

 

1) измерение тесноты связи признаков в уравнении регрессии;

2) определение по величине показателей тесноты связи практической

пригодности построенной модели связи;

3) проверка статистической значимости показателей тесноты связи.

 

Этап 1. Для построенной регрессионной модели измерение тесноты связи признаков X и Y осуществляется на основе следующих показателей:

R² - коэффициент (или индекс) детерминации, показывающий, какая часть общей вариации расчетных (теоретических) значений признака Y объясняется вариацией фактора X;

R – индекс корреляции (называемый также теоретическим корреляционным отношением и обозначаемый ), оценивающий степень тесноты связи между факторными значениями x_i и расчетными результативными значениями ;

r – линейный коэффициент корреляции, используемый для измерения тесноты связи признаков в регрессионной модели в случае линейной функции связи f(x).

Расчет этих показателей определяется следующими формулами:

, , ,

где a₁ – коэффициент регрессии в регрессионной модели связи.

Используя соотношения для вычисления соответствующих дисперсий, а также известное из математической статистики правило сложения дисперсий

,

для расчета показателей R² и R применяют формулы:

, (25)

 

В многофакторном корреляционном анализе для характеристики тесноты связи между изучаемым показателем и влияющими на него факторами, используются парные, частные и множественные коэффициенты корреляции, множественный коэффициент детерминации.

 

 

Этап 2. Практическая пригодность синтезированной регрессионной модели связи оценивается по величине показателей r (в случае линейности модели), R или R².

Значения всех трех показателей R², R и r зависят от того, какая типовая форма (8) уравнения регрессии была выбрана, поэтому эти показатели можно рассматривать как измерители степени близостивыбранной теоретической линии регрессии к фактическим данным.

Качественный характер такой близости может оцениваться в соответствии со шкалой Чэддока.

Индекс корреляции R принимает значения в пределах . При этом:

· близость R к единице означает, что связь между признаками достаточно хорошо описывается избранным уравнением корреляционной зависимости (8);

· если R равен или близок к нулю, это означает, что между фактором X и теоретическими значениями признака Y либо нет связи, либо если она и существует, то не может быть охарактеризована выбранным для моделирования типовым аналитическим выражением связи (8).

Аналогичные утверждения имеют место и для линейного коэффициента корреляции r, принимающего значения в пределах -1 r 1:

· близость к единице свидетельствует о хорошей аппроксимации фактических данных полученной линейной связи = a₀ + a₁x;

· близость к нулю, означает, что уравнение регрессии не может быть линейным.

Пригодность построенной регрессионной модели для практического использования можно оценить и по величине коэффициента детерминации R² :

· неравенство R² >0,5 позволяет считать, что построенная модель пригодна для практического применения, т.к. в ней достигается высокая степень тесноты связи признаков X и Y, при которой более 50% вариации признака Y объясняется влиянием фактора Х;

· неравенство означает, что построенная модель связи практического значения не имеет ввиду недостаточной тесноты связи между признаками X и Y, при которой менее 50% вариации признака Y объясняется влиянием фактора Х, и, следовательно, фактор Х влияет на вариацию Y в значительно меньшей степени, чем другие (неучтенные в модели) факторы.

 

В тех случаях, когда синтезированы различные (альтернативные) регрессионные модели, коэффициент детерминации R² используется в качестве критерия предпочтительности того или иного уравнения регрессии:

наилучшей считается модель с наибольшим значением R².

 

Этап 3. Так как показатели тесноты связи R или r рассчитываются на основе ограниченной совокупности наблюдаемых эмпирических данных (x_i, y_i), значения которых могли быть искажены влиянием случайных факторов, то найденные по уравнению регрессии показатели тесноты связи r, R проверяются на их неслучайность (значимость).

 

В начале проверки задается уровень значимости (чаще всего обозначаемый буквой греческого алфавита «альфа» — ), кото­рый показывает вероятность принятия ошибочного решения. Возможность ошибки вытекает из того факта, что для определе­ния взаимосвязи используются данные не всей совокупности, а лишь ее части. Обычно принимает следующие значения: 0,05; 0,02; 0,01; 0,001. Например, если = 0,05, то это означает, что в среднем в пяти случаях из ста принятое решение о значимости (или незначимости) найденных показателей тесноты связи будет ошибочным; при = 0,001 — в одном случае из тысячи и т.д.

 

Для оценки значимости линейного коэффициента корреляции r применяется t-критерий Стьюдента, фактическое значение которого рассчитывается по формуле:

.

Расчетное значение критерия t_r сравнивается с критическим t_табл, определяемым по таблице значений t-критерия Стьюдента с учетом заданного уровня значимости и числа степеней свободы k=n-2. Если t_расч>t_табл, то величина коэффициента корреляции признается значимой.

 

Для оценки значимости индекса корреляции R применяется F-критерий ФишераF_R, фактическое значение которого определяется по формуле

, (26)

где m – число параметров уравнения регрессии.

Расчетная величина F_R сравнивается с критическим табличным значением F_табл (), которое определяется по таблице F- критерия с учетом принятого уровня значимости и числа степеней свободы k₁=m-1 и k₂=n-m.

Если F_расч>F_табл, то величина найденного индекса корреляции R признается значимой.

Значимость показателя тесноты связи R или r означает, что зависимость между признаками Х и Y регрессионной модели является статистически существенной, т.е. построенная регрессионная модель в целом адекватна исследуемому процессу.

Следовательно, выводы, сделанные на основе регрессионной модели, построенной по данным ограниченной выборки, можно с достаточной вероятностью распространить на всю генеральную совокупность.

В качестве критериев адекватности регрессионной модели в практике экономико-статистического анализа, помимо показателя тесноты связи r, R и R²_, используются также показатели, оценивающие погрешность модели:

· средняя квадратическая ошибка уравнения регрессии , представляющая собой среднее квадратическое отклонение эмпирических значений признака Y от теоретических:

(27)

· средняя ошибка аппроксимации , выраженная в процентах:

(28)

В адекватных моделях ошибки и не должны превышать 12%-15%.

Показатели и могут использоваться при подборе наилучшей модели функции связи: более адекватной считается та модель, в которой меньше величина (или ).

 

 

5. Экономическая интерпретация построенной регрессионной модели

Взаимосвязи признаков

Экономическая интерпретация построенной регрессионной модели взаимосвязи признаков – это, по существу, перевод модели с языка статистики и математики на язык экономики.

Интерпретация уравнения регрессии используется для целей анализа и прогноза взаимосвязей признаков.

 

1. В процессе анализа прежде всего выясняются, как факторный признак влияет на величину результативного признака. Чем больше величина коэффициентов регрессии a_k при k-ой степени фактора Х, тем значительнее влияние данного признака на результативный.

В случае линейного уравнения регрессии = a₀ + a₁x величина коэффициента регрессии a₁ показывает, на сколько в среднем (в абсолютном выражении) изменяется значения результативного признака Y при изменении фактора Х на единицу его измерения. Знак при a₁ показывает направление этого изменения.

Если в соответствии с экономической теорией факторный признак должен иметь положительное значение, а коэффициент a_k имеет знак «-», то необходимо проверить расчёты параметров связи, а также возможность ошибок при сборе и обработке информации.

 

2. С целью расширения возможностей экономического анализа используется коэффициент эластичности , который показывает, на сколько процентов изменяется в среднем результативный признак при изменении факторного признака на 1%.

3. Анализируя остатки , характеризующие отклонения i -ых наблюдений от значений , которые следует ожидать в среднем, можно сделать ряд практических выводов об эффективности экономической деятельности рассматриваемых хозяйствующих субъектов и выявить скрытые резервы их развития и повышения деловой активности. При этом наиболее значительный экономический интерес представляют наибольшие и наименьшие положительные и отрицательные отклонения .

4. Уравнение регрессии может использоваться для краткосрочного прогнозирования ожидаемых значений результативного признака Y в зависимости от тех или иных значений фактора Х.

Прогнозные значения результативного признака Y получают путем подстановки в уравнение регрессии ожидаемых значений признака Х.

Уравнение позволяет определять в рассматриваемом временном периоде возможные значения признака Y при значениях Х из интервала [ x_min - x_max ].

 

При прогнозировании значений Y за рассмотренными пределами изменения фактора X необходимо соблюдать следующее ограничение: нельзя подставлять в уравнение регрессии факторные значения x_i, существенно отличающиеся от тех, на основе которых это уравнение было получено.

Для выполнения этого требования при выборе факторных производных значений х рекомендуется не выходить за пределы 1/3 размаха вариации фактора x (x_max-x_min) как сверх минимального (x_min), так и сверх максимального (x_max)значений, которые имеются в исходных эмпирических данных.

 

Пример (анализ взаимосвязей с проверкой возможности применения метода КРА).

 

Имеются следующие данные по отрасли:

 

№ предприятия	Ежегодная сумма амортизации, млн.руб.	Балансовая прибыль, млн.руб.	№ предприятия	Ежегодная сумма амортизации, млн.руб.	Балансовая прибыль, млн.руб.
			8.
			9.
			10.
			11.
			12.
			13.
			14.

 

Провести на основе приведенных данных исследование взаимосвязи балансовой прибыли и ежегодной суммы амортизации; проверить аналитическое выражение связи на достоверность.

 

Результативный признак – балансовая прибыль (y);

Факторный признак – ежегодная сумма амортизации (х).

 

1. Первичная информация проверяется на однородность по признаку-фактору с помощью коэффициента вариации. Для этого определим среднюю годовую сумму амортизации:

 

 

Результаты расчетов среднего квадратического отклонения приведены в табл. 1:

 

 

Таблица 1.

 

№ предприятия	Годовая сумма амортизации (х)		№ предприятия	Годовая сумма амортизации (х)
	–	–	Итого	–

 

Среднее квадратическое отклонение:

 

 

Коэффициент вариации:

Следовательно, совокупность можно считать однородной.

 

1. Проверка первичной информации на нормальность распределения с помощью правила «трех сигм».

Для этого определим интервалы для значений признака-фактора:

,

 

т.е. (37-4,1)¸(37+4,1); (37-2*4,1)¸(37+2*4,1) и (37-3*4,1)¸(37+3*4,1).

 

Первичная информация по признаку-фактору не подчиняется закону нормального распределения, однако это не является основанием для отказа использования корреляционно-регрессионного анализа (табл. 2).

Таблица 2.

 

Интервалы значений признака х, млн.руб.	Число единиц, входящих в интервал	Удельный вес единиц, входящих в интервал, в общем их числе, %	Удельный вес числа единиц, входящих в интервал при нормальном распределении, %
32,9 – 41,1		64,3	68,3
28,8 – 45,2		100,0	95,4
24,7 – 49,3		100,0	99,7

 

1. Исключение из первичной информации резко выделяющихся единиц, которые по признаку-фактору не попадают в интервал .

По данным таблицы, резко выделяющихся единиц в первичной информации нет.

 

1. Для установления факта наличия связи производится аналитическая группировка по признаку-фактору. Группировка выполняется при равных интервалах и числе групп 4 (табл. 3). Величина интервала определяется по формуле:

Таблица 3.

Пример.

По имеющимся данным о группе рабочих предприятия определите тесноту связи между признаками Среднемесячная производительность труда и Премия по итогам года с помощью:

1) коэффициента Фехнера;

2) коэффициента ранговой корреляции Спирмена;

3) коэффициента линейной корреляции.

Решение.

Построим вспомогательную таблицу для определения необходимых значений, где х – Среднемесячная производительность труда (факторный признак), у – Премия по итогам года (результативный признак).

Таблица

Таблица для расчетов показателей тесноты связи

	x	y						ранг x	ранг y	D ^2
	6,5	15,7	– 0,6	– 0,6	0,36	0,36	0,36
	6,6	15,5	– 0,5	– 0,8	0,40	0,25	0,64
	6,8	16,2	– 0,3	– 0,1	0,03	0,09	0,01
	6,9	16,1	– 0,2	– 0,2	0,04	0,04	0,04
	7,0	15,9	– 0,1	– 0,4	0,04	0,01	0,16
	7,0	15,8	– 0,1	– 0,5	0,05	0,01	0,25
	7,1	17,6	0,0	1,3	0,00	0,00	1,69
	7,1	16,4	0,0	0,1	0,00	0,00	0,01
	7,2	16,5	0,1	0,2	0,02	0,01	0,04
	7,3	16,4	0,2	0,1	0,02	0,04	0,01
	7,4	16,0	0,3	– 0,3	– 0,09	0,09	0,09
	7,5	16,7	0,4	0,4	0,16	0,16	0,16
	7,5	16,3	0,4	0,0	0,00	0,16	0,00
	7,6	17,2	0,5	0,9	0,45	0,25	0,81
å	99,5	228,3			1,49	1,47	4,27

; ;

 

1) = = 0,857 Þ связь тесная, прямая.

 

2) Þ связь тесная, прямая.

 

3) ; ;

 

Þ связь заметная, прямая.

3.Ранговый коэффициент корреляции Кендэла, как и ранговый коэффициент Спирмена, применяют с целью выявления взаимосвязей, как к количественным, так и к качественным признакам, если последние можно проранжировать. Формула для расчёта коэффициентов корреляции рангов Кендэла имеет вид:

 

, (31)

где S=P+Q;

P – число следующих наблюдений для переменной Y с большими значениями рангов;

Q – число следующих наблюдений для переменной Y с меньшими значениями рангов.

 

Значения рангового коэффициента корреляции Кендэла изменяются в пределах от –1 до +1. При >0,5 – связь существенная.

 

4.Коэффициент конкордации (множественный коэффициент ранговой корреляции) применяется для измерения связи между двумя и более признаками при условии, что их можно проранжировать; определяется по формуле:

, (32)

где m – количество признаков, между которыми устанавливается связь;

n – число наблюдений;

S – разность между суммой квадратов рангов и средним квадратом рангов.

Значения коэффициента конкордации изменяются в пределах от 0 до +1;

при значении больше 0,5 можно говорить о тесной связи между признаками.

5.Коэффициенты ассоциации и контингенции применяются для измерения связи между двумя качественными альтернативными признаками.

 

С этой целью строится четырёхклеточная таблица, которая показывает связь между двумя явлениями, каждое из которых обладает альтернативными признаками.

 

Например, обследуют группу населения одного из регионов России в отчётном периоде.

I вопрос – о месте проживания (следует выбрать правильный ответ):

1. Проживаю в городе.

2. Проживаю в сельской местности.

II вопрос – о принадлежности к полу (следует выбрать правильный ответ):

1. Мужчина.

2. Женщина.