Восстановление зависимости методом наименьших квадратов на основе непараметрической модели с периодической составляющей Текст научной статьи по специальности «Математика»

УДК 519.2:303.732.4

UDC 519.2:303.732.4

ВОССТАНОВЛЕНИЕ ЗАВИСИМОСТИ МЕТОДОМ НАИМЕНЬШИХ КВАДРАТОВ НА ОСНОВЕ НЕПАРАМЕТРИЧЕСКОЙ МОДЕЛИ С ПЕРИОДИЧЕСКОЙ СОСТАВЛЯЮЩЕЙ

Орлов Александр Иванович д.э.н., д.т.н., к.ф.-м.н., профессор

Московский государственный технический университет нм. Н.Э. Баумана, Россия, 105005, Москва, 2-я Бауманская ул., 5, рго{-ог!оу(сйта)I.ги

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

Ключевые слова: МЕТОД НАИМЕНЬШИХ КВАДРАТОВ, НЕПАРАМЕТРИЧЕСКИЕ МЕТОДЫ, ПЕРИОДИЧЕСКАЯ СОСТАВЛЯЮЩАЯ, ОЦЕНИВАНИЕ, ПРОГНОЗИРОВАНИЕ

RENEWAL DEPENDENCE METHOD OF LEAST SQUARES BASED NONPARAMETRIC MODEL WITH PERIODIC COMPONENT

Orlov Alexander Ivanovich

Dr.Sci.Econ., Dr.Sci.Tech., Cand.Phys-Math.Sci.,

professor

Bauman Moscow State Technical University, Moscow, Russia

We consider the nonparametric problem of reneval dependence, which is described by the sum of a linear trend and periodic function with a known period. We obtain the asymptotic distribution of the parameter estimates and the trend component. The methods of estimating the periodic component and designing interval forecast. In the model of the points of observation, natural for applications, justified by the conditions of use. In particular, we prove an asymptotically unbiased estimate of the coefficient of the linear term

Keywords: METHOD OF LEAST SQUARES, NONPARAMETRIC METHODS, PERIODIC COMPONENT, ESTIMATIOM, FORECASTING

Метод наименьших квадратов восстановления зависимости - один из наиболее распространенных статистических методов анализа данных. В статье рассмотрена непараметрическая постановка: восстанавливаемая зависимость - сумма линейной функции и периодической составляющей произвольного вида (с известным периодом), распределение случайных погрешностей (остатков, невязок) произвольно.

1. Задача восстановления линейной зависимости

Начнем с простейшего случая - задачи восстановления линейной зависимости. Пусть / - независимая переменная, ах- зависимая. Рассмотрим задачу восстановления зависимости х = х({) на основе набора п пар чисел (4, хк), к = 1,2,...,я, где t]i — значения независимой переменной, а Х/( - соответствующие им значения зависимой переменной.

Восстанавливать зависимость можно на основе различных моделей. Обычно

применяют модели временных рядов, включающие три составляющие: трендовую (7), периодическую (5) и случайную (/;’). Рассматривают, как в [1] и аналогичных изданиях, аддитивную модель Т + £ + Е и мультипликативную модель Тх 8 х Е. Простейшая аддитивная модель имеет вид

хк=а(гк- г)+ (1+ ек= а^к - 1) + й?+Х4) + Ек, к = 1,2,...,п. (1)

Здесь трендовая составляющая - линейная функция а (4 - г) + <3 (такая запись тренда предпочтительнее для облегчения выкладок); периодическая составляющая/(7) обычно описывает сезонность, т.е период известен (в зависимости от моделируемой ситуации он равен году, неделе, суткам и т.п.); случайная составляющая представлена слагаемыми Ек, которые являются реализациями независимых одинаково

распределенных случайных величин с нулевым математическим ожиданием и

2

дисперсией ° , неизвестной статистику. В рассматриваемой модели ек /(1к) + Ек, = 1,2,... ,/7, т.е. отклонения от линейного тренда ек не являются одинаково распределенными. Однако их распределения отличаются лишь сдвигами (на значения детерминированной периодической составляющей).

Соответствующая модели (1) мультипликативная модель имеет вид

У к =[^Лх/1(Ох[1 + 8Д к = \,2,...,п ф

В модели (2) сомножители имеют описанный выше смысл. При логарифмировании модель (2) переходит в аналог модели (1), следовательно, достаточно рассматривать модель (1).

Иногда принимают предположение о нормальности распределения погрешностей. Однако давно известно, что распределения реальных данных, как правило, отличаются от нормальных [2]. Поэтому далее рассматриваем непараметрическую модель, не предполагающую, что распределение погрешностей входит в то или иное параметрическое семейство. Отказ от задания распределения погрешностей в параметрическом виде - одно из оснований для того, чтобы именовать рассматриваемые модель и метод непараметрическими. Второе основание - отказ от выбора периодической составляющей из какого-либо параметрического

Научный журнал КубГАУ, №91(07), 2013 года семейства функций.

Практическая значимость модели (1) очевидна. Однако расчетные методы, описанные в [1] и аналогичных изданиях, являются эвристическими. Цель настоящей статьи - построить непараметрическую вероятностно-статистическую теорию прогноза временного ряда на базе линейного тренда с учетом аддитивной периодической составляющей.

Метод наименьших квадратов был разработан К. Гауссом в 1794 г. [2]. Согласно нему этому методу для расчета наилучшей функции, приближающей линейным образом зависимость х от ? в модели (1), следует рассмотреть функцию двух переменных

п _

/(а, </) = £ (хк - а{1к - г) - й)2

к=\

Оценки метода наименьших квадратов (кратко: оценки МНК) - это такие значения а* и й?*, при которых функция 1(а,с1) достигает минимума по всем значениям аргументов. Как известно (см., например, [2]), оценки МНК имеют вид

а* = х = -^~

(3)

а* _ и-----------------------^ = X = — у хк

Следуя эвристическому подходу [1], изучим асимптотическое поведение оценок МНК а* ий?*, заданных формулами (3), установим их асимптотическую нормальность в предположениях модели (19), а затем состоятельно оценим периодическую составляющую ДО и построим интервальный прогноз для х(1).

2. Асимптотические распределения оценок параметров

Из формулы (3) следует, что

п к=1 П к=1 п к=1 п к=1 п к=1 _ (4)

Согласно Центральной предельной теореме (для выполнения ее условий необходимо предположить, например, что погрешности ек, к = 1, 2, ...,«, финитны

или имеют конечный третий абсолютный момент; однако заострять внимание на этих внутриматематических «условиях регулярности» здесь нет необходимости) оценка й?* имеет асимптотически нормальное распределение с математическим ожиданием

<*+-£/('.) <т2/

п к=1 и дисперсией /и , ее оценка приводится ниже. Из формул (3) и (4)

вытекает, что

2 п 1 п

— х = а(ґк-ї) + сї + ек-сі---------^ек =Фк ~0 + ек--------------

п к=і п к=1

(■хк - х)(ь - о=а(ь - о2+ек (*к - о - ^ Е

п к=1

Последнее слагаемое во втором соотношении при суммировании по к обращается в 0, поэтому

(^-0

а* = а + ^скек = а + 1>ь/ХО + '

к=1 £=1 £=1

(5)

£=1

Формулы (5) показывают, что оценка а* является асимптотически нормальной с

математическим ожиданием ; и дисперсией

В(а*) = Х40(£,) =

о2

к=\

Отметим, что многомерная нормальность имеет быть, когда каждое слагаемое в

формуле (5) мало сравнительно со всей суммой, т.е.

Нш шах \ ік-і\ /{^(ік - і)2}1'2 = 0

к=1

(6)

Условие (6) выполнено, например, если 4 образуют (полную, т.е. без пропусков) арифметическую прогрессию, число членов которой безгранично растет.

Итак, дисперсии оценок МНК параметров а* и с/* линейного тренда - те же, что и при отсутствии сезонных искажений (см., например, [2]). А вот их математические ожидания зависят от периодической составляющей. Однако в случае

£ло=£(',-о/(<,)=о

г=1 ;=1 (7)

оценки а* и й* являются несмещенными.

Условия (7) являются необходимыми и достаточными для несмещенности и состоятельности оценок МНК коэффициентов линейной зависимости. Проверка условий (7) рассмотрена в конце статьи.

Несмещенность (в предположениях (7) и асимптотическая нормальность оценок метода наименьших квадратов позволяют легко указывать для них асимптотические доверительные границы и проверять статистические гипотезы, например, о равенстве определенным значениям, прежде всего 0.

3. Асимптотическое распределение трендовой составляющей

Из формул (4) и (5) следует, что при справедливости соотношений (7)

М{а *(t-t) + d*} = M(a*)(t -t) + M(d*) = a(t - t) + d, т.е. оценка y*(t) = a* (tk - f)+d* трендовой составляющей y(t) = a (t

рассматриваемой зависимости является несмещенной. Поэтому

)+ d

D(y * (?)) = D(a*)(t -1)2 + 2M{(a * -a)(d * -d)(t - ?)} + D(d*).

При этом, поскольку погрешности Ek независимы в совокупности иМ(Ек) = 0, то

М{(а * -a){d * -d)(t -t)} = -^ck(t- t)M(Ek) = -(t - t)ack = 0 .

n k=1 n k=1

Таким образом,

2

1,(1,-if

k=1

(8)

Итак, оценка у*(^) является несмещенной и асимптотически нормальной. Для ее практического использования (построения доверительных интервалов, проверки статистических гипотез) необходимо состоятельно уметь оценивать остаточную

дисперсию ~а .

В частности, не представляет труда выписывание нижней и верхней границ для

Научный журнал КубГАУ, №91(07), 2013 года трендовой составляющей прогностической функции:

УниЖн(1) = а*^-1) + ё*-5(0, УеерЛ*) = а*^-1) + с!*+8(*>

где полуширина доверительного интервала ^ ^ имеет вид

<5(0 = и(Г)^О*(у*(0) = Щур *

1

п

(9)

Здесь ^ - доверительная вероятность, ) - квантиль нормального

Цу ^(7) = ф-|1±11

распределения порядка 2 , т.е. ^ 2 ' , где - функция стандартного

нормального распределения с математическим ожиданием 0 и дисперсией 1. При У =

0,95 (наиболее применяемое значение) имеем = 1,96. В формуле (9) £)*(>’*(0) . состоятельная оценка дисперсии >'*(/). В соответствии с (8) она является произведением состоятельной оценки ° * среднего квадратического отклонения ° случайных погрешностей Е* на известную исследователю детерминированную функцию от /.

4. Математическое ожидание остаточной суммы квадратов

В точках 4, к = 1, 2, ..., п, имеются исходные значения зависимой переменной и восстановленные значения .у *(4)- Рассмотрим остаточную сумму квадратов

® = Ё (У • (!,)-*,У =Ё !(<>* -а)((, -<) + (</* -<0 - /((, )-Е,)г.

к-1 £=1

При отсутствии периодической составляющей используют [2] состоятельные оценки ° * среднего квадратического отклонения ° случайных погрешностей,

№ * Ж

а* = „— о* =

построенные на основе остаточной суммы квадратов » 11 или \п-2 Однако при наличии периодической составляющей так делать нельзя. Приходится использовать «обходный путь».

В соответствии с формулами (4) и (5) при справедливости условий (7)

п п

= Щ

і=і Ь=і

У=1

сД^-0+-

п

і=1

=1ж(.

£=1

Найдем математическое ожидание каждого из слагаемых:

М(Ж»)=М|Х

е,-Яч)-еЛ ==м|1

Е\ -

-2М|Х

О'; 'І' '

п

%\(Лік) + Ек)+М(Лік)-Ек)

Поскольку Ек независимы, одинаково распределены и имеют нулевое математическое ожидание, то

м5

с,(ч-Ь+-

п

1

Е-\ =£г'('*-?)+н

Далее,

-гм|2

‘-.І'

п

Еі (ДО + Ек) = -2\ск(ік -1) +

Наконец,

М(/^к) -Ек)2 = /2(?к) +<з

На основе трех последних равенств можно показать, что при выполнении условия асимптотической нормальности (6)

1ш1М(Щ = /2(4)+а2

Следовательно,

м\^]=о2+1-±/\Ч) п ) п к=1

(10)

В правой части (10) первое слагаемое соответствует вкладу случайной составляющей, второе - вкладу периодической составляющей.

В некоторых случаях второе слагаемое в правой части (10) может быть известно

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

5. Оценивание периодической составляющей

В литературе рассматривают как параметрические, так и непараметрические подходы. Популярный метод исходит из того, что достаточно гладкую функцию можно разложить в ряд Фурье и получить хорошее приближение с помощью небольшого числа гармоник. В простейшем случае - одна гармоника. Так, динамику индекса инфляции можно попытаться изучать с помощью модели

xk=a(tk- *) + d+f(tk) + Ek = a (tk - *) + d+gco^2nt^ +Ek,k= 1,2,... ,n (время t измеряется в годах). Тогда неизвестные параметры a, b, g оцениваются методом наименьших квадратов.

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

Пусть в согласии с предположениями (7) рассматривается целое число периодов, т.е. п = mq, где п - объем наблюдений, т - количество периодов, q - число наблюдений в одном периоде. Предполагается, что первые q моментов наблюдения при сдвиге на длину периода дают следующие q моментов времени, при сдвиге на две длины периода дают третий набор из q моментов наблюдения, и т.д. Тогда в соответствии с определением периодической составляющей справедливы равенства

ДО = /(W = /(*29+,) = ••• = /(*(«- 1)«+,X s =

Если наблюдения проводятся ежемесячно в течение т лет, то число наблюдений в одном периоде <7=12, общий объем наблюдений п = 12/77, далее s - номер месяца в году, s = 1, 2, ..., 12. Пусть gs - общее значение в (11). Для оценки периодической составляющей требуется оценить gb g2, ..., gq.

Естественный подход состоит в том. чтобы усреднить т значений Х/( -соответствующих моментам времени, отстоящим друг от друга на целое число периодов. Другими словами, усреднить «очищенные» от трендовой составляющей

исходные данные, соответствующие одноименным месяцам различных лет. Речь идет об оценках

^ т

£* = —5 = 1,2

т ]=1

(12)

Оценка периодической составляющей распространяется на весь интервал наблюдений очевидным образом:

/ * (О = / * (V.) = f * (*2«+*) = - = /* (^(«-1)9+*) = £*’ 5 = ? (13)

Сложив восстановленные значения трендовой и периодической оставляющей, получим оценку зависимости, «очищенную» от случайной составляющей

**(0 = 3'*(0 + /*(0 = - 0+ <**+/*(0 (14)

Здесь оценки а* и й?* находят по формулам (3), а оценки /*(?) - по формулам (12) -(13).

С помощью формулы (14) можно строить точечный прогноз, используя ее вне интервала наблюдений. Для этого достаточно распространить сезонную составляющую (?) вплоть до рассматриваемого момента времени по правилу (13) и суммировать ее с прогнозом трендовой составляющей у*({). Интерполяция и экстраполяция на моменты времени /, не входящие в исходное множество {4, к= 1,2, ..., п} и множества, полученные из него сдвигами на целое число периодов, может быть осуществлена путем линейной интерполяции ближайших значений или иным методом сглаживания.

Обсудим свойства оценок (12) - (14).

При безграничном росте объема данных и справедливости условий (6) и (7) оценки а* и (I* параметров трендовой составляющей являются состоятельными и несмещенными, а потому, как можно показать, в рассматриваемых в настоящей статье условиях суммы (12) оценивают периодическую составляющую состоятельно (при т 00) и несмещенно. Как следствие,

-±и*т2--±/2ы^о П к=1 П к=1 (15)

по вероятности при п х . В соответствии с (10) последнее соотношение дает

„2

возможность оценить ° , а затем построить интервальный прогноз для трендовой составляющей согласно (9).

Отметим, что в рассматриваемой ситуации, как правило, п растет, увеличиваясь на величины, кратные ц - числу наблюдений в одном периоде. Как следствие, уменьшаемое в (15) - константа, зависимости от п нет. Эти особенности связаны с тем, что выполнение условий (7) предполагает рассмотрение целого числа периодов. Рассмотрим оценки (12) подробнее. Как вытекает из (4.1.19), (11) и (12),

1 т 1 т

е*=до - (а * -«)—Е(о(,-1), - о - (<** +—2Ж+0--1),»5=

т~ц т~1

С учетом (4), (5) и (7) получаем, что

£=до- --ИЕк+-ИЕ^и-1)9, *=1,2,-,?

4^=1 у\т -/=1 у к=1 -/=1

Таким образом,

&=до+£л*а> ^=\,2,-,ч

к=1 (16)

7 1 1 7 1

К = -СЛ-- + — 7 Г , / • 14 1 К=-СЛ~-

где и т , если + 0 1)^7 1,2, и и при всех остальных

1

г.

значениях индекса суммирования к? и

да

У=1

Соотношение (16) означает, что рассматриваемые оценки есть суммы независимых случайных величин, а потому с помощью Центральной предельной теоремы можно построить доверительные интервалы для рассматриваемых значений периодической составляющей (в предположении справедливости условий (6)).

6. Интервальный прогноз

Точечный прогноз строят по формуле (11) на основе х*(/) - оценки зависимости, «очищенной» от случайной составляющей, но включающей трендовый и периодический компоненты. Если выполнены условия (7), то

Мх;*(0 = х(0 = - О + й? +Х0>

т.е. оценка х*(Т) является несмещенной.

При справедливости условий (7) с учетом (4), (5) и (16) получаем, что для момента времени входящего в исходное множество {4, к = 1, 2, ..., п) или в множества, полученные из него сдвигами на целое число периодов,

х * (о - х(о=(г - оХ скЕк+- ЕЕк+Е кьЕк

к=1 п к=1 к=\ _ (17)

В (17) при определении значений коэффициентов /?/„ в качестве л- следует взять

номер наименьшего из исходных моментов времени {4, к = 1,2, ...,«}, отстоящих от

рассматриваемого момента ? на целое число периодов. С помощью (16) заключаем,

что

х*(о-х(о=2уьЕк к=1

= с +

где Ь т ^ если ^ + О-= 1>2,...,«?} ^ и ~ск(*~* “О ПрИ всех

остальных значениях индекса суммирования &, и - то же, что и в формуле (16).

В правой части формулы (17) стоит сумма независимых случайных величин, поэтому оценка х*(() является асимптотически нормальной (при справедливости условий (6)) с математическим ожиданием х(/) и дисперсией

£>(*(0) = £"£/)(£*)= а 2£ і

(18)

к-1 к-1

у / А Д* А

Следовательно, нижняя !"<жп ' и верхняя ; доверительные границы для прогностической функции (с учетом как трендовой, так и периодической составляющих) имеют вид:

(0 = а * (I -1) + а * +/ * (0 - Д(г), Хверх (0 = а * (Г -1) + (I * +/ * (0 + Д(г)

где

Д(0 = Е/(у)л/Я*(**(0)=£'(У)а

(19)

1 £=1

доверительная вероятность, и' - квантиль нормального

1 + у

распределения порядка 2 . В формуле (19) ^*(х*(0) _ состоятельная оценка дисперсии точечного прогноза х*(г). В соответствии с (18) она является произведением состоятельной оценки ° * среднего квадратического отклонения ° случайных погрешностей Е* на известную статистику детерминированную функцию

от Величину ° * рассчитывают согласно (10) и (15).

7. Пример применения непараметрического метода наименьших квадратов в модели с периодической составляющей

Обработаем фактические данные ОАО «Магнитогорский металлургический комбинат» о закупочных ценах на лом черных металлов [3]. Как показано в [3], может быть использована модель (1) линейного тренда с периодической составляющей. Для облегчения расчетов оставим из каждого квартала данные только по одному месяцу. Введем условные моменты времени, а именно, будем измерять время в кварталах, начиная с первого квартала 2003 г. Исходные данные для демонстрации примера применения непараметрического метода наименьших квадратов в модели с периодической составляющей - пары чисел (4, хк), к = 1, 2, ..., 12, - представлены в табл.1 в столбцах (3) и (4) соответственно.

По формулам (3) найдем оценки параметров а* и^, что позволяет построить оценку трендовой составляющей

у*( 0 = а* (г -*) + с1* =212,26 (г - 6,5) + 3967,17 = 212,26 г + 2587,48. Численные значения трендовой составляющей приведены в столбце (5) табл.1. Рассчитав отклонения исходных значений закупочных цен от оценок трендовой составляющей (столбец (6) табл.1), возведя их в квадрат и сложив, получаем остаточную сумму квадратов Ж? = 4 539 214 и XV п = 88112 = 378 267,843.

Таблица 1

Построение модели прогнозирования цен на лом марки ЗА

№ п/п Периоды Условны Закупоч- Оценка Отклоне Восста-н Кажу-щи

времени е ные тренда ния от овлен-н еся

моменты времени цены, руб./т оценки тренда ые зна-чени я невяз-ки

к h хк Хк-у\ч) * ** хк-хк

(1) (2) (3) (4) (5) (6) (V) (8)

1 янв.ОЗ 1 2 750 2 800 -50 2 424 326

2 апр.ОЗ 2 3 800 3 012 788 3 545 255

3 июл.ОЗ 3 2 900 3 224 -324 2 655 245

4 окт.03 4 3 100 3 437 -337 3 848 -748

5 янв.04 5 2 761 3 649 -888 3 273 -512

6 апр.04 6 4 602 3 861 741 4394 208

7 июл.04 7 3 540 4 073 -533 3504 36

8 окт. 04 8 5 268 4 286 982 4 697 571

9 янв.05 9 4 307 4 498 - 191 4 122 185

10 апр.05 10 4 779 4 710 69 5 243 -464

11 июл.05 11 4 071 4 922 -851 4 353 -280

12 окт.05 12 5 723 5 135 588 5546 177

Сгруппировав отклонения исходных значений закупочных цен от оценок трендовой составляющей по месяцам (табл.2), наглядно убеждаемся в наличии периодической составляющей. Взяв среднее арифметическое отклонений от тренда за

конкретный месяц, рассчитываем оценку ^ периодической составляющей (в соответствии с формулой (12)). Результаты приведены в табл.2.

Рассчитав по формуле (13) оценки периодической составляющей на весь интервал времени и сложив их с оценками трендовой составляющей, получаем в соответствии с формулой (14) оценку зависимости, «очищенную» от случайной составляющей, т.е. восстановленные значения (столбец (7) табл.1). Кажущиеся невязки, т.е. отклонения исходных значений закупочных цен от восстановленных значений, приведены в столбце (8) табл.1. Сравнивая столбцы (6) и (8), убеждаемся в

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

Таблица 2

Оценивание периодической составляющей

Номер квартала я Месяц Отклонения от тренда Оценка £ = /*(0 периодической составляющей

В 2003 г. В 2004 г. В 2005 г.

1 Январь -50 - 888 -191 -376

2 Апрель 788 741 69 533

3 Июль -324 - 533 -851 -569

4 Октябрь -337 982 588 411

Возведя в квадрат оценки периодической составляющей (табл.2), сложив эти

1Ё(/*й))2

квадраты, умножив на число лет и поделив на п, получаем, что п к=1 = 229 537.

В соответствии с формулой (10) оценкой дисперсии случайной составляющей является

/ *42 . \2

(?*) =-------ЕС/ ('*))

п =378 267,83 -229 537 = 148 731,

а оценкой среднего квадратического отклонения

а* = д/148731 =385,7

В соответствии с формулами (4) и (5) оценим дисперсии оценок параметров

. *\2

” (п * г 1487^1

Д ■ («*) = >(£,)= ,’ =—— = 1040.

143

2> л2

к=\

(а*)2 148731

в *= = = 12394.

п 12

Средние квадратические отклонения а* и с1* оцениваются как 32,25 и 111,33

соответственно, а доверительные интервалы для доверительной вероятности 0,95 таковы: КЛх] = [149,05; 275,47], К1П;</_] = [3748,96; 4185,38]

Первое из условий (7) выполнено в силу построения оценок периодической составляющей по целому числу периодов. Действительно, согласно данным табл.2 сумма оценок периодической составляющей для 12 точек наблюдений равна (-3), незначительное отклонение от 0 вызвано ошибками округления.

В соответствии с формулой (5) смещение оценки а* оценивается как

!>»/•(/,) = ^-------------= = 38>94

1.(1,-О2 143

к = 1

Таким образом, смещение имеет тот же порядок, что и среднее квадратичное отклонение оценки а*, и заведомо меньше, чем полуширина доверительного интервала. Дальнейшее сравнение может быть проведено на основе оценки дисперсии смещения - случайной величины

-<)/*(<*)

1 =

к=1

На,-*)1

к=1

Алгоритм вычисления дисперсии Ъ аналогичен таковым для периодической составляющей и интервального прогноза (см. (16) и (18) соответственно), но более сложен, поэтому не включен в статью. Таким образом, можно считать, что предположения (7) модели (1) выполнены для данных табл.1.

Перейдем к оценке дисперсий значений периодической составляющей. Как следует из равенства (16),

г = 1,2,.„л

к=1

; 1 1 ; 1

К- = ~сл—+— 1 ( , < • 1\ • 1 о 1 К = ~сл —

где п , если { (-7 ’’-,т}, и п при иных значениях

^ т

~ 2(^+(у-1)д О

г..

индекса суммирования к, и

да

У=1

Начнем со значения ^ = 1 (периодическая составляющая для января). Тогда

г, =Ц(\- 6,5) + (5 - 6,5) + (9 - 6,5)) = -1,5

3 . Понадобятся значения

_ ^ ^ - 6,5 _ к - 6,5

±цк-1У 143 143

к=\

Расчет удобно проводить с помощью таблицы (табл.З).

г — 1 5

-^ =—— = -0,01049

В табл. 3 столбец (3) получен из столбца (2) умножением на 143 143 ,

каждый элемент столбца (6) равен сумма элементов столбцов (3), (4) и (5), стоящих в той же строке, а в столбце (7) стоят квадраты соседних элементов из столбца (6). Цель построения табл.З - расчет суммы элементов столбца (7). Эта сумма равна 0,28275. Следовательно,

т/в*(&)=с*.£1%1 = 385,7 х л/0,28275 = 204,8

V к = 1

Таблица 3

Расчет дисперсии периодической составляющей

к СЛ -1 !п + 1//72 \\ Ли

(1) (2) (3) (4) (5) (6) (V)

1 -5,5 0,0577 - 0,0833 0,3333 0,3077 0,09468

2 -4,5 0,0472 - 0,0833 - -0,0361 0,00130

3 -3,5 0,0367 - 0,0833 - - 0,0466 0,00217

4 -2,5 0,0262 - 0,0833 - - 0,0571 0,00326

5 -1,5 0,0157 - 0,0833 0,3333 0,2657 0,07060

6 -0,5 0,0052 - 0,0833 - -0,0781 0,00610

7 0,5 - 0,0052 - 0,0833 - - 0,0885 0,00783

8 1,5 -0,0157 - 0,0833 - - 0,0990 0,00980

9 2,5 - 0,0262 - 0,0833 0,3333 0,2238 0,05009

10 3,5 - 0,0367 -0,0833 - 0,1200 0,01440

11 4,5 - 0,0472 -0,0833 - 0,1305 0,01703

12 5,5 - 0,0577 -0,0833 - 0,1410 0,01988

Доверительный интервал для значения периодической составляющей в январе (-376 - 1,96 х 204,8; -376 + 1,96 х 204,8) захватывает 0 (при доверительной вероятности 0,95), отличие значения периодической составляющей от 0 не значимо (на уровне значимости 0,05).

Аналогичный расчет для значения ^ = 2 (периодическая составляющая для апреля) дает

Yjil, =0,25524, p*(g2)= а * f>,22 = 385,7x^0,25524 =194,86

к-1 V к-1

Доверительный интервал для значения периодической составляющей в апреле (533 - 1,96 х 194,86; 533 + 1,96 х 194,86) = (533 - 381,93; 533 + 381,93) не захватывает 0 (при доверительной вероятности 0,95), отличие значения периодической составляющей от 0 значимо (на уровне значимости 0,05).

Приступим к завершающему этапу анализа данных табл.1 - построению

, - ч 1

М!ь =ckV-t ~г*) + —

интервального прогноза. Необходимо рассчитать величины т , если

+ ^ и м>ь -с^-г -г^ ПрИ всех остальных значениях индекса

суммирования к, где г3 - то же, что и в формуле (16), поскольку точечный прогноз х*(() является несмещенным, асимптотически нормальным, а его дисперсия оценивается согласно (18) так:

к=1

Начнем с прогноза на январь 2006 г. (по данным за 2003 - 2005 гг.). Тогда t = 13,

s = 1 ri=~I,5 w«“8Ct + 3 есди к е {1 + 40 -1), j = 1,2,3} и wk\ = %ск

%ск = 8

при всех остальных

к - 6,5 U - 52

значениях индекса суммирования. При этом 143 143 Расчет удобно

проводить с помощью таблицы (табл.4).

Сумма значений, стоящих в последнем столбце табл.4, равна 0,61299. Согласно формуле (19)

Д(13) = и(0,95)^//) * (х * (13)) = 1,96 х 385,7 х ^/О,61299 = 591,88

Согласно (14) точечный прогноз прогностической функции таков:

х *(13) = а* (13-*) + </*+/*(13) = 212,26 х 13 + 2587,48 + (-376) = 4971

Нижняя и верхняя доверительные границы для прогностической функции (с учетом как трендовой, так и периодической составляющих) имеют вид:

*™(13) = 4971 - 592 = 4379, хверт(13) = 4971 + 592 = 5563 _

Таблица 4

Расчет дисперсии прогностической функции

к 8к — 52 143 \/т

1 - 0,3077 0,3333 0,0256 0,00066

2 -0,2517 - -0,2517 0,06336

3 -0,1958 - -0,1958 0,03834

4 -0,1399 - -0,1399 0,01957

5 - 0,0839 0,3333 0,2494 0,06220

6 - 0,0280 - - 0,0280 0,00078

7 0,0280 - 0,0280 0,00078

8 0,0839 - 0,0839 0,00700

9 0,1399 0,3333 0,4732 0,22392

10 0,1958 - 0,1958 0,03834

11 0,2517 - 0,2517 0,06336

12 0,3077 - 0,3077 0,09468

Реальное значение (табл.4.1.7) - 4336. Оно практически совпадает с нижней доверительной границей прогностической функции

Аналогичные расчеты для апреля 2006 г. (? = 14, £ = 2, г2=-0>5) дают Д(14) = 1,96 х 385,7 х ^0,72480 = 643,60 Хочечный прогноз равен х*(14) = 6092, а нижняя и

х (14) = 5448 х (14) = 6736 верхняя доверительные границы таковы: ’ ’ ':срх . Реальное

значение (табл.4.1.7) - 5430. Оно практически совпадает с нижней доверительной

Границей ПрОГНОСТИЧеСКОЙ фуНКЦИИ Хшжн0-4)

8. Интервальный прогноз индивидуальных значений

Формула (19) позволяет строить интервальный прогноз для прогностической функции, т.е. для математического ожидания временного ряда. Наблюдаемое значение отличается от него на величину невязки. Распределение невязки можно оценить по значениям кажущихся невязок (см. столбец (8) в табл.1). Напомним, что это распределение не является нормальным, не описывается элементом какого-либо параметрического семейства. Интервальный прогноз индивидуального значения построить, скорректировав интервальный прогноз для прогностической функции с помощью выборочных квантилей кажущихся невязок.

Для рассмотренного выше примера вариационный ряд п = 12 кажущихся невязок таков: -748, - 512, - 464, - 280, 36, 177, 185, 208, 245, 255, 326, 571. Нижний дециль оценим как второй член вариационного ряда (-512), верхний - как предпоследний (одиннадцатый) член вариационного ряда 326. Для расчета нижней доверительной границы индивидуального значения надо взять нижнюю доверительную границу прогностической функции и отнять 512. Для расчета верхней доверительной границы индивидуального значения надо взять верхнюю доверительную границу

прогностической функции и прибавить 326.

Итак, для данных табл.1 индивидуальные значения лежат «глубоко внутри» доверительных интервалов. Прогнозы полностью оправдались.

9. О проверке условий (7)

Рассмотрим три вопроса. Верны ли условия (7) в моделях, соответствующих

реальным ситуациям? Как проверять справедливость условий по результатам наблюдений? Каковы свойства оценок, если эти условия оказываются невыполненными?

В условиях (7) важную роль играет система точек наблюдения ^ к = 1, 2, , п.

Более тщательно рассмотрим ранее принятую модель с целым числом периодов, для которой справедливо соотношение (11). При этом объем наблюдений п = mq, где т -количество периодов, д - число наблюдений в одном периоде. Предполагается, что первые д моментов наблюдения при сдвиге на длину периода дают следующие д моментов времени, при сдвиге на две длины периода дают третий набор из д моментов наблюдения, и т.д. Для значений периодической составляющей выше построены точечные оценки и доверительные интервалы (в предположении, что количество периодов т безгранично растет), в чем и состоит оценивание периодической составляющей. (Для гладкой функции ДО при безграничном росте числа наблюдений д в одном периоде можно получить сходимость оценок периодической составляющей не только в д точках, но и на всем периоде. При этом от оценок в д точках придется перейти к оценкам на всем периоде, например, кусочно-линейным, соединив соседние точки графика отрезками прямых.)

Описанная модель справедлива, когда, например, в течение некоторого числа лет имеются поквартальные или помесячные данные бухгалтерского учета. При изучения посещений сайта или торгового заведения - почасовые данные за целое число недель. Если в ряду наблюдений есть пропуски (временной ряд не является полным) -предпосылки модели не выполняются. Если система точек наблюдения не образует арифметическую прогрессию,

В рассматриваемой модели естественно принять, что

Ж) = о

(20)

суммарное отклонение значений восстанавливаемой функции от линейного тренда за один период является нулевым. Тогда первое из условий (7) выполнено:

Ело=даЕ/(^)=0

z=l к=1

В реальных ситуациях система точек наблюдения может включать в себя, кроме целого числа периодов, еще несколько начальных точек следующего периода. Можно априори принять первое условие (7), для этого изменив - при необходимости -величину свободного члена d в модели тренда (та же логика рассуждений, что и при принятии условий М(вк) = 0 - в модели без периодической составляющей - и М(Еь) = 0 в общем случае). Однако возникает противоречие между первым условием (7) и условием (20). Условие первое условие (7) автоматически обеспечивается методом наименьших квадратов, а условие (20) соответствует логике моделирования. Однако поскольку рассматриваем асимптотическую теорию при безграничном росте числа периодов, указанное различие исчезает при т 00 . Таким образом, первое из условий (7) вытекает из свойств рассматриваемой модели и потому вообще не требует проверки по экспериментальным данным, в отличие от второго условия (7), которое выполнено не всегда.

Добавим к модели с целым числом периодов два предположения -

симметричности множества {4, к = 1, 2, ..., п} относительно 1 и четности периодической составляющей f{t) относительно той же точки. Эти предположения выполнены, если, например, график f(t) симметричен относительно середины года. Тогда второе условие (7) выполнено. Ясно, что обычно нет оснований априори считать, что реальные данные описываются такой моделью.

10. Проверка второго условия (7) по экспериментальным данным

Естественно использовать статистику

у=Ъ,-<)/*(',)

у=1

f*(t)

где v j> - ранее построенная оценка периодической составляющей f(t). Оценка /*(0

v 3' является несмещенной, а потому

м<У) = І,«, -ОАО

У=1

При справедливости (6) распределение 7 является асимптотически нормальным (при безграничном росте количества периодов т). Для проверки второго условия (7), т.е. для проверки нулевой гипотезы Но: М{ У) = 0 при альтернативной гипотезе о неравенстве математического ожидания 0 достаточно оценить дисперсию 7.

В соответствии с (11) формулу (16) можно записать для любого_/ = 1, 2, ... , и, если под к = к([) понимать к{[) = / aq при максимально возможном а, при котором к(/) остается положительным, т.е. к([) - это остаток от деления / на q, если этот остаток ненулевой, и к(/) = q при нулевом остатке. Таким образом,

/*(*;■) = Д0) + Т,Ки)Ег’ 1

2=1

(21)

где - те же, что и в формуле (16). В соответствии с определением 7 из (21) следует, что

у=1

у=1

у=1

(22)

Изменим порядок суммирования во втором слагаемом в (22):

I

у=1

(/-02]к(])Ег

2=1

2=1

ЩЛ

.1=і

Следовательно, поскольку Е1 - независимые одинаково распределенные

2

случайные величины с математическим ожиданием 0 и дисперсией ° , то

ЩЛ

у=1

(23)

2 ^

Величину ° оцениваем по формулам (10) и (15), величины ти) описаны после

формулы (16). Подставив оценку ° в (23), получаем оценку Л*(7) дисперсии 7.

В соответствии с асимптотической нормальностью 7 правило принятия решений при проверке гипотезы Н0: М(7) = 0 таково: если

Y

л/D*{Y)

<С{а) = Ф_1|1-|

где - критическое значение, соответствующее уровню значимости а , то нулевая гипотеза принимается (второе условие (7) выполнено), если же неравенство (24) не выполнено, то принимается альтернативная гипотеза (второе условие (7) не выполнено).

11. Асимптотическая несмещенность оценки параметра а

Приведем пример, когда второе условие (7) не выполнено. Измерять время будем в месяцах. Пусть данные берутся на середину квартала. Тогда последовательность моментов времени такова: 2, 5, 8, 11, 14, 17, 20, 23, ... Задан период - год. Периодическая составляющая задается четырьмя числами: g\ = - 1, g2 = -2, g3 = -3, g4 =

£Ж) = -1-2-3 + 6 = 0

6. Для таких данных выполнено равенство (20), т.е. k=l

Следовательно, выполнено первое условие (7). Используя это условие, можно

упростить второе условие (7):

£(', -Ь/м=Ё<,/(<()-'Ё/(0=!<,/(<,)=о

2=1 2=1 2=1 2=1

Для простоты расчетов ограничимся двумя годами. Тогда

п

Еtjfc) = 2 х (-1) + 5 х (-2) + 8 х (-3) +11 х 6 +14 х (-1) +17 х (-2) + 20 х (-3) +

2 = 1

+ 23 х 6 = (-2) + (-10) + (-24) + 66 + (-14) + (-34) + (-60) +138 = 60

Второе условие (7) не выполнено. Оно не будет выполнено и для любого иного числа лет. Действительно, если х - начало года (для первого года х = 0, для второго х = 12, и т.д.), то вклад этого года в рассматриваемую сумму будет равен

(х + 2) х (-1) + (х + 5) х (-2) + (х + 8)(-3) + (х +11) х 6 =

= 2 х (-1) + 5 х (-2) + 8 х (-3) + 11x6 = 30

Причина нарушения второго условия (7) ясна - периодическая составляющая не симметрична в течение года. Такое поведение периодической составляющей естественно для сельскохозяйственных предприятий. Противоположную ситуацию

демонстрирует периодическая составляющая для временного ряда цен на лом черных металлов (по данным Магнитогорского металлургического комбината), проанализированного выше.

Смещение оценки параметра а равно

!/(*,)(',- о 1/(0*,

М(а*) -а = —----------- —--------

2>,-о2 Ей-о2

(25)

В рассматриваемом примере числитель за т лет равен 30т. А знаменатель, очевидно, имеет порядок т3. Смещение имеет порядок т~2, т.е. быстро убывает с ростом числа периодов. Оценка а* параметра а является асимптотически несмещенной.

Нетрудно показать, что для модели с целым числом периодов всегда имеет асимптотическая несмещенность оценки а* параметра а. Если второе условие (7) выполнено - эта оценка является несмещенной, если не выполнено - смещенной, но смещение стремится к 0 при росте числа периодов. Таким образом, выполнение второго условия (7) не является необходимым для применения рассматриваемых методов. Тем не менее проверка второго условия (7) по экспериментальным данным является полезным для решения о том, можно ли пользоваться асимптотической несмещенностью оценки при имеющемся объеме данных.

12. Заключительные замечания

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

1) дать общее обоснование этим алгоритмам в рамках асимптотических методов математической статистики и указать условия их применимости (формула (6));

2) выявить принципиально важные условия (7), необходимые и достаточные для несмещенности и состоятельности рассматриваемых оценок;

3) построить доверительные интервалы для зависимости (прогностической

функции), трендовой и периодической составляющих, индивидуальных значений временного ряда.

Обсуждение отдельных сторон рассматриваемой проблемы проведено в работах [2,4, 5].

В рамках математической статистики удается провести анализ не всех распространенных эвристических алгоритмов. Так, довольно часто рекомендуют вначале провести сглаживание («выравнивание») временного ряда, например, методом скользящих средних [1, с. 137]. При этом периодическая (сезонная) составляющая меняется (также сглаживается), а погрешности (отклонения от суммы трендовой и периодической составляющих) становятся зависимыми случайными величинами, что делает невозможным применение описанных в настоящей статье методов.

Теория устойчивости [6] отвергает идею поиска оптимального метода, поскольку зачастую оказывается, что для любого выбранного для рассмотрения метода анализа данных можно подобрать такое понимание оптимальности, что именно этот метод является оптимальным. Например, метод наименьших квадратов в определенном смысле оптимален, если погрешности имеют нормальное распределение, в то время как метод наименьших модулей оптимален, если погрешности имеют распределение Лапласа. В задаче проверки однородности двух независимых выборок установлено [V], что для любого из обычно используемых критериев однородности существует такое распределение на множестве альтернативных гипотез, что рассматриваемый критерий является оптимальным (в том смысле, который определен в [7]).

Работа выполнена в рамках новой парадигмы прикладной (математической) статистики [8, 9]. Изучена непараметрическую модель, не предполагающая, что распределение погрешностей (ошибок, невязок) входит в то или иное параметрическое семейство. Второе основание для того, чтобы именовать рассматриваемые модель и метод непараметрическими - оценивание периодической составляющей произвольного вида, т.е. отказ от выбора периодической составляющей из какого-либо параметрического семейства функций.

Полученные в статье [10] научные результаты, касающиеся средних величин и законов больших чисел в пространствах произвольной природы, могут быть применены для анализа данных в различных научных и прикладных областях. В отличие от них результаты настоящей работы нацелены прежде всего на анализ динамических рядов экономических показателей (временных рядов), необходимость которого часто возникает при организационно-экономическом моделировании с целью решения задач управления хозяйственными единицами [11]. Именно потребности экономики и управления ставят во главу угла модели с одной независимой переменной - временем. Длина периода задается существом рассматриваемой прикладной задачи (для оценки длины периода по статистическим данным нужен другой математический аппарат, разработанный в [12]). Рассмотренные в настоящей работе постановки можно относить к эконометрике [2], т.е. статистическим методам в экономике. Полученные результаты могут быть применены для прогнозирования и построения экономико-математических моделей, в частности, в рамках неформальной информационной экономики будущего [13].

Литература

1.Практикум по эконометрике: Учеб. пособие / И.И. Елисеева, С.В. Курышева, Н.М. Гордеенко и др.; Под ред. И.И. Елисеевой. - М.: Финансы и статистика. 2001. - 192 с.

2.Орлов А.И. Эконометрика. Изд. 4-е, доп. и перераб. Учебник для вузов. - Ростов-на-Дону: Феникс, 2009. - 572 с.

3.Крюкова Е.М. Применение методов организационно-экономического прогнозирования в отрасли лома черных металлов // Заводская лаборатория. Диагностика материалов. - 2008. - Т.74. -№7. - С.67 - 72.

4.Орлов А.И. Непараметрический метод наименьших квадратов: учет сезонности // Статистические методы оценивания и проверки гипотез: межвуз. сб. науч. тр. Вып. XXI. - Пермь: Перм. ун-т, 2008. - С. 135-148.

5. Орлов А. И. Непараметрический метод наименьших квадратов с периодической составляющей: условия применимости // Статистические методы оценивания и проверки гипотез: межвуз. сб. науч. тр. Вып. XXII. - Пермь: Перм. ун-т, 2010. - С.96-108.

6.Орлов А.И. Устойчивые экономико-математические методы и модели. Разработка и развитие устойчивых экономико-математических методов и моделей для модернизации управления предприятиями. - Saarbrucken: LAP, 2011. - 436 с.

7.Никитин Я.Ю. Асимптотическая эффективность непараметрических критериев. - М.: Наука, 1995.-240 с.

8.Орлов А.И. Новая парадигма прикладной статистики // Заводская лаборатория. Диагностика материалов. - 2012. - Том 78. - №1, часть I. - С.87-93.

9.Орлов А.И. Основные черты новой парадигмы математической статистики / А.И. Орлов // Политематический сетевой электронный научный журнал Кубанского государственного аграрного

университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2013. -№06(090). С. 188-214. - ША [article ГО]: 0901306013. - Режим доступа:

http://ei.kubagro.ru/2013/06/pdf/13.pdf, 1,688 у.п.л.

10.Орлов А.И. Средние величины и законы больших чисел в пространствах произвольной природы / А.И. Орлов // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. -Краснодар: КубГАУ, 2013. - №05(89). - С. 554 - 584. ГОА [article ГО]: 0891304038. - Режим доступа: http://ei.kubagro.ru/2013/05/pdf/38.pdf, 1,938 у.п.л., импакт-фактор РИНЦ=0,577

11.Орлов А.И. Организационно-экономическое моделирование при решении задач управления хозяйственными единицами / А.И. Орлов // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2013. - №03(87). - С. 679-705. - Режим доступа: http://ej.kubagro.ru/2013/03/pdf/53.pdf

12.Орлов А.И. Метод оценивания длины периода и периодической составляющей сигнала // Статистические методы оценивания и проверки гипотез. Межвузовский сборник научных трудов. -Пермь: Изд-во Пермского государственного университета, 1999. - С.38-49.

13. Орлов А.И. Проблемы методологии государственной политики и управления в неформальной информационной экономике будущего / А.И. Орлов // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2013. - №04(88). - С. 592 - 618. -ГОА [article ГО]: 0881304041. - Режим доступа: http://ei.kubagro.ru/2013/04/pdf/41.pdf

References

1. Praktikum ро jekonometrike: Ucheb. posobie / I.I. Eliseeva, S.V. Kurysheva, N.M. Gordeenko i dr.; Pod red. I.I. Eliseevoj. - М.: Finansy i statistika. 2001. - 192 s.

2. Orlov A.I. Jekonometrika. Izd. 4-e, dop. i pererab. Uchebnik dlja vuzov. - Rostov-na-Donu: Feniks, 2009. - 572 s.

3. Krjukova E.M. Primenenie metodov organizacionno-jekonomicheskogo prognozirovanija v otrasli loma chemyh metallov // Zavodskaja laboratorija. Diagnostika materialov. - 2008. - T.74. - №7. -S.67 - 72.

4. Orlov A.I. Neparametricheskij metod naimen'shih kvadratov: uchet se-zonnosti // Statisticheskie metody ocenivanija i proverki gipotez: mezhvuz. sb. nauch. tr. Vyp. XXI. - Perm': Perm, un-t, 2008.-S. 135-148.

5. Orlov A.I. Neparametricheskij metod naimen'shih kvadratov s periodicheskoj sostavljajushhej: uslovija primenimosti // Statisticheskie metody ocenivanija i proverki gipotez: mezhvuz. sb. nauch. tr. Vyp. XXII. - Perm': Perm, un-t, 2010. - S.96-108.

6. Orlov A.I. Ustojchivye jekonomiko-matematicheskie metody i modeli. Razrabotka i razvitie ustoichivyh iekonomiko-matematicheskih metodov i modelei dlja modemizacii upravleniia predpriiatiiami.

- Saarbrucken: LAP, 2011. - 436 s.

7. Nikitin Ja.Ju. Asimptoticheskaja jeffektivnost' neparametricheskih kriteriev. - М.: Nauka, 1995.-240 s.

8. Orlov A.I. Novaja paradigma prikladnoj statistiki // Zavodskaja laboratorija. Diagnostika

materialov. - 2012. - Tom 78. - №1, chast' I. - S.87-93.

9. Orlov A.I. Osnovnye cherty novoj paradigmy matematicheskoj statistiki / A.I. Orlov //

Politematicheskij setevoj jelektronnyj nauchnyj zhurnal Kubanskogo gosudarstvennogo agramogo universiteta (Nauchnyj zhurnal KubGAU) [Jelektronnyj resurs], - Krasnodar: KubGAU, 2013. - №06(090). S.188-214. - IDA [article ID]: 0901306013. - Rezhim dostupa: http://ej.kubagro.ru/2013/06/pdf/13.pdf,

1,688 u.p.l.

10. Orlov A.I. Srednie velichiny i zakony bol'shih chisel v prostranstvah proizvol'noj prirody /

A.I. Orlov // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal Kubanskogo gosudarstvennogo agramogo universiteta (Nauchnyj zhurnal KubGAU) [Jelektronnyj resurs], - Krasnodar: KubGAU, 2013. -

№05(89). - S. 554 - 584. ША [article ID]: 0891304038. - Rezhim dostupa:

http://ej.kubagro.ru/2013/05/pdf/38.pdf, 1,938 u.p.l., impakt-faktor RINC=0,577

11. Orlov A.I. Organizacionno-jekonomicheskoe modelirovanie pri reshenii zadach upravlenija

hozjajstvennymi edinicami / A.I. Orlov // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal Kubanskogo gosudarstvennogo agramogo universiteta (Nauchnyj zhurnal KubGAU) [Jelektronnyj resurs], -Krasnodar: KubGAU, 2013. - №03(87). - S. 679-705. - Rezhim dostupa:

http://ej.kubagro.ru/2013/03/pdf/53.pdf

12. Orlov A.I. Metod ocenivanija dliny perioda i periodicheskoj sostav-ljajushhej signala //

Statisticheskie metody ocenivanija i proverki gipotez. Mezhvuzov-skij sbornik nauchnyh trudov. - Perm': Izd-vo Permskogo gosudarstvennogo universi-teta, 1999. - S.38-49.

13. Orlov A.I. Problemy metodologii gosudarstvennoj politiki i uprav-lenija v neformal'noj

informacionnoj jekonomike budushhego / A.I. Orlov // Polite-maticheskij setevoj jelektronnyj nauchnyj zhurnal Kubanskogo gosudarstvennogo ag-rarnogo universiteta (Nauchnyj zhurnal KubGAU) [Jelektronnyj resurs], - Krasnodar: KubGAU, 2013. - №04(88). - S. 592 - 618. - ША [article ID]: 0881304041. -Rezhim dos-tupa: http://ej.kubagro.ru/2013/04/pdf/41.pdf