Ольков_С_Г_Аналитическая юриспруденция
.pdfАвтокорреляция первого порядка (first-order autocorrelation) – оценивает степень зависимости между последовательными значениями временного ряда. Автокорреляция второго порядка
(second-order autocorrelation) – оценивает степень зависимости между значениями временного ряда, разделенными двумя временными интервалами. Автокорреляция p-го порядка (pth- order autocorrelation) – оценивает степень зависимости между значениями временного ряда, разделенными n временными интервалами.
Аналитические характеристики (показатели) временного ряда – абсолютный прирост (цепной и к базе), темпы роста (цепные и к базе), темпы прироста (цепные и к базе), коэффициент опережения, коэффициент вариации, коэффициент осцилляции, размах, стандартное отклонение и другие.
Барометрические методы прогнозирования – методы, учитывающие эффекты совпадения, опережения и запаздывания различных факторов, связанных с прогнозируемым.
Блочная диаграмма (box plot) – диаграмма, представляющая квартили и среднее (характеристики распределения).
Веменной ряд (time series) – это набор числовых данных, полученных в течение последовательных периодов времени.
Доверительный интервал тренда (confidence interval of the trend) – рассчитываемый по формуле диапазон, в пределах которого ожидается прогнозируемое по тренду значение переменной величины (постоянная величина).
Доверительный интервал прогноза (confidence interval of the forecast) – рассчитываемый по формуле диапазон, в пределах которого ожидается значение прогнозируемой переменной. Ширина доверительного интервала зависит от N – длины временного ряда (зависимость обратная – чем длиннее временной ряд, тем ýже доверительная трубка) и периода упреждения l (зависимость прямая – чем дальше горизонт прогнозирования, тем шире доверительный интервал).
202
Методы анализа временных рядов (time-series forecasting methods) – позволяют предсказывать значение численной переменной на основе её прошлых и настоящих значений.
Методы качественного прогнозирования (qualitative forecasting methods) – неколичественные методы прогнозирования в условиях более или менее высокой степени неопределенности, например, «мозговой штурм», «экспертные оценки», «интуитивные догадки» следует рассматривать как предварительные и поисковые, носят весьма субъективный характер.
Методы количественного прогнозирования (quantitative forecasting methods): 1) методы анализа временных рядов (timeseries forecasting methods); 2) методы анализа причинноследственных зависимостей (causal forecasting methods).
Принцип экономии (parsimony) – из нескольких равноправных моделей прогнозирования следует выбирать простейшую.
Прогноз (от греч. рrognōsis – предвидение, предсказание) – научно-обоснованное предсказание, суждение о состоянии какоголибо явления в будущем.
Прогнозирование (forecasting) – разработка прогноза, а «прогностика» – наука о законах и способах разработки прогнозов.
Прогнозирование по переменным, зависящим от сходных факторов – это прогнозирование основанное на предположении, что прогнозируемая и прогнозирующая переменные являются следствиями общего причинного комплекса.
Сезонность юридического процесса (seasonal component) – периодические колебания юридического процесса, вызванные сезонными (связанными с временем года), климатическими и иными процессами.
Скользящее среднее (moving average, also called rolling average, rolling mean or running average) – последовательность средних значений для выбранного временного периода, имеющего длину L, например, 3 года, 5 лет и т.п., вычисленных для последовательностей длины L. В англоязычной литературе сокращенно обозначают MA(L). Вычисляется по формуле:
203
MA(L) = y1 + y2 +... + yL . Используется для сглаживания временных
L
рядов и прогнозирования.
Средняя ошибка аппроксимации вычисляется по формуле:
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
= |
1 |
å |
yi |
|
- yi |
×100 или |
||||||
A |
|||||||||||||||
|
N |
|
y |
||||||||||||
|
|
|
|
|
|
|
i=1 |
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
2 |
|
|
|
|
|
|
|
|
|
100 |
|
|
å(yi - yi ) |
|
|
|
|||
|
А = |
× |
|
i=1 |
|
|
. |
||||||||
|
|
у |
|
|
N |
|
|||||||||
Применяется для оценки качества аппроксимирующих моделей. Аппроксимация считается приемлемой, если средняя ошибка аппроксимации не превышает 8%.
Тенденция (от лат. tendentia – направленность (англ. tendency)) – это отличное от стационарного (стабильного) течение какого-либо процесса.
Тренд (trend) – уравнение, выражающее тенденцию.
Фиктивная (искусственная или индикаторная) переменная (dummy variable) – нули и единицы, присваиваемые каким-либо качественным показателям (категориальным переменным), которые нужно включить в количественную математическую модель, в частности для проведения множественного регрессионного анализа.
F-статистика Фишера (F-statistic) – показатель качества регрессионного уравнения в целом (эмпирический коэффициент сравнивается с табличным).
Экспоненциальное сглаживание (exponential smoothing) – метод прогнозирования, получивший свое название от последовательности экспоненциально взвешенных скользящих средних. Позволяет делать, как долгосрочные, так и краткосрочные прогнозы. Вычисление экспоненциально сглаженного значения в i-
м |
периоде |
осуществляется |
по |
формуле: |
E1 =Y1 , |
Ei =WY i +(1 -W )Ei−1 , i = 2,3,4... , |
|
|
|
||
где Ei – значение экспоненциально сглаженного ряда, вычисленное для i-го периода, Ei-1 – значение экспоненциально сглаженного ряда, вычисленное для i-1 периода, Yi – наблюдаемое значение временного ряда в i-м периоде, W – субъективный вес,
204
или сглаживающий коэффициент, принимающий значение в пределах от нуля до единицы, не включая крайние значения.
ОСНОВНЫЕ ТЕРМИНЫ:
блочная диаграмма (box plot); влияние сезонности (seasonal component), влияние случайности (irregular, or random component), временной ряд (time series), горизонт прогнозирования (период упреждения), доверительный интервал прогноза (confidence interval of the forecast), доверительный интервал тренда (confidence interval of the trend), множественный регрессионный анализ (multiple regression), прогноз (forecast), тенденция (tendency), тренд (trend), сезонный индекс (seasonal index), сезонная поправка (seasonal adjustment), статистика Фишера (F-statistic), скользящая средняя (moving average), фиктивная (искусственная, индикаторная) переменная (dummy variable), циклическая составляющая (cyclical component), экспоненциальное сглаживание (exponential smoothing).
РЕШЕНИЕ ТИПОВЫХ ЗАДАЧ
Задача №1.
Дано: поквартальные первичные статистические данные о числе зарегистрированных грабежей в городе Сургуте с 1998 по 2004 годы.
Таблица №1. Временной поквартальный ряд грабежей в городе Сургуте с 1998 по 2004 годы.
|
t |
|
|
|
кварталы |
|
y |
t |
по |
t |
грабежи, |
годы |
порядку |
кварталы |
шт. |
|
|
|
|
|
|
|
|
1998 |
1 |
1 |
268 |
|
|
|
|
1998 |
2 |
2 |
164 |
|
|
|
|
1998 |
3 |
3 |
147 |
|
|
|
|
1998 |
4 |
4 |
412 |
|
|
|
|
205
1999 |
5 |
1 |
261 |
1999 |
6 |
2 |
122 |
1999 |
7 |
3 |
115 |
1999 |
8 |
4 |
185 |
2000 |
9 |
1 |
99 |
2000 |
10 |
2 |
76 |
2000 |
11 |
3 |
69 |
2000 |
12 |
4 |
125 |
2001 |
13 |
1 |
126 |
2001 |
14 |
2 |
128 |
2001 |
15 |
3 |
122 |
2001 |
16 |
4 |
178 |
2002 |
17 |
1 |
204 |
2002 |
18 |
2 |
188 |
2002 |
19 |
3 |
206 |
2002 |
20 |
4 |
278 |
2003 |
21 |
1 |
250 |
2003 |
22 |
2 |
207 |
2003 |
23 |
3 |
297 |
2003 |
24 |
4 |
312 |
2004 |
25 |
1 |
317 |
2004 |
26 |
2 |
253 |
2004 |
27 |
3 |
298 |
2004 |
28 |
4 |
427 |
Требуется найти: 1) скользящие средние грабежей; 2) частные индексы сезонности, как частное от деления числа грабежей на их скользящие средние; 3) средние сезонные индексы для первого, второго, третьего и четвертого кварталов; 4) построить график сезонной составляющей временного ряда грабежей; 5) число грабежей с поправкой на сезонность; 6) трендовую составляющую в виде линейного уравнения. Сделать: 7) прогноз грабежей в городе Сургуте на 2005-2006 годы (29-36 кварталы) с учетом трендовой и сезонной составляющих временного ряда.
Ответ:
Таблица №2. Расчет скользящих средних и частных
индексов сезонности.
206
|
|
w - |
Частные |
|
|
скользящее |
индексы |
t, |
y, грабежи, |
среднее, |
сезонности |
кварталы |
шт. |
шт. |
y/w |
|
|
|
|
|
|
|
|
1 |
268 |
нет |
нет |
|
|
|
|
2 |
164 |
нет |
нет |
|
|
|
|
3 |
147 |
246,88 |
0,595 |
|
|
|
|
4 |
412 |
240,75 |
1,711 |
|
|
|
|
1 |
261 |
231,50 |
1,127 |
|
|
|
|
2 |
122 |
199,13 |
0,613 |
|
|
|
|
3 |
115 |
150,50 |
0,764 |
|
|
|
|
4 |
185 |
124,50 |
1,486 |
|
|
|
|
1 |
99 |
113,00 |
0,876 |
|
|
|
|
2 |
76 |
99,75 |
0,762 |
|
|
|
|
3 |
69 |
95,63 |
0,722 |
|
|
|
|
4 |
125 |
105,50 |
1,185 |
|
|
|
|
1 |
126 |
118,63 |
1,062 |
|
|
|
|
2 |
128 |
131,88 |
0,971 |
|
|
|
|
3 |
122 |
148,25 |
0,823 |
|
|
|
|
4 |
178 |
165,50 |
1,076 |
|
|
|
|
1 |
204 |
183,50 |
1,112 |
|
|
|
|
2 |
188 |
206,50 |
0,910 |
|
|
|
|
3 |
206 |
224,75 |
0,917 |
|
|
|
|
4 |
278 |
232,88 |
1,194 |
|
|
|
|
1 |
250 |
246,63 |
1,014 |
|
|
|
|
2 |
207 |
262,25 |
0,789 |
|
|
|
|
3 |
297 |
274,88 |
1,080 |
|
|
|
|
4 |
312 |
289,00 |
1,080 |
|
|
|
|
1 |
317 |
294,88 |
1,075 |
|
|
|
|
2 |
253 |
309,38 |
0,818 |
|
|
|
|
3 |
298 |
нет |
нет |
|
|
|
|
4 |
427 |
нет |
нет |
|
|
|
|
Расчет скользящих средних производится следующим
образом: |
w1 |
= |
(0,5× 268 +164 +147 + 412 + 261 ×0,5) |
= 246,88 |
||
|
||||||
|
|
|
4 |
|
|
|
w2 = |
(0,5×164 +147 + 412 + 261 +122 ×0,5) |
= 240,75 |
и т.д. |
|||
|
|
4 |
||||
|
|
|
|
|
|
|
Таблица №3. Расчет средних сезонных индексов для
первого, второго, третьего и четвертого кварталов.
207
|
|
|
|
|
|
s, |
Частные |
|
|
|
|
|
средний |
индексы |
Кварталы |
1кв. |
2 кв. |
3кв. |
4 кв. |
сезонный |
сезонности |
|
|
|
|
|
индекс по |
y/w |
|
|
|
|
|
кварталам |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Нет |
1 |
нет |
|
|
|
1,044 |
|
|
|
|
|
|
|
Нет |
2 |
|
нет |
|
|
0,811 |
|
|
|
|
|
|
|
0,595 |
3 |
|
|
0,595 |
|
0,817 |
|
|
|
|
|
|
|
1,711 |
4 |
|
|
|
1,711 |
1,289 |
|
|
|
|
|
|
|
1,127 |
1 |
1,127 |
|
|
|
1,044 |
|
|
|
|
|
|
|
0,613 |
2 |
|
0,613 |
|
|
0,811 |
|
|
|
|
|
|
|
0,764 |
3 |
|
|
0,764 |
|
0,817 |
|
|
|
|
|
|
|
1,486 |
4 |
|
|
|
1,486 |
1,289 |
|
|
|
|
|
|
|
0,876 |
1 |
0,876 |
|
|
|
1,044 |
|
|
|
|
|
|
|
0,762 |
2 |
|
0,762 |
|
|
0,811 |
|
|
|
|
|
|
|
0,722 |
3 |
|
|
0,722 |
|
0,817 |
|
|
|
|
|
|
|
1,185 |
4 |
|
|
|
1,185 |
1,289 |
|
|
|
|
|
|
|
1,062 |
1 |
1,062 |
|
|
|
1,044 |
|
|
|
|
|
|
|
0,971 |
2 |
|
0,971 |
|
|
0,811 |
|
|
|
|
|
|
|
0,823 |
3 |
|
|
0,823 |
|
0,817 |
|
|
|
|
|
|
|
1,076 |
4 |
|
|
|
1,076 |
1,289 |
|
|
|
|
|
|
|
1,112 |
1 |
1,112 |
|
|
|
1,044 |
|
|
|
|
|
|
|
0,910 |
2 |
|
0,91 |
|
|
0,811 |
|
|
|
|
|
|
|
0,917 |
3 |
|
|
0,917 |
|
0,817 |
|
|
|
|
|
|
|
1,194 |
4 |
|
|
|
1,194 |
1,289 |
|
|
|
|
|
|
|
1,014 |
1 |
1,014 |
|
|
|
1,044 |
|
|
|
|
|
|
|
0,789 |
2 |
|
0,789 |
|
|
0,811 |
|
|
|
|
|
|
|
1,080 |
3 |
|
|
1,08 |
|
0,817 |
|
|
|
|
|
|
|
1,080 |
4 |
|
|
|
1,08 |
1,289 |
|
|
|
|
|
|
|
1,075 |
1 |
1,075 |
|
|
|
1,044 |
|
|
|
|
|
|
|
0,818 |
2 |
|
0,818 |
|
|
0,811 |
|
|
|
|
|
|
|
Нет |
3 |
|
|
нет |
|
0,817 |
|
|
|
|
|
|
|
Нет |
4 |
|
|
|
нет |
1,289 |
|
|
|
|
|
|
|
|
Среднее |
1,044 |
0,811 |
0,817 |
1,289 |
|
|
|
|
|
|
|
|
208
Рис. №1. Сезонная составляющая временного ряда грабежей.
Таблица №4. Рассчет прогнозных значений числа
грабежей в г.Сургуте с учетом трендовой и сезонной составляющих.
|
|
|
|
|
|
|
Грабежи с |
|
s, |
|
|
|
поправкой |
|
средний |
|
С |
t, |
на |
|
сезонный |
|
поправко |
квартал |
сезонность: |
z=q(t)= |
индекс по |
t, кварталы |
й на |
ы по |
|
5,232∙t+138,8 |
квартала |
по порядку |
сезонност |
порядку |
q=y∙s |
|
м |
|
ь Yоц=z∙s |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
279,79 |
144,032 |
1,044 |
1 |
150,37 |
|
|
|
|
|
|
2 |
133,00 |
149,264 |
0,811 |
2 |
121,05 |
|
|
|
|
|
|
3 |
120,10 |
154,496 |
0,817 |
3 |
126,22 |
|
|
|
|
|
|
4 |
531,07 |
159,728 |
1,289 |
4 |
205,89 |
|
|
|
|
|
|
5 |
272,48 |
164,96 |
1,044 |
5 |
172,22 |
|
|
|
|
|
|
6 |
98,94 |
170,192 |
0,811 |
6 |
138,03 |
|
|
|
|
|
|
7 |
93,96 |
175,424 |
0,817 |
7 |
143,32 |
|
|
|
|
|
|
8 |
238,47 |
180,656 |
1,289 |
8 |
232,87 |
|
|
|
|
|
|
9 |
103,36 |
185,888 |
1,044 |
9 |
194,07 |
|
|
|
|
|
|
10 |
61,64 |
191,12 |
0,811 |
10 |
155,00 |
|
|
|
|
|
|
11 |
56,37 |
196,352 |
0,817 |
11 |
160,42 |
|
|
|
|
|
|
12 |
161,13 |
201,584 |
1,289 |
12 |
259,84 |
|
|
|
|
|
|
13 |
131,54 |
206,816 |
1,044 |
13 |
215,92 |
|
|
|
|
|
|
14 |
103,81 |
212,048 |
0,811 |
14 |
171,97 |
|
|
|
|
|
|
15 |
99,67 |
217,28 |
0,817 |
15 |
177,52 |
|
|
|
|
|
|
16 |
229,44 |
222,512 |
1,289 |
16 |
286,82 |
|
|
|
|
|
|
209
17 |
212,98 |
227,744 |
1,044 |
17 |
237,76 |
18 |
152,47 |
232,976 |
0,811 |
18 |
188,94 |
19 |
168,30 |
238,208 |
0,817 |
19 |
194,62 |
20 |
358,34 |
243,44 |
1,289 |
20 |
313,79 |
21 |
261,00 |
248,672 |
1,044 |
21 |
259,61 |
22 |
167,88 |
253,904 |
0,811 |
22 |
205,92 |
23 |
242,65 |
259,136 |
0,817 |
23 |
211,71 |
24 |
402,17 |
264,368 |
1,289 |
24 |
340,77 |
25 |
330,95 |
269,6 |
1,044 |
25 |
281,46 |
26 |
205,18 |
274,832 |
0,811 |
26 |
222,89 |
27 |
243,47 |
280,064 |
0,817 |
27 |
228,81 |
28 |
550,40 |
285,296 |
1,289 |
28 |
367,75 |
29 |
|
290,528 |
1,044 |
29 |
303,31 |
30 |
|
295,76 |
0,811 |
30 |
239,86 |
31 |
|
300,992 |
0,817 |
31 |
245,91 |
32 |
|
306,224 |
1,289 |
32 |
394,72 |
33 |
|
311,456 |
1,044 |
33 |
325,16 |
34 |
|
316,688 |
0,811 |
34 |
256,83 |
35 |
|
321,92 |
0,817 |
35 |
263,01 |
36 |
|
327,152 |
1,289 |
36 |
421,70 |
Число грабежей с поправкой на сезонность (q) находится как произведение значений столбца грабежей (у) на столбец средних сезонных индексов (s) попарно (y1∙s1; y2∙s2 и т.д.): q=y∙s.
Трендовая составляющая временного ряда (z) получается методом наименьших квадратов, как линейное уравнение, где независимой переменной выступает время (t) - в данном случае берутся кварталы по порядку: 1, 2, 3, 4, 5, 6, 7…28, а зависимой – число грабежей с поправкой на сезонность (q): z=q(t).
В нашем случае z=q(t)=5,232∙t+138,8. Откуда получаем трендовые значения, включая прогнозные, например, пятое: 5,232∙5+138=164,96.
210
Рис. №2. Уравнение тренда: z=5,232∙t+138,8
Чтобы получить окончательный прогноз с учетом трендовой и сезонной составляющей, попарно умножаем трендовые значения (z) на средние сезонные индексы (s): z∙s, а далее строим соответствующий график, включающий прогнозные значения полученные с учетом трендовой и сезонной составляющих временного ряда грабежей.
Рис. №3. Прогноз грабежей в г.Сургуте на 2005-2006 годы (29-36 кварталы), учитывающий тенденцию: z=5,232∙t+138,8, а также сезонные колебания.
211
