Контрольная работа Линейные уравнения парной и множественной регрессии

Работа добавлена на сайт bukvasha.net: 2015-10-25

Поможем написать учебную работу

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

Предоплата всего

от 25%

Подписываем

договор

Имя

Выберите тип работы:

Принимаю Политику конфиденциальности

Скидка 25% при заказе до 25.11.2024

НЕГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ

ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

"ВОЛГОГРАДСКИЙ ИНСТИТУТ БИЗНЕСА"

Кафедра

Математики и естественных наук
Домашняя контрольная работа

Дисциплина

Эконометрика

Тема: Линейные уравнения парной регрессии
Студента (ки)

Иванова Ивана Ивановича
Волгоград 2010

Задача№ 1
По данным приведенным в таблице:

1)                построить линейное уравнение парной регрессии y на x;

2)                рассчитать линейный коэффициент парной корреляции и оценить тесноту связи;

3)                оценить статистическую значимость параметров регрессии и корреляции, используя F-статистику, t-статистику Стьюдента и путем расчета доверительных интервалов каждого из показателей;

4)                вычислить прогнозное значение y при прогнозном значении x, составляющем 108% от среднего уровня.

5)                оценить точность прогноза, рассчитав ошибку прогноза и его доверительный интервал;

6)                полученные результаты изобразить графически и привести экономическое обоснование.
Таблица №1

По территориям Центрального района известны данные за 1995 г.

Район	Средний размер назначенных ежемесячных пенсий, тыс.руб., y	Прожиточный минимум в среднем на одного пенсионера в месяц, тыс.руб., х
Брянская обл.	240	178
Владимирская обл.	226	202
Ивановская обл.	221	197
Калужская обл.	226	201
Костромская обл.	220	189
Московская обл.	237	215
Орловская обл.	232	166
Рязанская обл.	215	199
Смоленская обл.	220	180
Тульская обл.	231	186
Ярославская обл.	229	250

x_i	178	202	197	201	189	215	166	199	180	186	250
y_i	240	226	221	226	220	237	232	215	220	231	229

Х	Y
178	240
202	226
197	221
201	226
189	220
215	237
166	232
199	215
180	220
186	231
250	229

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

Таблица № 4 Параметры (коэффициенты) уравнения регрессии

Коэффициенты

Y-пересечение

227,7117993

Переменная X 1

-0,003619876

На основании этих данных запишем уравнение регрессии: .

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

Столбец 1

Столбец 2

Столбец 1

1

Столбец 2

-0,010473453

1

Для оценки качества уравнения регрессии в целом необходимо проверить статистическую значимость индекса детерминации: проверяется нулевая гипотеза , используется .
Таблица №6

Регрессионная статистика

R-квадрат

0,000109693

.

Т.к. Значение детерминации R-квадрат имеет малое значение, которое менее 1%, то дальнейшее решение не имеет смысла, т.к. вероятность того что прогноз будет верным меньше 1%.

Задача №2

Используя данные, приведенные в таблице: построить линейное уравнение множественной регрессии;

1)                оценить значимость параметров данного уравнения и построить доверительные интервалы для каждого из параметров, оценить значимость уравнения в целом, пояснить экономический смысл полученных результатов;

2)                рассчитать линейные коэффициенты частной корреляции и коэффициент множественной детерминации, сравнить их с линейными коэффициентами парной корреляции, пояснить различия между ними;

3)                вычислить прогнозное значение y при уменьшении вектора x на 6 % от максимального уровня, оценить ошибку прогноза и построить доверительный интервал прогноза;
Таблица №5

номер наблюдения, i	Накопления семьи, Y (y.e.)	Доход семьи, X_{1 (} _y. _e.)	Расходы на питание, X₂ ( y. e.)
1	2	20	5
2	6	27	6
3	7	26	7
4	5	19	5
5	4	15	5
6	2	15	5
7	7	28	10
8	6	24	7
9	4	14	6
10	5	21	7
11	5	20	10
12	3	18	6

Таблица №6 Параметры (коэффициенты) уравнения регрессии

	Коэффициенты
Y-пересечение	-1,767785782
x1	0,232792618
x2	0,24953991

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

На основании этих данных запишем уравнение регрессии:

.
Таблица №7 Регрессионная статистика

R-квадрат	0,663668925
Нормированный R-квадрат	0,588928686

! Параметр R-квадрат, представляет собой квадрат коэффициента корреляции r_xy2 и называется коэффициентом детерминации. Величина данного коэффициента характеризует долю дисперсии зависимой переменной y, объясненную регрессией (объясняющей переменной x). Соответственно величина 1 - r_xy2 характеризует долю дисперсии переменной y, вызванную влиянием всех остальных, неучтенных в эконометрической модели объясняющих переменных. Доля всех неучтенных в полученной эконометрической модели объясняющих переменных приблизительно составляет: 0,663668, или 66,3%.

Находим, что численное значение

, а скорректированный (нормированный, исправленный) коэффициент детерминации равен

1) Для оценки качества уравнения регрессии в целом необходимо проверить статистическую значимость индекса детерминации

: проверяется нулевая гипотеза

, используется

.

Наблюдаемое значение критерия

и оценку его значимости находим в Таблице №8
Таблица №8 Дисперсионный анализ:

F	Значимость F
8,87967358	0,007420813

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

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

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

2) Для статистической оценки значимости коэффициентов регрессии () используем статистику Стьюдента.

Проверяется нулевая гипотеза .

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

Таблица №9

t-статистика

P-Значение

-1,127971079

0,28850322

2,838964459

0,01943598

1,130728736

0,28740002

В этой же таблице находим границы доверительных интервалов для каждого из параметров:

Нижние 95%

Верхние 95%

-5,313097658

1,777526094

0,047297697

0,418287538

-0,249694323

0,748774142

3. Значения парных коэффициентов корреляции найдем из соответствующей матрицы.
Таблица №10 Корреляционная матрица

y

x1

x2

y

1

x1

0,784786247

1

x2

0,60206001

0,531178469

1

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

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

Частные коэффициенты корреляции найдем по формулам
,

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

4. Рассчитаем прогнозное значение результата, если прогнозные значения факторов составляют 110% их максимального значения. Найдем прогнозные значения факторов и подставим их в полученное уравнение регрессии.

По условию прогнозные значения составляют 110% их максимального значения.

Таблица №11

maxX1

maxX2

28

10

Далее вычисляем прогнозные значения факторов: . Затем, подставив эти значения в уравнение регрессии, получим прогнозное (предсказанное) значение фактора . Доверительный интервал прогноза оценивается формулой: , где - ошибка прогноза,стандартная ошибка регрессии.
Таблица №12

Стандартная ошибка

1,104878833

;

- коэффициент Стьюдента, который в данном случае имеет смысл кратности случайной (стандартной) ошибки прогноза ;

- число, которое получим в результате операций над матрицами:
-
матрица значений факторных переменных ,

транспонированная матрица ;

- произведение матриц ;

- матрица, обратная к матрице ;

- матрица прогнозных значений факторов;

- транспонированная матрица прогнозов.

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

Максимальную ошибку прогноза =11,07714043: 1) нижняя граница прогноза =44,92285957, 2) верхнюю границу прогноза =67,07714043. Интервал прогнозных значений результативного признака
=>

Задача № 3

Используя данные, представленные в таблице проверить наличие гетероскедастичности, применяя тест Голдфельда-Квандта.

Таблица№13. Данные

Страна

Индекс человеческого развития, У

Расходы на конечное потребление в текущих ценах, % к ВВП, Х

Австрия

0,904

75,5

Австралия

0,922

78,5

Англия

0,918

84,4

Белоруссия

0,763

78,4

Бельгия

0,923

77,7

Германия

0,906

75,9

Дания

0,905

76,0

Индия

0,545

67,5

Испания

0,894

78,2

Италия

0,900

78,1

Канада

0,932

78,6

Казахстан

0,740

84,0

Китай

0,701

59,2

Латвия

0,744

90,2

Нидерланды

0,921

72,8

Норвегия

0,927

67,7

Польша

0,802

82,6

Россия

0,747

74,4

США

0,927

83,3

Украина

0,721

83,7

Финляндия

0,913

73,8

Франция

0,918

79,2

Чехия

0,833

71,5

Швейцария

0,914

75,3

Швеция

0,923

79,0

1) Найдем параметры линейного уравнения множественной регрессии и значения остатков.

Определим остаточные суммы квадратов и , то есть суммы квадратов остатков регрессии по "урезанным выборкам".
Таблица№14

№

Y

X

Yp

ei

(ei) ^2

1

0,932

78,6

77,90431365

0,695686352

0,483979501

2

0,927

67,7

77,85057558

-10,15057558

103,0341846

3

0,927

83,3

77,85057558

5,44942442

29,69622651

4

0,923

77,7

77,80758513

-0,107585125

0,011574559

5

0,923

79,0

77,80758513

1, 192414875

1,421853234

6

0,922

78,5

77,79683751

0,703162488

0,494437485

7

0,921

72,8

77,7860899

-4,986089898

24,86109247

8

0,918

84,4

77,75384706

6,646152943

44,17134894

S1

9

0,918

79,2

77,75384706

1,446152943

2,091358334

206,2660556

10

0,914

75,3

77,7108566

-2,410856603

5,812229559

11

0,913

73,8

77,70010899

-3,900108989

15,21085013

12

0,906

75,9

77,62487569

-1,724875694

2,975196159

13

0,905

76,0

77,61412808

-1,61412808

2,60540946

14

0,904

75,5

77,60338047

-2,103380467

4,424209388

15

0,900

78,1

77,56039001

0,539609988

0,291178939

16

0,894

78,2

77,49590433

0,704095669

0,495750712

17

0,833

71,5

76,8402999

-5,3402999

28,51880303

18

0,802

82,6

76,50712388

6,092876121

37,12313943

19

0,763

78,4

76,08796695

2,312033052

5,345496834

20

0,747

74,4

75,91600513

-1,51600513

2,298271555

21

0,744

90,2

75,88376229

14,31623771

204,9546622

22

0,740

84,0

75,84077183

8,159228165

66,57300425

23

0,721

83,7

75,63656718

8,063432824

65,0189489

24

0,701

59,2

75,4216149

-16,2216149

263,1407901

S2

25

0,545

67,5

73,74498718

-6,244987181

38,99986489

743,7878055

1) Находим наблюдаемое значение критерия . По условию задачи . Из таблицы значений Фишера находим, что

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

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

Задача № 4

По данным таблицы построить уравнение регрессии, выявить наличие автокорреляции остатков, используя критерий Дарбина - Уотсона, и проанализировать пригодность полученного уравнения для построения прогнозов.
Таблица №15

Год	Выпуск продукции в США в среднем за 1 час, % к уровню 1982 г., Х	Среднечасовая заработная плата в экономике США, в сопоставимых ценах 1982 г., Y
1960	65,6	6,79
1961	68,1	6,88
1962	73,3	7,07
1963	76,5	7,17
1964	78,6	7,33
1965	81,0	7,52
1966	83,0	7,62
1967	85,4	7,72
1968	85,9	7,89
1969	85,9	7,98
1970	87,0	8,03
1971	90,2	8,21
1972	92,6	8,53
1973	95,0	8,55
1974	93,3	8,28
1975	95,5	8,12

Найдем параметры линейного уравнения множественной регрессии и значения остатков.

Дополним таблицу данных столбцами "", "Квадрат разности остатков " и "Квадрат остатка " и заполним их.
Таблица №16

Y

X

Yi

et

et-1

(et-et-1) ^2

et^2

6,79

65,6

6,667235239

0,122765

0,015071

6,88

68,1

6,815288112

0,064712

0,122765

0,003370136

0,004188

7,07

73,3

7,123238088

-0,05324

0,064712

0,013912197

0,002834

7,17

76,5

7,312745766

-0,14275

-0,05324

0,008011624

0,020376

7,33

78,6

7,437110179

-0,10711

-0,14275

0,001269895

0,011473

7,52

81,0

7,579240937

-0,05924

-0,10711

0,002291464

0,003509

7,62

83,0

7,697683236

-0,07768

-0,05924

0,000340118

0,006035

7,72

85,4

7,839813994

-0,11981

-0,07768

0,001775001

0,014355

7,89

85,9

7,869424568

0,020575

-0,11981

0,019709191

0,000423

7,98

85,9

7,869424568

0,110575

0,020575

0,008100000

0,012227

8,03

87,0

7,934567833

0,095432

0,110575

0,000229318

0,009107

8,21

90,2

8,12407551

0,085924

0,095432

0,000090396

0,007383

8,53

92,6

8,266206268

0,263794

0,085924

0,031637467

0,069587

8,55

95,0

8,408337026

0,141663

0,263794

0,014915922

0,020068

8,28

93,3

8,307661073

-0,02766

0,141663

0,028670633

0,000765

8,12

95,5

8,437947601

-0,31795

-0,02766

0,084266268

0,101091

Суммы

0,218589631

0,298494

По формуле вычислим значение статистики :

Так как , то значение статистики

равно .

По таблице критических точек Дарбина Уотсона определим значения критерия Дарбина-Уотсона (нижнее) и (верхнее) для заданного числа наблюдений , числа независимых переменных модели и уровня значимости . Итак, находим, что , .

По этим значениям числовой промежуток разбиваем на пять отрезков:
,

,

,

,

.
На основании выполненных расчетов находим, что наблюдаемое значение статистики принадлежит первому интервалу.

Вывод: существует отрицательная автокорреляция, то есть гипотеза отклоняется и с вероятностью принимается гипотеза .

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

Задача № 5

В таблице приводятся данные о динамике выпуска продукции Финляндии (млн. долл.).
Таблица №17

Год	Выпуск продукции, yt млн.долл.
1989	23 298
1990	26 570
1991	23 080
1992	29 800
1993	28 440
1994	29 658
1995	39 573
1996	38 435
1997	39 002
1998	39 020
1999	40 012
2000	41 005
2001	39 080
2002	42 680

Задание:

1.                 Постройте график временного ряда.

2.                 Сделайте вывод о присутствии или отсутствии тренда при доверительной вероятности 0,95.

3.                 Найдите среднее значение, среднеквадратическое отклонение и коэффициенты автокорреляции (для лагов

) заданного ВР.

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

;

5. Найдите уравнение тренда ВР

, предполагая, что он линейный, и проверьте его значимость на уровне

.

6. Дайте точечный и интервальный (с надежностью 0,95) прогнозы индивидуального значения выпуска продукции на 2003 год.
Таблица №18

Год	t	Выпуск продукции, yt млн.долл.
1989	1	23 298
1990	2	26 570
1991	3	23 080
1992	4	29 800
1993	5	28 440
1994	6	29 658
1995	7	39 573
1996	8	38 435
1997	9	39 002
1998	10	39 020
1999	11	40 012
2000	12	41 005
2001	13	39 080
2002	14	42 680

2. Для обнаружения тенденции в данном ВР воспользуемся критерием "восходящих и нисходящих" серий.

Критерий "восходящих и нисходящих" серий

1) Для исследуемого ВР определяется последовательность знаков, исходя из условий: (+), если

, (-), если

.

При этом, если последующее наблюдение равно предыдущему, то учитывается только одно наблюдение.

2) Подсчитывается число серий . Под серией понимается последовательность подряд расположенных плюсов или минусов, причем один плюс или один минус считается серией.

3) Определяется протяженность самой длинной серии .

4) Значение находят из следующей таблицы:
Таблица №25

Длина ряда,
Значение	5	6	7

5) Если нарушается хотя бы одно из следующих неравенств, то гипотеза об отсутствии тренда отвергается с доверительной вероятностью 0,95

Определим последовательность знаков:

Таблица №19

t	Выпуск продукции, yt млн.долл.
1	23 298
2	26 570	+
3	23 080	-
4	29 800	+
5	28 440	-
6	29 658	+
7	39 573	+
8	38 435	-
9	39 002	+
10	39 020	+
11	40 012	+
12	41 005	+
13	39 080	-
14	42 680	+

Определим число серий

. Определим протяженность самой длинной серии : .

, так как

. Проверим выполнение неравенств:

Вывод: второе неравенство не выполняются, следовательно, тренд (тенденция) в динамике выпуска продукции имеется на уровне значимости 0,05. Среднее значение

. Среднее значение

. Вычислим коэффициенты автокорреляции первого и второго порядков, то есть для лагов

. Подготовим данные для вычисления коэффициентов автокорреляции первого и второго порядков. Дополним таблицу данных двумя столбцами

.

Таблица №20

t	Yt	Yt-1	Yt-2
1	23 298
2	26 570	23 298
3	23 080	26 570	23 298
4	29 800	23 080	26 570
5	28 440	29 800	23 080
6	29 658	28 440	29 800
7	39 573	29 658	28 440
8	38 435	39 573	29 658
9	39 002	38 435	39 573
10	39 020	39 002	38 435
11	40 012	39 020	39 002
12	41 005	40 012	39 020
13	39 080	41 005	40 012
14	42 680	39 080	41 005

.

.
Вывод:

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

2) исследуемый ряд содержит только тенденцию, так как наиболее высоким оказался коэффициент автокорреляции первого порядка (0,85>0,83).

Скользящие средние найдем по формуле: , здесь . При

Вычисляем:

и так далее.

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

t

yi

yt

1

23 298

2

26 570

24 315,76

3

23 080

26 483,07

4

29 800

27 106,40

5

28 440

29 299,04

6

29 658

32 556,67

7

39 573

35 888,31

8

38 435

39 002,94

9

39 002

38 818,61

10

39 020

39 344,27

11

40 012

40 011,93

12

41 005

40 031,93

13

39 080

40 921,26

14

42 680

Таблица № Параметры (коэффициенты) уравнения тренда.
Таблица №22

Коэффициенты

Y-пересечение

22686,54945

t

1543,250549

Анализ данных таблицы Дисперсионного анализа показывает, что получено статистически значимое уравнение, так как наблюдаемое значение

, равное 52,785, превышает его табличное значение

. Вывод: Таким образом, параметры уравнения тренда статистически значимы на уровне

: уравнение тренда можно использовать для прогноза.

Сделаем точечный и интервальный (с надежностью 0,95) прогнозы среднего и индивидуального значений прогнозов на 2003 год.

Определим точечный прогноз

Вычислим интервальный прогноз:

Так как тренд является прямой, то доверительный интервал можно представить в виде:

.

Здесь стандартная ошибка предсказания по линии тренда

вычисляется по формуле:

,
здесь величина

является стандартной ошибкой регрессии, и ее значение находится в таблице Регрессионная статистика
Таблица №23

Стандартная ошибка

1637,180026

кратность ошибки (надежность) находят по таблице значений критерия Стьюдента;

уровень значимости;

число степеней свободы.

Итак, по условию задачи имеем:

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

необходимо вычислить

и сумму

.
Таблица № 24

t	yt	(t1-tcr) ^2
1	23 298	42,25
2	26 570	30,25
3	23 080	20,25
4	29 800	12,25
5	28 440	6,25
6	29 658	2,25
7	39 573	0,25
8	38 435	0,25
9	39 002	2,25
10	39 020	6,25
11	40 012	12,25
12	41 005	20,25
13	39 080	30,25
14	42 680	42,25
7,5	Сумма	227,5

Вычисляем

(млн. долл.)

По таблице значений критерия Стьюдента найдем

Максимальная ошибка прогноза будет равна:

(млн. долл.).
Нижняя граница прогноза имеет значение

(млн. долл.)

Верхняя граница прогноза имеет значение

(млн. долл.)

Вывод:

1) значение выпуска продукции Финляндии в 2003 составит 20111,2 млн. долл.

2) с надежностью 0,95 данное значение будет находиться в интервале

Bukvasha