Автор работы: Пользователь скрыл имя, 07 Апреля 2013 в 14:12, лабораторная работа
Строим линейную модель Y=a + b *X с наиболее информативным фактором Х4.
Для удобства выполнения расчетов мы предварительно упорядочили всю таблицу исходных данных по возрастанию факторной переменной Х4 (Данные → Сортировка).
Используя программу РЕГРЕССИЯ, нашли коэффициенты модели.
Коэффициенты |
Стандартная ошибка |
t-статистика |
P-Значение |
Нижние 95% |
Верхние 95% |
Нижние 95,0% |
Верхние 95,0% | |
Y-пересечение |
13,50586 |
13,12787 |
1,028793 |
0,323868 |
-15,0973 |
42,10904 |
-15,0973 |
42,10904 |
Х4 |
1,924704 |
0,597716 |
3,220101 |
0,007353 |
0,622394 |
3,227015 |
0,622394 |
3,227015 |
Для этой модели остаточная сумма квадратов SS1 = 2376,18.
С помощью программы РЕГРЕССИЯ построим модель по последним 14-ти наблюдениям (регрессия-2):
ВЫВОД ИТОГОВ | |
Регрессионная статистика | |
Множественный R |
0,840751 |
R-квадрат |
0,706862 |
Нормированный R-квадрат |
0,682434 |
Стандартная ошибка |
30,2291 |
Наблюдения |
14 |
Дисперсионный анализ |
|||||
df |
SS |
MS |
F |
Значимость F | |
Регрессия |
1 |
26442,04 |
26442,04 |
28,93639 |
0,000166 |
Остаток |
12 |
10965,59 |
913,7988 |
||
Итого |
13 |
37407,62 |
Коэффициенты |
Стандартная ошибка |
t-статистика |
P-Значение |
Нижние 95% |
Верхние 95% |
Нижние 95,0% |
Верхние 95,0% | |
Y-пересечение |
-65,2289 |
37,77644 |
-1,72671 |
0,109847 |
-147,537 |
17,07886 |
-147,537 |
17,07886 |
Переменная X 1 |
3,283442 |
0,61039 |
5,379256 |
0,000166 |
1,953517 |
4,613366 |
1,953517 |
4,613366 |
Для этой модели остаточная сумма квадратов SS2 =10965,59.
Рассчитаем статистику критерия: F= SSmax/ SSmin=10965,59/2376,18= =4,61.
Критическое значение при уровне значимости α = 5% и числах степеней свободы k1 = k2 =14−1−1=12 составляет Fкр=2,69 (функция FРАСПОБР).
Сравним: F=4,61>Fкр=2,69 , следовательно, свойство постоянства дисперсии остатков не выполняется, модель является гетероскедастичной.
В) Для проверки независимости уровней ряда остатков используем критерий Дарбина-Уотсона:
1) Определяем d-статистику по формуле:
По столбцу остатков с помощью функции СУММКВРАЗН можно определить , а с помощью функции СУММКВ – .
Таким образом, d = СУММКВРАЗН/ СУММКВ=1,998.
По таблице критических точек Дарбина-Уотсона определяем значения d1 и d2: d1= 1,425, а d2=1,54.
Т.к. d2=1,54<d =1,998<2, следовательно, свойство независимости остатков выполняется.
В учебных целях проверим выполнение свойства независимости ряда остатков по первому коэффициенту автокорреляции:
С помощью функции СУММПРОИЗВ найдем для остатков:
.
Т.о., r(1)= СУММПРОИЗВ/СУММКВ=-0,006
Критическое значение для коэффициента автокорреляции определяется как отношение:
и составляет для данной задачи r кр=0,31.
Т.к. | r(1)|= 0,006< r кр=0,31, то свойство независимости остатков выполняется.
Г) Соответствие ряда остатков нормальному закону распределения проверим с помощью R/S – критерия.
С помощью функций МАКС и МИН для ряда остатков определим E max = 60,97 и E min = −49,69.
Стандартная ошибка модели найдена программой РЕГРЕССИЯ и составляет SЕ =25,53 (таблица «Регрессионная статистика»).
Тогда R/S=(60,97-(−49,69))/ 25,53 =4,33.
Определим верхнюю и нижнюю границы по таблице критических границ: НГ=3,625, ВГ=5,1.
4,33∈(3,625;5,1), значит, для построенной модели свойство нормального распределения остаточной компоненты выполняется.
ВЫВОД: Проведенная проверка показала, что для модели выполняются не все условия Гаусса-Маркова. Для улучшения модели необходимо взять более однородную выборку. Далее модель можно использовать в учебных целях, т.к. она не является адекватной.