Автор работы: Пользователь скрыл имя, 02 Декабря 2014 в 12:10, контрольная работа
Задание по эконометрическому моделированию стоимости квартир в Московской области
1. Рассчитайте матрицу парных  коэффициентов корреляции; оцените  статистическую значимость коэффициентов корреляции.
2. Постройте поле корреляции  результативного признака и наиболее  тесно связанного с ним фактора.
3. Рассчитайте параметры линейной  парной регрессии для всех факторов X.
4. Оцените качество каждой модели  через коэффициент детерминации, среднюю ошибку аппроксимации  и F-критерий Фишера. Выберите лучшую модель.
Таблица 5.
| Коэффициенты | Стандартная ошибка | t-статистика | P-Значение | Нижние 95% | Верхние 95% | Нижние 95,0% | Верхние 95,0% | |
| Y-пересечение | 1,944 | 0,249 | 7,809 | 0,00010 | 1,356 | 2,533 | 1,356 | 2,533 | 
| t (наблюдение) | 2,633 | 0,044 | 59,516 | 9,923E-11 | 2,529 | 2,738 | 2,529 | 2,738 | 
В соответствии с полученными расчетными данными модель регрессии в линейной форме будет выглядеть следующим образом:
Таким образом, в течение недели спрос на кредитные ресурсы увеличится на 2,633 млн. руб.
В соответствии с полученными результатами оценим качество линейной парной регрессии через коэффициент детерминации, среднюю ошибку аппроксимации и F-критерия Фишера.
Коэффициент детерминации в соответствии с расчетом равен . Вариация спроса на кредитные ресурсы (Y) на 99,8% объясняется изменением времени (t). Значение близко к единице, поэтому качество модели можно признать удовлетворительным.
Оценку значимости уравнения регрессии проведем с помощью F-критерия Фишера. В результате расчетов .
Определим табличное значение F-критерия с помощью функции MS Excel FPACПOБP (для , k1=1, k2=7). .
, следовательно, уравнение регрессии признается статистически значимым.
3. Произведем оценку 
3.1 Оценка адекватности построенно
где р – количество поворотных точек определяемых из графика остаточной компоненты, N – длина временного ряда; квадратные скобки означают, что от результата вычисления следует взять целую часть.
Для определения числа поворотных точек необходимо построить график остаточной компоненты. Построим с использованием имеющихся результатов расчета график остаточной компоненты (рис.2)
Рис. 2. График остаточной компоненты.
Исходя из построенного графика, определяем количество поворотных точек: р=5. При выполнении расчетов по приведенной выше формуле получаем, что 5>2, следовательно, неравенство соблюдается, ряд остатков можно считать случайным, а значит, свойство случайности возникновения отдельных отклонений от тренда выполняется.
3.2. Оценка адекватности 
Результаты расчета представлены в таблице 6.
Таблица 6.
| t (номер наблюдения) | Y (спрос, млн. руб.) | Предсказанное Y (спрос, млн. р) | E(t) | E(t)-E(t-1) | (E(t)-E(t-1))2 | E2(t) | 
| 1 | 5 | 4,578 | 0,422 | – | – | 0,18 | 
| 2 | 7 | 7,211 | -0,211 | -0,63 | 0,40 | 0,04 | 
| 3 | 10 | 9,844 | 0,156 | 0,37 | 0,13 | 0,02 | 
| 4 | 12 | 12,478 | -0,478 | -0,63 | 0,40 | 0,23 | 
| 5 | 15 | 15,111 | -0,111 | 0,37 | 0,13 | 0,01 | 
| 6 | 18 | 17,744 | 0,256 | 0,37 | 0,13 | 0,07 | 
| 7 | 20 | 20,378 | -0,378 | -0,63 | 0,40 | 0,14 | 
| 8 | 23 | 23,011 | -0,011 | 0,37 | 0,13 | 0,00 | 
| 9 | 26 | 25,644 | 0,356 | 0,37 | 0,13 | 0,13 | 
| Сумма | 1,88 | 0,82 | ||||
| d | 2,28 | |||||
При отсутствии автокорреляции d~2, а при полной автокорреляции равно 0 или 4. Следовательно, оценки, получаемые по критерию, являются не точечными, а интервальными. Верхние (d2) и нижние (d1) критические значения, позволяющие принять или отвергнуть гипотезу об отсутствии автокорреляции, зависят от количества уровней динамического ряда и числа независимых переменных модели. Значения этих границ для уровня значимости α=0,05: d1=l,08, d2=l,36. Так как d>d2 принимается нулевая гипотеза о равенстве нулю серийных корреляций и делается вывод об адекватности построенной модели. Свойство независимости выполняется.
3.3. Оценка адекватности 
,
где , в соответствии с результатами таблицы 3 имеем S=0,343.
В соответствии с расчетом (табл. 6) Еmах= 0,422, Еmin = –0,478. Тогда . Расчетное значение RS-критерия не попадает в табличный интервал [2,7; 3,7], следовательно, нормальный закон распределения выполняется. Модель по этому критерию адекватна.
4. Для оценки точности модели 
определим среднюю ошибку 
Результаты расчета представлены в таблице 8.
Таблица 8.
| Наблюдение | Y (спрос, млн. руб.) | E(t) | ABS(E(t)) | ABS(E(t)/Y) | 
| 1 | 5 | 0,422 | 0,422 | 0,084 | 
| 2 | 7 | -0,211 | 0,211 | 0,030 | 
| 3 | 10 | 0,156 | 0,156 | 0,016 | 
| 4 | 12 | -0,478 | 0,478 | 0,040 | 
| 5 | 15 | -0,111 | 0,111 | 0,007 | 
| 6 | 18 | 0,256 | 0,256 | 0,014 | 
| 7 | 20 | -0,378 | 0,378 | 0,019 | 
| 8 | 23 | -0,011 | 0,011 | 0,000 | 
| 9 | 26 | 0,356 | 0,356 | 0,014 | 
| Сумма | 0,225 | |||
| Eотн | 2,5 | |||
Средняя относительная ошибка построенной модели равна 2,5%, следовательно, модель имеет удовлетворительный уровень точности.
5. Осуществим точечный прогноз спроса на следующие две недели.
В течение первой недели (k=1, t=10) спрос будет равен Y(10)=28,274. В течение второй недели (k=2, t=11) спрос будет равен Y(11)=30,907.
На базе точечных прогнозов разрабатываем интервальные прогнозы. С этой целью рассчитывается ширина доверительного интервала:
где S – стандартная ошибка оценки, которая определяется по формуле:
В соответствии с результатами таблицы 3 имеем S = 0,343. Произведем расчет интервального прогноза для первой недели, для этого определим ширину доверительного интервала. При расчете используем Кр=1,08, n=9, m=l, k=l. Результаты расчета представлены в таблице 9.
Таблица 9.
| t (номер наблюдения) | |||
| 1 | –4 | 16 | |
| 2 | –3 | 9 | |
| 3 | –2 | 4 | |
| 4 | –1 | 1 | |
| 5 | 0 | 0 | |
| 6 | 1 | 1 | |
| 7 | 2 | 4 | |
| 8 | 3 | 9 | |
| 9 | 4 | 16 | |
| Сумма | 45 | 60 | |
| 
 | 5 | ||
| U(k) | 2,54 | ||
В результате расчета имеем U(k) = 2,54.
Таким образом, прогнозное значение Y(10)=28,274, будет находиться между верхней границей, равной 28,274+2,54=30,814, и нижней границей, равной 28,274–2,54=25,734.
Произведем расчет интервального прогноза для второй недели, для этого определим ширину доверительного интервала. При расчете используем Кр=1,08, n=9, m=l, k=2. Результаты расчета представлены в таблице 10.
Таблица 10.
| t (номер наблюдения) | ||||||
| 1 | –4 | 16 | ||||
| 2 | –3 | 9 | ||||
| 3 | –2 | 4 | ||||
| 4 | –1 | 1 | ||||
| 5 | 0 | 0 | ||||
| 6 | 1 | 1 | ||||
| 7 | 2 | 4 | ||||
| 8 | 3 | 9 | ||||
| 9 | 4 | 16 | ||||
| Сумма | 45 | 60 | ||||
| 
 | 5 | |||||
| U(k) | 2,68 | |||||
В результате расчета имеем U(k) = 2,68.
Таким образом, прогнозное значение Y(ll)= 30,907, будет находиться между верхней границей, равной 30,907+2,68=33,587, и нижней границей, равной 30,907–2,68=28,227.
Таблица прогнозов
| Неделя наблюдения | Прогноз | Нижняя граница | Верхняя граница | 
| 10 | 28,274 | 25,734 | 30,814 | 
| 11 | 30,907 | 28,227 | 33,587 | 
6. На графике (рис.3) представлены графически фактические значения показателя, результаты моделирования, результаты прогнозирования.