Автор работы: Пользователь скрыл имя, 21 Ноября 2011 в 14:09, контрольная работа
1. Построим поле корреляции, для чего отложим на плоскости в прямоугольной системе координат точки (xi и yi)
ОДНОФАКТОРНЫЙ РЕГРЕССИОННО-КОРРЕЛЯЦИОННЫЙ АНАЛИЗ 3
АНАЛИЗ И ПРОГНОЗИРОВАНИЕ ВРЕМЕННЫХ РЯДОВ 18
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 35
Непосредственное вычисление (вычисление «вручную») вектора оценок параметров регрессии а весьма громоздко, т.к. матрица независимых переменных Х имеет довольно высокую размерность (20 х 3), матрица Y – размерности (20 х 1). В табл. 8 приведены размерности матриц – результатов промежуточных действий.
Таблица 8
ХТ | (3 х 20) |
ХТХ | (3 х 3) |
(ХТХ)–1 | (3 х 3) |
(ХТХ)–1ХТ | (3 х 20) |
(ХТХ)–1ХТY | (3 х 1) |
Задача существенно упрощается при использовании средств Excel. Операции, предписанные формулой целесообразно проводить с помощью следующих встроенных в Excel функций:
Для вычисления вектора оценок параметров регрессии а в Excel необходимо выполнить следующие действия:
1. Ввести данные .
2. Выделить диапазон ячеек для записи вектора а, соответствующий его размерности (3 х 1) .
3. Используя встроенные в Excel функции, ввести формулу, определяющую вектор а.
4. Нажать одновременно клавиши CTRL + SHIFT + ENTER. Появится результат (рисунок 14).
Таким образом, имеем
Уравнение регрессии зависимости объема реализации от затрат на рекламу и индекса потребительских расходов можно записать в виде:
Рисунок 14 - Результат вычислений – вектор оценок
параметров
регрессии а
Расчетные значения Y определяются путем последовательной подстановки в эту модель значений факторов, взятых для каждого момента времени t.
Применение инструмента Регрессия
Для проведения регрессионного анализа с помощью Excel выполните следующие действия:
1) выберите команду Сервис => Анализ данных;
2) в диалоговом окне Анализ данных выберите инструмент Регрессия. Щелкните по кнопке ОК;
3) в диалоговом окне Регрессия в поле «Входной интервал Y» введите адрес диапазона ячеек, который представляет зависимую переменную Y. В поле «Входной интервал X» введите адреса одного или нескольких диапазонов, которые содержат значения независимых переменных (переменные Х3, Х5). Если выделены заголовки столбцов, то установить флажок «Метки в первой строке»;
5) выберите параметры вывода. В данном примере – установите переключатель «Новая рабочая книга»;
6) в поле «Остатки» поставьте необходимые флажки;
Результаты
представлены на рисунке 15 и заключены
в таблицах.
Пояснения к таблице «Регрессионная статистика» (рисунок 15)
Таблица 9
Регрессионная статистика | |||
Наименования
в отчете Excel |
Принятые
наименования |
Формула | |
Множественный R |
Коэффициент множественной корреляции, индекс корреляции | | |
R – квадрат |
Коэффициент детерминации, R2 |
||
Нормированный R2 | Скорректированный R2 | ||
Стандартная ошибка | Стандартная ошибка оценки | ||
Наблюдения | Количество наблюдений, n | n |
Рисунок 15 -Результаты регрессионного анализа, проведенного
с
помощью Excel
Пояснения к таблице «Дисперсионный анализ» (рисунок 15)
Таблица 10
D f – число степеней свободы | SS – сумма квадратов | MS |
F-критерий Фишера | |
Регрессия |
k |
|||
Остаток | n– k –1 | |||
Итого | n – 1 |
Во
втором столбце таблицы дисперсионного
анализа (рисунок 15) содержатся коэффициенты
уравнения регрессии а0,
а1, а2, в третьем столбце
содержатся стандартные ошибки коэффициентов
уравнения регрессии, в четвертом –
f-статистика, используемая для проверки
значимости коэффициентов уравнения регрессии.
Рисунок
16 - График остатков
Оценка качества модели
В таблице «вывод остатка» (рисунок 15) приведены вычисленные по модели значения и значения остаточной компоненты е.
Исследование
на наличие автокорреляции остатков
проведем с помощью d-критерия
Дарбина – Уотсона. Для определения величины
d-критерия воспользуемся расчетной
таблицей 11.
Таблица 11
Наблюю
дение |
Y | Предсказанное Объем реализ. | ε(t) | ε(t)2 | (ε(t)-ε(t-1))2 | ε(t)-ε(t-1) | (Y-Yср)2 |
1 | 328 | 300,3189493 | 27,6810507 | 766,240569 | 290,7025 | ||
2 | 336 | 303,2864106 | 32,7135894 | 1070,17893 | 25,32644612 | 905,5465286 | 81,9025 |
3 | 342 | 306,7243707 | 35,2756293 | 1244,37002 | 6,564048351 | 1153,992455 | 9,3025 |
4 | 276 | 306,6521366 | -30,6521366 | 939,553478 | 4346,470321 | -1081,27341 | 4767,9025 |
5 | 296 | 315,0840215 | -19,0840215 | 364,199878 | 133,8212861 | 584,9660352 | 2405,9025 |
6 | 296 | 315,1020801 | -19,1020801 | 364,889463 | 0,00032611 | 364,5445073 | 2405,9025 |
7 | 307 | 334,2370809 | -27,2370809 | 741,858575 | 66,17823841 | 520,2848996 | 1447,8025 |
8 | 312 | 329,8581231 | -17,8581231 | 318,912559 | 87,96484995 | 486,403142 | 1092,3025 |
9 | 400 | 380,9750845 | 19,0249155 | 361,94741 | 1360,358533 | -339,749282 | 3019,5025 |
10 | 395 | 345,5822515 | 49,4177485 | 2442,11386 | 923,7242959 | 940,1684885 | 2495,0025 |
11 | 391 | 339,647329 | 51,352671 | 2637,09682 | 3,743925194 | 2537,733378 | 2111,4025 |
12 | 391 | 345,0395186 | 45,9604814 | 2112,36585 | 29,07570875 | 2360,193479 | 2111,4025 |
13 | 250 | 325,4520774 | -75,4520774 | 5693,01599 | 14741,00944 | -3467,8138 | 9034,5025 |
14 | 283 | 345,0846649 | -62,0846649 | 3854,50562 | 178,6877182 | 4684,416945 | 3850,2025 |
15 | 358 | 350,4768545 | 7,52314549 | 56,5977181 | 4845,247268 | -467,071967 | 167,7025 |
16 | 363 | 380,2878834 | -17,2878834 | 298,870914 | 615,5871567 | -130,059262 | 322,2025 |
17 | 374 | 368,4180384 | 5,58196161 | 31,1582955 | 523,029813 | -96,5003018 | 838,1025 |
18 | 395 | 392,9632873 | 2,03671267 | 4,14819849 | 12,5687901 | 11,36885192 | 2495,0025 |
19 | 402 | 404,7247813 | -2,72478129 | 7,42443307 | 22,67182469 | -5,54959657 | 3243,3025 |
20 | 406 | 411,0850564 | -5,08505638 | 25,8577984 | 5,570898505 | 13,85566648 | 3714,9025 |
Сумма | 6901 | 6901 | 6,2528E-13 | 23335,3064 | 27927,60089 | 8975,456758 | 45904,95 |
Имеем:
В качестве критических табличных уровней при n = 20, двух объясняющих факторах при уровне значимости γ = 0,05 возьмем величины dL = 1,10 и dU =1,54.
Расчетное значение попало в интервал от dL = 1,1 до dU =1,54 (рисунок 17)
Есть
положительная автокорреляция остатков.
Н0 отклоняется. С вероятностью (1 – γ) принимается Н1. |
Зона неопределенности |
Нет оснований
отклонять Н0
(автокорреляция остатков отсутствует) |
Зона неопределенности |
Есть отрицательная
автокорреляция остатков.
Н0 отклоняется. С вероятностью (1 – γ) принимается Н1*. |
0
dL d dU 2
4–dU
4– dL 4
Рисунок 17 - Сравнение расчетного значения d-критерия
Дарбина
– Уотсона с критическими значениями
dL и dU
Так как расчетное значение d-критерия Дарбина–Уотсона попало в зону неопределенности, то нельзя сделать окончательный вывод об автокорреляции остатков по этому критерию.
Для определения степени автокорреляции вычислим коэффициент автокорреляции и проверим его значимость при помощи критерия стандартной ошибки. Стандартная ошибка коэффициента корреляции рассчитывается по формуле:
Коэффициенты автокорреляции случайных данных должны обладать выборочным распределением, приближающимся к нормальному с нулевым математическим ожиданием и средним квадратическим отклонением, равным
Если коэффициент автокорреляции первого порядка r1
находится в интервале
–1,96 · 0,224 < r1 < 1,96· 0,224,
то можно считать, что данные не показывают наличие автокорреляции первого порядка.
Используя расчетную таблицу 6, получаем:
Так как
–0,4782 < r1 = 0,3846 < 0,4782,
то свойство независимости остатков выполняется.
Вычислим для построенной модели множественный коэффициент детерминации
= 1 –
23335,3/45904,95 = 0,4917.
множественный коэффициент детерминации показывает долю вариации результативного признака под воздействием включенных в модель факторов Х2 и Х5. Т.о., около 50 % вариации зависимой переменной (объема реализации) в построенной модели обусловлено влиянием включенных факторов Х2 (расходы на рекламу) и Х5 (индекс потребительских расходов).
Проверку значимости уравнения регрессии проведем на основе F-критерия Фишера
Табличное значение F-критерия при доверительной вероятности 0,95, степенями свободы
v 1 = k = 2
и
v2 = (п – k – 1) = 20 – 2 – 1 = 17
составляет Fтабл = 3,5915.
Поскольку
то уравнение регрессии следует признать адекватным.
Значимость
коэффициентов уравнения
Табличное значение, t-критерия Стьюдента при уровне значимости 0,05 и степенях свободы (20 – 2 –1) = 17 составляет . Так как
Информация о работе Анализ и прогнозирование временных рядов