Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
6
Добавлен:
19.01.2025
Размер:
7 Мб
Скачать
☆

Построим и сравнив аналогичные графики для пола Female и в среднем, без учета пола. Результат представлен на рисунке 4.2.

Рис. 4.1. Построение диаграммы разброса по Male

Рис. 4.2. Сравнение диаграмм разброса

Чем ближе точки на графике разброса расположены к воображаемой прямой, тем теснее линейная взаимосвязь. MS Excelпозволяет построить на графике разброса линию регрессии: «Конструктор диаграмм» → «Добавить элемент диаграммы» → «Линейный».

51

Судя по результатам, для пола Male значения Х и Y связаны слабо или не связаны вообще; для пола Female, так и для третьего, усредненного случая, – связь отрицательная.

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

являться то,

что здесь рассматриваются только

первые 70 стран (менее

половины).

 

 

По графикам очевидно, что связь между показателями значительно теснее

для пола

Female. Убедимся в этом, рассчитав

значение коэффициента

корреляции Пирсона и проверив его на значимость. Коэффициент корреляции для пары показателей можно рассчитать с помощью функции «КОРРЕЛ». Тем не

менее, в качестве иллюстрации рассчитаем значения

 

вручную по формуле

(4.1).

 

 

 

 

 

∑

̅∙

,

 

(4.1)

 

∑

̅∙∑

 

 

 

где – объем выборки, ̅– выборочное среднее Х, – выборочное среднее Y. Расчеты на примере выборки со значением Male отображены на рисунке

4.3.

Сначала необходимо рассчитать выборочное среднее ̅ 9,87 и 22,62. Затем для каждой страны были найдены отклонения текущих значений и от соответствующих средних. Суммы их преобразований были подставлены в формулу (4.1). Получим значение коэффициента корреляции, равное -0,17.

Рис. 4.3. Расчеты коэффициента парной корреляции Пирсона

52

Положительное значение коэффициента говорит о прямой связи между показателями, отрицательное – об обратной. Чем ближе значение отдельного коэффициента к 1 или -1, тем теснее взаимосвязь между соответствующими переменными. Для интерпретации можно использовать шкалу Чеддока, представленную в таблице 4.1. Общее потребление алкоголя на душу населения (15+) слабо связано с вероятностью смерти в возрасте от 30 до 70 лет от любого из сердечно-сосудистых заболеваний, рака, диабета или хронических респираторных заболеваний (при условии пола Male), направление связи – обратное.

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

Н0:

 

,

 

 

 

 

 

 

 

 

Н1:

 

0.

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

Таблица 4.1 - Шкала Чеддока

 

 

 

 

 

 

 

Значения

 

0,1 … 0,3

0,3 … 0,5

 

0,5 … 0,7

0,7 … 0,9

0,9 … 0,99

 

Оценка

 

слабая

умеренная

 

заметная

высокая

весьма высокая

Для проверки значимости коэффициента корреляции необходимо найти

расчётное .(

̂) (формула 4.2) и критическое ( крит) значения критерия:

крит

 

 

 

 

 

 

 

 

 

; 2

 

 

̂

∙

 

 

 

 

 

(4.2)

 

 

,

 

 

 

 

 

 

 

 

 

 

 

Результаты расчетов – на рисунках 4.4, 4.5. Так как расчетное значение ̂по модулю меньше критического, то нулевую гипотезу следует принять – коэффициент корреляции незначим, связь между показателями (при условии Male) в реальности отсутствует.

Рис. 4.4. Оценка значимости коэффициента корреляции (при условии

Male)

53

Рис. 4.5. Оценка значимости коэффициента корреляции при всех вариантах

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

Построим на этой выборке модель парной линейной регрессии вида (4.3), отражающую зависимость Y от Х.

где

– случайная компонента модели

, ,

–,коэффициенты

модели.

(4.3)

 

Традиционно коэффициенты

модели оцениваются методом наименьших

 

 

 

 

 

 

квадратов (МНК). Посмотрим, как применять его в MS

Excelс помощью

надстройки «Анализ данных».

 

 

 

 

 

 

 

Выбираем «Данные» → «Анализ данных» → «Регрессия». Заполнение

полей в открывшемся окне настройки процедуры представлено на рисунке 4.6 – в основном, оно аналогично предыдущим задачам. Выделяем в качестве входных интервалов столбцы вместе с первой строкой, содержащей заголовки, и отмечаем опцию «Метки».

Рис. 4.6. Настройка расчетов по построению парной линейной регрессии

54

По умолчанию константа в модели присутствует, поэтому опцию «Константа – ноль» оставляем без изменений. Отчет по вычислениям представлен на рисунке 4.7.

Рис. 4.7. Отчет по построению парной линейной регрессии

Значения оценок находятся в столбце «Коэффициенты»: «Y-пересечение»

-это , Х - . Воспользуемся формулой (4.3) и запишем построенную модель

вявном виде:

Можно сделать вывод,

22,39

.

2,11 ∙

 

что в

среднем, без учета влияния прочих

детерминированных факторов, исследуемая вероятность смерти в возрасте от 30 до 70 лет (для женщин) составляет 22,39%. Если общее потребление алкоголя среди женщин станет больше на 1 литр (чистого алкоголя), то вероятность смерти сократится в среднем на 2,11%. Заметим, что вывод об обратной связи мы уже получали в рамках корреляционного анализа.

Столбцы «Нижние 95%» и «Верхние 95%» показывают соответственно нижнюю и верхнюю границы доверительного интервала оценок параметров. Эти значения говорят о том, что, например, реальное значение константы модели с вероятностью 0,95 лежит в интервале от 19,98 до 24,79.

Значения столбца «Стандартная ошибка» для коэффициентов модели равны величине saj, стандартной ошибке параметра регрессии.

В связи с необходимостью получения надежных результатов анализа встает вопрос о критериях оценки качества модели. Эта проблема решается в

55

рамках дисперсионного анализа модели линейной регрессии. Оценим значимость модели с помощью таблицы «Дисперсионный анализ».

С помощью F-критерия Фишера проверяется гипотеза о незначимости регрессии (модели) в целом:

Н0 : а0 = а1 = … = аm = 0

m

Н1 : aj 0

j 0

Расчетное значение статистики Фишера, согласно отчету, равно 33,1 («F»). Критическое значение статистики Фишера равно F ( m; n – m – 1) = F0,05(1; 70 – 1 – 1) =F0,05(1; 68), где m – количество факторов модели. В MS Excelесть встроенная функция для определения этого критического значения. Это функция =F.ОБР.ПХ( ; m; n – m – 1). Для решения примера необходимо задать параметры =F.ОБР.ПХ(0,05;1;68), что даст результат 3,98. Расчетное значение больше критического, следовательно, нулевую гипотезу о незначимости модели в целом принимать не следует.

Альтернативный способ проверки гипотезы – по р-значению F-критерия («Значимость F»). Согласно рисунку 34, р-значение меньше уровня значимости, следовательно, модель значима в целом.

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

Н0 : аj = 0 Н1 : аj ≠ 0

Критическое значение статистики Стьюдента равно t /2(n – m – 1) = t0,025(70

– 1 – 1) =t0,025(68). В MS Excelесть встроенная функция для определения критического двухстороннего значения. Это функция =СТЬЮДЕНТ.ОБР.2Х( ; n – m – 1). Для решения примера необходимо задать параметры =СТЬЮДЕНТ.ОБР.2Х(0,05;68), что даст результат 1,996.

Теперь необходимо сравнить каждое значение столбца «t-статистика», чтобы оценить значимость полученных оценок. Например, расчетное значение t- статистики константы («Yпересечение- », а0) равно 18,56, что больше критического значения. Следовательно, нулевую гипотезу о незначимости константы следует отвергнуть, то есть константа значима.

Отметим, что здесь также можно использовать альтернативный способ проверки гипотез – по р-значению критерия Стьюдента (столбец «Р-значение»).

Врезультате мы показали, что построенная модель имеет смысл (значима

вцелом) и содержит только значимые коэффициенты. Теперь оценим (с

56

помощью таблицы «Регрессионная статистика»), насколько качественно эта модель описывает реальные данные.

Общепринятым показателем качества построения модели множественной

регрессии является коэффициент детерминации R2 . Он показывает долю вариации зависимой переменной, которая объясняется влиянием совокупности независимых переменных модели. Значения коэффициента детерминации

изменяются от 0 до 1. Чем ближе R2 к 1, тем большая доля изменений независимой переменной объясняется влиянием включенных в модель независимых переменных, то есть тем выше качество модели. Согласно рисунку 34, коэффициент детерминации составляет 33%, что говорит об умеренном качестве модели

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

В конце рассчитаем коэффициент эластичности по формуле (4.4). Коэффициент эластичности показывает, на сколько процентов изменится зависимая переменная, если соответствующая независимая изменится на 1%.

Э

∙

 

̅

(4.4)

 

Подставляем в формулу (4.4) значения:

Э

2,11 ∙

, ,

0,92.

Иными словами, на основе построенной модели получаем, что при изменении общего потребления алкоголя на 1%, вероятность смерти от 30 до 70 лет (для женщин) снизится на 0,92%.

Пример реализации лабораторной работы №4

В качестве исходных данных к лабораторной работы 4 рассматриваются два показателя «ожидаемая продолжительность жизни» и «уровень младенческих смертей» по двум группам стран «Развитые страны» и «Развивающиеся страны». Число наблюдений первой группы стран составляет 28 стран. Вторая группа стран содержит 116 наблюдений.

Фрагмент исходных данных представлен на рисунке 4.8.

57

Рис. 4.8. Фрагмент исходных данных для лабораторной работы 4.

Задание 1. Графики разброса по двум группам стран представлены на рисунках 4.9, 4.10. Из диаграмм видно, что можно выделить группы стран, со схожими характеристиками.

Рис. 4.9. Диаграмма разброса по развитым странам

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

58

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

Рис. 4.10. Диаграмма разброса по развивающимся странам

Задание 2.

Расчеты коэффициента корреляции Пирсона представлены в таблице 4.2. Таблица 4.2 -Расчет коэффициента корреляции по формулам вручную

задания 2. Лабораторная работа 4.

 

 

 

 

Развитые страны

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,38

 

-67,46

 

 

-93,00

 

 

1,90

 

4551,43

0,08

 

-7,46

 

-0,59

 

 

 

0,01

 

55,72

7,68

 

2,54

 

19,47

 

 

58,96

 

6,43

-7,02

 

64,54

 

 

-453,13

 

49,30

 

4164,86

-3,52

 

23,54

 

-82,88

 

 

12,40

 

553,93

1,68

 

-20,46

 

 

-34,35

 

 

2,82

 

418,79

2,68

 

-0,46

 

-1,24

 

 

 

7,17

 

0,22

7,68

 

-4,46

 

-34,28

 

 

58,96

 

19,93

-5,72

 

63,54

 

 

-363,52

 

32,73

 

4036,79

1,18

 

-24,46

 

 

-28,83

 

 

1,39

 

598,50

-0,12

 

-7,46

 

0,91

 

 

 

0,01

 

55,72

59

Продолжение таблицы 4.2

1,18

 

-16,46

 

 

-19,40

 

1,39

 

271,07

2,18

 

-16,46

 

 

-35,87

 

4,75

 

271,07

-6,92

 

82,54

 

 

-571,27

47,91

 

6812,14

-7,92

 

95,54

 

 

-756,78

62,75

 

9127,07

0,38

 

-8,46

 

-3,20

 

 

0,14

 

71,64

0,08

 

-18,46

 

-1,45

 

 

0,01

 

340,93

0,38

 

-15,46

 

-5,85

 

 

0,14

 

239,14

0,28

 

-12,46

 

-3,47

 

 

0,08

 

155,36

-4,02

 

-61,46

 

 

247,17

16,17

 

3777,86

7,68

 

4,54

 

34,83

 

58,96

 

20,57

-6,52

 

61,54

 

 

-401,30

42,53

 

3786,64

5,68

 

2,54

 

14,40

 

32,25

 

6,43

1,28

 

-15,46

 

 

-19,77

 

1,63

 

239,14

0,98

 

-19,46

 

 

-19,05

 

0,96

 

378,86

1,88

 

-22,46

 

 

-42,20

 

3,53

 

504,64

-0,32

 

-2,46

 

0,79

 

 

0,10

 

6,07

-2,22

 

-59,46

 

 

132,10

 

4,93

 

3536,00

Сумма

 

 

 

 

-2522

 

 

504

44007

Коэффициент корреляции Пирсона примет значение

 

2522

0,536

√

504 44007

Задание 3. Матрица

корреляций,

полученная с помощью Excel,

представлена для развитых стран на рисунке 4.11, для развивающихся стран на рисунке 4.12.

Рис. 4.11. Матрица корреляций по группе развитых стран

Рис. 4.12. Матрица корреляций по группе развивающихся стран

60