Тема 2. Классическая и обобщенная модели множественной регрессии.

Экономические явления определяются, как правило, большим числом совокупно действующих факторов. В связи с этим часто возникает задача исследования зависимости одной переменной Y от нескольких объясняющих переменных X₁, X₂, X_n. Эта задача решается с помощью множественного регрессионного анализа.

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

- они должны быть количественно измеримыми (качественным факторам необходимо придать количественную определенность);

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

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

∙ затрудняется интерпретация параметров множественной регрессии как характеристик действия факторов в «чистом виде», поскольку факторы связаны между собой; параметры линейной регрессии теряют экономический смысл;

∙ оценки параметров ненадежны, имеют большие стандартные ошибки и меняются с изменением объема наблюдений.

 

Пусть - матрица-столбец значений зависимости переменной размера п,

 

 

- матрица значений объясняющих переменных;

 

- матрица-столбец (вектор) параметров размера т+1;

- матрица-столбец (вектор) остатков размера п;

Тогда в матричной форме модель множественной линейной регрессии запишется следующим образом:

(1)

При оценке параметров уравнения регрессии (вектора b ) применяется метод наименьших квадратов (МНК). При этом делаются определенные предпосылки.

1. В модели (1) ε – случайный вектор, Х – неслучайная (детерминированная) матрица.

2. Математическое ожидание величины остатков равно нулю: М(ε ) = 0_п.

3. Дисперсия остатков ε _i постоянна для любого i (условие гомоскедастичности), остатки ε _i и ε _j при i ≠ j не коррелированны:

4. ε – нормальное распределенный случайный вектор, т.е. ε ~ N(0_n; σ ²E_n).

5. r (X) = m+1 < n. Столбцы матрицы Х должны быть линейно независимыми (ранг матрицы Х максимальный, а число наблюдений п превосходит ранг матрицы).

Модель (1), в которой зависимая переменная, остатки и объясняющие переменные удовлетворяют предпосылкам 1-5 называется классической нормальной линейной моделью множественной регрессии. Если не выполняется только предпосылка 4, то модель называется классической линейной моделью множественной регрессии (КЛММР).

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

 

 

Решением этой задачи является вектор

Одной из наиболее эффективных оценок адекватности модели является коэффициент детерминации R², определяемый формулой:

 

 

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

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

 

,

где:

п – число наблюдений;

т – число параметров при переменных х.

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

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

Средняя относительная ошибка аппроксимации рассчитывается по формуле:

 

 

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

 

 

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

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

	=		+
Общая сумма квадратов отклонений		Сумма квадратов отклонений, объясненная регрессией		Остаточная сумма квадратов отклонений

 

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

 

Компоненты дисперсии	Сумма квадратов	Число степеней свободы	Оценка дисперсии на одну степень свободы
Общая			-------
Объясненная
Остаточная

 

В данной таблице п – число наблюдений, т – число параметров при переменных х.

Сравнивая полученные оценки объясненной и остаточной дисперсии на одну степень свободы, определяют значение F -критерия Фишера, используемого для оценки значимости уравнения регрессии:

 

С помощью F -критерия проверяется нулевая гипотеза о равенстве дисперсий Н₀: S_R² = S².

Если нулевая гипотеза справедлива, то объясненная и остаточная дисперсии не отличаются друг от друга. Для того, чтобы уравнение регрессии было значимо в целом (гипотеза Н₀ была отвергнута) необходимо, чтобы объясненная дисперсия превышала остаточную в несколько раз. Критическое значение F -критерия определяется по таблице Фишера - Снедекора.

Расчетное значение сравнивается с табличным, и если оно превышает табличное ( F_расч > F_табл ), то гипотеза Н₀ отвергается, и уравнение регрессии признается значимым.

Если F_расч < F_табл, то уравнение регрессии считается статистически незначимым. Нулевая гипотеза Н₀ не может быть отклонена.

Расчетное значение F -критерия связано с коэффициентом детерминации R² следующим соотношением:

где:

т – число параметров при переменных х;

п – число наблюдений.

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

 

 

Проверка гипотезы проводится с помощью t -критерия Стьюдента. Для этого расчетное значение t -критерия:

 

,

где:

- коэффициент регрессии при ;

- средняя квадратическая ошибка коэффициента регрессии .

сравнивается с табличным t_табл при заданном уровне значимости α и числе степеней свободы

( п-2 ).

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

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

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

Уравнение регрессии в стандартизованной форме:

 

 

где: . – стандартизованные переменные.

 

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

 

Коэффициенты «чистой» регрессии связаны со стандартизованными коэффициентами следующим соотношением:

 

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

Средние коэффициенты эластичности вычисляются по формуле:

 

 

Коэффициент эластичности показывает, на сколько процентов (от средней) изменится в среднем Y при увеличении только фактора X_i на 1%.

Рассмотрим пример построения модели множественной регрессии с помощью средств приложения Microsoft Excel.

 

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

Х₁ – объем товарных запасов, тыс. руб.;

Х₂ – фонд оплаты труда, тыс. руб.;

Х₃ – издержки обращения, тыс. руб.;

Х₄ – объем продаж по безналичному расчету, тыс. руб.

Таблица 2

Месяц	Y	Х1	Х2	Х3	Х4
1	2	3	4	5	6
1.	41321, 57	300284, 10	19321, 80	42344, 92	100340, 02
2.	40404, 27	494107, 21	20577, 92	49000, 43	90001, 35
3.	37222, 12	928388, 75	24824, 91	50314, 52	29301, 98
4.	37000, 80	724949, 11	28324, 87	48216, 41	11577, 42
5.	29424, 84	730855, 33	21984, 07	3301, 30	34209, 84
6.	20348, 19	2799881, 13	11000, 02	21284, 21	29300, 00
7.	11847, 11	1824351, 20	4328, 94	28407, 82	19513, 92
8.	14320, 64	1624500, 80	7779, 41	40116, 00	17343, 20
9.	18239, 46	1115300, 93	18344, 11	32204, 98	4391, 00
10.	22901, 52	1200947, 52	20937, 31	30105, 29	14993, 25
11.	27391, 92	1117850, 93	27344, 30	40294, 40	104300, 00
12.	44808, 37	1379590, 02	31939, 52	42239, 79	119804, 33
13.	40629, 28	588365, 77	29428, 60	55584, 35	155515, 15
14.	31324, 80	434281, 91	30375, 82	49888, 17	60763, 19
15.	34847, 92	1428243, 59	33000, 94	59866, 55	8763, 25
16.	33241, 32	1412181, 59	31322, 60	49975, 79	4345, 42
1	2	3	4	5	6
17.	29971, 34	1448274, 10	20971, 82	3669, 92	48382, 15
18.	17114, 90	4074616, 71	11324, 93	26032, 95	10168, 00
19.	8944, 94	1874298, 99	8341, 52	29327, 21	22874, 40
20.	17499, 58	1525436, 47	10481, 14	40510, 01	29603, 05
21.	19244, 80	1212238, 89	18329, 90	37444, 69	16605, 16
22.	34958, 32	1154327, 22	29881, 52	36427, 22	32124, 63
23.	44900, 83	1173125, 03	34928, 60	51485, 62	200485, 00
24.	57300, 25	1435664, 93	41824, 92	49959, 92	88558, 62

 

Задание:

1. Для заданного набора данных постройте линейную модель множественной регрессии.

2. Оцените точность и адекватность построенного уравнения регрессии.

3. Выделите значимые и незначимые факторы в модели.

4. Постройте уравнение регрессии со статистически значимыми факторами. Дайте экономическую интерпретацию параметров модели.

 

Решение:

Для получения отчета по построению модели в среде EXCEL необходимо выполнить следующие действия:

1. В меню Сервис вбираем строку Анализ данных. На экране появится окно

Рис. 1.

2. В появившемся окне выбираем пункт Регрессия. Появляется диалоговое окно, в котором задаем необходимые параметры (рис. 2).

Рис. 2.

3. Диалоговое окно рис. 2 заполняется следующим образом:

Входной интервал Y – диапазон (столбец), содержащий данные со значениями объясняемой переменной;

Входной интервал X – диапазон (столбцы), содержащий данные со значениями, объясняющих переменных.

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

Константа-ноль – флажок, указывающий на наличие или отсутствие свободного члена в уравнении регрессии ( β ₀ );

Выходной интервал – достаточно указать левую верхнюю ячейку будущего диапазона, в котором будет сохранен отчет по построению модели;

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

Если необходимо получить значения и графики остатков ( e_i ), установите соответствующие флажки в диалоговом окне. Нажмите кнопку ОК.

Вид отчета о результатах регрессионного анализа представлен на рис. 3.

Рис. 3.

Рассмотрим таблицу «Регрессионная статистика».

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

R-квадрат – это R². В нашем примере значение R² = 0, 8178 свидетельствует о том, сто изменения зависимой переменной Y (балансовой прибыли) в основном (на 81, 78%) можно объяснить изменениями включенных в модель объясняющих переменных – X₁, X₂, X₃, X₄. такое значение свидетельствует об адекватности модели.

Нормированный R-квадрат – поправленный (скорректированный по числу степеней свободы) коэффициент детерминации.

Стандартная ошибка регрессии - необъясненная дисперсия (мера разброса зависимой переменной вокруг линии регрессии); п – число наблюдений (в нашем примере равно 24), т – число объясняющих переменных (в нашем примере равно 4).

Наблюдения – число наблюдений п.

Рассмотрим таблицу с результатами дисперсионного анализа.

df – degrees of freedom – число степеней свободы связано с числом единиц совокупности п и с числом определяемых по ней констант (т +1).

SS – sum of squares – сумма квадратов (регрессионная (RSS – regression sum of squares), остаточная (ESS – error sum of squares) и общая (TSS – total sum of squares), соответственно).

MS – mean sum – сумма квадратов на одну степень свободы.

F – расчетное значение F–критерия Фишера. Если нет табличного значения, то для проверки значимости уравнения регрессии в целом можно посмотреть Значимость F. На уровне значимости α = 0.05 уравнение регрессии признается значимым в целом, если Значимость F < 0.05, и незначимым, если Значимость F ≥ 0.05.

Для нашего примера имеем следующие значения:

 

	df	SS	MS	F	Значимость F
Регрессия	т = 4	RSS = 2, 82Е+09	7, 04Е+08	21, 32	8, 28Е-07
Остаток	п – т – 1 = 19	ESS = 6, 27Е+08	3, 30Е+07
Итого	п - 1 = 23	TSS = 3, 44Е+09

 

В нашем случае значение F–критерия Фишера составляет 21, 32. Значимость F = 8, 28Е-07, что меньше 0, 05. Таким образом, полученное уравнение в целом значимо.

В последней таблице приведены значения параметров (коэффициентов) модели, их стандартные ошибки и расчетные значения t–критерия Стьюдента для оценки значимости отдельных параметров модели.

 

	Коэффициенты	Стандартная ошибка	t – статистика	Р – Значение	Нижние 95%	Верхние 95%
Y	b0 = 7825, 51	тb0 = 5350, 78	tb0 = 7825, 51/5350, 78 = 1, 4625	0, 1599	-3373, 80≤ b0 ≤ 19024, 83
X1	b1 = -0, 00098	тb1 = 0, 00172	tb1 = -0, 569	0, 5762	-0, 0046≤ b1 ≤ 0, 0026
X2	b2 = 0, 8806	тb2 = 0, 15891	tb2 = 5, 5417	0, 00002	0, 5480≤ b2 ≤ 1, 2132
X3	b3 = 0, 0094	тb3 = 0, 09754	tb3 = 0, 0961	0, 9244	-0, 1948≤ b3 ≤ 0, 2135
X4	b4 = 0, 0617	тb4 = 0, 02647	tb4 = 2, 3312	0, 0309	0, 0063≤ b4 ≤ 0, 1171

 

Анализ таблицы для рассматриваемого примера позволяет сделать вывод о том, что на уровне значимости α = 0.05 значимыми оказываются лишь коэффициенты при факторах X₂ и X₄, так как только для них Р – Значение меньше 0.05. Таким образом, факторы не существенны X₁ и X₃, и их включение в модель нецелесообразно.

Поскольку коэффициент регрессии в эконометрических исследования имеет четкую экономическую интерпретацию, то границы доверительного интервала для коэффициента регрессии не должны содержать противоречивых результатов, как например, – -0, 1948 ≤ b₃ ≤ 0, 2135. такого рода запись указывает, что истинное значение коэффициента регрессии одновременно содержит положительные и отрицательные величины и даже ноль, чего не может быть. Это также подтверждает вывод о статистической незначимости коэффициентов регрессии при факторах X₁ и X₃.

Исключим несущественные факторы X₁ и X₃ и построим уравнение зависимости Y (балансовой прибыли) от объясняющих переменных X₂ и X₄. результаты регрессионного анализа приведены в таблице 3.

Таблица 3.

ВЫВОД ИТОГОВ

Регрессионная статистика
Множественный R	0, 9024465
R - квадрат	0, 8144098
Нормированный R - квадрат	0, 7967345
Стандартная ошибка	5515, 53984
Наблюдения

 

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

	df	SS	MS	F	Значимость F
Регрессия				46, 076253	2, 08847Е-08
Остаток		638844774, 1	30421179, 72
Итого

 

	Коэффициенты	Стандартная ошибка	t – статистика	Р – Значение	Нижние 95%	Верхние 95%
Y	5933, 1025	2844, 611998	2, 085733487	0, 0493883	17, 40698	11848, 798
X2	0, 9162546	0, 132496978	6, 915286693	7, 384Е-07	0, 640712	1, 1917972
X4	0, 0645183	0, 024940789	2, 58686011	0, 0172036	0, 012651	0, 1163856

 

Оценим точность и адекватность модели.

Значение R² = 0, 8144 свидетельствует о том, что вариация зависимой переменной Y (балансовой прибыли) по-прежнему в основном (на 81, 44%) можно объяснить вариацией включенных в модель объясняющих – X₂, и X₄. это свидетельствует об адекватности модели.

Значение поправленного коэффициента детерминации (0, 7967) возросло по сравнению с первой моделью, в которую были включены все объясняющие переменные (0, 7794).

Стандартная ошибка регрессии во втором случае меньше, чем в первом (5515 < 5745).

Расчетное значение F–критерия Фишера составляет 46, 08. Значимость F = 2, 08847Е-08, что меньше 0.05. таким образом, полученное уравнение в целом значимо.

Далее оценим значимость отдельных параметров построения модели. Из таблицы 3 видно, что теперь на уровне значимости α = 0.05 все включенные в модель факторы являются значимыми: Р – Значение < 0.05.

Границы доверительного интервала для коэффициентов регрессии не содержат противоречивых результатов:

- с надежностью 0.95 (с вероятностью 95 %) коэффициент b₁ лежит в интервале

0, 64 ≤ b₁ ≤ 1, 19;

- с надежностью 0.95 (с вероятностью 95 %) коэффициент b₂ лежит в интервале

0, 01 ≤ b₂ ≤ 0, 12.

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

Рассмотрим теперь экономическую интерпретацию параметров модели.

Коэффициент b₁ = 0, 916, означает, что при увеличении только фонда оплаты труда (X₂) на 1 тыс. рублей балансовая прибыль в среднем возрастет на 0, 916 тыс. рублей, а то, что коэффициент b₂ = 0, 065, означает, что увеличение только объема продаж по безналичному расчету (X₄) на 1 тыс. рублей приводит в среднем к увеличению балансовой прибыли на 0, 065 тыс. рублей. Как было отмечено выше, анализ Р–значений показывает, что оба коэффициента значимы.

При эконометрическом моделировании реальных экономических процессов предпосылки КЛММР нередко оказываются нарушенными: дисперсии остатков модели не одинаковы (гетероскедастичность остатков), или наблюдается корреляция между остатками в разные моменты времени (автокоррелированные остатки). Тогда предпосылка 3 запишется следующим образом:

Принимая, что дисперсии объясняющих переменных могут быть произвольными, мы получаем обобщенную линейную модель множественной регрессии (ОЛММР).

В этом случае оценка параметров модели осуществляется обобщенным методом наименьших квадратов (ОМНК):

Если модель гетероскедастична, то матрица Ω -диагональная. Тогда имеем:

В этом случае обобщенный метод наименьших квадратов называется взвешенным методом наименьших квадратов, поскольку мы «взвешиваем» каждое наблюдение с помощь коэффициента 1/σ _i.

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

где f(x_i) – квадратичная функция.

Далее по полученному уравнению рассчитывают теоретическое значения и определяют набор весов Затем вводят новые переменные и находят уравнение . Полученная оценка и есть оценка взвешенного метода наименьших квадратов.

Проверить модель на гетероскедастичность можно с помощью следующих тестов: ранговой корреляции Спирмена; Голдфельда-Квандта; Уайта; Глейзера.

Рассмотрим тест на гетероскедастичность, применяемый в случае, если ошибки регрессии можно считать нормально распределенными случайными величинами, - тест Голдфельда-Квандта.

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

Гипотеза о гомоскедастичности равносильна тому, что значение остатков е₁, …, е_т и

е_п-т+1, … е_п представляют совой выборочные наблюдения нормально распределенных случайных величин, имеющих одинаковые дисперсии.

Гипотеза о равенстве дисперсий двух нормально распределенных совокупностей проверяется с помощью F–критерия Фишера.

Расчетное значение вычисляется по формуле (в числителе всегда большая сумма квадратов):

.

Гипотеза о равенстве дисперсий двух наборов по т наблюдений (т.е. гипотеза об отсутствии гетероскедастичности остатков) отвергается, если расчетное значение превышает табличное F > F_{α; m-p; m-p}, где р – число регрессоров.

Мощность теста (вероятность отвергнуть гипотезу об отсутствии гетероскедастичности, когда гетероскедастичность действительно нет) максимальна, если выбирать т порядка п/3.

Тест Голдфельда-Квандта позволяет выявить факт наличия гетероскедастичности, но не позволяет описать характер зависимостей дисперсий ошибок регрессии количественно.

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

Как правило, если автокорреляция присутствует, то наибольшее влияние на последующее наблюдение оказывает результат предыдущего наблюдения. Наличие автокорреляции между соседними уровнями ряда можно определить с помощью теста Дарбина-Уотсона. Расчетное значение определяется по следующей формуле:

.

Затем по таблицам находятся пороговые значения d_в и d_н. если расчетное значение:

- d_в < d < 4–d_в, то гипотеза об отсутствии автокорреляции не отвергается (принимается);

- d_н < d < d_в или 4-d_в < d < 4-d_н, то вопрос об отвержении или принятии гипотезы остается открытым (расчетное значение попадает в зону неопределенности);

- 0 < d < d_н, то принимается альтернативная гипотеза о наличии отрицательной автокорреляции.

- 4-d_н < d < 4, то принимается альтернативная гипотеза о наличии отрицательной автокорреляции.

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

Поэтому наряду с тестом Дарбина-Уотсона для проверки наличия автокорреляции используются тест серии (Бреуша-Годфри), Q-тест Льюинга-Бокса и другие.

Наиболее распространенным приемом устранения автокорреляции во временных рядах является построение авторегрессионных моделей.

 

Пример 2. Рассмотрим полученную в предыдущем примере модель зависимости балансовой прибыли предприятия торговли Y (тыс. руб.) от следующих переменных:

Х₂ – фонд оплаты труда, тыс. руб.

Х₄ – объем продаж по безналичному расчету, тыс. руб.

Задание: Для полученной модели проверить выполнение условия гомоскедастичности остатков, применив тест Голдфельда-Квандта.

 

Решение.

Для выполнения этого задания снова воспользуемся «Пакетом анализа», встроенным в EXCEL.

В соответствии со схемой теста Голдфельда-Квандта упорядочим данные по возрастанию переменной Х₄, предполагая, что дисперсии ошибок зависят от величины этой переменной.

В нашем примере т = п/3 = 8.

Результаты дисперсионного анализа модели множественной регрессии, построенной по первым 8 наблюдениям (после ранжирования по возрастанию переменной Х₄), приведенные в таблице 4.

Таблица 4.

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

	df	SS	MS	F	Значимость F
Регрессия		5, 07Е+08	2, 53Е+08	20, 95996	0, 003707
Остаток		ЕSS1 = = 6, 04Е+07	1, 21Е+07
Итого		5, 67Е+08

 

Результаты дисперсионного анализа модели, построенной по последним 8 наблюдениям, приведены в таблице 5.

 

Таблица 5.

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

	df	SS	MS	F	Значимость F
Регрессия		1, 77Е+08		1, 111617	0, 398654
Остаток		ЕSS2 = = 3, 98Е+07
Итого		5, 75Е+08

 

Рассчитаем статистику F_расч = ЕSS₂ / ЕSS₁(т.к. ЕSS₂ > ЕSS₁). Для нашего примера получаем: F = 3, 987Е+08 / 6, 04Е+07 = 6, 58.

Для того, чтобы узнать табличное значение, воспользуемся встроенной в EXCEL функцией FРАСПОБР(0, 05; 6; 6) с параметрами 0, 05 – заданная вероятность ошибки гипотезы Н₀; т-р = 8-2 = 6; т-р = 6 – параметры распределения Фишера. Данная функция находится в категории «статистических» функций.

Статистика F_расч больше табличного значения FРАСПОБР(0, 05; 6; 6) = 4, 28. Следовательно, модель гетороскедастична.

 

Пример 3. Рассмотрим полученную в примере 1 модель зависимости балансовой прибыли предприятия Y (тыс. руб.) от следующих переменных:

Х₂ – фонд оплаты труда, тыс. руб.; Х₄ – объем продаж по безналичному расчету, тыс. руб.

 

Задание: проверить полученную модель на наличие автокорреляции остатков с помощью теста Дарбина-Уотсона.

 

Решение.

Прежде всего, по эмпирическим данным необходимо методом наименьших квадратов построить уравнение регрессии и определить значения отклонений для каждого наблюдения i (i = 1, 2, …, п).

Для этого в диалоговом окне Регрессия в группе Остатки следует установить одноименный флажок Остатки.

Зачем рассчитываем статистику Дарбина-Уотсона по формуле:

.

Результаты расчетов представлены в таблице 6.

Таблица 6

		^2	^2
11211, 00896			1, 3Е+08
9809, 816986	11211, 01	1963338, 9	9, 6Е+07
6652, 565001	9809, 817	9968240, 1	4, 4Е+07
4367, 949639	6652, 565	5219467, 4	1, 9Е+07
1141, 570741	4367, 95
2445, 881613	1141, 571	1701226, 8
687, 4294812	2445, 882	3092153, 9
140, 6630821	687, 4295	298953, 5	19786, 1
-4784, 81741	140, 6631		2, 3Е+07
-3182, 828283	-4784, 82	2566369, 2	1Е+07
-10324, 78476	-3182, 83		1, 1Е+0, 8
1880, 960336	-10324, 8
-2301, 490224	1880, 96
-6360, 626521	-2301, 49		4Е+07
-1887, 83539	-6360, 63
-1671, 617647	-1887, 84	46750, 112
1701, 17565	-1671, 62
149, 2560547	1701, 176	2408454, 4	22277, 4
-6106, 936579	149, 2561		3, 7Е+07
53, 14551195	-6106, 94		2824, 45
-4554, 494657	53, 14551		2, 1Е+07
-426, 4897698	-4554, 49
-5970, 720141	-426, 49		3, 6Е+07
7331, 218328	-5970, 72		5, 4Е+07
СУММА:	6, 5Е+08	6, 4Е+08

 

Таким образом, расчетное значение равно d = 6, 5Е+08 / 6, 4Е+08 = 1, 02.

По таблице критических точек распределение Дарбина-Уотсона для заданного уровня значимости α , числа наблюдений п и количества объясняющих переменных т определить два значения: d_н – нижняя граница и d_в – верхняя граница (таблица 7).

Таблица 7

⇐ Предыдущая12345678Следующая ⇒

Поделиться:

Популярное:

1. I. Понятие и система криминалистического исследования оружия, взрывных устройств, взрывчатых веществ и следов их применения.
2. III. РАСПРЕДЕЛЕНИЕ ЧАСОВ КУРСА ПО ТЕМАМ И ВИДАМ РАБОТ
3. V1: Понятие, объект, предмет и система криминологии
4. V2: Линейная модель множественной регрессии
5. V6: Нелинейные модели регрессии
6. V7: Система линейных одновременных уравнений
7. VI. Немецкая классическая философия.
8. Автоматизированная система телемеханического управления (АСТМУ)
9. Административная реформа и система органов исполнительно власти.
10. Административное право - публичное право. Административное право как отрасль права и система правового регулирования государственного управления.
11. Аксиологическое «Я» педагога как система ценностных ориентаций
12. Алмазно-расточный станок модели 2А78