Автор работы: Пользователь скрыл имя, 29 Мая 2013 в 00:25, контрольная работа
Задача 1 По предприятиям легкой промышленности региона получена информация, характеризующая зависимость объема выпуска продукции (y, млн. руб.) от объема капиталовложений (x, млн. руб.). Требуется: Для характеристики y от x построить следующие модели: линейную, степенную, показательную, гиперболическую. Оценить каждую модель, определив: индекс корреляции, среднюю относительную ошибку, коэффициент детерминации, F- критерий Фишера.
Составить сводную таблицу вычислений, выбрав лучшую модель, дать интерпретацию рассчитанных характеристик.
Рассчитать прогнозное значение результативного признака, если прогнозное значение фактора увеличится на 10% относительно среднего уровня. Результаты расчетов отобразить на графике.
- коэффициент детерминации R2. Расчетное значение этого коэффициента приводится в табл.4.1 (R-квадрат) R2 = 0,712. Он показывает долю вариации результативного признака под воздействием изучаемых факторов х1 и х2. Следовательно, около 71,2% вариации зависимой переменной учтено в модели и обусловлено влиянием включенных факторов.
- средние коэффициенты эластичности Э1 и Э2. Учитывая, что коэффициенты регрессии невозможно использовать для непосредственной оценки влияния факторов на зависимую переменную из-за различия единиц измерения, рассчитаем средние коэффициенты эластичности Э1 и Э2 для факторов х1 и х2 соответственно по формулам:
где, - коэффициент регрессии для фактора хj, - среднее значение для фактора хj (табл.1), - среднее значение для результирующей переменной y:
Э1 = 0,089*59,4/78,8 = 0,067
Э2 = 0,833*47,6/78,8 = 0,503.
Коэффициент эластичности показывает, насколько процентов изменяется зависимая переменная у от своего среднего значения при изменении фактора хj на один процент от своего среднего значения .
- бетта-коэффициенты ( ). Учитывая, что коэффициент регрессии невозможно использовать для непосредственной оценки влияния факторов на зависимую переменную из-за различия единиц измерения, используем бетта-коэффициенты, значения которых рассчитываются по формулам:
где - коэффициент регрессии для фактора хj, - среднеквадратическое отклонение фактора хj; - среднеквадратическое отклонение зависимой переменной у. Расчет среднеквадратических отклонений для факторов хj и зависимой переменной у представлен в таблице 5.
Таблица 5.
(Х1-Хср)2 |
(Х2-Хср)2 |
(У-Уср)2 |
11,56 |
309,76 |
353,44 |
129,96 |
57,76 |
116,64 |
54,76 |
12,96 |
219,04 |
1,96 |
384,16 |
46,24 |
43,56 |
5,76 |
0,64 |
6,76 |
70,56 |
84,64 |
129,96 |
5,76 |
125,44 |
43,56 |
70,56 |
10,24 |
112,36 |
153,76 |
174,24 |
73,96 |
207,36 |
231,04 |
608,4 |
1278,4 |
1361,6 |
Корень из S | ||
Sx12 |
67,60 |
8,22 |
Sx22 |
142,04 |
11,92 |
Sy2 |
151,29 |
12,30 |
Значения бетта-коэффициентов составят:
Бета-коэффициент с
- дельта-коэффициенты ( ) оценивают долю влияния фактора в суммарном влиянии всех факторов и рассчитываются по формулам:
где - коэффициент парной корреляции между фактором хj и зависимой переменной у (табл.2), - бетта-коэффициент фактора хj; R2 - коэффициент детерминации (табл. 4.1).
4. Оценка надежности и качества уравнения регрессии
Качество модели оценивается по адекватности и точности на основе анализа остатков регрессии . Расчетные значения ŷi получаются путем подстановки в модель фактических значений всех включенных в модель факторов. В таблице 4.4 приведены вычисленные по модели (1) расчетные значения ŷi и значения остаточной компоненты , представляющие собой разности между фактическими уi и расчетными ŷi значениями, т.е. =yi - ŷi.
Анализ остатков позволяет получить представление о правильности выбора вида модели и метода оценки коэффициентов. Наглядное представление об остатках можно получить с помощью графика остатков.
График остатков хорошо показывает
резко отклоняющиеся от модели наблюдения
– выбросы. По данным графика остатков
выбросами можно считать
Проверка значимости полученного уравнения множественной линейной регрессии в целом осуществлялась с помощью F-критерия Фишера.
Расчетное значение F-критерия Фишера определялось по формуле:
Табличное значение F-критерия Фишера при доверительной вероятности 0,05 при 1 = k = 2 и 2 = n – k – 1 = 10 – 2 – 1 = 7 составляет 4,74.
Поскольку Fрасч. > Fтабл., уравнение регрессии следует признать адекватным.
5. Оценить с помощью t-критерия Стьюдента статистическую значимость уравнения множественной регрессии
Значимость коэффициентов уравнения регрессии a0, a1, a2 оценим с использованием t-критерия Стьюдента:
ta0 = 1,859; ta1 = 0,238; ta2 = 3,206.
Расчетные значения t-критерия Стьюдента для коэффициентов уравнения регрессии а0, a1, a2 приведены в четвертом столбце таблицы 4.3 (использована программа Excel). Табличное значение t-критерия Стьюдента t(0,05; 7) равно tтабл = 2,36. Так как |tрасч.| < tтабл., то коэффициенты а0, a1 не значимы, а коэффициент a2 значим, так как |tрасч.| > tтабл.
Прогнозируемые точечные значения результирующего показателя упрогн по модели (1) можно определить с помощью методов экспертных оценок, с помощью средних абсолютных приростов или вычислить на основе экстраполяционных методов. В соответствии с методом экстраполяции необходимо определить прогнозируемые значения по факторам X1 и X2 на первом и на втором шагах. То есть необходимо определить значения факторов X1 и X2 в одиннадцатом X1,11, X2,11 и двенадцатом X1,12 и X2,12 наблюдениях по наилучшим подобранным моделям. Наилучшей считается та модель, для которой значение индекса (коэффициента) детерминации R2 максимально. В качестве анализируемых рассмотрим линейную, квадратичную и степенную зависимости значений факторов X1 и X2 от номера наблюдения (времени t). Исходные данные представлены в табл.1.
Для фактора X1 Ставки по кредитам выбрана модель вида:
X1 = 0,1212*t2 + 0,4727*t + 52,133,
имеющая наибольшее значение индекса детерминации R2 = 0,4551. По данной модели при t = 11 и t = 12 получен прогноз на 2 шага вперед. График модели временного ряда Ставки по кредитам приведен на рис. 6.1. Прогнозируемые значения по фактору X1 составили:
X1,11 = 0,1212*112 + 0,4727*11+ 52,133 = 71,9979
X1,12 = 0,1212*122 + 0,4727*12 + 52,133 = 75,2582.
Для фактора X2 Ставки по депозитам выбрана модель вида:
X1 = - 0,0227*t2 + 3,6439*t + 28,433,
которая имеет наибольшее значение индекса детерминации R2 = 0,7436
Рис. 6.1 Ставки по кредитам
. По данной модели при t = 11 и t = 12 получен прогноз на 2 шага вперед. График модели временного ряда Ставки по депозитам приведен на рис. 6.2. Прогнозируемые значения по фактору X2 составили:
X2,11 = - 0,0227*112 + 3,6439*11 + 28,433 = 65,9735
X2,12 = - 0,0227*122 + 3,6439*12 + 28,433 = 68,891 .
Рис. 6.2. Ставки по депозитам
Для получения точечных оценок прогноза зависимой переменной упрогн в одиннадцатом упрогн., 11 и двенадцатом упрогн., 12 наблюдениях по модели:
у = 33,846 + 0,089*x1 + 0,833*x2
подставим в нее найденные прогнозные значения факторов X1, X2 соответственно в одиннадцатом X1,11, X2,11 и двенадцатом X1,12, X2,12 наблюдениях:
упрогн., 11 = 33,846 + 0,089*71,9979 + 0,833*65,9735= 95,2097,
упрогн., 12 = 33,846 + 0,089*75,2582 + 0,833*68,891= 97,9302.
Предельная ошибка прогноза, которая в 95% случаев не будет превышена, составит:
где
- среднеквадратическое
= 7,488 (табл.4.1. стандартная ошибка);
tкр – критическое или табличное значение критерия Стьюдента при уровне значимости = 0,05, т.е. при доверительной вероятности 95%, и числе степеней свободы = n - m – 1 = 10 - 2 - 1 = 7; tкр (0,05; 7) = 2,36;
X - матрица значений факторов X0, X1, X2, включенных в модель (табл.3); размерность матрицы X – (10 3);
XT - транспонированная матрица X; размерность матрицы XT – (3 10);
(XT*X)-1 - матрица обратная к матрице XT*X; размерность матриц XT*X и (XT*X)-1 – (3 3);
Xпр - матрица-столбец прогнозируемых значений факторов; размерность матрицы Xпр – (3 1);
XпрT - транспонированная матрица прогнозируемых значений факторов Xпр; размерность матрицы XпрT – (1 3).
Расчеты интервального прогноза представлены в приложениях (в таблицах Excel).
Предельная ошибка прогноза на первом шаге, т.е. в одиннадцатом наблюдении составит:
Предельная ошибка прогноза на втором шаге, т.е. в двенадцатом наблюдении составит:
Доверительный интервал прогноза будет иметь следующие границы:
- верхняя граница прогноза: упрогн. +
- нижняя граница прогноза: упрогн. -
Результаты расчетов прогнозных оценок модели регрессии представлены в таблице 6.
Результаты прогнозирования
Шаг |
Точечн. прогноз |
Предел. ошибка |
Нижняя граница |
Верхняя граница |
1 |
95,2097 |
11,586 |
83,6236 |
106,7959 |
2 |
97,9302 |
13,512 |
84,4187 |
111,4417 |
1. Построение линейной модели |
t |
y |
x |
y*x |
x*x |
yi-yср |
(yi-yср)2 |
xi-xcр |
(xi-xcр)2 |
yрасч |
ei=yi-yрасч |
abs(ei/yi) *100% |
(yi-yср)* (xi-xср) |
(y-yрасч)/y |
1 |
40 |
22 |
880 |
484 |
-13,429 |
180,327 |
-15,857 |
251,449 |
40,792 |
-0,792 |
1,979 |
212,939 |
-0,020 |
2 |
44 |
28 |
1232 |
784 |
-9,429 |
88,898 |
-9,857 |
97,163 |
45,573 |
-1,573 |
3,575 |
92,939 |
-0,036 |
3 |
48 |
30 |
1440 |
900 |
-5,429 |
29,469 |
-7,857 |
61,735 |
47,167 |
0,833 |
1,735 |
42,653 |
0,017 |
4 |
52 |
32 |
1664 |
1024 |
-1,429 |
2,041 |
-5,857 |
34,306 |
48,761 |
3,239 |
6,229 |
8,367 |
0,062 |
5 |
56 |
44 |
2464 |
1936 |
2,571 |
6,612 |
6,143 |
37,735 |
58,324 |
-2,324 |
4,150 |
15,796 |
-0,041 |
6 |
64 |
51 |
3264 |
2601 |
10,571 |
111,755 |
13,143 |
172,735 |
63,903 |
0,097 |
0,152 |
138,939 |
0,002 |
7 |
70 |
58 |
4060 |
3364 |
16,571 |
274,612 |
20,143 |
405,735 |
69,481 |
0,519 |
0,741 |
333,796 |
0,007 |
итого |
374 |
265 |
693,714 |
1060,857 |
0,000 |
18,562 |
845,429 |
-0,008 | |||||
ср. знач |
53,429 |
37,857 |
2143,429 |
1584,714 |
2,652 |
||||||||
дисперия |
|||||||||||||
Урасч.=а+bx |
b |
0,797 |
|||||||||||
a |
23,259 |
||||||||||||
ry,x |
0,986 |
||||||||||||
линейный коэф парной регрессии |
|||||||||||||
R2=ry,x2 |
0,971 |
||||||||||||
коэф детерминации |
|||||||||||||
F |
168,714 |
||||||||||||
оценка значимости ур-я регрессии |
|||||||||||||
Еср.отн |
2,652 |
||||||||||||
средняя относительная ошибка |
|||||||||||||