По
данным табл. 2.15 уравнение тренда
=779506,87+9318,27t,
где ti отсчитываются от начала ряда. Взвешенные индексы
сезонности составляют:
для I квартала:
для II квартала:
для III квартала:
для IV
квартала:
Умножив
уровни тренда соответствующих кварталов на эти средние индексы сезонности (в
разах), получаем уровни, учитывающие и тренд, и сезонные колебания,
обозначенные как
.
Далее
можно вычислить средние за 3 года показатели силы сезонных колебаний:
=55195,01
=
что говорит об умеренной колеблемости. При этом учтена потеря двух степеней свободы колебаний (по числу параметров линейного тренда).
С учетом наличия сезонных колебаний общую сумму квадратов отклонений фактических уровней динамического ряда от среднего уровня за весь изучаемый период можно разложить на составляющие элементы:
.
- сумму квадратов отклонений за счет тренда (за счет
фактора времени);
.
- сумму квадратов отклонений за сезонности.
.
- сумму квадратов отклонений за счет случайных
колебаний.
Таким
образом (формула (2.21)):
(2.21)
Рассчитаем
эти суммы для данного примера (табл. 2.14).
Таблица
2.14 Расчет величин для разложения общей суммы квадратов отклонений фактических
уровней от их средней (
= 16523,75 )
|
Год |
Квартал |
Общая площадь квартир введенных в действие ЖСК, м2, y |
Выровненные уровни yt=779506,87 +9318,27t |
Тренд с учетом сезонности |
|
|
|
|
|
2002 |
I |
15986,33 |
46626,94 |
12910,03 |
938846981,2 |
906199037,9 |
1136830020 |
9463621,69 |
|
|
II |
16023,78 |
50058,01 |
13158,87 |
1158328812 |
1124543240 |
1361546533 |
8207709,31 |
|
|
III |
16550,21 |
53489,08 |
13712,06 |
1364480117 |
1366431925 |
1582211320 |
8055095,42 |
|
|
IV |
16482,68 |
56920,15 |
13684,98 |
1635188980 |
1631865093 |
1869279925 |
7827125,29 |
|
2003 |
I |
15832,4 |
60351,22 |
16710,00 |
1981925334 |
1920842744 |
1904556083 |
770181,76 |
|
|
II |
16005,9 |
63782,29 |
16766,61 |
2282583441 |
2233364877 |
2210474166 |
578679,704 |
|
|
III |
16455,14 |
67213,36 |
17230,31 |
2576396898 |
2569431493 |
2498305287 |
600888,529 |
|
|
IV |
16176,56 |
70644,43 |
16984,63 |
2966748862 |
2929042592 |
2879374136 |
652977,125 |
|
2004 |
I |
16897,5 |
74075,5 |
20509,97 |
3269323684 |
3312198173 |
2869266004 |
13049939,5 |
|
|
II |
17045,6 |
77506,57 |
20374,34 |
3655528893 |
3264091705 |
11080510 |
|
|
|
III |
17561,9 |
80937,64 |
20748,56 |
4016484421 |
4149142784 |
3622725351 |
10154802 |
|
|
IV |
17267 |
84368,71 |
20284,28 |
4502639485 |
4602931813 |
4106814168 |
9103978,6 |
|
Итого |
12 |
198285 |
785973,9 |
203074,66 |
30348475908 |
30464892009 |
29305474699 |
79545508,9 |
Используя формулу (2.21), разложим общую сумму квадратов отклонений уровня от средней величины за весь период на составляющие элементы. При этом ввиду округления всех уровней в табл. 2.16 точного равенства не получится, но расхождение сумм квадратов менее одной тысячной совершенно не существенно для выводов о процессе.
Общая сумма квадратов или 100%
в том числе за счет тренда
за счет сезонности
за счет случайной колеблемости
На основе полученных данных можно сделать вывод, что
случайные колебания в исходном ряду были весьма незначительны. Основные факторы
колеблемости уровней исследуемого ряда - тренд и сезонность.
.3 Автокорреляция в рядах
динамики. Построение моделей авторегрессии
Во
многих рядах динамики можно наблюдать зависимость t-го уровня
от предшествующих
.
Например, урожайность сельскохозяйственных культур в отдельные годы также может
быть связана с урожайностью в предшествующие периоды.
Зависимость между последовательными (соседними) уровнями ряда динамики называется в статистике автокорреляцией. Исследование рядов на автокорреляцию - одна из частных, но важных задач при статистическом изучении рядов динамики. В частности, если установлено наличие автокорреляции, то эту зависимость можно выразить уравнением авторегрессии. В отдельных случаях приходится устранять влияние автокорреляции на взаимосвязь между исследуемыми показателями. Так возникает необходимость измерения автокорреляции.
Измерить
автокорреляцию между уровнями ряда можно с помощью коэффициента автокорреляции
, исчисленный по формуле парного линейного
коэффициента корреляции:
(2.22)
Коэффициент
автокорреляции
можно рассчитывать либо между соседними уровнями,
либо между уровнями, сдвинутыми на любое число единиц времени т. Этот сдвиг,
именуемый временном лагом, определяет порядок коэффициента автокорреляции: 1-го
порядка при т == 1, т.е. между соседними уровнями; 2-го порядка при т == 2,
т.е. при сдвиге уровней на 2 периода, и т.д.
Коэффициент автокорреляции 1-го порядка.
Если
исходные фактические уровни ряда, относящиеся к определенному моменту времени
(или периоду) t обозначить через
, то
сдвинутые уровни (в зависимости от направления сдвига) соответственно
обозначают
или
Тогда
формулу коэффициента автокорреляции можно записать в двух вариантах:
(2.23)
Предпочтение
обычно отдается второй формуле. При достаточно большом числе уровней ряда
значения средних уровней и средних квадратических
отклонений у исходного и сдвинутого рядов практически совпадают, т. е.
и
Используя
эти равенства и отдавая предпочтение средней
и
дисперсии
, рассчитанным для всех
членов
исходного ряда, получим приближенную формулу коэффициента автокорреляции:
(2.25)
или
тождественную ей
(2.26)
Чтобы
иметь возможность пользоваться данными формулами для коротких рядов, у которых
первый и последний уровни отличаются незначительно, сдвинутый ряд условно
дополняют, принимая
(чтобы сдвинутый ряд не укорачивался и чтобы средний
уровень и дисперсия одного ряда были соответственно равны среднему уровню и
дисперсии второго ряда).
Расчет коэффициента автокорреляции приведен на условном примере, известны данные о размерах посевной площади по Тюменской области Таблица 2.15
|
Год |
Посевная площадь, млн. га.
(фактические уровни) |
|
|
|
|
1990 |
1,6 |
(0,78) |
1,248 |
2,56 |
|
1991 |
1,57 |
1,6 |
2,512 |
2,46 |
|
1992 |
1,49 |
1,57 |
2,339 |
2,22 |
|
1993 |
1,38 |
1,49 |
2,056 |
1,91 |
|
1994 |
1,23 |
1,38 |
1,697 |
1,51 |
|
1995 |
1,19 |
1,23 |
1,464 |
1,42 |
|
1996 |
1,18 |
1,19 |
1,404 |
1,39 |
|
1997 |
1,01 |
1,18 |
1,192 |
1,02 |
|
1998 |
0,96 |
1,01 |
0,969 |
0,92 |
|
1999 |
0,94 |
0,96 |
0,902 |
0,88 |
|
2000 |
0,93 |
0,94 |
0,874 |
0,86 |
|
2001 |
0,91 |
0,93 |
0,846 |
0,83 |
|
2002 |
0,85 |
0,91 |
0,774 |
0,72 |
|
2003 |
0,78 |
0,85 |
0,663 |
0,61 |
|
|
|
|
|
|
По итоговым данным табл. 2.1.5. найдены:
;
;
;
Полученные
значения необходимо подставить в формулу (2.25)
Тат
же результат будет получен и при использовании формулы (2.26)
Найденное значение коэффициента автокорреляции само по себе еще не говорит о наличие или отсутствии автокорреляции. Его необходимо сравнить с критическим.
Фактическое
значение коэффициента автокорреляции
необходимо
сравнить с табличным (критическим).
Рассчитанное
в данном примере значение коэффициента автокорреляции
=0,67 сравнивается с табличным при 5-процентном уровне
значимости, для
=14 при а = 0,05 критическое значение коэффициента
автокорреляции равно 0,335. Так как рассчитанное значение
(0,67) больше табличного (0,335), то с вероятностью Р
=0,95 (Р =
) можно сделать вывод о наличии автокорреляции в
исследуемом ряду .
Часто приходится решать вопрос о наличии или отсутствии автокорреляции не между уровнями ряда, а между их отклонениями и от тренда или от среднего уровня, т. е. между так называемыми остаточными величинами. Сумма таких остаточных величии и средняя из них равны нулю.
Из
формулы (2.26) коэффициент автокорреляции, для рядов со средним уровнем равным
нулю
, равен: (2.27)
Формула
коэффициента корреляции для остаточных автокорреляции для остаточных величин
примет вид:
(2.28)
Кроме
показателя
, для обнаружения автокорреляции между соседними
остаточными величинами часто используется критерий, разработанный Дурбиным и
Ватсоном.
Критерии
Дурбина-Ватсона, обозначаемый как
(иногда DW),
рассчитывается по формуле
(2.29)
Этот
показатель можно связать с формулой (2.28) коэффициента автокорреляции для
остаточных величин. Так, если предположить, что
то,
возведя в квадрат числитель критерия
, можно
записать
(2.30)
Очевидно,
что если автокорреляция отсутствует, т. е.
= 0, то
= 2. Если же имеет место полная автокорреляция, т. е.
равен 1 или -1, то значение
будет соответственно 0 или 4.
Для
того чтобы иметь возможность более точно судить об отсутствии автокорреляции в
остаточных величинах
, составлена таблица, в которой для разного числа
наблюдений
и числа независимых переменных в уравнении регрессии
определены верхнее
и нижнее
критические границы критерия
.
Для
проверки нулевой гипотезы об устранении автокорреляции в остаточных величинах
рассчитанное фактическое значение
сравнивается
с табличными
и
:
)
если
(до 4-
),
гипотеза об отсутствии автокорреляции принимается;
)
если
, гипотеза об отсутствии автокорреляции отвергается;
)
если
или (4 -
)
, ничего определенного сказать нельзя и требуется
дальнейшее исследование;
)
если
, имеет место отрицательная автокорреляция.
Для
иллюстрации расчета критерия Дурбина-Ватсона
, а также
используются данные табл. 2.1.4.
На
основе фактических и выровненных уровней рассчитываются остаточные величины
их необходимо проверить на автокорреляцию. Расчеты
представлены в табл. 2.1.6.
При
аналитическом выравнивании рядов динамики остаточные величины проверяются на
автокорреляцию. Цель проверки - определить адекватность подобранной функции
(линии тренда), используемой для отражения тенденций развития в исследуемый
период. Если в остаточных величинах обнаруживается автокорреляция, это признак
неадекватности выбранного уравнения тренда.