Прогнозирование на основе сглаженного временного ряда

Задача 5.

Динамика удельного расхода условного топлива на производство теплоэнергии (, кг/Гкал) на ТЭЦ по годам представлена в таблице.

Требуется:

1) произвести сглаживание ряда методом трехлетней скользящей средней;

2) выровнять ряд по прямой – т.е. оценить параметры линейного тренда методом наименьших квадратов;

3) начертить графики первичного и сглаженных рядов;

4) на уровне значимости проверить согласованность линейной трендовой модели с результатами наблюдений;

5) методом экстраполяции найти точечные и интервальные (с доверительной вероятностью ) оценки прогноза экономического показателя на 2002 и 2003г.г.

	169,2	168,1	168,6	168,4	167,9	167,6	167,8	166,9	167,1

Временным рядом называется последовательность значений (уровней) некоторого экономического показателя , расположенных в порядке возрастания времени. Уровни ряда должны отражать значения экономического показателя за одинаковые или через одинаковые промежутки времени.

Одной из важнейших задач исследования временного ряда является задача выявления основной тенденции развития (тренда) изучаемого процесса.

Решение этой задачи необходимо для прогнозирования. При этом исходят из того, что тенденция развития, установленная в прошлом, может быть распространена (экстраполирована) на будущий период.

Наиболее простыми и часто применяемыми способами выявления основной тенденции развития являются сглаживание временного ряда методом скользящей средней или выравнивание по прямой методом наименьших квадратов.

1) Метод скользящей средней основан на переходе от начальных значений членов ряда к их средним значениям на интервале времени, длина которого определена заранее. При этом сам выбранный интервал времени "скользит" вдоль ряда, получаемый таким образом ряд скользящих средних ведет себя более гладко, чем исходный ряд.

Для нашего примера скользящие средние находим по формуле

, .

Например, при

,

при .

По результатам получим сглаженный ряд:

	-	168,6	168,4	168,3	168,0	167,8	167,8	167,3	-

2) По статистическим данным найдем оценки и параметров линейного тренда методом наименьших квадратов. Для этого применим известные формулы [1]:

, ,

где , , , ,

, .

Здесь и в дальнейшем – номер уровня ряда: 1993 г. соответствует номер 1, …2001 году – номер 9.

Вычисление средних значений организуем в форме расчетной таблицы.

	169,2		28628,64	169,2
	168,1		28257,61	336,2
	168,6		28425,96	505,8
	168,4		28358,56	673,6
	167,9		28190,41	839,5
	167,6		28089,76	1005,6
	167,8		28156,84	1174,6
	166,9		27855,61	1335,2
	167,1		27922,41	1503,9
		1511,6		253885,8	7543,6
		167,955	31,67	28209,53	838,18

; ;

; .

Таким образом искомые оценки параметров линейного тренда равны:

, . Уравнение линейного тренда имеет вид:

.

3) На рисунке цифрой (1) отмечен первичный ряд, цифрой (2) – скользящая трехлетняя средняя, цифрой (3) помечен ряд, выровненный по прямой.

4) Проверка согласованности линейной трендовой модели с результатами наблюдений выполняется как решение задачи проверки статистической гипотезы об отсутствии линейной статистической связи переменных и на заданном уровне значимости . Для проверки гипотезы используется коэффициент детерминации и применяется статистика Фишера с и степенями свободы.

В рассматриваемом случае , , .

Критическое значение статистики Фишера равно

.

Так как , то выдвинутая гипотеза отвергается, что свидетельствует о согласии линейной трендовой модели с результатами наблюдений.

5) По полученному уравнению линейного тренда найдем точечные (индивидуальные) прогнозы показателя на 2002 и 2003 г.г.

Для 2002г.

.

Для 2003г.

.

Дать интервальную оценку тренда – значит указать границы интервала, в который попадет возможное значение переменной с заданной доверительной вероятностью (в нашем примере ).

Этот интервал определяется по известным формулам [3]

,

где - точность прогноза ,

здесь - число степеней свободы, , ищется по таблице критических точек распределения Стьюдента для двусторонней критической области (см., например [4]); в нашем случае ; ;

. (Можно воспользоваться так же таблицами [3]).

- исправленное среднеквадратическое отклонение (С.К.О.) индивидуальных значений зависимой переменной

.

Из этой формулы видно, чем больше , тем меньше точность прогноза.

- исправленное С.К.О. ошибок линейной регрессии

.

Вычисление доверительных интервалов прогнозов организуем в виде таблицы

	169,2	168,9266	0,2734	0,07475
	168,1	168,6837	-0,5837	0,34071
	168,6	168,4408	0,1592	0,02534
	168,4	168,1979	0,2021	0,04084
	167,9	167,9550	-0,055	0,00303
	167,6	167,7121	-0,1121	0,01257
	167,8	167,4692	0,3308	0,10943
	166,9	167,2263	-0,3263	0,10647
	167,1	166,9834	0,1166	0,01360
	–	–	–	0,72674

.

Дальнейшие вычисления проводим отдельно для (2002 г.) и (2003 г.)

Для

.

,

.

Итак, с вероятностью , удельный расход условного топлива в 2002 г. будет принадлежать интервалу (кг/Гкал)

.

Аналогично для 2003 г. , получим

;

,

,

, .

Рекомендуемая литература

1. Ефимова М.Р., Петрова Е.В., Румянцев В.Н. Общая теория статистики: Учебник. – М.:ИНФРА-М, 2001. – 416 с.

2. Ефимова М.Р., Ганченко О.И., Петрова Е.В. Практикум по общей теории статистики: Учеб. пособие. – М.: Финансы и статистика, 2001. – 280 с.

3. Кремер Н.Ш., Путко Б.А. Эконометрика: Учебник для вузов / Под ред. проф. Н.Ш. Кремера. – М.:ЮНИТИ-ДАНА, 2002. – 311 с.

4. Гмурман В.Е. Теория вероятностей и математическая статистика: Учеб. пособие для вузов. – М.: Высш. школа, 1972. – 368 с.

5. Сборник задач по математике для вузов. Специальные курсы. / Под ред. А.В. Ефимова. – М.: Наука, 1984. – 608 с.

Подписано в печать 06.06.2005. Формат 60х84/16. Бумага для множ.аппаратов.

Печать плоская. Усл. печ. л. 1,0. Уч.-изд. л.1,1. Тираж 150 экз. Заказ № 333-р.

––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––-

Ризограф РГППУ. Екатеринбург, ул. Машиностроителей, 11

Дата добавления: 2014-12-16