АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомДругоеЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция

Пример обнаружения гетероскедастичности в Gretl

Читайте также:
  1. Exercises for Lesson 3. Requests and offers / Просьбы и предложения. Способы выражения, лексика, примеры.
  2. Exercises for Lesson 3. Requests and offers / Просьбы и предложения. Способы выражения, лексика, примеры.
  3. Exercises for Lesson 3. Requests and offers / Просьбы и предложения. Способы выражения, лексика, примеры.
  4. IХ. Примерный перечень вопросов к итоговой аттестации
  5. Автоматическая Система Обнаружения и Тушения Пожаров (АСОТП) «Игла-М.5К-Т»
  6. Алгоритм оценки и проверки адекватности нелинейной по параметрам модели (на примере функции Кобба-Дугласа).
  7. Алгоритм теста Глейзера на наличие или отсутствие гетероскедастичности случайных возмущений.
  8. Алгоритм управления запасами. Пример алгоритма с критическим уровнем.
  9. Анализ нестабильности условий деятельности фирм на примере «Apple»
  10. АСПЕКТЫ ИСТОРИЧЕСКОЙ КОМПЛЕМЕНТАРНОСТИ НА ПРИМЕРЕ КАВКАЗСКОЙ ВОЙНЫ (1817-1864)
  11. б/ Примерное рассуждение
  12. В 2. Дефекты по месту расположения: методы их обнаружения и предупреждения.

 

Шаг 1. Выберем команду меню File\Open Data\Sample File …и на закладке Verbeek двойным щелчком мыши откроем встроенный разработчиками набор данных bwages.gdt ( рисунок 3).

 

Рисунок 3 – Открытие встроенного набора данных bwages.gdt

Шаг 2. Просмотрим текстовую информацию о наборе данных, выбрав команду Data\Print Description. В открывшемся окне (рисунок 4) появится информация о наборе данных bwages.gdt в целом и о каждой переменной.

Файл содержит 1472 наблюдения, относящихся к 1994, группе бельгийских семей Европейского Сообщества. Тип данных undated (срез данных для фиксированного момента времени – cross-sectional). Переменные:

wage – заработная плата в час, до выплаты налогов (Евро),

educ – уровень образования от 1 (низкий) до 5 (высокий),

exper – опыт работы (лет)

male – фиктивная переменная, принимает значение «1», если мужчина и «0», если женщина. Также представлены логарифмы переменных wage, educ, exper.

Требуется оценить влияние уровня образования (educ) в бельгийских семьях на величину заработной платы (wage), используя фактические данные 1994 года по 1472 семьям.

Рисунок 4 – Текстовое описание набора данных bwages.gdt

 

Шаг 3. Проведём предварительный графический анализ данных, построив диаграммы рассеяния пар переменных educ -wage, exper-wage, male –wage.

Выберем команду меню View\Multiple graphs Vars\ X-Y scatters… (рисунок 5). В открывшемся окне выберем переменную wage (по оси ОУ), а переменные educ, exper, male (по оси ОХ) нажатием кнопок «Choose» и «Add» соответственно.

Как показывает рисунок 6, дисперсия wage (которую будем рассматривать как зависимую переменную Y) увеличивается при росте каждого из факторов educ, exper, male(которые будем рассматривать как независимые переменные, Хi), что свидетельствует о гетероскедастичности.

 

 

 

Рисунок 5 – Построение диаграмм рассеяния educ-wage, exper-wage, male-wage

 

Рисунок 6 – Диаграммы рассеяния пар переменных educ -wage, exper-wage, male –wage

Шаг 4. Проведём оценивание модели обычным методом наименьших квадратов 1МНК (OLS).Для этого выберем команду меню Model\Ordinary Least Squares… и в открывшемся окне спецификации модели выберем зависимую переменную wage при помощи кнопки «Choose» и независимую переменную educ при помощи кнопки «Add». После нажатия кнопки ОК появится окно с результатами моделирования (рисунок 7).

 

 

Рисунок 7 – Построение модели методом 1МНК (OLS)

 

По результатам 1МНК в полученной модели параметры являются значимыми при уровнях значимости 5% и 1%, поскольку P-VALUE=0,001%<1%< 5% (принимаем альтернативную гипотезу Н1: ).

Хотя по данным 1МНК можно установить существенность влияния переменной educ (образование) на переменную wage (з.пл.), необходимо сделать поправку на гетероскедастичность в связи с ошибочностью расчёта величин COEFFICIENT, STDERROR и T-STAT и ненадёжностью результатов данного оценивания -возможностью принятия неправильной гипотезы.

Сохраним величины квадратов остатковв отдельную переменную usq1набора данных, обратившись к команде Save\Squared Residuals меню окна результатов моделирования (рисунок 7) и нажав кнопку ОК.

Шаг 5. Подтвердим наличие гетероскедастичности, используя формальный тест Уайта. Обратимся к команде Tests\Heteroskedasticity меню окна результатов моделирования (рисунок 8).

 

Рисунок 8 – Проведение теста Уайта на гетероскедастичность

 

В результате теста получим вспомогательную модель регрессии остатков относительно переменной educ и её квадрата (рисунок 9).

 

Рисунок 9 – Окно результатов теста Уайта на гетероскедастичность

 

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

В окне результатов теста (рисунок 9) значение p-value для статистики теста n*R2 (TR^2) составило 0,000%, что меньше уровней значимости 1% и 5% и свидетельствует о наличии гетероскедастичности (адекватности вспомогательной модели, неравенстве нулю всех её параметров).

Т.о. отвергаем гипотезу Ho о гомоскедастичности (равенстве нулю всех её параметров), поскольку расчётное значение статистики n*R2 =TR^2=56,429 больше критического ().

Найдём критическое для уровня значимости 1%, обратившись к команде основного меню Tools\Statistical Tables и введя в открывшемся окне на закладке Chi-Square число степеней свободы (2, равное числу регрессоров вспомогательной модели) и уровень значимости 0.001 (рисунок 10). Получим =13,8155 (Critical value).

 

Рисунок 10 – Нахождение критического значения для df=2, p=1%

 

Построим график данного распределения, выбрав команду основного меню Tools\Distribution Graphs, закладку Chi-square и df=2 в открывшемся окне (рисунок 11).

Рисунок 11 - График распределения, df=2

 


1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 |

Поиск по сайту:



Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Студалл.Орг (0.004 сек.)