Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

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

.pdf
Скачиваний:
125
Добавлен:
12.03.2016
Размер:
1.94 Mб
Скачать

Решение с помощью Excel

1.Ввод исходных данных по объему продаж цветных карандашей.

2.Расчет абсолютных приростов y по исходным данным y

иоформление результатов расчетов в виде табл. 4.1.2.2.

Та б л и ц а 4.1.2.2

t

y

y

t

y

y

1

801,13

 

9

1582,65

137,46

2

859,23

58,1

10

1722,41

139,76

3

938,27

79,04

11

1880,44

158,03

4

1015,27

77

12

2045,77

165,33

5

1106,56

91,29

13

2219,39

173,62

6

1211,15

104,59

14

2404,3

184,91

7

1326,03

114,88

15

2589,6

185,3

8

1445,19

119,16

 

 

 

3. Определение типа роста по «Линейчатой» диаграмме, построенной для абсолютных приростов (рис. 4.1.2.1).

Рис. 4.1.2.1. Абсолютные приросты продаж цветных карандашей

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

y

t

b b t b t2

и

y

t

b bt.

 

0 1

2

 

 

0 1

4. Подготовка исходных данных для построения указанных моделей и оформление их в виде табл. 4.1.2.3.

61

Т а б л и ц а 4.1.2. 3

t

 

t2

 

y

ln y

t

 

t2

 

y

ln y

 

1

 

1

801,13

6,69

 

9

 

81

1582,65

7,37

 

2

 

4

859,23

6,76

 

10

 

100

1722,41

7,45

 

3

 

9

938,27

6,84

 

11

 

121

1880,44

7,54

 

4

 

16

1015,27

6,92

 

12

 

144

2045,77

7,62

 

5

 

25

1106,56

7,01

 

13

 

169

2219,39

7,7

 

6

 

36

1211,15

7,1

 

14

 

196

2404,3

7,79

 

7

 

49

1326,03

7,19

 

15

 

225

2589,6

7,86

 

8

 

64

1445,19

7,28

 

 

 

 

 

 

5. Нахождение коэффициентов трендовых моделей с помощью «Пакета анализа» Excel (см. Вывод итогов 4.1.2.1 и Вывод итогов 4.1.2.2).

Таким образом, в рассматриваемом случае парабола имеет вид

yt 747,62 46,86t 5,088t2.

Поскольку

b elnb0

2,7186,58 727,42;

b elnb1

2,7180,08 1,08,

0

 

1

 

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

yt 727,42 1,08t .

ВЫВОД ИТОГОВ 4.1.2.1

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

Множественный R

0,9999875

 

 

 

 

 

R-квадрат

0,999975

 

 

 

 

 

Нормированный R-

 

 

 

 

 

 

квадрат

0,9999708

 

 

 

 

 

Стандартная

 

 

 

 

 

 

ошибка

3,1338086

 

 

 

 

 

Наблюдения

15

 

 

 

 

 

 

 

 

 

 

 

 

Дисперсионный анализ

 

 

 

 

 

 

df

SS

MS

F

Значимость F

Регрессия

2

4714355

2357177

240019,9

2,44E-28

 

Остаток

12

117,8491

9,820757

 

 

 

Итого

14

4714473

 

 

 

 

 

 

 

 

 

 

 

 

Коэффи-

Стандарт-

t-статис-

P-Зна-

Нижние

Верхние

 

циенты

ная ошибка

тика

чение

95%

95%

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

747,61277

2,791389

267,8282

4,94E-24

741,5309

753,6947

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

46,860806

0,802812

58,37084

4,22E-16

45,11163

48,60998

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

5,0886304

0,048791

104,2937

4,04E-19

4,982323

5,194938

62

ВЫВОД ИТОГОВ 4.1.2.2

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

Множественный R

 

0,9998213

 

 

 

 

 

R-квадрат

 

0,9996426

 

 

 

 

 

Нормированный R-

 

 

 

 

 

 

 

квадрат

 

0,9996151

 

 

 

 

 

Стандартная

 

 

 

 

 

 

 

ошибка

 

0,0075111

 

 

 

 

 

Наблюдения

 

15

 

 

 

 

 

 

 

 

 

 

 

 

 

Дисперсионный анализ

 

 

 

 

 

 

 

 

df

SS

MS

F

Значимость F

Регрессия

 

1

2,051269

2,051269

36358,94

8,56E-24

 

Остаток

 

13

0,000733

5,64E-05

 

 

 

Итого

 

14

2,052002

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффи-

Стандарт-

t-статис-

 

Нижние

Верхние

 

 

циенты

ная ошибка

тика

P-Зна-чение

95%

95%

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

 

6,5895102

0,004081

1614,586

7,45E-36

6,580693

6,598327

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

 

0,0855918

0,000449

190,6802

8,56E-24

0,084622

0,086562

6. Вычисление расчетных значений объема продаж по построенным моделям и оформление результатов в виде табл. 4.1.2.4.

7.Расчет отклонений расчетных значений и фактических, их квадратов и средних квадратических отклонений. Оформление результатов в виде табл. 4.1.2.5.

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

Т а б л и ц а 4.1.2.4

t

y

yˆпарабола

yˆпоказательная

t

y

yˆпарабола

yˆпоказательная

1

801,13

799,56

792,43

9

1582,65

1581,54

1571,58

2

859,23

861,69

863,24

10

1722,41

1725,08

1712,02

3

938,27

933,99

940,38

11

1880,44

1878,81

1865,01

4

1015,27

1016,47

1024,42

12

2045,77

2042,71

2031,67

5

1106,56

1109,13

1115,96

13

2219,39

2216,78

2213,22

6

1211,15

1211,97

1215,68

14

2404,3

2401,04

2411,00

7

1326,03

1324,98

1324,32

15

2589,6

2595,47

2626,45

8

1445,19

1448,17

1442,66

 

 

 

 

8.Подготовка данных для расчета числителя критерия Дарбина – Уотсона в виде табл. 4.1.2.6.

9.Окончательный расчет критерия Дарбина – Уотсона

d252,00 2,138. 117,84

63

Т а б л и ц а 4.1.2.5

 

 

 

 

 

yˆ

 

2

 

 

 

2

y

yˆ

 

 

 

y yˆ

 

 

ˆ

парабола

 

показательная

парабола

 

 

 

y

y показательная

 

 

 

 

 

 

 

801,13

 

799,56

 

 

792,43

2,4580

 

 

 

 

75,7191

859,23

 

861,69

 

 

863,24

6,0462

 

 

 

 

16,0879

938,27

 

933,99

 

 

940,38

18,2939

 

 

 

 

4,4587

1015,27

1016,47

 

 

1024,42

1,4498

 

 

 

 

83,6409

1106,56

1109,13

 

 

1115,96

6,6181

 

 

 

 

88,3400

1211,15

1211,97

 

 

1215,68

0,6696

 

 

 

 

20,5462

1326,03

1324,98

 

 

1324,32

1,0998

 

 

 

 

2,9305

1445,19

1448,17

 

 

1442,66

8,8897

 

 

 

 

6,3944

1582,65

1581,54

 

 

1571,58

1,2341

 

 

 

 

122,5496

1722,41

1725,08

 

 

1712,02

7,1496

 

 

 

 

107,9803

1880,44

1878,81

 

 

1865,01

2,6702

 

 

 

 

238,1671

2045,77

2042,71

 

 

2031,67

9,3929

 

 

 

 

198,8850

2219,39

2216,78

 

 

2213,22

6,8028

 

 

 

 

38,0646

2404,30

2401,04

 

 

2411,00

10,6562

 

 

 

 

44,8528

2589,60

2595,47

 

 

2626,45

34,4182

 

 

 

 

1357,7565

Сумма

квадратов отклонений

117,8491

 

 

 

 

2406,3736

Средний квадрат отклонений

7,8566

 

 

 

 

160,4249

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

2,8030

 

 

 

 

12,6659

 

 

 

 

 

 

 

Т а б л и ц а 4.1.2.6

t

 

y

 

 

yˆпарабола

y yˆ парабола

 

 

 

et et 1 2

1

 

 

801,13

 

 

799,56

1,5678

 

 

 

 

2

 

 

859,23

 

 

861,69

-2,4589

 

 

 

16,2143

3

 

 

938,27

 

 

933,99

4,2771

 

 

 

45,3743

4

 

 

1015,27

 

 

1016,47

-1,2041

 

 

 

30,0438

5

 

 

1106,56

 

 

1109,13

-2,5726

 

 

 

1,8727

6

 

 

1211,15

 

 

1211,97

-0,8183

 

 

 

3,0774

7

 

 

1326,03

 

 

1324,98

1,0487

 

 

 

3,4857

8

 

 

1445,19

 

 

1448,17

-2,9816

 

 

 

16,2430

9

 

 

1582,65

 

 

1581,54

1,1109

 

 

 

16,7484

10

 

 

1722,41

 

 

1725,08

-2,6739

 

 

 

14,3246

11

 

 

1880,44

 

 

1878,81

1,6341

 

 

 

18,5585

12

 

 

2045,77

 

 

2042,71

3,0648

 

 

 

2,0469

13

 

 

2219,39

 

 

2216,78

2,6082

 

 

 

0,2085

14

 

 

2404,3

 

 

2401,04

3,2644

 

 

 

0,4306

15

 

 

2589,6

 

 

2595,47

-5,8667

 

 

 

83,3768

 

 

 

 

Числитель критерия

Дарбина – Уотсона

 

 

 

252,0053

При 5%-м уровне значимости для 15 наблюдений и двух переменных в модели нижняя граница критерия d1 0,95, а верхняя

d2 1,54. Так как d 2, то с критическим значениям сравнивается не сам коэффициент d , а 4 d , равный 1,862. Таким обра-

64

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

10.Расчет прогнозных оценок и их доверительных границ, учитывая, что t0,95 12 2,179. Оформление результатов в виде табл. 4.1.2.7.

 

 

 

Т а б л и ц а 4.1.2.7

 

 

 

 

Месяц

Прогнозные оценки

Нижняя граница

Верхняя граница

объема продаж

прогнозной оценки

прогнозной оценки

 

16

2800,08

2792,18

2807,97

17

3014,86

3006,53

3023,19

18

3239,82

3230,98

3248,67

4.1.3. Задание для самостоятельной работы

Задание 4.1.3.1. По данным табл. 4.1.2.1 для каждого товара, кроме цветных карандашей, определить тип роста временного ряда, отражающего динамику соответствующего объема продаж. Применяя среднеквадратический критерий, определить среди функций, используемых для моделирования данного типа роста, наиболее подходящую для прогнозных расчетов. С помощью критерия Дарбина – Уотсона проверить адекватность прогнозной модели и получить точечные и интервальные прогнозы на четыре периода. Построить «точечный» график для фактических и расчетных значений, включая прогнозные.

4.2. МОДЕЛИРОВАНИЕ СЕЗОННЫХ КОЛЕБАНИЙ

4.2.1. Расчетные формулы

4.2.1.1. Модель с аддитивной сезонной составляющей yt f(t) S(t) t ,

где f t – тренд; S t – сезонная составляющая; t – случайная компонента.

4.2.1.2. Модель с мультипликативной сезонной составляющей:

yt f(t) S(t) t .

4.2.1.3. Модель временного ряда с циклическими колебаниями периодичностью k:

yt b0 b1t c1x1 c2x2 ck 1xk 1 t ,

65

где

1для j nk,

n 0,1,2,

xj

 

.

 

0

для всех остальных случаев

4.2.2. Решение типовых задач

Задание 4.2.2.1. Некий предприниматель, оптовый поставщик фруктов на рынки г. Воронежа, желает спланировать свою деятельность на 2004г. таким образом, чтобы получить максимум прибыли. Большую часть в общем объеме его продаж занимают груши. Данные по этому виду фруктов представлены в табл. 4.2.2.1. Для того чтобы получить прогнозные оценки объема продаж груш на требуемый период, ему посоветовали построить модель с аддитивной сезонной компонентой.

Решение с помощью Excel

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

 

 

 

 

 

Т а б л и ц а 4.2.2.1

Год

Сезон

Объем продаж, т.

Год

Сезон

 

Объем продаж, т.

 

 

 

 

 

 

 

 

зима

3836

 

зима

 

5007

2000

весна

1831

2002

весна

 

2096

лето

11588

лето

 

14744

 

 

 

 

осень

30945

 

осень

 

39328

 

зима

4447

 

зима

 

5437

2001

весна

2067

2003

весна

 

2358

лето

12938

лето

 

16644

 

 

 

 

осень

34892

 

осень

 

44381

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

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

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

2.3.Вычисление сезонной компоненты в виде разницы фактических значений и центрированных скользящих средних.

2.4.Оформление результатов расчетов в виде табл. 4.2.2.2.

66

Т а б л и ц а 4.2.2.2

Год

Сезон

Объем

Скользящее

Центрированное

Оценка

среднее

скользящее

сезонной

продаж

 

 

за 4 периода

среднее

компоненты

 

 

 

 

зима

3836

-

-

-

2000

весна

1831

12050,22

-

-

лето

11588

12202,94

12126,58

-538,48

 

 

осень

30945

12261,76

12232,35

18713,00

 

зима

4447

12599,16

12430,46

-7983,58

2001

весна

2067

13585,73

13092,45

-11025,72

лето

12938

13725,87

13655,80

-718,14

 

 

осень

34892

13733,26

13729,57

21162,10

 

зима

5007

14184,86

13959,06

-8951,61

2002

весна

2096

15293,95

14739,40

-12643,15

лето

14744

15401,30

15347,62

-603,56

 

 

осень

39328

15466,75

15434,03

23894,00

 

зима

5437

15941,81

15704,28

-10267,41

2003

весна

2358

17205,05

16573,43

-14215,40

лето

16644

-

-

-

 

 

осень

44381

-

-

-

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

3.1.Формирование из оценок сезонной компоненты, полученных в предыдущем пункте, табл. 4.2.2.3, удобной для расчета средних значений этой же компоненты.

3.2.Расчет итоговых значений сезонной компоненты.

3.3.Определение средних значений итоговой компоненты.

Та б л и ц а 4.2.2.3

Показатели

Год

 

Сезон

 

Зима

Весна

Лето

Осень

 

 

 

2000

-

-

-538,48

18713,00

 

2001

-7983,58

-11025,72

-718,14

21162,10

 

2002

-8951,61

-12643,15

-603,56

23894,00

 

2003

-10267,41

-14215,40

-

-

Итого за сезон

 

-27202,59

-37884,28

-1321,69

45056,09

Средняя оценка

 

-9067,53

-12628,09

-440,56

15018,70

сезонной компоненты

 

 

 

 

 

 

Скорректированная

 

-7288,16

-10848,72

1338,81

16798,07

сезонная компонента

 

 

 

 

 

 

3.4. Определение корректирующего коэффициента

k( 9067,53 12628,09 440,56 15018,70):4 1779,37.

3.5.Расчет скорректированных значений сезонной компоненты путем вычитания корректирующего коэффици-

67

ента из средних оценок сезонной компоненты (сумма скорректированных значений равна 0).

4. Вычисление основных составляющих сезонной модели

(f,S, ).

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

4.2.Построение по данным элиминированного временного ряда трендовой модели с помощью МНК. Оформление результатов моделирования в виде табл. 4.2.2.4.

Та б л и ц а 4.2.2.4

Результаты моделирования

Коэффициент регрессии

554,41

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

209,77

Множественный R

0,58

Константа

9821,28

Число наблюдений

16

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

14

F - критерий

6,99

4.3.Получение расчетных значений по трендовой модели.

4.4.Расчет значений уровня ряда по аддитивной модели.

4.5.Вычисление отклонений расчетных значений от фак-

тических.

4.6.Оформление результатов расчетов в виде табл. 4.2.2.5.

Та б л и ц а 4.2.2.5

Год

Сезон

yi

Si

yi Si

f

f S

yi f S

 

зима

3835,98

-7288,16

11124,14

10375,68

3087,53

748,46

2000

весна

1831,44

-10848,72

12680,16

10930,09

81,37

1750,07

лето

11588,10

1338,81

10249,29

11484,50

12823,31

-1235,21

 

 

осень

30945,35

16798,07

14147,28

12038,91

28836,98

2108,38

 

зима

4446,88

-7288,16

11735,04

12593,31

5305,15

-858,27

2001

весна

2066,72

-10848,72

12915,44

13147,72

2299,00

-232,28

лето

12937,67

1338,81

11598,86

13702,13

15040,94

-2103,27

 

 

осень

34891,66

16798,07

18093,59

14256,53

31054,60

3837,06

 

зима

5007,45

-7288,16

12295,60

14810,94

7522,78

-2515,34

2002

весна

2096,25

-10848,72

12944,97

15365,35

4516,63

-2420,38

лето

14744,07

1338,81

13405,26

15919,76

17258,57

-2514,50

 

 

осень

39328,02

16798,07

22529,95

16474,16

33272,23

6055,79

68

 

зима

5436,87

-7288,16

12725,03

17028,57

9740,41

-4303,54

2003

весна

2358,03

-10848,72

13206,75

17582,98

6734,26

-4376,23

лето

16644,32

1338,81

15305,51

18137,39

19476,19

-2831,87

 

 

осень

44380,99

16798,07

27582,92

18691,79

35489,86

8891,13

5.Построение с помощью «Мастера диаграмм» графика (рис. 4.2.2.1), отражающего динамику объема продаж груш на рынках г. Воронежа (фактических, рассчитанных по трендовой и аддитивной моделям).

6.Расчет прогнозных значений для каждого сезонного периода

2004г.

yзима 9821,28 554,41 17 7288,16 11958,04; yвесна 9821,28 554,41 18 10848,72 8951,88; yлето 9821,28 554,41 19 1338,81 21693,82; yосень 9821,28 554,41 20 16798,07 37707,49.

Рис. 4.2.2.1. Динамика фактического и расчетных объемов продаж груш, т.

Задание 4.2.2.2. Известно, что стоимость репетиторских услуг зависит от спроса на такие услуги, который распределен по периодам подготовки к вступительным экзаменам. Условно можно выделить четыре периода: 1) август –октябрь (низкая стоимость); 2) ноябрь – декабрь (средняя стоимость); 3) январь – март (стоимость выше средней); апрель – июль (высокая стоимость). Усредненные значения стоимости репетиторских услуг в г. Воронеже за четыре года с разбивкой по периодам приведены в табл. 4.2.2.6. Абитуриенты решили построить модель сезонных

69

колебаний для расчета ожидаемой стоимости репетиторских услуг в 2004 г.

 

 

 

 

 

Таблица 4.2.2.6

Год

Период

Цена занятия,

Год

Период

Цена занятия,

руб.

руб.

 

 

 

 

 

1

65

 

1

145

2000

2

95

2002

2

190

3

140

3

250

 

 

 

4

190

 

4

310

 

1

110

 

1

175

2001

2

155

2003

2

240

3

185

3

310

 

 

 

4

250

 

4

375

Решение с помощью Excel

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

2.Построение графика временного ряда, характеризующего стоимость репетиторских услуг (см. рис. 4.2.2.2).

.

400

 

 

 

 

, руб

300

 

 

 

 

занятия

200

 

 

 

цена

100

 

 

 

 

Цена

0

 

 

 

 

0

5

10

15

20

 

 

 

 

Время, период

 

 

 

Р и с. 4.2.2.2. Динамика стоимости репетиторских услуг

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

3.Расчет оценок сезонной компоненты для мультипликативной модели.

70