45

Лабораторна робота № 1

Тема. Застосування електронних таблиць EXCEL та пакетів прикладних програм для розв'язування економетричних задач

Мета роботи: ознайомитися з порядком застосування електоронних таблиць та пакетів прикладних програм у сатистичних та економетричних розрахунках.

Завдання

Ознайомитися з прийомами використання електронних таблиць EXCEL для економетричних розрахунків.

Ознайомитися з функціональними можливостями професійних пакетів прикладних програм статистичної обробки даних STATGRAFICS, SPSS.

Хід роботи

1) Для ознайомлення з можливостями застосування електронних таблиць в економетричних розрахунках скласти в оболонці EXCEL розрахункову табл. 1.1. Занести вихідні дані - ряди даних для змінних X, Y. Розрахувати значення граф 4 - 7, а також значення параметрів A, B за формулами:

,

.

Таблиця 1.1

Макет розрахункової таблиці для виконання завдання 1

№ спотереження
1
2
…
n
Сума
Середнє значення

Для розрахунків використати функції СУММ, СРЗНАЧ, СТЕПЕНЬ, КОРЕНЬ, СУММПРОИЗВ, ЧСТРОК.

2) Ознайомитися з можливостями EXCEL при виконанні операцій з матрицями. Використовуючи функції ТРАНСП, МУМНОЖ, МОБОР, виконати дії з матрицями (завдання 1.2).

Вихідні дані для розрахунків:

матриця D = (12 х 4)

y - вектор розмірністю (12 х 1)

Для виконання завдання необхідно згадати елементи матричного обчислювання.

Елементи лінійної алгебри

1. Матриці

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

Основні визначення

Матриці - це прямокутні таблиці елементів, розташованих по рядках та стовпцях:

.

Матриця називається прямокутною матрицею порядку m на n або (m x n) (m - число рядків, n - число стовпців).

Елемент, який знаходиться в i-му рядку та в j-му стовпці, позначається через (перший індекс - номер рядка, другий - стовпця).

Матриця, в якій число рядків дорівнює числу стовпців, називається квадратною. Порядок квадратної матриці визначається одним числом - кількістю рядків (стовпців).

Матриця, яка складається з одного рядка елементів, називається вектором-рядком.

Матриця, яка складається з одного стовпця елементів, називається вектором-стовпцем.

Квадратна матриця називається діагональною, якщо елементи, які не належать головній діагоналі, дорівнюють нулю.

Діагональна матриця, в якої кожен елемент головної діагоналі дорівнює одиниці, називається одиничною і позначається буквою .

Одинична матриця має вигляд:

.

Матриця, в якої всі елементи дорівнюють нулю, називається нульовою.

Квадратна матриця порядку n називається симетричною, якщо виконується умова для всіх елементів цієї матриці.

Рівність двох матриць. Матриця дорівнює матриці , якщо вони однакових розмірів, наприклад (m x n), і мають однакові відповідні елементи:

.

Дії над матрицями

Додавання матриць

Додавання матриць вводиться тільки для матриць одного порядку. Сумою двох матриць і порядку (m x n) називається матриця , яка має такий самий порядок (m x n), причому кожен елемент матриці дорівнює сумі відповідних елементів матриць і :

.

Множення числа на матрицю

Добутком числа на матрицю порядку (m x n) називається матриця порядку (m x n), кожний елемент якої дорівнює добутку числа на відповідний елемент матриці :

.

Для додавання і множення матриць на число справедливі такі операції:

а)

- комутативний закон додавання матриць;

б)

- асоціативний закон додавання матриць;

в)

- асоціативний закон множення чисел на матрицю;

г)

- дистрибутивний закон множення числа на суму матриць;

ґ)

- дистрибутивний закон множення суми чисел на матрицю.

Добуток матриць

Добуток двох матриць вводиться лише для узгоджених матриць. Дві матриці і називаються узгодженими, якщо кількість стовпців першої матриці дорівнює кількості рядків другої матриці .

Якщо матриці порядку (m x n) і порядку (n x p) узгоджені, то добутком цих матриць називається матриця порядку (m x p), для якої елемент дорівнює добутку кожного елемента і-го рядка матриці на j-й стовпець матриці .

Взагалі операція множення матриць не комутативна:

.

Квадратну матрицю можна помножити саму на себе, тобто піднести до квадрата.

Для дій над матрицями справедливі такі властивості:

а)

- асоціативний закон множення матриць;

б)

- дистрибутивний закон множення матриці на суму матриць;

в)

- комутативний закон множення квадратної матриці на одиничну матрицю такого ж порядку.

Транспонування матриць

Матриця ' називається транспонованою відносно матриці , якщо кожен стовпець матриці ' є відповідним рядком матриці , тобто перший стовпець матриці 'є першим рядком матриці , відповідно другий стовпець матриці ' є другим рядком матриці і т.д.

Для елементів транспонованих матриць виконується умова

.

Якщо квадратна матриця симетрична, то виконується умова .

Властивості транспонованих матриць:

Інвертування матриць

Розглянемо невироджену матрицю n-го порядку:

.

Квадратна матриця називається невиродженою, якщо її визначник не дорівнює нулю, тобто , і виродженою, якщо її визначник дорівнює нулю, тобто .

Квадратна матриця називається оберненою до квадратної матриці того ж порядку, якщо їх добуток дорівнює одиничній матриці:

Визначення рангу матриці

Якщо у будь-якій матриці виділити r довільних столбців та r довільних рядків, то з елементів матриці, які містяться на перетині цих рядків і стовпців, можна скласти визначник r-го порядку. Його називають мінором r-го порядку.

Рангом матриці називають число, яке дорівнює найвищому порядку її мінора, відмінного від нуля (rang [A]).

Диференціальне обчислювання в матричній формі

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

1.Похідна від скалярного добутку векторів () по одному з них дорівнює другому:

.

2.Розглянемо добуток , де А - квадратна симетрична матриця порядку n, x - вектор розмірністю n.

або

.

.

Друга частинна похідна по вектору х :

.

Для побудови та аналізу економетричних моделей, а також для прогнозування економічних процесів застосовується ряд професійних пакетів прикладних програм. Такими є пакет STATGRAFICS, SPSS. В рамках лабораторної роботи необхідно поверхньо ознайомитися з призначенням цих пакетів, їх функціональними можливостями та особливостями, а також послідовністю операцій, які виконуються з їх застосуванням.

Завдання для самостійної роботи студентів

Завдання 1.1

Згадати правила виконання операцій з матрицями (додавання, множення, транспонування, інвертування, диференціювання).

Завдання 1.2

Виконати дії над матрицями:

,

,

,

,

 (E - одинична матриця).

Вихідні дані для розрахунків:

, abc - три останні цифри шифру студента,

.

Лабораторна робота № 2

Тема. Парна лінійна регресія

Мета роботи: навчитися будувати парну лінійну регресійну модель економічних процесів.

Завдання

1. На основі спостережених даних показника Y і фактора X знайти оцінки:

коефіцієнтів кореляції і детермінації;

параметрів лінії регресії .

2. Побудувати ANOVA-таблицю для парної регресії.

3. Використовуючи критерій Фішера, з надійністю P=0.95 оцінити адекватність прийнятої моделі статистичним даним.

4. Розрахувати інші показники якості моделі.

5. Використовуючи t-статистику, з надійністю Р=0.95 оцінити значущість коефіцієнта кореляції.

6. Використовуючи t-статитстику, з надійністю Р=0.95 оцінити значущість параметрів моделі та визначити інтервали довіри для параметрів.

7. Якщо модель адекватна статистичним даним, то знайти:

з надійністю Р=0.95 надійні зони базисних даних;

точковий прогноз показника;

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

8. На основі одержаної економетричної моделі зробити висновки.

Хід роботи

1) Коєфіцієнт кореляції є мірою щільності зв'язку між змінними.

Коєфіцієнт кореляції між двома рядами спостережуваних змінних X та Y розраховується за формулою:



Коефіцієнт детермінації дорівнює квадрату коефіцієнта кореляції.

Вводиться гіпотеза, що між фактором Х та показником Y існує лінійна стохастична залежність

.

Оцінки параметрів та парної регресіїобчислюються методом 1МНК за формулами:

,

(або

)

,

де n - кількість спостережень.

Для роботи використовується пакет EXCEL. Складається розрахункова таблиця за макетом (табл.2.1) і розраховуються оцінки параметрів:

Таблиця 2.1

Розрахункова таблиця для оцінки параметрів парної лінійної моделі (за формулами (2.1), (2.3))

№ спостереження	X	Y	XY	X2
1	2	3	4	5
1
2
…
n
Сума	x	х
Середнє значення			х	х
Прогнозне значення

Результат розрахунків - вектор параметрів .

2. Для проведення дисперсійного аналізу складається ANOVA-таблиця (табл. 2.2):

Таблиця 2.2

ANOVA-таблиця

Джерело варіації	Кількість ступенів вільності	Сума квадратів	Середні квадрати
Зумовлене регресією (модель)	К-1
Не пояснюване за допомогою регресії (помилка)	n-K
Загальне	n-1		-

У разі парної регресії К=2 - кількість оцінюваних параметрів.

Для розрахунку ANOVA-таблиці розрахункова табл. 2.1 додається такими графами :

Продовження табл. 2.1

№ спостереження		()2		()2	()2
1	6	7	8	9	10
1
2
…
n
Сума			0
Середнє значення	х	Х	0	х	х
Прогнозне значення

3. Перевірка моделі на адекватність за допомогою критерія Фішера здійснюється за 6-ти-кроковою схемою.

КРОК 1. Формулюються нульова та альтернативна гіпотези:

- незалежна змінна Х не впливає на значення залежної Y.

- значення Х впливає на значення Y.

КРОК 2. Задається рівень значущості : .

КРОК 3. Обчислюється F-відношення:

.

КРОК 4. Знаходиться критичне значення F-розподілу Фішера при заданому рівні значущості та з (К-1), (n-K) ступенями вільності (функція FРАСПОБР в EXCEL) - .

КРОК 5. Порівнюється розрахункове та критичне значення функції F-розподілу.

КРОК 6. Робиться висновок. Якщо , тоді гіпотеза відхиляється, якщо , то приймається.

4. Розраховуються інші показники адекватності моделі:

Середня помилка прогнозу ME:

;

Дисперсія помилок VAR:

та стандартне відхилення:

;

Середній квадрат помилки MSE (з ANOVA-таблиці):

або сума квадратів помилок SSE:

.

Абсолютна середня процентна помилка MAPE:

 ()

Якщо MAPE<10% - існує висока точність прогнозу;

10%< MAPE<20% - добра точність;

20%< MAPE<50% - задовільна точність;

MAPE>50% - незадовільна точність.

Середня процентна помилка MPE:

(MPE<|5%|)

Середня абсолютна помилка MAE:

.

5. Оцінка значущості коефіцієнта кореляції здійснюється за допомогою t-теста (6 кроків).

КРОК 1. Формулюються нульова та альтернативна гіпотези:

- в генеральній сукупності немає зв'язку між X та Y

- коефіцієнт кореляції статистично значущий

КРОК 2. Обирається рівень значущості: .

КРОК 3. Знаходиться розрахункове значення t-статистики:

,

де R - вибірковий коефіцієнт кореляції.

КРОК 4. За таблицями t-розподілу Ст'юдента знаходиться критичне значення функції розподілу (функція СТЬЮРАСПОБР в EXCEL).

КРОК 5. Розрахункове значення t-статистики порівнюється з табличним. Знаходиться критична зона (рис. 2.1).

КРОК 6. Якщо розрахункове значення t-статистики потрапляє в критичну зону, то відхиляється, у ішшому випадку - приймається.

Рис.2.1. Графічне зображення критичної зони для розрахункового значення t-статистики.

6. Етапи тестування за критерієм Ст'юдента на значимість параметрів моделі та .

КРОК 1. Формулюються нульова та альтернативна гіпотези:

- оцінка параметру у генеральній сукупності статистично не значимий,

- оцінка параметру статистично значимий

КРОК 2. Обирається рівень значущості .

КРОК 3. Будується t-статистика для кожного параметру:

,

де - 1МНК оцінка дисперсії параметру ,

;

.

КРОК 4. За таблицями t-розподілу Ст'юдента знаходиться критичне значення функції розподілу (функція СТЬЮРАСПОБР в EXCEL).

КРОК 5. Розрахункове значення t-статистики порівнюється з табличним. Знаходиться критична зона.

КРОК 6. Якщо значення не потрапляє в критичну зону, то можна стверджувати з ймовірністю 95%, що оцінка є статистично незначимою - приймається гіпотеза . Інакше - гіпотеза відхиляється.

Для того, щоб визначити, як параметри та пов'язані з дійсними параметрами та, будуються шнтервали довіри для параметрів моделі за формулою:

.

7. Прогнозування за моделлю простої лінійної регресії.

Точковий прогноз дає значення залежної змінної для відповідного значення , виходячи з побудованої моделі:

.

Надійні зони для базисних даних та прогнозні інтервали знаходяться за формулами:

а) для окремоно значення y:

;

б) для математичного сподівання y:

.

Для розрахунків доповнити розрахункову табл. 2.1 графами:

Продовження табл.2.1

№ спостереження	(розрах. за (2.19))	(розрах.за (2.19))
1	11	12
1
2
…
n
Сума	х	х
Середнє значення	х	х
Прогнозне значення

8. Зробити висновки щодо:

економічної інтерпретації параметрів моделі;

інтерпретації коефіцієнта кореляції;

адекватності побудованої моделі;

прогнозу показника.

Оформити звіт про виконання лабораторної роботи.

Завдання для самостійної роботи студентів

Завдання 2.1

Виконати завдання лабораторної роботи № 2 на основі даних спостереження (табл 2.3):

Таблиця 2.3

Вихідні дані для побудови простої лінійної регресійної моделі

№ спостереження	Незалежна змінна Х	Залежна зміннаY
1	0,12	625
2	0,15	624
3	0,19	362
4	0,12	580
5	0,25	425
6	0,60	аbc
Прогнозне значення	0.90

(abc - три останні цифри шифру студента)

Завдання 1.2

Припустимо, що Ви збираєте дані про річний продаж фірмою продукції (y) і суми, які витрачено на наукові дослідження (x). Ви маєте таку статистику:

cov(x,y)=300;

var(y)=125;

var(x)=880.

Середній річний продаж ()=1200.

Середня сума витрат на наукові дослідження ()=895.

Підрахуйте коефіцієнт кореляції між продажем і сумою, використаною на наукові дослідження. Визначте коефіцієнт детермінації. Знайдіть параметри регресії та .

Завдання 1.3

Проведено оцінку регресії та розраховані SSE та SSR:

SSE=53.27

SSR=202.91.

Розрахуйте SST, R², r.

Завдання 1.4

Вивчаючи зміну попиту на товар залежно від його ціни, отримано такі результати:

;

;

;

;

;

n=4.

Фірма встановлює на товар ціну: 1,75 грн. Спрогнозуйте попит і побудуйте 95%-й інтервал довіри для математичного сподівання прогнозу.

Завдання 1.5

Ви оцінюєте таку регресію:

;

;

n=28;

.

Перевірте значимість нахилу при 95%-ному рівні довіри.

Побудуйте 90%-ний інтервал довіри для нахилу.

Завдання 1.6

На яку додаткову оплату може очікувати особа, яка навчалась додатково 1 рік, якщо співвідношення між заробітною платою (в грн.) - y і освітою (в роках) - x має вигляд:

.

Завдання 1.7

Припустимо, що Ви підрахували кореляцію між двома випадковими змінними, яка дорівнює 0.62. Якщо для оцінки коефіцієнта кореляції було використано 25 спостережень, використайте 5%-ний рівень значимочті, щоб перевірити значимість коефіцієнта кореляції.

Завдання 1.8

Припустимо,що Ви оцінюєте залежність доходу відповідно до кількості років навчання, використовуючи 30 спостережень. Середньоквадратичні відхилення параметрів подано в дужках.

(4,8) (127)

а) перевірте значимість нахилу при 5%-ному рівні значимості;

б) побудуйте 95%-ний інтервал довіри для нахилу.

Завдання 1.9

Припустимо, що в регресії із завдання 1.8 SSE=75, SSR=81. Використайте F-тест для перевірки адекватності регресії.

Лабораторна робота № 3

Тема. Парна нелінійна регресія

Мета роботи: навчитися будувати парну нелінійну регресійну модель економічних процесів.

Завдання

Виконати завдання для самостійної роботи 3.1.

На основі статистичних даних показника Y і фактора X (вихідні дані з л.р.№1) знайти оцінки параметрів лінії регресії, якщо припустити, що стохастична залежність між фактором X і показникомY має визляд однієї з вищерозраховуваних функцій (завдання 1)

Використовуючи критерій Фішера з надійністю P=0.95, оцінити адекватність прийнятої моделі статистичним даним.

Побудувати ANOVA-таблицю для нелінійної моделі.

Якщо з заданою надійністю прийнята модель адекватна експерементальним даним, то знайти:

а) з надійністю Р=0.95 довірчу зону базисних даних;

б) точкову оцінку прогнозу;

в) з надійністю Р=0.95 інтервальну оцінку прогнозу

Зробити висновки щодо ступеню апроксимації вихідних даних лініями простої лінійної та нелінійної регресії.

Побудувати графіки:

а) фактичних даних;

б) лінії прямої регресії та її довірчу зону (л.р.№2);

в) нелінійної функції та її довірчу зону.

Хід роботи

Зведення кривих зростання до лінійної функції дає змогу оцінити параметри методом 1МНК та використовувати подальший аналіз моделі.

Приклад 3.1

Зведення нелінійної парної регресії до лінійної виконується заміною , . У результаті маємо лінійну регресію: .

Вводиться гіпотеза, що між фактором X та показником Y (вихідні ряди даних беруться з л.р.№1) існує нелінійна залежність (завдання 3.1, за вибором). Заміною приводиться нелінійна парна регресія до парної лінійної вигляду:

.

Параметри оцінюються за формулами, аналогічними (2.1) - (2.3) із застосуванням пакету EXCEL.

Перевірка моделі на адекватність за критерієм Фішера проводиться за 6-ти кроковою схемлю (п.3, л.р.№2).

Розраховується ANOVA-таблиця (табл.2.2).

За формулами (2.18) - (2.20) у разі адекватності моделі розраховуються з надійністю Р=0.95 довірчі зони базисних даних, точкова оцінка прогнозу, з надійністю Р=0.95 інтервальна оцінка прогнозу індивідуального заначення показника Y та його математичного сподівання.

На основі результатів дисперсійного аналізу зробити висновки щодо порівняної якості двох побудованих моделей - парної лінійної та нелінійної (критерій - min SSE).

Проводиться графічний аналіз даних:

а) будується діаграма розсіювання (на координатну площину наносяться сточки спостережуваних даних);

б) зображуються графіки оцінених лінійної та нелінійної функцій регресії та їх довірчі зони

в) робляться висновки.

Завдання для самостійної роботи студентів

Завдання 3.1

Шляхом необхідних перетворень та заміни змінних звести наведені нелінійні функції до лінійного вигляду.

; ; ; ; ; ; ; ;	; ; ; ; ; ; ; .

Завдання 3.2

Чи можна параметри модифікованої експоненти розрахувати за методом найменших квадратів? Поясніть відповідь.

Завдання 3.3

Наведені такі дані (табл.3.1):

Таблиця 3.1. Вихідні дані для побудови моделі

86	3
79	7
76	12
69	17
65	25
ab	bc

(abc - три останні цифри шифру студента)

Побудуйте за наведеними даними модель вигляду

,

оцініть її параметри.

Лабораторна робота № 4

Тема. Багатофакторна модель лінійної регресії

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

Завдання

Підприємство має велику кількість філіалів, і керівництво цього підприємства хотіло б знати, як Y (річний товарообіг одного філіалу, млн.грош.од.) функціонально залежить від X2 - торговельної площі, тис. м², та X3 - середньоденної інтенсивності потоку покупців, тис.чол/день. Конкретно необхідно визначити, яке значення має кожний коефіцієнт такого регресійного рівняння:

Для дванадцяти філіалів за певний рік маємо фіксовані значення показників Y, X2 та X3 (табл. 4.1).

Таблиця 4.1

Просторові дані за філіалами підприємства

№ філіалу	Значення Y	Значення X2	Значення X3
1	2.93	0.31	10.24
2	5.27	0.98	7.51
3	6.85	1.21	10.81
4	7.01	1.29	9.89
5	7.02	1.12	13.72
6	8.35	1.49	13.92
7	4.33	0.78	8.54
8	5.77	0.94	12.36
9	7.68	1.29	12.27
10	3.16	0.48	11.01
11	1.52	0.24	8.25
12	3.15	0.55	9.31

Оцінити параметри моделі за методом 1МНК (у матричній формі). Інтерпретувати отримані оцінки.

Оцінити стандартизовані регресійні коефіцієнти ("бета-коефіцієнти"). Інтерпретувати оцінені стандартизовані коефіцієнти регресії.

Скласти до числового прикладу вектори , ,, , а також матриці X та D.

Розрахувати значення величин , , .

Оцінити еластичність товарообігу відносно торговельної площі та відносно середньоденної частоти потоку покупців, обчисливши коефіцієнти еластичності.

Перевірити значимість окремих коефіцієнтів регресії (провести t-тестування), визначити їх інтервали довіри.

7. Розрахувати та інтерпретувати коефіцієнт детермінацї, частинний коефіцієнт детермінації та зкоректований коефіцієнт детермінації.

8. Перевірити модель на адекватність за допомогою F-критерію Фішера.

9.У разі адекватності моделі обчислити та інтерпретувати для регресії:

точковий прогноз товарообігу t+1-го філіалу;

99%-ний прогнозний інтервал математичного сподівання товарообігу цього філіалу;

99%-ний прогнозний інтервал безпосередньо самого товарообігу y_t+1 цього філіалу, якщо задані такі значення регресорів:

Хід роботи

Для знаходження вектора оцінок параметрів багатофакторної лінійної моделі застосовується метод 1МНК у матричній формі:

Параметри лінійної регресії інтерпретуються так: зміна величини к-го регресора на одиницю свого виміру за інших рівних умов призведе до зміни оціненої величини на число одиниць свого виміру, яке дорівнює значенню .

Стандартизовані коефіцієнти регресії обчислюються за формулою:

, (k=2,…,k),де

1МНК-оцінка регресійного коефіцієнта ;

- емпіричне стандартне (середньоквадратичне) відхилення k-го регресора x_k

- емпіричне стандартне (середньоквадратичне) відхилення регресанда y.

Емпіричний стандартизований регресійний коефіцієнт вказує на те, який великий за інших рівних умов типовий ефект впливу k-го регресора у порівнянні з типовим ефектом зміни регресанда.

Складаються такі вектори і матриці:

- вектор спостережуваних значень показника Y;

- вектор оцінених значень регресанда;

- вектор дійсних, але невідомих параметрів регресії;

- вектор 1МНК-оцінок параметрів моделі;

 - 1МНК-оцінка вектору помилок;

- матриця регресорів;

- матриця даних.

Розраховуються величини: - суму помилок регресії (має дорівнювати 0), - суму квадратев помилок, - дисперсію помилок.

Коефіцієнти еластичності розраховуються за формулою:

,

де - значення регресанда і к-го регресора, що визначають точку регресійної функції, для якої обчислюється коефіцієнт еластичності. Можна використовувати та - середні значення.

t-тест для перевірки гіпотези про числові значення окремих коефіцієнтів регресії проводиться за 6-тикроковою схемою (схема наведена у л.р.№ 2, п. 6).

Інтервал довіри для регресійного коефіцієнта при рівні довіри (1-) є інтервалом з випадково залежними межами.

Довірчий інтервал для дійсного значення регрессійного коефіцієнта :

Коефіцієнт детермінації дрівнює квадрату емпіричного множинного коефіцієнта кореляції між двома рядами спостережень: емпіричних значень регресанда() та його розрахунковим значенням ().

Три формули для розрахунку коефіцієнта детермінації:

;

;

.

Регресійне рівняння оцінене тим краще, чим більше за інших рівних умов .

Частинний коефіцієнт детермінації називається граничним вкладом k-го регресора в і показує, на яку величину зменшується коефіцієнт детермінації, якщо k-й регресор (і тільки він) буде виключений із групи K регресорів.

,

де - коефіцієнт детермінації, одержаний при включенні усіх К регресорів;

- квадрат обчисленого значення t-статистики для k-го регресійного коефіцієнта-

;

Т - довжина ряду спостережень;

К - кількість регресорів;

Т-К - кількість ступенів вільності.

З двох варіантів регресійних рівнянь, які відрізняються на величину зкоректованого коефіцієнта детермінації, але мають однаково гарні інші критерії якості, обирають варіант з більшим значенням зкоректованого коефіцієнта детермінації.

Зкоректований коефіцієнт детермінації за Тейлом:

.

Зкоректований коефіцієнт детермінації за Амемієй:

.

F-тест може бути проведений за схемою, яка складається з 6-ти кроків.

КРОК 1. Формулюється пара гіпотез:

- жоден регресор не впливає на регресанд

існує хоча б один регресор, який впливає на регресанд: .

КРОК 2. Обирається рівень значимості.

КРОК 3. Визначається табличне значення F-критерію (функція FРАСПОБР в EXCEL).

КРОК 4. Числове значення F-статистики може бути розраховане за спрощеною формулою із застосуванням коефіцієнта детермінації :

.

КРОК 5. Порівнюється розрахована величина F з її табличним значенням та приймається рішення у відповідності з правилом застосування F-тесту:

відхиляється, якщо .

КРОК 6. Інтерпретуються результати тесту.

У разі адекватності моделі вона застосовується для прогнозування економічного показника.

Точковий прогноз регресанда одержують, виходячи з оціненого регресійного рівняння:

.

Прогнозні інтервали для математичного сподівання та індивідуального значення регресанда визначається при рівні довіри (1-) таким чином:

,

де - оцінена стандартна помилка прогнозу:

помилка прогнозу при оцінці математичного сподівання регресанда

(),

де .

помилка прогнозу при оцінці індивідуального значення регресанда

Як висновок наводиться геометрична інтерпретація прогнозних інтервалів. Обидва прогнозних інтервали є найменшими при .

Завдання для самостійної роботи студентів

Завдання 4.1

Є такі дані (табл. 4.2)

Таблиця 4.2

Вихідні дані для побудови класичної регресійної моделі

Відстань від фінішу, Y	Вага спортсмена, X2	Час, X3
925	210	4,8
850	185	4,7
1622	225	4,7
1121	215	4,6
658	180	4,9
Abс	212	4,6

(abc - три останні цифри шифру студента)

1. Оцініть параметри лінійної регресійної моделі . Інтерпретуйте коефіцієнти регресії.

2. Перевірте на значимість оцінені параметри з рівнем значимості 5% за t-тестом Ст'юдента.

3. Перевірте адекватність моделі за F-тестом при 5%-ному рівні значимості.

4. У разі адекватності моделі зробіть 95%-ний прогноз індивідуального значення регресанда та його математичного сподівання.

Завдання 4.2

Ви оцінюєте таку модель: , що базується на 20 спостереженнях, і отримали

,

Підрахуйте коефіцієнт детермінації . Чи “добре” модель пояснює виявлену закономірність?

Завдання 4.3

Скільки ступенів вільності мають чисельник та знаменник F-статистики в регресії, що складається з 50 спостережень та 5 незалежних змінних?

Лабораторна робота № 5. Тема. Виробничі функції

Мета роботи: навчитися аналізувати виробничі процеси за допомогою виробничої функції.

Завдання

Побудуйте виробничу функцію Кобба-Дугласа, використовуючи дані про випуск продукції, витрати основних виробничих фондів, витрати праці за 10 років, наведені в табл. 5.1. Розрахуйте характеристики:

середню продуктивність ресурсів;

граничну ефективність ресурсів;

еластичність випуску за ресурсами;

потребу в ресурсах;

фондоозброєність праці;

граничну норму заміщення витрат праці виробничими фондами;

еластичність заміщення ресурсів.

Таблиця 5.1

Вихідні дані для побудови виробничої функції

№ року	Випуск продукції, грш.од. (Y)	Витрати виробничих фондів, грош.од. (X1)	Витрати праці, люд-год., (X2)
1	6.2	3.9	2.3
2	6.4	4.8	2.9
3	7.2	5.6	3.2
4	8.2	7.3	3.6
5	9.5	8.4	4.2
6	10.3	9.5	4.5
7	11.3	11.4	5.2
8	12.7	12.6	5.4
9	13.9	13.4	6
10	14.5	14.5	6.7

Хід роботи

1. Для розрахунку параметрів функції використовуємо рівняння:

.

Прологарифмуємо вихідне рівняння:

.

Введемо заміну змінних:

,

,

,

.

Одержуємо лінійну форму виробничої функції:

,

параметри якої оцінемо за методом 1МНК.

2. Розрахуємо основні характеристики цієї функції.

Середня продуктивність праці:

.

Із зростанням витрат праці (величины Х2) середня продуктивність праці знижується (показник степеню від'ємний). Це зумовлено тим, что кількість засобів праці залишається незмінною, і тому фондоозброєність праці спадає, а разом із нею і продуктивність. Із збільшенням же основних фондів (величины Х1) середня продуктивність праці зростає.

Середня фондовіддача:

Із збільшенням фондів середня фондовіддача спадає (при незмінних трудовых ресурсах), а зростання залученої до виробництва робочої сили (за фіксованої величини фондів) призводить до зростання фондовіддачі.

Гранична продуктивність праці:

.

Із зростанням витрат праці за незмінних фондів гранична продуктивність праці знижується. Із зменшенням обсягу фондів за незмінних трудових ресурсах (тобто при збільшенні фондоозброєності праці), гранична продуктивність праці зростає.

Гранична фондовіддача:

.

При збільшенні обсягу виробничих фондів за незмінних трудовых ресурсах гранична фондовіддача знижується. При збільшенні обсягу трудових ресурсів при незмінних фондах гранична фондовіддача зростає.

Еластичність випуску продукції за витратами праці:



Таким чином, показник ступеню функції Кобба-Дугласа є коефіцієнтом еластичності випуску продукції за витратами праці і показує, що при збільшенні витрат праці на 1% и незмінних основних фондах випуск продукції зростає на %.

Еластичність випуску продукції за витратами засобів праці:

Показник ступеню функції Кобба-Дугласа є коефіцієнтом еластичності випуску продукції за витратами засобів праці. Тобто, збільшення використованих у виробництві основних фондів на 1% сприяє (за інших рівних умов) зростанню випуску продукції на %.

Потреба в ресурсах

Виробнича функція дозволяє розрахувати потребу в одному з виробничих факторів при заданому обсязі випуску та величині іншого фактора:

;

.

Фондоозброєність праці:

Виробнича функція дозволяє досліджувати питання співвідношення, взаємодії, заміщення ресурсів. На основі отриманих вище співвідношень визначається важливий економічний показник - фондоозброєність праці:.

Гранична норма заміщення праці основними фондами дорівнює:

.

Еластичність заміщення ресурсів

Функція Кобба-Дугласа має еластичність заміщення праці капіталом, яка дорівнює 1.

Завдання для самостійної роботи студентів

Завдання 5.1

Побудуйте виробничу функцію Кобба-Дугласа, використовуючи дані про випуск продукції, витрати основних виробничих фондів, витрати праці за 6 років, наведені в табл. 5.2. Розрахуйте характеристики:

середню продуктивність ресурсів;

граничну ефективність ресурсів;

еластичність випуску за ресурсами;

потребу в ресурсах;

фондоозброєність праці;

граничну норму заміщення витрат праці виробничими фондами;

еластичність заміщення ресурсів.

Таблиця 5.2

Вихідні дані для побудови виробничої функції

№ року	Випуск продукції, грш.од. (Y)	Витрати виробничих фондів, грош.од. (X1)	Витрати праці, люд-год., (X2)
1	6.2	3.9	2.3
2	6.4	4.8	2.9
3	7.2	5.6	3.2
4	8.2	7.3	3.6
5	9.5	8.4	4.2
6	ab	b	c

(abc - три останні цифри шифру студента).

Інтерпретувати одержані результати та зробити економічні висновки.

Лабораторна робота № 6

Тема. Взаємозалежні економетричні моделі

Мета роботи: навчитися оцінювати параметри системи взаємозалежних рівнянь.

Завдання

1. На основі статистичних даних за 10 періодів

ендогенних величин:

Y1 - експорт,

Y2 - імпорт

і екзогенних величин:

Х1 - національний дохід;

Х2 - оборот зовнішньої торгівлі,

використовуючи метод 2МНК, оцінити параметри структурної системи регресій:

2. Для даних значень екзогенних величин знайти точкові оцінки прогнозу ендогенних величин та зробити аналіз взаємного впливу величин.

Вихідні дані для розрахунків:

Таблиця 6.1

Вихідні дані для оцінки параметрів структурної системи регресії

Період часу	Х1	Х2	Y1	Y2
1	606,2	371,0	70,2	50,3
2	617,1	396,8	74,6	55,7
3	623,0	400,8	81,4	59,9
4	626,9	411,3	87,9	64,9
5	635,5	432,7	92,5	70,0
6	636,4	442,9	99,6	45,2
7	644,5	444,8	106,5	80,7
8	647,2	472,4	112,2	85,8
9	654,2	478,9	116,5	89,9
10	659,4	491,0	121,3	95,8
Прогнозні заначення екзогенних величин	667,3	507,6

Хід роботи

1. Для вибору методу оцінки параметрів сумісної системи регресій перевіримо її ідентифікацію. Умова ідентифікації регресії має вигляд

,

де n - число ендогенних величин у системі регресій,

m - число екзогенних величин у системі регресій,

n_i і m_i - відповідно число ендогенних і екзогенних величин в і-й регресії. Для першої регресії ця умова запишеться так:

2+2 - (2+2)=0

Оскільки 0<1, то умова ідентифікації невиконана. Тому для оцінки параметрів структурної системи регресій (6.1) використовуємо метод 2МНК.

ПЕРШИЙ КРОК МНК. Будемо вважати, що визначник матриці не дорівнює нулю, тоді можна перетворити структурну систему регресій у приведену :

,

або в розгорнутій формі:



Оцінки параметрів для прогнозної системи регресій (6.4) знаходимо за методом 1МНК за формулою в матричній формі:

.

ДРУГИЙ КРОК МНК. Величини , , які знаходяться справа, вважаємо передвизначеними і, використовуючи МНК, знаходимо оцінки параметрів структурної системи регресій за такою формою:

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

Зробити висновок щодо прогнозних значень ендогенних величин та взаємного впливу ендогенних змінних (інтерпретувати коефіцієнти ).

Завдання для самостійної роботи студентів

Завдання 6.1

На основі статистичних даних за 6 періодів, використовуючи метод 2МНК, оцінити параметри структурної системи регресій:

Для даних значень екзогенних величин знайти точкові оцінки прогнозу ендогенних величин та зробити аналіз взаємного впливу величин.

Література

Грубер Й. Эконометрия. Т.1. Введение в эконометрию. - К., 1996.

Замков О.О., Толстопятенко А.В.,Черемных Ю.Н. Математические методы в экономике. - М.:МГУ им. М.В. Ломоносова, 1997.

Лук'яненко І.Г., Краснікова Л.І. Економетрика:Підручник. - К.: Тов. “Знання”, 1998.

Лук'яненко І.Г., Краснікова Л.І. Економетрика:Практикум з використанням комп'ютера. - К.: Тов. “Знання”, 1998.

Толбатов Ю.А. Економетрика:Підручник для студентів екон.спец.вищ.навч.закл. - К.:Четверта хвиля,1997.