markdown:
# «Самарский научный центр Российской академии наук В.К. СЕМЁНЫЧЕВ, В.Н. КОЖУХОВА АНАЛИЗ И ПРЕДЛОЖЕНИЯ МОДЕЛЕЙ ЭКОНОМИЧЕСКОЙ ДИНАМИКИ С КУМУЛЯТИВНЫМ ...»
Представляется оправданным инструментарий, способный решать подобные сложные задачи идентификации на основе обоб щенных параметрических моделей авторегрессии-скользящего среднего (ARMA-моделей) [58], современного аппарата вычисли тельной математики, нейросетей [47], методов генетического мо делирования [57].
Для многокомпонентных моделей с логистическим трендом можно воспользоваться методом итерационной параметрической тренд-сезонной декомпозиции, предложенным в [62], который по зволяет получить аналитические выражения для отдельных ком понент.
Суть метода заключается в последовательном выделении не скольких трендов либо тренда и колебательной компоненты. На первой итерации находятся параметры логистического тренда временного ряда, затем на второй итерации после вычитания най денной модели из Yk (в случае аддитивной структуры) выделяется дополнительный тренд либо сезонная составляющая.
Их параметры находятся с помощью метода наименьших квадратов (например для полиномиальных трендов) либо методом конструирования обобщенных параметрических ARMA-моделей при помощи Z-преобразования (например для обобщенного экспо ненциального тренда или колебательной компоненты в виде сум мы нескольких гармоник).
Процесс последовательного выделения компонент ряда по вторяется до тех пор, пока наблюдается прирост в точности по вы бранным критериям моделирования и прогнозирования.
ВЫВОДЫ ПО ПЕРВОЙ ГЛАВЕ 1. Логистическая динамика носит довольно общий характер, часто встречается в объектах анализа в силу имеющей обычно ме сто ограниченности ресурсов этих объектов. Логистическими кри выми моделируются процессы в социологии и экономике;
как пра вило, это макропоказатели. Эволюционный характер логистиче ской кривой заложен в последовательной смене трех стадий роста.
Логистические функции отражают кумулятивный характер изме нения показателей во времени (например суммарные продажи то вара за весь период времени на рынке), в то время как их первые производные отражают текущие значения показателя, носят им пульсный характер (например текущие продажи товара). Модели рование логистической динамики является важной составляющей анализа глобальных показателей экономического развития.
2. Проведенный обзор известных приложений логистических траекторий и известных логистических моделей позволил сделать предложения по обобщению известных моделей Рамсея, Гомпертца и по феноменологическим мультилогистическим моделям. Практи ческой ценностью обладает приведенный автором расчет точек пе региба и точек начала стадий логистического процесса, который представляет интерес при анализе логистической динамики.
3. Все известные и предложенные модели логистических функций существенно нелинейны, порой многокомпонентны, ино гда являются суммой или мультипликацией указанных моделей. Их сложность позволила определить математический аппарат, при помощи которого целесообразно создавать инструментарий, соот ветствующий достижению цели исследований.
4. Сформулирована актуальная задача разработки и реализа ции методики оценки точности моделей и методов идентификации для моделирования и прогнозирования логистической динамики в диапазонах соотношений до 20-30% дисперсий помехи и детерми нированной компоненты, а также в широких диапазонах изменения параметров из комплекса логистических моделей.
5. В рамках данной методики необходимо сравнить точность различных методов и обосновать выбор лучшего метода. В качест ве меры точности следует, видимо, принять точечные характери стики оценок модели, прогноза, а также параметров моделей и до верительные интервалы для логистических траекторий.
6. Условием решения предыдущей задачи является разработ ка вычислительного комплекса, реализующего методы идентифи кации, методику оценки точности на тестовых и реальных выбор ках.
ГЛАВА 2. ИССЛЕДОВАНИЕ ТОЧНОСТИ МОДЕЛЕЙ РЯДОВ ДИНАМИКИ С ЛОГИСТИЧЕСКИМ ТРЕНДОМ И ОБОСНОВАНИЕ ВЫБОРА МЕТОДОВ ИХ ИДЕНТИФИКАЦИИ 2.1. Методики исследования точности идентификации временных рядов и выбор критериев оценки точности моделирования и прогнозирования В первой главе монографии была обоснована необходимость разработки и реализации методики оценки точности моделирова ния и прогнозирования временных рядов логистической динамики при больших соотношениях дисперсий стохастической и детерми нированной компонент в различных диапазонах изменения пара метров для комплекса логистических моделей.
Возникает задача обоснования выбора наиболее точного ме тода идентификации основных предложенных в монографии моде лей в рамках принятой методики.
Обычно в известных публикациях методы идентификации эконометрических моделей проверяются на реальных выборках, при этом редко сравниваются показатели точности их моделирова ния и прогнозирования.
Известна попытка сравнения точности нескольких методов идентификации на одной тестовой детерминированной выборке Dk логистического тренда Ферхюльста при дисперсии аддитивной помехи k в 5% от дисперсии этой детерминированной выборки [80].
Существенным развитием этого подхода можно считать пред ложенную в [61] методику исследования точности в широких диа пазонах изменения значений параметров моделей со стохастиче ской компонентой, дисперсия которой меняется в диапазоне зна чений коэффициента шум/сигнал (2.1), определяемого отношени ем дисперсий стохастической и детерминированной компонент:
Kn/ s 2 / 2 D. (2.1) Применим этот подход для комплекса анализируемых логи стических моделей. Анализу будет подвергаться тестовая (генери руемая компьютером) выборка наблюдений ряда динамики M Yk Dk k, k k.
При применении тестовых выборок необходимо разбить вы борку на две части: рабочую (n наблюдений) и прогнозную (l на блюдений).
По рабочей части выборки (наблюдения k 0, n 1 ) строится модель, рассчитываются критерии точности модели.
По прогнозной части выборки рассчитываются критерии точ ности прогноза (наблюдения k n, n l 1 ).
Приведем основные характеристики, используемые при вы боре моделей, построенных на основе реальных данных, и оценки их точности.
Необъясненная (остаточная) сумма квадратов возмущений k отражает меру близости значений Yko, полученных по предложен ной модели, к реальным значениям наблюдений Yk :
n USS Yk Yko.
k Для обоснования выбора лучшей по адекватности модели из комплекса используют и информационные критерии Акаике и Шварца.
Оба критерия «штрафуют» за увеличение числа параметров модели, критерий Шварца является более строгим [83, 130]. Дан ные критерии рекомендуется использовать при n / m 40.
Для задач, где n / m 40, существует скорректированный кри терий Акаике AICc AIC 2m(m 1) / (n m 1).
При наличии больших выборок можно обращаться к более точным (и, естественно, более сложным) критериям Диболда, Ма риано [1].
Наиболее распространенной до настоящего времени характе ристикой точности модели является доля дисперсии Y, объяснен ная выбранной моделью, которая описывается безразмерным ко эффициентом детерминации R 2, изменяемым от 0 до 1.
Для сравнения моделей с разным числом параметров приме няют скорректированный коэффициент детерминации Radj. Для оценки точности прогнозирования используются не сколько характеристик. Наиболее известной можно считать MAPE оценку (среднюю относительную ошибку прогноза) [81].
Хорошим считается прогноз, если значение MAPE-оценки не превышает 10%, а удовлетворительным – если находится в преде лах 20%. MAPE-оценка малочувствительна к ошибкам прогноза больших значений. В то же время ее расчет затруднителен, если значения наблюдений ряда близки к нулю. Она чувствительна к случайному выбросу в выборке, поскольку при расчете осуществ ляется деление на абсолютное значение ряда.
Поэтому вместо MAPE-оценки часто используют и второй ко эффициент Тейла (2.2) [1, 8, 81, 136]:
n l Y Yko k kT 2 100%.
k n (2.2) n l 1 n l Y Y o k k k n k n Второй коэффициент Тейла (в сравнении с коэффициентом детерминации) более устойчив к случайным выбросам, поскольку при его расчете используются сглаженные значения: суммы квад ратов наблюдений исходного и модельного рядов.
Тем не менее и второй коэффициент Тейла также рассчитыва ется по абсолютным значениям ряда. Следовательно, при наличии растущего логистического тренда в знаменателе получим его большую величину, чем при наличии падающего логистического тренда, в то время как числитель останется примерно одинаковым (поскольку он рассчитывается как сумма квадратов отклонений модельного ряда от исходной реальной выборки).
Данную проблему сравнения точности прогнозирования рас тущих и падающих логист позволяет исправить критерий, предло женный Н.Г. Загоруйко в [18]:
1 nl Z Yk Yko / Ymax Ymin 100%.
l k n В отличие от MAPE-оценки и второго коэффициента Тейла данный критерий использует для расчета не только прогнозную часть выборки: минимальное и максимальное значения берутся по всей выборке.
При условии большого диапазона изменения наблюдений ис ходного ряда данный критерий позволяет устранить разницу в оценивании прогноза для растущих и падающих трендов.
К расчету критериев оценки точности моделирования и про гнозирования при использовании тестовых данных можно подхо дить по-разному.
Приведем две методики предложенного расчета для коэффи циента детерминации и второго коэффициента Тейла.
В первой методике осуществляется сравнение рассчитывае мых модельных значений ряда Yko с зашумленными наблюдениями Yk. Тогда критерии точности определятся как в формуле (2.2). Тем самым рассчитаем, какими будут показатели точности, если ука занный метод будет применяться на реальных выборках с таким же соотношением шум/сигнал.
Можно сравнивать модельные значения ряда Yko и с исходны ми (детерминированными, генерированными) уровнями Dk (вто рая методика):
n l D Yko k kT 2 100%.
k n n l 1 n l Y D o k k k n k n Тем самым определяется, была ли найдена (и насколько точ но) истинная, предложенная исследователем модель по зашумлен ной выборке.
В данной монографии в качестве основной методики при тес тировании моделей принята первая, поскольку она дает прибли женные к реальным данным результаты – такие, как если бы иден тифицировался ряд динамики социально-экономического показа теля (подверженный влиянию множества случайных факторов, ошибок наблюдений и т.д.).
Вторая методика, отражающая точность самого примененного метода идентификации, будет продемонстрирована в п. 2.2. и 2. для симметричной логистической модели Ферхюльста и асиммет ричной модели Гомпертца соответственно.
Известно, что для моделей без присутствия в них свободного члена свойства коэффициента детерминации могут нарушаться, поэтому в модели Ферхюльста, Гомпертца и Рамсея введены кон станты.
Сравнение точности идентификации одной и той же модели разными методами на тестовых выборках осуществляется при по мощи коэффициента детерминации, а для разных моделей (приме ры идентификации которых на тестовых выборках приведены в третьей главе) применяется скорректированный Radj.
Помимо исследования точности идентификации модели и про гноза по ней, поставлена и решена задача оценки точности опреде ления параметров модели io, которые должны обладать свойствами несмещенности, эффективности, состоятельности [1, 8, 37, 42].
Несмещенные, эффективные и состоятельные оценки пара метров моделей парной и множественной линейной регрессии при применении метода наименьших квадратов (МНК) можно полу чить, если стохастическая компонента отвечает условиям Гаусса Маркова [66, 77, 81, 82].
В настоящих исследованиях их можно применить лишь при расчете параметров моделей методом итеративной тренд-сезонной декомпозиции с применением метода конструирования обобщен ных параметрических ARMA-моделей (п. 2.4, 3.1.2, 3.2.2, 3.5).
Для нелинейной регрессии можно предположить свойства центрированности (для обеспечения отсутствия систематической погрешности) и гомоскедастичности стохастической компоненты.
Можно, видимо, говорить и о законе распределения стохастической компоненты, близком к нормальному, в силу того что стохастиче скую компоненту формирует большое число независимых факто ров. Во всяком случае распределение, наверное, можно считать симметричным.
Разработанная методика оценки точности моделирования и прогнозирования временных рядов позволяет рассчитать точеч ные оценки точности при помощи множества сгенерированных выборок и расчета числовых характеристик этих оценок: матема тического ожидания и смещения оценок, дисперсии и среднеквад ратического отклонения, коэффициента вариации.
Сравнивая дисперсии оценок параметров для нескольких ме тодов идентификации одной модели, можно определить метод, дающий более эффективные оценки параметров в классе рассмот ренных методов.
Данная методика позволяет также рассчитать доверительные интервалы для математического ожидания оценок, определить, на крывает ли доверительный интервал известное истинное значение параметра, а также рассчитать доверительный интервал для про гноза.
Для известного ряда статистических реальных данных стро ится модель, рассчитываются оценки параметров, затем вычисля ется ряд остатков, вычисляется эмпирический коэффициент шум/сигнал согласно формуле (2.1).
Для полученных оценок параметров и рассчитанного Kn/ s ге нерируются тестовые выборки, вычисляются точечные оценки точности оценок параметров модели.
Затем можно перейти от точечных оценок точности к интер вальным, т.е. рассчитать доверительные интервалы для оценок па раметров модели с доверительной вероятностью (2.3):
S S M t,n1, M t,n1 (2.3) n n где t,n1 – квантиль распределения Стьюдента.
Расчет доверительного интервала прогноза осуществляется следующим образом (2.4) [80]:
Ykol t,nm SY Yk l Ykol t,nm SY, (2.4) Y Yko / n m – где m – число параметров модели, SY k средняя квадратическая ошибка.
В третьей главе монографии будут приведены примеры рас чета точечных и интервальных оценок точности с использованием тестовых выборок для параметров и прогноза модели, построенной по реальным социально-экономическим данным.
2.2. Обоснование метода идентификации симметричной трехпараметрической модели Ферхюльста Рассмотрим наиболее распространенную, во всяком случае в отечественной эконометрической практике, модель логистическо го тренда Ферхюльста (1.2) (здесь и далее с аддитивной («канони ческой») структурой стохастической компоненты k ) (2.5):
A k.
Yk (2.5) 1 A1e k В известной научной литературе до настоящего времени от сутствует сравнение точности моделирования и прогнозирования известных методов идентификации модели (2.5) в диапазоне соче таний параметров и при различной дисперсии помехи, а приводят ся лишь результаты реализации этих методов на тех или иных кон кретных выборках.
Все известные методы идентификации модели Ферхюльста основываются на использовании тех или иных преобразований мо дели тренда, приводящих к моделям, линейным по параметрам.
При этом для стохастической компоненты принимают «удоб ные» (для применения МНК) предположения о ее свойствах и/или о месте ее вхождения в структуру модели.
Наиболее распространенный в настоящее время метод иден тификации модели (2.5) заключается в переходе к обратным зна чениям уровней ряда.
При этом предполагаются «удобные» для реализации метода свойства модели: стохастическая компонента k находится в зна менателе модели тренда, и, в отличие от принятой канонической структуры, она мультипликативна по отношению к экспоненци альной функции, и для нее постулируют логарифмически нормаль ное распределение (2.6). Кроме того, уровень насыщения A0 счита ется априорно известным (в то время как он обычно представляет самостоятельный интерес для определения).
Такую «преобразованную» модель ряда можно логарифмиро вать, а затем применить МНК к полученной парной линейной рег рессии с параметрами ln A1, и нормальным законом распределе ния помехи ln k [82]:
A A A, 0 1 A1e k k ln 0 1 ln A1 k ln k. (2.6) Yk 1 A1e k k Yk Yk В этом случае стохастическая компонента анализируемого ря да динамики имеет нормальный закон распределения, соответст вуя по этому показателю условию Гаусса-Маркова для применения МНК.
Тем не менее МНК не обеспечивает оптимальность оценок па раметров анализируемой логисты. Дело в том, что МНК минимизи рует в этом случае другую среднеквадратическую невязку в (2.6) на логарифмах значений ln A0 / Yk 1 и ln A1 :
A n, o arg min ln 0 1 ln A1 k.
ln A o Yk ln A1, k Одним из условий Гаусса-Маркова, обеспечивающих опти мальность получения оценок методом наименьших квадратов, яв ляется нормальный закон распределения помехи [1, 8, 37, 42, 61] или, по крайней мере, соблюдение условия симметричности закона распределения и центрированности стохастической компоненты, которые означают отсутствие систематической ошибки в оценках параметров модели.
Поэтому постулирование асимметричного логнормального распределения невязки «удобно», но неоправданно [16, 81]. Тем не менее указанный метод идентификации широко распространен, в том числе и в зарубежной научной литературе, а преобразование вида (2.6) (без учета присутствия стохастической компоненты в модели) носит название преобразования Фишера и Прая («Fisher Pry transform») [95, 114]. Оно применяется и к большинству других моделей логистических функций [103].
Принятие указанных условий можно считать существенными недостатками известного метода идентификации, ограничиваю щими область его применения и снижающими его точность, поэто му перейдем к обзору и сравнению других методов идентификации модели (2.5).
В известном методе трех сумм исходный ряд динамики раз бивается на три равных отрезка, затем вычисляются суммы значе ний ряда внутри каждого отрезка, и определяются разности этих сумм, по которым однозначно выявляются значения оценок пара метров [80].
Методы Фишера и Готеллинга предполагают исследование дифференциального уравнения логистической функции dY / dt Y Y / A0. В обоих методах используются приближения для расчета производной.
В методе Фишера [80] приближенно вычисляются темпы при роста:
n Y, o, Bo arg min 0,5ln Yk 1 / Yk 1 BYk 1, A0 o / Bo.
Yk 1 2 o 0,5ln k Yk Yk 1, B k В методе Готеллинга [74] производная рассчитывается как разность текущего и предшествующего значений показателя Yk :
Y / k Yk Yk 1.
Данные методы, как и приведенные ниже, предполагают ад дитивное «включение» в модель стохастической компоненты уже после проведения линеаризующих преобразований.
Обычно такое искусственное введение в модель стохастиче ской компоненты после приведения ее к виду, удобному для при менения МНК, хотя и признается некорректным, но оправдывается нелинейностью модели и, как следствие, сложностью или невоз можностью других методов идентификации [1, 2, 9, 10, 11, 81].
В частности, для метода Готеллинга будем иметь Yk Yk 1 Yk21 / A0 k, где k – «новая» (не равная k в (2.5)) стохастическая компо нента. Оценки параметров, B находятся с помощью МНК:
n, B arg min Yk Yk 1 Yk 1 BYk 12, A0o o / Bo.
o o, B k Метод Юла сводит задачу идентификации параметров модели (2.5) к идентификации параметров регрессии темпов прироста Yk 1 Yk / Yk на Yk 1 :
Yk 1 Yk / Yk e 1 Yk 1 e 1 / A0 k, где k – «новая» стохастическая компонента.
Нелинейно входящие параметры в уравнении заменяют ли нейными e 1 Q, e 1 / A0 P и применяют МНК:
n Q, P arg min Yk 1 Yk / Yk Q PYk 1.
o o k Q, P Метод Родса [80] основан на взятии разности между соседни ми обратными значениями ряда, при этом идентификация осуще ствляется относительно параметров 1 e / A0 и e :
1 e e 1, k Yk 1 A0 Yk где k – «новая» стохастическая компонента.
Метод Нейра использует регрессию разности соседних обрат ных значений ряда на их сумму. Решение осуществляется относи тельно параметров (e 1) / (e 1) и (e 1) / (e 1) (2 / A0 ) :
1 e 1 2 e 1 1 1 k, Yk 1 Yk e 1 A0 e 1 Yk 1 Yk где k – «новая» стохастическая компонента.
Все перечисленные методы предполагают последовательное вычисление оценок параметров модели: на первом этапе произво дится расчет оценок параметров o, A0 ;
на втором – оценки пара o метра A1o.
В работах [74, 80] предлагается оценивать значение парамет ра A1 методом моментов:
A0 / Yk 1 A1 ek, ln A1 k ln A0 / Yk 1, а затем находить средние значения:
1 n1 1 n1 1 n ln A1 k ln A0 / Yk 1.
n k 0 n k 0 n k n k 0,5n(n 1), получаем оценку Учитывая, что k 1 n A exp 0,5 n(n 1) ln A0 / Yk 1.
o o o n k Однако при наличии в выборке значений, превышающих най денную оценку уровня насыщения логистической кривой A0, метод o моментов становится неработоспособным, поскольку возникает отрицательное число под знаком логарифма.
Казалось бы, оценка параметра A1o может быть найдена с по мощью перехода к обратным значениям уровней ряда и МНК:
1 e k o n A arg min o A1.
o A 1 o k 0 Yk A0 A Однако нужно учесть, что в этом случае получим гетероскеда стическую стохастическую компоненту:
A0 k 1 A1e k 1 1 A1e k A k Yk,, Yk A0 k 1 A1e k 1 A1e k 1 A1e k k 1 A1e k 1 1 A1e k k, k.
Yk A0 A0Yk Для компенсации гетероскедастичности (уменьшения тем са мым неэффективности оценок параметров модели [8]) можно рас смотреть возможность развития в этом направлении предложенного в работе [58] метода обобщенных параметрических моделей авто регрессии-скользящего среднего (ARMA-моделей) для модели (2.5).
Для реализации этого метода осуществим замену переменных модели A0 1/ A00, A1 A10 / A00, а затем с помощью прямого и обрат ного Z-преобразования [15] сконструируем разностную схему (для k 2 ):
Dk ( Dk 1 Dk 2 ) Dk 1, (2.7) e, где – детерминированная часть модели, Dk k Dk A00 A10e.
Первый способ реализации этого метода (ARMA I) учитывает соотношение Dk 1/ Yk k и приводит к следующей модели авто регрессии:
1/ Yk Gk Gk 1 Gk 2 Gk 1 k, k k k 1 (k 2 k 1 ), где k – гетероскедастическая стохастическая компонента.
Оценка параметра (и с учетом обозначений в (2.7) парамет ра ) находится с помощью взвешенного МНК для компенсации ге тероскедастичности. В качестве оценок весов wk можно использо вать:
1) уровни ряда Gk2, если помехи малы по сравнению с уровня ми ряда;
2) МНК в два этапа – на первом этапе модель идентифициру ется с помощью обычного МНК, на втором этапе полученные оцен ки Dk применяются в качестве весов;
3) один из методов непараметрического сглаживания: вна чале оцениваются значения k Yk Yko, где Yko – полученные сгла женные значения ряда, а затем по соседним r точкам строится ряд оценок дисперсии k, при этом теряются r-1 значения:
1 r 1 1 r k i k i.
Sk r 1 i 0 r i Полученные оценки дисперсии используются в качестве весов wk :
n 2 k G Gk 1 Gk 2 Gk 1.
o arg min k 2 wk Второй способ (ARMA II) учитывает соотношение Dk 1/ Yk k :
Yk 1 Yk Yk 2 Yk 2 Yk 1 Yk Yk 1 k 2 Yk 2 k 1 k 1 k k k 1 k 2, а затем с помощью МНК определяется:
Y Y Y n arg min 2 Yk 1 Yk k 2 k 1 k.
o k 2 wk Yk Второй этап идентификации для обоих способов реализации ARMA-метода одинаковый: находятся оценки параметров A00, A10 с помощью взвешенного МНК, а затем вычисляются МНК-оценки па раметров A0, A1 :
n 2 k D A00 A10e k.
A00o, A10o arg min k 0 wk A00, A Для улучшения результатов идентификации (уменьшения ав токорреляции наблюдений в ARMA II) целесообразно использовать приемы прореживания и сглаживания исходной выборки.
Под приемом прореживания будем понимать выбор каждого i-го наблюдения при идентификации модели, в результате чего по лучим i прореженных выборок. При этом из выборки исключаются наблюдения, обладающие тесной взаимосвязью, что должно уменьшить смещение оценки. Шаг прореживания ограничен объ емом выборки, необходимым для построения ARMA-модели (в слу чае модели Ферхюльста минимальный объем выборки равен 3).
В итоге из всех полученных результатов меры точности моде лирования и/или прогнозирования при различных шагах прорежи вания выбирается тот, для которого точность окажется выше. В ра боте [71] показано, что данный прием снижает дисперсию оценок параметров и случайной компоненты. Оптимальный шаг прорежи вания равен обычно порядку авторегрессии.
При использовании приема сглаживания по исходным данным (консолидации) берут для выполнения идентификации средние значения, получаемые из 2, 3-х и т.д. исходных наблюдений выбор ки. При этом, с одной стороны, уменьшается объем используемой выборки, но с другой стороны - уменьшается дисперсия и автокор реляция остатков. Вновь выбираются тот шаг сглаживания и те оценки параметров, при которых получаются более высокие ре зультаты оценивания точности модели.
Результат действия обоих приемов практически одинаков, но второй несколько проще в программной реализации.
Отметим, что модель (2.5) также может быть идентифициро вана с помощью численного решения МНК. Наиболее известным алгоритмом, который сводит задачу минимизации нелинейной функции МНК к итерационной минимизации линейных функций, является метод Гаусса-Ньютона1 (нелинейный итерационный МНК).
Постановка задачи идентификации в этом случае выглядит так: найти такое значение вектора параметров модели, которое бы определило локальный минимум функции ошибки n E Yk f, k min, k где f, k – регрессионная модель.
Функция f ( z ), k раскладывается в ряд Тейлора до первого члена разложения в окрестности точки ( z ) :
f ( z ), k m f, k f, k i, (z) i i где z – номер итерации, m – число параметров модели.
Перед началом работы алгоритма следует задать начальный вектор параметров.
На каждом шаге итерации минимизируется функция:
n arg min Yk f ( z ), k, ( z ) ( z ) t и предыдущее значение заменяется вектором ( z ) ( z ).
Чтобы найти значение ( z ), необходимо решать систему ли нейных уравнений ( z ) J T J J T Y f ( z ), k, где J – якобиан функции f ( z ), k в точке ( z ). Алгоритм за вершает работу, когда приращения параметров ( z ) будут малы.
Существует и улучшенный метод – метод Левенберга– Марквардта2, являющийся комбинацией градиентного метода и Wikipedia, the free encyclopedia: Gauss–Newton algorithm. URL: http: // en.wikipedia.org / wiki/Gauss–Newton_algorithm.
2 Библиотека алгоритмов ALGLIB. Алгоритм Левенберга-Марквардта. URL:
http://alglib.sources.ru.
метода Гаусса-Ньютона. В данном случае при нахождении искомых значений ( z ) ( z ) J T J diag J T J J T Y f ( z ), k вводится параметр регуляризации 0, назначаемый на каждой итерации алгоритма.
2.2.1. Исследование точности идентификации модели Ферхюльста с растущим логистическим трендом Проведем исследование точности моделирования и прогнози рования для модели Ферхюльста с растущим логистическим трен дом при аддитивной помехе девятью методами: Фишера, Готел линга, Юла, Родса, Нейра, трех сумм, ARMA I, ARMA II и решением МНК методом Левенберга-Марквардта.
Данную задачу целесообразно осуществить лишь на одних и тех же тестовых детерминированных выборках функции Ферхюль ста с регулированием ее параметров в различных диапазонах изме нения, с наложением стохастической компоненты, дисперсия кото рой меняется в широком диапазоне по отношению к дисперсии де терминированной компоненты.
В отечественной литературе известно лишь одно подобное исследование [80] оценки параметров симметричной логистиче ской функции Ферхюльста (2.5) шестью методами (трех точек, трех сумм, Фишера, Юла, Родса, Нейра), но только на одной единствен ной тестовой выборке.
Усложним задачу как по расширению диапазона значений па раметров, так и по обеспечению репрезентативности результатов.
Для реализации предложенной методики генерировались тес товые выборки логистического тренда (2.5) объемом 24, 36, 48 на блюдений.
Для характеристики области возможного применения мето дов целесообразно назначать не отдельные значения параметров тренда, а варьировать значения параметров в достаточно широком диапазоне.
На каждую выборку аддитивно накладывалась генерируемая помеха, дисперсия которой задавалась с помощью коэффициента шум/сигнал Kn / s, который варьировался от 0 до 0,3 с шагом 0,05.
Для каждого из девяти сравниваемых методов генерирова лись 12 600 выборок, а результаты усреднялись по 1 800 выборкам.
Точность моделирования оценивалась помощью коэффициен та детерминации R 2, а точность прогнозирования – с помощью MAPE-оценки.
Рассмотрим и точность оценивания отдельных параметров модели: смещение, среднеквадратическое отклонение, коэффици ент вариации оценок параметров. Для этого было сгенерировано 10 000 выборок объемом 36 наблюдений, глубиной прогноза 12 на блюдений и коэффициентом шум/сигнал в 10%.
Значения параметров модели, использованные при генерации тестовых выборок, а также найденные с помощью различных мето дов идентификации оценки параметров представлены в таблице А.1 приложения А.
На рисунках А.1-А.4 приложения А представлены результаты оценивания точности моделирования и прогнозирования для че тырех лучших по точности методов. В результате исследования по лучены близкие результаты для методов Фишера и Готеллинга, ARMA I и ARMA II.
Методы Фишера и Готеллинга позволяют приближенно вы числять производные функции, обладающие низкой точностью в условиях «зашумленности» выборки, а точность идентификации модели для всех аналитических методов сильно зависит от наборов исходных данных.
Применение взвешенного МНК в методах ARMA I и ARMA II существенно не улучшает точность идентификации.
Методы Родса, Юла и Нейра не дали удовлетворительного ре зультата даже в случае добавления «шума», дисперсия которого не превышает 5% дисперсии детерминированной компоненты (т.е.
коэффициенты детерминации объясняют менее 10% исходных данных).
Анализ показал, что только численное решение МНК методом Левенберга-Марквардта дает наиболее приемлемые результаты по точности моделирования и прогнозирования, а также несмещен ные и эффективные в классе рассмотренных методов оценки пара метров модели (2.5).
Проведем сравнение двух методик исследования точности идентификации временных рядов, представленных в п. 2.1, для ло гистической модели Ферхюльста (2.5) методом Левенберга Марквардта.
В качестве критерия точности моделирования примем коэф фициент детерминации, как и в предыдущем исследовании, а в ка честве критерия прогнозирования будем вместо MAPE-оценки применять второй коэффициент Тейла, а также критерий Н.Г. Заго руйко. Исходные параметры для генерации тестовых выборок со ответствуют значениям из таблицы А.1 приложения А.
Зависимость R2, kT2 и Z от Kn/s для сравниваемых методик рас чета оценок точности представлена на рисунке А.5 приложения А.
Для каждого значения коэффициента Kn/s (каждая точка на графи ке) генерировалось 1 800 выборок и результаты усреднялись.
Проанализированный метод идентификации по второй мето дике позволяет восстановить заданную исходную модель с точно стью 98% при дисперсии шума в 30% дисперсии детерминирован ной компоненты. В случае реальной зашумленной выборки (по первой методике) этот показатель равняется 78%.
Точность прогноза по первой методике не превышает 16%, а по второй – составляет не более 5%. Стоит отметить, что в сравне нии с предыдущим исследованием (см. рис. А.4 б приложения А) второй коэффициент Тейла и критерий Н.Г. Загоруйко дают более высокие значения точности прогноза.
Для выборки в 24 наблюдения MAPE-оценка не превышает 24%, а kT2 и Z –16% и 11% соответственно при шуме в 30%.
Видим, что оценки точности, рассчитанные по первой методи ке, всегда будут хуже полученных по второй, поскольку разброс значений зашумленного ряда больше.
Для растущего тренда в применении критерия Н.Г. Загоруйко нет видимых преимуществ по отношению ко второму коэффициен ту Тейла, поэтому в дальнейшем для растущих трендов будет при меняться критерий оценки точности прогноза по второму коэффи циенту Тейла.
Продолжим исследования точности идентификации модели Ферхюльста теперь в зависимости от объема используемой выбор ки.
2.2.2. Исследование зависимости точности идентификации модели Ферхюльста с растущим логистическим трендом от объема выборки и положения точки перегиба Как было отмечено в п. 1.2.1, наиболее важными характери стиками функции Ферхюльста являются уровень насыщения A (горизонтальная асимптота кривой), определяемый согласно фор муле (1.3), и точка перегиба - (1.4).
Исследуем точность идентификации модели Ферхюльста и точность прогноза в зависимости от положения точки перегиба.
Будем считать, что кривая Ферхюльста достигла уровня насыщения при Dk 0,95 A0. Тогда абсцисса данной точки определится следую щим образом:
0, A /.
0,95 A0 k ** ln 1 A1e k ** A Определим относительную длину выборки, необходимую для получения оптимальной оценки параметра A0 (уровня насыщения) модели Ферхюльста.
Для этого рассмотрим три случая:
1. Число наблюдений в выборке достигло значения точки пе региба k *.
2. Число наблюдений в выборке достигло половины расстоя ния между точкой перегиба и 95% уровня насыщения k * k ** / 2.
3. Число наблюдений в выборке достигло 95% уровня насы щения k **.
Исследование проводилось на тестовых выборках. Исходные значения параметров представлены в таблице 2.1.
Дисперсия стохастической компоненты задается с помощью коэффициента шум/сигнал Kn/ s. При исследовании коэффициент шум/сигнал варьируется от 0 до 0,3 с шагом 0,05.
Таблица 2. Истинные значения параметров модели Параметр Значение A 10;
50;
A 0,2;
0,4;
0,6;
0, Минимальное число наблюдений в выборке с учетом порядка авторегрессионных моделей составляют три наблюдения. Расчет длин выборок для исследования представлен в таблице А.2 прило жения А.
Для каждого сочетания параметров генерировалась 1000 вы борок. Глубина прогноза в каждом случае определялась одной тре тью длины выборки.
Выше было рассмотрено восемь аналитических методов иден тификации модели Ферхюльста с аддитивной структурой стохас тической компоненты и один численный метод решения МНК – ме тод Левенберга-Марквардта. Таким образом, для идентификации моделей использовалось девять методов.
Для повышения точности идентификации для методов ARMA I и ARMA II использовались прореживание и сглаживание исходной выборки. Точность полученных оценок определялась с помощью коэффициента вариации.
Для каждой выборки выбирался наилучший результат иден тификации по критерию коэффициента детерминации R 2. Точ ность прогноза оценивалась с помощью второго коэффициента Тейла kT2.
Результаты исследования точности идентификации и прогно зирования при различном положении точки перегиба представле ны на рисунке 2.1.
По критерию R 2 качество идентификации остается высоким даже при дисперсии шума, равной 30% от дисперсии детерминиро ванной компоненты. Более высокая точность идентификации по лучена на самых коротких выборках, равных расстоянию до точки перегиба. Точность прогноза улучшается с приближением к уровню насыщения. Результаты точности оценивания параметров модели при различном положении точки перегиба представлены на рисун ке А.6 приложения А.
40% 0, 30% 0, 20% 0, 10% 0, 0,75 0% 0 0,05 0,1 0,15 0,2 0,25 0,3 0 0,1 0,2 0, Точка перегиба Точка перегиба 1/2 расстояния до уровня насыщения 1/2 расстояния до уровня насыщения Уровень насыщения Уровень насыщения а) б) Рис. 2.1. Зависимость R (а) и kT2 (б) от Kn / s при различном положении точки перегиба Разброс получаемых оценок параметра относительно сред него значения остается приемлемым при коэффициенте шум/сигнал ниже 0,15. Точность остальных оценок крайне низка даже при шуме в 5%. Параметр A1 оценивается хуже остальных па раметров модели, в том числе при длине выборки, приближающей ся к 95% уровня насыщения.
Оценка уровня насыщения A0 возможна лишь при длине вы борки, приближающейся к уровню насыщения при уровне шума 10%. Исследование также показало, что на выборках в 16 наблюде ний и менее разброс получаемых оценок велик.
Рассмотрены также выборки длиной 25, 35 и 45 наблюдений.
Зависимость значений коэффициентов вариации параметров от Kn / s представлена на рисунке А.7 приложения А.
Показано, что, как и в предыдущем случае, оценки параметра A1 обладают наиболее низкой точностью.
Разброс оценок параметра находится в пределах 40% даже при шуме в 30%, а оценка уровня насыщения практически посто янна на выборках в 25 и в 35 наблюдений при дисперсии шума в пределах 25% от дисперсии детерминированной компоненты.
Таким образом, применяемые аналитические и численные ме тоды дают высокие результаты по критерию точности идентифи кации R 2 и точности прогноза по критерию kT2, что, однако, не оз начает, что получаемые оценки будут обладать свойствами несме щенности и эффективности.
С другой стороны, в большинстве практически важных случа ев истинные значения параметров неизвестны, в результате чего имеется возможность лишь ориентироваться на экспертные выво ды, существующие ограничения ресурсов и физиче ский/экономический смысл задач.
Методы не работоспособны на коротких выборках (менее 16-20 наблюдений), если не достигнута хотя бы половина расстоя ния между точкой перегиба и уровнем насыщения.
2.2.3. Исследование точности идентификации обобщенной модели Ферхюльста с падающим логистическим трендом Модели с растущим логистическим трендом на практике при меняются до настоящего времени чаще, а применение падающих логист сдерживается в определенной мере трудностью их иденти фикации.
Рассмотренная в п. 1.2 монографии запись тренда Ферхюльста (1.2) является общепринятой. Однако следует учитывать, что ниж няя горизонтальная асимптота данной функции равна нулю.
При моделировании условие ненулевого нижнего предела за частую не является обязательным. Например, при описании продаж товара в начальный момент времени (t = 0) они могут быть ненуле выми. Это значение определяется точкой пересечения логистиче ской кривой (1.2) с осью ординат T (0) A0 / 1 A1.
Совершенно иной представляется ситуация с падающей логи стической кривой: как правило, динамика подобных показателей (примеры были приведены в п. 1.1) стремится к некоторому нену левому уровню спада.
В связи с этим запись функции (1.2) представляется неполной.
Добавим к модели функции Ферхюльста константу, обеспечив тем самым ненулевую нижнюю горизонтальную асимптоту (или уро вень спада при рассмотрении падающей кривой):
A T (t ) C.
1 A1et Для удобства выразим точку перегиба в качестве параметра t0 ln A1 /, тогда обобщенный тренд Ферхюльста и его точка пере гиба определятся следующими выражениями:
A, t t0, T (t ) C A0 / 2.
T (t ) C t t 1 e Исследуем точность идентификации обобщенной модели Ферхюльста с аддитивной структурой стохастической компоненты (2.8):
A k.
Yk C (2.8) k k 1 e Идентификацию моделей с падающим логистическим трен дом будем проводить методом Левенберга-Марквардта. Объем ис ходных выборок зададим в 24 наблюдения, а горизонт прогноза ус тановим в восемь наблюдений. Истинные значения параметров представлены в таблице 2.2.
Таблица 2. Истинные значения параметров обобщенной модели Ферхюльста Параметр Минимальное значение Максимальное значение Шаг 10 10 C A0 50 100 –0,8 –0,2 0, k0 5 15 Коэффициент шум/сигнал варьируется от 0 до 0,3 с шагом 0,05. Расчеты производились как по первой, так и по второй мето дике исследования точности идентификации моделей.
В качестве критериев оценки точности были выбраны коэф фициент детерминации – для оценки точности моделирования, второй коэффициент Тейла и критерий Н.Г. Загоруйко – для оценки точности прогноза.
Результаты усреднялись по 1 800 выборкам (каждая точка на графиках), всего генерировалось 12 600 выборок. Результаты рас четов представлены на рисунках А.8-А.10 приложения А.
На каждом графике для сравнения представлены зависимости критериев точности от коэффициента шум/сигнал.
Результаты исследования хорошо отражают разницу как в ме тодиках оценки точности моделей, так и в точности прогноза моде лей с падающим трендом. Очевидно, что вторая методика, как и в предыдущем случае, дает более высокие результаты по критериям точности моделирования.
По критерию точности прогноза, напротив, результаты отли чаются как для методик (первая дает более низкие результаты примерно в два раза), так и для показателей точности прогноза.
Критерий второго коэффициента Тейла чувствителен к сни жению тенденции ряда, в то время как критерий Н.Г. Загоруйко, учитывающий всю выборку при расчетах, показывает стабильно высокую точность даже в широком диапазоне изменения исходно го показателя (точность прогноза не превышает 11% по первой ме тодике и 7% по второй при дисперсии шума в 30%).
Рассчитаем теперь точечные оценки точности получаемых в результате идентификации оценок параметров обобщенной моде ли Ферхюльста. Для этого будем генерировать 10 000 выборок с ко эффициентом шум/сигнал 0,2, объемом в 24 наблюдения и гори зонтом прогноза в восемь наблюдений. Истинные значения пара метров C 50;
A0 50;
0,3;
k0 14.
Рассмотрим точность получаемых оценок, сравнивая резуль таты для двух случаев:
1) расчет производится по всем выборкам, даже если иденти фикация некоторых выборок не удалась (т.е. R 2 0,5, либо оценка уровня насыщения/спада меньше нуля, либо 1 );
2) расчет производился только по удачно идентифицирован ным выборкам.
Результаты расчетов сведены в таблицы А.3, А.4 приложения А. Цифрами (1) и (2) обозначены первая и вторая методики расчета критериев точности соответственно.
Для падающих логистических кривых идентификация не уда лась для 12,35% выборок. В связи с этим заметно различие в точеч ных оценках при исключении неудачных выборок из расчетов: не сколько неудачных выборок сильно искажают общую картину ре зультатов идентификации.
При исключении неудачно идентифицированных выборок из расчетов имеем несколько заниженную оценку параметра C и не сколько завышенные оценки параметров и k0.
Причина получения таких результатов для падающей кривой заключается в том, что оценка уровня спада (по сравнению с оце ниванием уровня насыщения) чаще может оказаться отрицатель ной величиной, и в среднем для удачных выборок вариация пара метров выше, чем для моделей с растущим логистическим трендом.
2.3. Обоснование метода идентификации асимметричной модели Гомпертца Рассмотрим известные методы идентификации моделей ад дитивной структуры, включающих функцию Гомпертца в записях (1.5), (1.7), (1.8) в качестве тренда и стохастическую компоненту.
Прежде всего покажем, что традиционные методы не учиты вают изначальное присутствие стохастической компоненты в мо дели.
Так, наиболее известный способ идентификации модели с ле вой асимметрией с трендом в записи (1.7) заключается в следую щем. Параметр K данной функции не подлежит оцениванию, а предполагается известным – назначается априори экспертами.
Затем с помощью двойного логарифмирования (преобразова ние Фишера и Прая) получают уравнение в линеаризованном виде [10, 80, 95, 109]: ln ln k ln ln A k ln B.
T K Только на следующем шаге предполагается, что стохастиче ская компонента входит в модель аддитивно, а неизвестные значе ния тренда Tk заменяются реальными наблюдениями Yk :
Y ln ln k ln ln A k ln B k. (2.9) K После этого, получив линейную модель, применяют МНК:
n Y A, B arg min ln ln k ln ln A k ln B.
o o K k 0 A, B Но тогда в исходной модели стохастическая компонента должна быть аргументом экспоненциальной функции и умножать ся на B k в степени функции Гомпертца, что значительно отличает ся от принятых канонических структур:
k k Yk KAB.
e Если же изначально учитывать вхождение в модель стохасти ческой компоненты, например, аддитивно:
k Yk KAB k, то такую модель нельзя преобразовать в линейную по параметрам и применить МНК.
Для оценивания параметров линеаризованной модели (2.10) (заметим, без учета стохастической компоненты) ln Yk ln K Bk ln A (2.10) применяется описанный выше метод трех сумм, основанный на разбиении исходной выборки на три равных отрезка.
Для каждого отрезка рассчитывается сумма значений ряда как аналитически, так и эмпирически, и по полученным трем уравне ниям находятся оценки параметров модели.
Отмечено, что данный метод работает только в небольшом диапазоне изменения параметров модели, а тестирование данного метода в широком диапазоне изменения параметров, в частности на логисте Ферхюльста, дало смещенные и неэффективные оценки параметров уже при дисперсии шума в 5% [63].
Также для уточнения оценок параметров применяется итера ционная процедура П. Стонера [80]. Для линеаризованной модели (2.10) составляется система нормальных уравнений, которая сво дится к одному уравнению так, чтобы в нем присутствовал только параметр B.
Затем в полученное уравнение подставляется найденная ка ким-либо образом оценка параметра B o. Если левая часть такого уравнения больше правой, то оценку уменьшают, и наоборот. Такая процедура повторяется до тех пор, пока соотношение частей урав нения не сменится на противоположное. Уточнив оценку B o, пере считывают оценки других параметров.
Недостатком данной процедуры является то, что нормальные уравнения изначально строятся по линеаризованной модели, а не по исходной, следовательно получаемые оценки будут смещенны ми.
Описанные известные методы либо вовсе не учитывают при сутствие в модели стохастической компоненты, либо «удобно» до бавляют ее искусственно после некоторых преобразований функ ции тренда.
Применение аналитических методов, например описанного в [61] метода конструирования обобщенных параметрических ARMA моделей, невозможно из-за нелинейности модели по параметрам.
Таким образом, для идентификации модели Гомпертца требуются численные процедуры.
Запись функции (1.5) с левой асимметрией можно считать бо лее удобной, так как она включает в качестве параметра абсциссу точки перегиба. Прибавим к данной функции константу C, обеспе чив этим смещение функции вдоль оси ординат:
CA.
t t0 t t lim C A0 e e C, lim C A0 e e t t Точка перегиба функции определяется следующим образом:
t* t0, T (t*) C A0 / e.
Подобным же способом преобразуем предложенную функцию с правой асимметрией (1.8):
C A.
k k0 k k lim C A0 1 e e C, lim C A0 1 e e t t Точка перегиба функции может быть найдена следующим об разом:
t* t0, T (t*) C A0 / 1 1/ e.
Аналогично модели Ферхюльста для идентификации нели нейных моделей Гомпертца с левой и правой асимметрией и с ад дитивной структурой стохастической компоненты k k k, Yk C A0e e (2.11), k k Yk C A0 1 e e (2.12) k можно применять алгоритм Левенберга-Марквардта.
Возможно также применение эвристического алгоритма RPROP, разработанного в теории нейронных сетей [47]. При приме нении этого метода также минимизируется нелинейная функция ошибки E (). Аналогично перед началом работы алгоритма следу ет задать начальный вектор параметров, а на каждом шаге ите рации z предыдущее значение этого вектора заменяется вектором ( z ) ( z ).
Для вычисления ( z ) используется вектор градиента функ ции E (), но при этом учитывается лишь знак производной:
E (( z ) ) S, (z) ( z) ( z) ( z) sgn (z) где ( z ) – вектор коэффициентов минимизации;
S ( z ) – вектор знаков производной;
sgn x – функция, определяющая знак x:
1, x sgn x 0, x 0.
1, x Каждый элемент вектора ( z ) корректируется на каждом шаге по следующему правилу:
min ( ai ( z 1),max ), Si ( z ) Si ( z 1) i ( z ) max (bi ( z 1),min ), Si ( z ) Si ( z 1) 0, ( z 1), S ( z ) S ( z 1) i i i где a 1 и 0 b 1 – константы;
min – минимальное допустимое значение;
max – максимальное допустимое значение.
Если после предыдущего шага производная не изменила свой знак, то коэффициент минимизации увеличивается, а если изме нила – то уменьшается. Другими словами, если на предыдущем ша ге точка минимума не была достигнута, то шаг необходимо увели чить, а если алгоритм «проскочил» точку минимума, то к ней необ ходимо вернуться с меньшим шагом.
Для идентификации модели Гомпертца можно применять и генетический алгоритм, который решает задачу моделирования путём случайного подбора, комбинирования и вариации искомых параметров с использованием механизмов, напоминающих биоло гическую эволюцию1 [57].
Постановка задачи выглядит следующим образом: имеется целевая функция от многих переменных, у которой необходимо найти глобальный экстремум: f x1, x2,..., xn.
Например, это функция МНК, n 1 Q(C, A0,, k0 ) Yk C A0e k k e k которую необходимо минимизировать.
Независимые переменные необходимо представить в виде хромосом. Преобразование независимых переменных в хромосомы осуществляют с помощью кодирования, которое задается в двоич ном формате. Используются N бит для каждого параметра, причем N может быть различным для каждого параметра. Хромосома пред ставляется вектором длины L Nm, где m – число независимых пе ременных (параметров) функции.
Wikipedia, the free encyclopedia: Genetic algorithm. URL: http:// en.wikipedia.org / wiki/Genetic_algorithm.
Каждая особь состоит из массива X (0,...,L–1) и значения функ ции f на переменных, извлеченных из этого массива.
В общем случае генетический алгоритм состоит из следующих шагов:
1. Генерация начальной популяции – заполнение популяции особями, в которых независимые переменные (закодированные битами) заполнены случайным образом.
2. Выбор родительской пары: берутся K особей с минималь ными значениями целевой функции и из них составляются все воз можные пары K ( K 1) / 2.
3. Кроссинговер: случайным образом выбирается точка t на массиве X (0,...,L–1).
4. Все элементы массива с индексами 0-t новой особи (потом ка) заполняются элементами с теми же индексами, но из массива X первой родительской особи. Остальные элементы заполняются из массива второй родительской особи. Для второго потомка действия производятся наоборот: элементы 0-t берутся от второй родитель ской особи, а остальные – от первой.
5. Новые особи с некоторой вероятностью мутируют: инвер тируется случайный бит массива X этой особи. Вероятность мута ции обычно полагают порядка 1%.
6. Полученные особи-потомки добавляются в популяцию по сле переоценки. Обычно новую особь добавляют взамен самой не удачной старой особи при условии, что значение функции на новой особи меньше значения функции на старой особи.
7. Если самое лучшее решение в популяции не удовлетворяет, то осуществляется переход на шаг 2.
Критерием окончания процесса может служить заданное ко личество поколений или схождение популяции.
Более сложный генетический алгоритм содержит такие шаги, как отбор особей для размножения и генерация пар из отобранных особей. При этом каждая особь может быть задействована в одной и более паре в зависимости от алгоритма.
2.3.1. Исследование точности идентификации моделей Гомпертца при левой и правой асимметрии с растущим логистическим трендом Проведем исследование точности идентификации «классиче ской» модели Гомпертца с левой асимметрией с растущим трендом (2.11), а затем – предложенной модели с правой асимметрией (2.12) на тестовых выборках.
Тестовые выборки логистического тренда Гомпертца генери ровались объемом в 24 наблюдения и прогнозом в восемь наблю дений. Исходные значения параметров приведены в таблице 2.3.
Дисперсия генерируемой помехи задавалась с помощью коэффици ента Kn/ s, который варьировался от 0 до 0,3.
Для модели с левой асимметрией рассмотрим три метода идентификации: метод Левенберга-Марквардта, алгоритм RPROP и генетический алгоритм.
Для каждого Kn/ s результаты усреднялись по 1 800 выборкам, всего было сгенерировано по 12 600 выборок для каждого метода.
Таблица 2. Исходные значения параметров для генерации тестовых выборок Параметр Минимальное значение Максимальное значение С 10 50 A 0,2 0, 5 k Рассчитывались значения оценок параметров и меры точно сти по двум методикам, представленным в п. 2.1: точность модели рования оценивалась с помощью коэффициента детерминации R 2, а точность прогнозирования – с помощью второго коэффициента Тейла kT2.
Результаты оценки точности моделирования и прогнозирова ния для логистической модели Гомпертца с левой асимметрией представлены на рисунках Б.1, Б.2 приложения Б.
Заметим, что метод Левенберга-Марквардта значительно ус тупил по точности модели и прогноза двум другим методам. Тем не менее точность прогноза, достигаемая всеми методами идентифи кации, остается высокой (в пределах 20%) даже при дисперсии шу ма в 30% от дисперсии детерминированной компоненты.
Результаты алгоритма RPROP и генетического алгоритма практически совпадают по точности в обеих методиках, поэтому в дальнейшем не будем применять генетический алгоритм, требую щий больших временных затрат на расчеты.
Проведем аналогичное исследование с теми же исходными данными параметров для модели с правой асимметрией двумя ме тодами – Левенберга-Марквардта и RPROP. Результаты оценки точ ности моделирования и прогнозирования для логистической моде ли Гомпертца с правой асимметрией представлены на рисунках Б. и Б.4 приложения Б.
По результатам первой методики оценки точности моделиро вания и прогнозирования обоими методами практически совпада ют.
По результатам второй методики видно, что алгоритм RPROP дает несущественно более точные результаты по отысканию изна чально заданной модели.
Таким образом, для идентификации асимметричной справа модели Гомпертца с растущим логистическим трендом возможно применение любого из двух методов.
2.3.2. Исследование точности идентификации моделей Гомпертца при левой и правой асимметрии с падающим логистическим трендом Интерес может представлять и оценка точности идентифика ции падающей логистической кривой при 0, когда логистиче ская кривая будет стремиться не к уровню насыщения, а к уровню спада – нижней горизонтальной асимптоте.
В таблице 2.4 приведены исходные значения параметров мо делей Гомпертца (с левой и правой асимметрией) с падающим трендом.
Таблица 2. Исходные значения параметров для генерации тестовых выборок Параметр Минимальное значение Максимальное значение С 10 50 A –0,8 –0, 5 k Все остальные исходные данные повторяют исследование, проведенное в п. 2.3.1. В качестве критерия оценки точности про гноза, помимо второго коэффициента Тейла, примем и критерий Н.Г. Загоруйко.
Из приведенных результатов для падающей логистической кривой с левой асимметрией (рис. Б.5, Б.6 приложения Б) видно, что точность идентификации рядов с падающим логистическим трен дом Гомпертца хуже, чем для рядов с растущим трендом.
Интересен тот факт, что при использовании первой методики получаемый прогноз является недостоверным.
Уже при дисперсии шума в 5% значения второго коэффициента Тейла намного превышают обычно рекомендуемый уровень в 20%.
Вместе с тем прогноз, соотнесенный с истинными (заданны ми) выборками, для всех трех методов является достоверным в пределах 20% соотношения шум/сигнал.
Применение критерия Н.Г. Загоруйко, в отличие от второго коэффициента Тейла, позволяет соотнести результаты оценки точ ности идентификации для растущего и падающего трендов (рис. Б.7, Б.8 приложения Б).
Полученные результаты говорят о том, что точность иденти фикации падающей кривой примерно совпадает с точностью иден тификации растущей кривой. Точность прогноза для падающей кривой по значениям второго коэффициента Тейла определять не корректно, поскольку он по-разному реагирует на рост и спад в ди намике показателя.
На основе полученных данных для модели Гомпертца и с уче том исследования точности идентификации обобщенной модели Ферхюльста с падающим трендом можно предположить, что точ ность хорошего прогноза по критерию Н.Г. Загоруйко должна нахо диться в пределах 10%, что примерно соответствует точности в 20% по значениям второго коэффициента Тейла для растущих ло гистических кривых.
Повторим данное исследование для падающей логистической кривой Гомпертца с правой асимметрией. На рисунках Б.9, Б. приложения Б представлены результаты оценки точности модели рования и прогнозирования.
Точность идентификации падающей логистической кривой Гомпертца с правой асимметрией несколько ниже, чем для падаю щей кривой с левой асимметрией.
Метод RPROP дает более высокие результаты и по точности моделирования, и по точности прогнозирования.
2.4. Обоснование метода идентификации асимметричной модели Рамсея и ее обобщений В работе [48] приведен метод идентификации классической двухпараметрической модели Рамсея с аддитивной структурой стохастической компоненты.
В данной монографии будет рассматриваться трехпараметри ческая модель Рамсея с дополнительным параметром C и аддитив ной структурой стохастической компоненты (2.13):
Yk B0 (1 (1 k )e k ) C k. (2.13) Рассмотрим метод идентификации этой модели на основе конструирования обобщенной параметрической ARMA-модели при помощи Z-преобразования.
Получим ARMA-модель следующего вида:
Yk (2 1)Yk 1 (2 2 )Yk 2 2Yk 3 k, (2.14) k k (2 1) k 1 (2 2 ) k 2 2 k 3, где k – гомоскедастическая стохастическая компонента, e.
Оценку параметра можно найти с помощью метода наи меньших квадратов:
n arg min Yk (2 1)Yk 1 ( 2 2 )Yk 2 2Yk 3.
o k Зная оценку параметра o, вычислим оценку параметра :
o ln o /.
Оценки параметров C и B0 находятся с помощью МНК:
n C, B arg min Yk B0 (1 (1 o k )e k ) C.
o o o k C, B Аналогично идентифицируются модели с модифицированны ми трендами Рамсея и аддитивным вхождением в модель стохасти ческой компоненты:
Yk C ( B0 k )ek k, (2.15) Yk C ( B0 B1k )ek k. (2.16) ARMA-модель (2.14) одинакова для всех трех моделей Рамсея.
Различие возникает лишь на этапе отыскания оценок параметров, входящих в модели линейно, для этого применим обычный МНК:
n C, B arg min Yk C ( B0 o k )e k, o o o k C, B n C o, B0, B1o arg min Yk C ( B0 B1k )e k.
o o k C, B0, B Также модель Рамсея можно идентифицировать и численно, аналогично модели Ферхюльста, методом Левенберга-Марквардта, что будет рассмотрено в п. 2.4.1.
2.4.1. Исследование точности идентификации модели Рамсея с растущим логистическим трендом Оценим точность идентификации модели Рамсея (2.13) на тестовых выборках двумя методами: аналитическим методом с по мощью конструирования обобщенной параметрической ARMA модели и численным методом Левенберга-Марквардта.
Точность моделирования оценивается с помощью коэффици ента детерминации R2, а точность прогноза – с помощью второго коэффициента Тейла kT2.
Будем рассматривать выборки объемом в 24, 36, 48 наблюде ний. Дисперсия стохастической компоненты, как и в предыдущем тестировании, задается с помощью коэффициента шум/сигнал Kn/ s.
При исследовании будем изменять коэффициент шум/сигнал от до 0,35 с шагом 0,05.
Всего генерируется 134 400 выборок для каждого метода, по 19 200 выборок с каждым набором параметров (для каждого значе ния Kn/ s ). Для улучшения результатов идентификации методом ARMA будем использовать приемы прореживания и сглаживания исходной выборки.
Рассмотрим также и точность оценивания отдельных пара метров модели: смещение, среднеквадратическое отклонение, ко эффициент вариации оценок параметров. Для этого сгенерируем 10 000 выборок для каждого метода объемом в 24 наблюдения, глубиной прогноза в восемь наблюдений и коэффициентом шум/сигнал в 20%.
Значения параметров модели, использованные при генерации тестовых выборок, представлены в таблице 2.5, а в таблице В. приложения В – найденные с помощью двух методов идентифика ции их оценки.
Таблица 2. Заданные параметры модели для исследования точности моделирования и прогнозирования Параметр Минимальное значение Максимальное значение Шаг 10 25 B 0,1 0,8 0, C 50 200 При идентификации методом ARMA выбирался лучший вари ант оценивания параметров: либо по исходной выборке, либо по прореженной, либо по сглаженной.
В итоге из 134 440 выборок, идентифицированных методом ARMA, лучшие результаты были получены по исходным данным на 19 200 выборках (14,29%), с прореживанием данных – на 84 выборках (63,05%), со сглаживанием – на 30 462 выборках (22,67%).
Из рисунков В.1, В.2 и таблицы В.1 приложения В видно, что результаты идентификации модели Рамсея методом обобщенных параметрических ARMA-моделей в среднем несколько выше по сравнению с численным методом Левенберга-Марквардта.
И в том, и в другом случае с увеличением объема выборки возрастает точность прогнозирования.
2.4.2. Исследование зависимости точности идентификации модели Рамсея с растущим логистическим трендом от объема выборки и отсутствия в выборке точки перегиба Проведем исследование, позволяют ли предложенные методы спрогнозировать уровень насыщения логисты, если в выборке от сутствует точка перегиба логисты?
Оценим точность идентификации трехпараметрической мо дели Рамсея (2.13) на тестовых выборках, для которых максималь ное значение номера наблюдения не превышает значения точки перегиба функции при заданном значении параметра.
Поскольку порядок авторегрессии модели (2.14) равен трем, то минимальный объем выборки составляет четыре наблюдения.
Таким образом, максимальное значение параметра для исследо вания составляет:
k * 1/ 4 0,25.
Поскольку параметры модели C и B0 обеспечивают только ее смещение вдоль оси ординат, то для исследования примем их по стоянными, а значение параметра будем изменять от 0,05 до 0, с шагом 0,05 (табл. 2.6):
Таблица 2. Истинные значения параметров модели Параметр Минимальное значение Максимальное значение Шаг 100 100 B 10 10 C 0,05 0,25 0, Максимальный объем выборки для исследования определяет ся значениями точек перегиба. Будем рассматривать выборки с максимальным объемом, равным точке перегиба, уменьшая затем каждый раз длину выборки на одно наблюдение до тех пор, пока не достигнем минимального объема в четыре наблюдения. Обозначим через s расстояние до точки перегиба.
В таблице 2.7 приведен расчет для относительно коротких выборок: объемом от 4 до 20 наблюдений.
Выборки объемом менее 10 наблюдений (в частности 4 и 5) выбраны для того, чтобы выяснить, можно ли получить прогноз уровня насыщения на таких коротких выборках (поскольку суще ствует известный метод трех точек, когда предлагается однозначно определить три параметра логистической функции всего лишь по трем наблюдениям [80]).
Таблица 2. Определение максимальной длины выборки Точка перегиба k Макс. объем выборки s 0,05 20 20 0- 0,1 10 10 0- 0,15 6,67 6 0, 1, 0,2 5 5 0, 0,25 4 4 При исследовании коэффициент шум/сигнал варьируется от до 0,35 с шагом 0,05. Для каждого значения коэффициента шум/сигнал будем генерировать 1000 выборок. Таким образом, всего генерируется 240 000 выборок. Глубина прогноза принимает ся равной одной третьей от длины выборки. На рисунке В.3 прило жения В приведены зависимости R2 и kT2 от Kn/s для разных рас стояний до точки перегиба.
Приведем и некоторые результаты оценивания параметров модели. Оценки параметра B0 для различных значений коэффици ента шум/сигнал и разных расстояний до точки перегиба пред ставлены на рисунке В.4 а приложения В. Значения оценок пара метра по сравнению с истинными значениями для различных значений коэффициента шум/сигнал представлены на рисунке В.4 б приложения В.
Несмотря на высокие показатели точности моделирования и прогноза, определяемые коэффициентом детерминации и вторым коэффициентом Тейла, полученные оценки параметра С не соот ветствуют истинному значению.
Таким образом, невозможно правильно определить уровень насыщения для выборки, длина которой не превышает точку пере гиба. Точность определения параметра также невысока – оценки параметра завышены, имеется смещение.
Таким образом, численные эксперименты на репрезентатив ных тестовых выборках показали, что ни устойчивые численные методы, ни аналитические методы не дают достоверного прогноза уровня насыщения логисты, если длина выборки составляет менее 16 наблюдений.
Необходимым, но не достаточным условием получения досто верного прогноза является наличие точки перегиба в исходной вы борке [21].
2.5. Идентификация многопараметрических моделей логистической динамики с произвольной асимметрией Известно, что с увеличением числа параметров модели суще ственно усложняется ее идентификация, уменьшается вычисли тельная устойчивость, уменьшается возможный диапазон измене ния параметров. Поэтому ограничимся рассмотрением вопросов идентификации «разумно» сложных моделей - содержащих четыре и пять параметров.
Начнем с моделей Ричардса (2.17) и Скиадаса (2.18) при адди тивном вхождении в структуру ряда стохастической компоненты:
k k M k, Yk C A0 1 e (2.17) Tk 1 1 Tk 1 / A k.
Yk Tk 1 (2.18) 1 Tk 1 / A0 1 Данные модели являются достаточно характерными предста вителями различных классов, по которым можно сформировать суждение о трудностях их идентификации.
Модель Ричардса относится к классу асимметричных моделей, записываемых в аналитическом виде, а модель Скиадаса принад лежит к классу моделей на основе дифференциальных уравнений, не имеющих аналитического решения.
Модель Скиадаса было предложено идентифицировать мето дом Гаусса-Ньютона, так как возможность его применения была убедительно показана зарубежными исследователями на реальных данных динамики доли производства стали кислородно конвертерным способом в нескольких странах [97, 133, 134].
Вместе с тем отметим, что приведенное в п. 2.1.1 описание дан ного метода требует знания производных по параметрам функции детерминированной части модели [126, 133, 134].
Аналогичным методом можно идентифицировать и логистиче скую модель Шарифа-Кабира (2.19):
Tk 1 1 Tk 1 / A k.
Yk Tk 1 (2.19) 1 Tk 1 / A0 1 Были проведены исследования точности идентификации мо делей Скиадаса и Шарифа-Кабира методом Гаусса-Ньютона на тес товых выборках.
Они показали, что данный метод неработоспособен в широком диапазоне сочетаний параметров и величины помехи. Он может быть применен только на узком динамическом диапазоне («подхо дящих») параметров при дисперсии шума не более 5% от диспер сии изучаемого тренда. Необходимо предложить другой метод, ко торый обеспечит высокую точность при различных сочетаниях па раметров и значительной величине помехи.
Модель Ричардса широко используется на реальных данных [89, 96, 99-102, 110, 104, 129, 141]. Ее идентификация производится, как правило, либо методом Левенберга-Марквардта, либо с помо щью его упрощенного варианта – метода Гаусса-Ньютона. Тестиро вание модели при различных сочетаниях параметров и помехи также в известной литературе не проводилось.
Для модели Басса (2.20) существует известный и достаточно простой метод идентификации, именно поэтому она очень часто применяется на практике зарубежными исследователями [85, 87, 101, 103, 102, 112, 132, 141]:
Yk Tk 1 Tk 1 C 1 Tk 1 / A0 k. (2.20) Этот метод предполагает раскрытие скобок в (2.20) и введе ние новых обозначений для параметров:
C C Yk C 1 Tk 1 Tk21 k, 1 Q, P.
A0 A0 A0 A Получим выражение Yk C QTk 1 PTk21 k, в которое все параметры входят линейно, т.е. для идентификации можно применить классический МНК:
C o, Q o, Po arg min Yk C QTk 1 PTk21.
C,Q,P Затем по найденным МНК-оценкам параметров можно вер нуться к идентифицируемым параметрам, решая квадратное урав нение:
A0 A0 C o Q o A P A0 1 Q A C 0, o2 o o P A o Qo 1 D o D 1 Q o, P o A0.
4P C, A o o o o 0o 2P Причем из двух найденных при этом значений A0 выбираем o то, которое удовлетворяет условию A0 0. В данном методе имеет o ся существенный недостаток: начальное значение логистической кривой T0 также является параметром.
Метод предусматривает либо априорное назначение данного параметра экспертами, исходя из содержания рассматриваемого приложения, либо его поиск проведением вариации этого парамет ра с малым шагом в заданном исследователем диапазоне и приня тие того значения T0, которое обеспечит лучшие точностные харак теристики моделирования.
Указанные обстоятельства усложняют идентификацию моде ли, увеличивают время расчета, требуют обоснования шага вариа ции, поэтому при построении модели обычно игнорируют присут ствие стохастической компоненты в исходной модели и искусст венно «добавляют» ее на конечном этапе преобразований, т.е. рас сматривают модель в виде Yk C QYk 1 PYk21 k.
Тогда начальное значение Y0 не является параметром, а опре деляется по имеющейся выборке.
Оценки параметров C o, Qo, Po также находятся при помощи МНК: C o, Q o, Po arg min Yk C QYk 1 PYk21.
C,Q,P Исследования точности идентификации модели Басса на тес товых выборках в широком диапазоне сочетаний параметров и дисперсии помехи в известной зарубежной литературе не приво дятся.
Покажем, что модели Скиадаса (2.18) и Басса (2.20) предпоч тительнее идентифицировать предложенным методом RPROP.
Для этого рассмотрим достаточно короткие выборки объемом в 24 наблюдения, но содержащие точку перегиба. Выборки в 36, и более наблюдений рассматривать не будем, так как с увеличени ем объема выборки точность идентификации логистических моде лей естественно повышается.
Точность моделирования оценим с помощью коэффициента детерминации, а точность прогнозирования с помощью второго ко эффициента Тейла. Коэффициент шум/сигнал изменяется в диапа зоне от 0 до 0,3 с шагом 0,05. Параметры T0 и A0 (нижняя и верхняя горизонтальные асимптоты) фиксированы: T0 0,01, A0 10. Пара метры и варьируются: от 0,5 до 0,9 с шагом 0,1;
- от 1 до 1, с шагом 0,1. Выбор параметров обусловлен условием наличия точ ки перегиба в выборке.
Помимо этого, по рекомендациям авторов модели наложим следующее ограничение на параметры: при / 2 возникает хао тическая осцилляция на уровне насыщения при приближенном расчете решения дифференциального уравнения методом Эйлера (который и используется в модели).
Выполнена генерация 21 000 выборок, по 100 с каждым набо ром параметров. Не удалось осуществить идентификацию (R2 0,5) только для 226 (1,08%) выборок. Результаты оценки точности идентификации даны в таблице 2.8.
Таблица 2. Оценки точности моделирования и прогнозирования модели Скиадаса методом RPROP Kn/ s 0 0,05 0,1 0,15 0,2 0,25 0, 1 0,957 0,915 0,873 0,827 0,786 0, R kT 2 0,01% 8,55% 14,97% 18,53% 25,65% 29,30% 32,05% Из таблицы 2.8 видно, что точность прогнозирования остается удовлетворительной при шуме в пределах 15% ( kT 2 20% ). Как пра вило, для реальных данных такая точность является достаточной.
Близкие по точности результаты данный метод идентифика ции дает и для модели Шарифа-Кабира.
Теперь исследуем точность идентификации модели Ричардса (2.17) методами Левенберга-Марквардта и RPROP. Генерируются тестовые выборки объемом в 36 наблюдений, содержащие точку перегиба (кривая Ричардса требует больший объем выборки, по скольку это пятипараметрическая функция). Коэффициент шум/сигнал варьируется в диапазоне от 0 до 0,3 с шагом 0,05. Точ ность моделирования также оценивается с помощью коэффициента детерминации, точность прогнозирования - с помощью второго ко эффициента Тейла.
Исходные значения параметров представлены в таблице 2.9.
Для каждого метода выполнена генерация 25 200 выборок, по 100 с каждым набором параметров. Методом Левенберга-Марквардта не удалось осуществить идентификацию (R2 0,5) для 638 (2,53%) выборок, а методом RPROP не удалось осуществить идентифика цию для 713 (2,83%) выборок. Результаты оценки точности иден тификации представлены на рисунке Г.1 приложения Г.
Из полученных результатов видно, что метод RPROP дает бо лее точный прогноз, в то время как по точности моделирования оба метода сравнимы.
Необходимо отметить, что применение метода Левенберга Марквардта, обычно используемого для построения модели Ричар дса, также оправдано, поскольку он дает приемлемый прогноз (kT2 находится в пределах 20% при шуме в 30%).
Таблица 2. Значения параметров кривой Ричардса для генерации тестовых выборок Параметр Минимальное значение Максимальное значение Шаг C 20 20 A0 30 30 0,2 0,8 0, k0 5 10 M –5 –1 Сравним и два метода идентификации модели Басса: извест ный метод, применяющий МНК при линейном вхождении всех па раметров, и метод RPROP. Будем генерировать выборки объемом в 36 наблюдений, изменяя коэффициент шум/сигнал от 0 до 0,3 с ша гом 0,05.
Параметры модели Басса, соответствующие начальному зна чению и уровню насыщения, примем постоянными: T0 0,01;
A0 10.
Остальные параметры варьируются: - от 0,2 до 0,8 с шагом 0,2;
C - от 0 до 0,1 с шагом 0,05.
Заметим, что при нулевом значении С модель Басса соответст вует модели Ферхюльста. Таким образом, МНК для модели Басса при С = 0 есть метод Готеллинга для модели Ферхюльста. Анало гично предыдущим исследованиям результаты оцениваются с по мощью коэффициента детерминации и второго коэффициента Тейла.
На рисунке Г.2 приложения Г представлены графики зависи мости R2 и kT2 от коэффициента шум/сигнал. Видим, что МНК нера ботоспособен в широком диапазоне изменения параметров модели.
Метод RPROP дает лучшие результаты как по точности моделиро вания, так и по точности прогнозирования, шире по возможному диапазону изменения значений параметров модели Басса.
Таким образом, для моделей с произвольной асимметрией лучший результат по точности обеспечил алгоритм RPROP. Он яв ляется альтернативой методу Гаусса-Ньютона, поскольку послед ний работоспособен лишь на «подходящих» практически незашум ленных выборках (с дисперсией помехи в пределах 5%) и лишь при определенных наборах параметров.
Алгоритм RPROP, напротив, показывает стабильные результа ты при увеличении шума до 15% (а для модели Басса – до 30%) при различных сочетаниях параметров модели. В случае модели Ричар дса оправдано применение как метода Левенберга-Марквардта, так и алгоритма RPROP, причем для получения более точного прогноза рекомендуется применять второй алгоритм.
ВЫВОДЫ ПО ВТОРОЙ ГЛАВЕ 1. Для обоснования выбора методов идентификации комплек са анализируемых логистических моделей была предложена мето дика, основанная на генерации тестовых выборок и последующем расчете критериев точности моделирования и прогнозирования. В нее вошли две методики расчета критериев точности в зависимо сти от целей исследования: первая отражает точность моделей по отношению к исходным зашумленным данным, вторая – точность самого метода идентификации. Методика также позволяет полу чать точечные и интервальные оценки точности оценок парамет ров моделей. Получение интервальных оценок будет рассмотрено в третьей главе с использованием реальных выборок временных ря дов.
2. Рассмотрен комплекс моделей логистических трендов: Фер хюльста, Рамсея, Гомпертца, моделей с произвольной асимметрией, полученных численным решением дифференциальных уравнений.
Для симметричной логистической модели Ферхюльста прове дено исследование девяти методов идентификации в широком диапазоне сочетаний параметров и при различной дисперсии по мехи. Показана некорректность преобразования Фишера и Прая, а также некорректность введения в модель логистического тренда стохастической компоненты после проведения линеаризующих преобразований. Обосновано применение метода Левенберга Марквардта для идентификации модели Ферхюльста.
3. Модель Ферхюльста в обобщенном виде позволяет иденти фицировать падающую симметричную логистическую кривую с оценкой уровня спада. Методом Левенберга-Марквардта проведе ны расчеты критериев качества модели и прогноза по тестовым выборкам, а также проанализированы точечные оценки точности получаемых оценок параметров моделей.
4. Для модели Рамсея сравнены два метода идентификации:
численный метод Левенберга-Марквардта и аналитический метод конструирования обобщенных параметрических ARMA-моделей.
Показано, что метод ARMA-моделей имеет преимущество по точно сти.
Доказана необходимость присутствия в анализируемой вы борке точки перегиба для идентификации логистической модели и получения несмещенных и эффективных оценок параметров. Про демонстрирована и необходимость ограничения, накладываемого на длину исходной выборки, на симметричной и асимметричной логистах методами Левенберга-Марквардта и конструирования обобщенных параметрических ARMA-моделей. Доказана необходи мость иметь выборку объемом не менее 16 наблюдений для полу чения точных прогнозов уровня насыщения и включать в выборку точку перегиба.
5. Предложена новая модификация модели Гомпертца с фик сированной правой асимметрией. При идентификации моделей ме тодами Левенберга-Марквардта, алгоритма RPROP и генетического алгоритма продемонстрированы различия в результатах примене ния методик расчета точности. Модели Гомпертца с высокой точ ностью идентифицируются всеми методами, но лучшие результаты для растущего тренда с левой асимметрией показали алгоритм RPROP и генетический алгоритм. Поскольку последний требует больших временных затрат на проведение расчетов, то рекоменду ется применять алгоритм RPROP. Для модели с правой асимметри ей возможно применение как метода Левенберга-Марквардта, так и алгоритма RPROP.
При анализе моделей с падающим трендом была продемонст рирована некорректность применения второго коэффициента Тей ла для оценки точности прогноза. Большей универсальностью и точностью обладает критерий Н.Г. Загоруйко, который дает сопос тавимые результаты как для растущего, так и для падающего трен дов. Показано, что для достоверных прогнозов значение критерия Н.Г. Загоруйко для падающих кривых не должно превышать 10%, что соответствует рекомендованному значению второго коэффи циента Тейла в 20% для растущих кривых.
Точность идентификации падающей логистической кривой Гомпертца с правой асимметрией несколько ниже, чем падающей логисты с левой асимметрией. В обоих случаях алгоритм RPROP да ет лучшие результаты и по точности моделирования, и по точности прогнозирования.
6. Для моделей, полученных по дифференциальным уравнени ям (моделей Скиадаса, Шарифа-Кабира и Басса), рекомендуется применять алгоритм RPROP. Часто используемый эконометристами метод Гаусса-Ньютона практически неработоспособен в широких диапазонах изменения параметров и при высокой дисперсии поме хи так же, как и классический МНК для модели Басса. Лишь метод RPROP дает точные прогнозы в пределах мощности шума 15% для модели Скиадаса и 30% - для модели Басса.
7. Модель Ричардса оправдано идентифицировать как мето дом Левенберга-Марквардта, так и алгоритмом RPROP, поскольку они дают сравнимые по точности результаты.
ГЛАВА 3. ПРИМЕНЕНИЕ КОМПЛЕКСА МОДЕЛЕЙ С ЛОГИСТИЧЕСКИМ ТРЕНДОМ И МЕТОДОВ ИХ ИДЕНТИФИКАЦИИ НА РЕАЛЬНЫХ ВЫБОРКАХ ЭКОНОМИЧЕСКОЙ ДИНАМИКИ Комплекс моделей логистической динамики и методов их идентификации реализован в виде программных модулей, напи санных на языке Delphi в среде CodeGear RAD Studio 2009, которые подключаются к четырем различным приложениям:
- «Эконометрическое моделирование динамических рядов»;

Content-Type: text/html; charset=WINDOWS-1251


markdown:
Đăng ký
Trợ giúp
Ghi nhớ?
Diễn đàn
Bài gửi hôm nay
Hỏi Đáp
Lịch
Cộng đồng
Danh sách thành viên
Hoạt động
Đánh dấu chuyên mục đã đọc
Liên kết nhanh
Danh sách BQT diễn đàn
Có gì mới?
ban si quan jean
Tìm nâng cao
Tin nhắn hệ thống
Bạn có thể
đăng ký để trở thành thành viên
và bắt đầu đăng bán sản phẩm của mình.
Tin nhắn hệ thống
Không có Chủ đề. Nếu bạn theo đúng địa chỉ, xin hãy báo cho
BQT
Chọn nhanh
Khu vực riêng
Tùy chỉnh
Tin nhắn riêng
Theo dõi
Đang lên mạng
Tìm các chuyên mục
Trang chính
Diễn đàn
Mua bán thời trang
Thời Trang Nam
Thời Trang Nữ
Giày dép, túi xách, balô
Nước Hoa - Mỹ Phẩm - Trang Sức - Làm đẹp
Mẹ và Bé
Mua Bán Điện Thoại
Chợ Mobile
Chợ Apple
Linh kiện điện thoại
Chợ sim số - Thẻ cào
Máy tính - Laptop
Máy tính bảng
Máy tính - Linh phụ kiện - Thiết bị mạng
Chợ Laptop
Thiết Bị Kỹ Thuật Số & Điện Tử
Thiết Bị Kỹ Thuật Số
Chợ thiết bị Điện Tử-Điện Lạnh-Gia Dụng
Đồ Chơi Mô Hình & Games
Dịch vụ - Sửa chữa - Cài đặt
Chợ ve chai điện tử
Mua bán tổng hợp
Rao vặt đồ linh tinh
Xây dựng - Nội thất
Chợ Xe máy, xe đạp
Phụ tùng Xe Máy
Mua bán - Cho thuê Ô tô
Mua bán - Thuê nhà đất
Hợp tác kinh doanh - Kêu gọi đầu tư
Website - Phần mềm - Mạng và ADSL
Ẩm thực & Giải trí
Siêu thị Thú nuôi - Cây cảnh
Thực phẩm - Đồ uống
Hoa - Quà tặng - Đồ trang trí
Sách - Học hành - Tuyển dụng
Tuyển dụng - Tìm việc
Chợ Sách
Học hành
Thảo Luận
Chuyện trò
Chuyện trò Xe máy
Chuyện trò Ô Tô
Chuyện trò Mobile
Chuyện trò Computer
Âm thanh
Sức Khoẻ
Giới thiệu Mualoi
Nội qui và hướng dẫn sử dụng
Quảng cáo Mualoi
Thông báo và góp ý
Quảng cáo của Mualoi
Mualoi.Com Trên Facebook
Liên hệ quảng cáo
Hổ trợ trực tuyến
-- vB4 Default Style
-- Default Mobile Style
-- English (US)
-- VietNam
Liên hệ
Mua Ban - Dien Dan Mua Ban
Lưu trữ
Lên trên
Theo giờ GMT +7. Bây giờ là
07:51 PM
.
Powered by
vBulletin®
Version 4.2.0
Copyright © 2013 vBulletin Solutions, Inc. All rights reserved.
Xem website tốt nhất ở chế độ
màn hình 1024 x 768
Website : Mualoi.com
Tel: 0908.511.316
Email:
vantien2050@gmail.com
Copyright © 2010 by MuaLoi -
Diễn Đàn Mua Bán
Liên hệ quảng cáo
:
Mualoi.com
Tel: (+84) 908.511.316
Email:
vantien2050@gmail.com
Hotline
:
0908.511.316 (Mr. Tiến)
Lưu ý
: MuaLoi.com là diễn đàn miễn phí dành cho quảng cáo mua bán các sản phẩm được pháp luật Việt Nam cho phép kinh doanh.
Chúng tôi không chịu trách nhiệm trên các quảng cáo của người dùng đăng tải trên diễn đàn. Các trường hợp vi phạm, chúng tôi sẽ tiến hành xử lý xóa bỏ ngay lập tức. Nếu bạn nhận thấy quảng cáo vi phạm, xin vui lòng
gởi thông báo cho chúng tôi
, xin cám ơn!
quan jean nam

WARC-Type: response

Content-Type: application/http; msgtype=response

WARC-Truncated: length

WARC-Identified-Payload-Type: application/xhtml+xml

Server: LiteSpeed

Connection: close

Pragma: private

Accept-Ranges: bytes

Cache-Control: private

Set-Cookie: bb_lastactivity=0; expires=Mon, 20-May-2019 12:51:50 GMT; path=/

Vary: Accept-Encoding

Content-Type: text/html; charset=UTF-8


markdown:
[Kinh nghiệm quốc tế](http://www.nclp.org.vn/kinh_nghiem_quoc_te)
|
|
Xác định người đại diện của pháp nhân trong tố tụng hình sự Pháp và một vài gợi ý trong tố tụng hình sự của Việt Nam
|
|
Xác định người đại diện của pháp nhân trong tố tụng hình sự (TTHS) đối với pháp nhân là một trong những vấn đề quan trọng trong quá trình truy cứu trách nhiệm hình sự (TNHS) pháp nhân. Bởi vì trong TTHS, xuất phát từ đặc thù của pháp nhân - là một chủ thể pháp luật được tạo thành từ các thể nhân, dẫn tới mọi hoạt động tố tụng của pháp nhân bị truy cứu TNHS phải được tiến hành thông qua một thể nhân - người đại diện. Việc thực hiện các hoạt động tố tụng của pháp nhân thông qua người đại diện không phải là đặc thù của TTHS mà cả trong tố tụng dân sự, thương mại. Nghiên cứu thực tiễn TTHS của các quốc gia khác mang lại cho chúng ta những kinh nghiệm hữu ích.
TNHS đối với pháp nhân được áp dụng tại Pháp từ năm 1994 với Bộ luật Hình sự (BLHS) mới ngày 22/7/1992
|TS. Nguyễn Văn Quân - Khoa Luật - Đại học Quốc gia Hà Nội.|

Server: Zope/(Zope 2.9.7-final, python 2.4.4, win32) ZServer/1.1 Plone/2.5.3-final

Via: 1.0 www.nclp.org.vn

Cache-Control: max-age=1800, public, must-revalidate

Vary: Accept-Encoding,User-Agent

Content-Language: vi


markdown: Hyperdimension Neptunia - Hello New World [Cập nhật lúc: 08:53 12/06/2016] Tác giả Idea Factory & Compile Heart Tình trạng Đang tiến hành Thể loại Action - Comedy - Fantasy - Slice of Life Lượt xem 5.108 Hyperdimension Neptunia - Hello New World Xếp hạng: 5/5 - 3 Lượt đánh giá. 13 Người Đã Theo Dõi Đọc từ đầu Đọc mới nhất Nội dung Đây là một bản spin-off của series game cùng tên.Vùng đất Gamindustri được chia làm 4 quốc gia: Planeptune, Lastation, Lowee, và Leanbox, mỗi quốc gia được đặt dưới sự cai quản của một nữ thần. Đây là những mẩu chuyện về cuộc sống của 4 nữ thần, sau khi hiệp định hòa bình được kí kết. Danh sách chương Số chương Cập nhật Lượt xem Chapter 14.5 12/06/16 3.100 Chapter 14 12/06/16 3.098 Chapter 13 24/05/16 3.087 Chapter 12 02/05/16 3.085 Chapter 11 13/04/16 3.109 Chapter 10 28/03/16 3.091 Chapter 9 16/03/16 3.101 Chapter 8 29/02/16 3.120 Chapter 7 20/02/16 3.113 Chapter 6 12/02/16 3.109 Chapter 5 31/01/16 3.116 Chapter 4.5 21/01/16 3.124 Chapter 4 14/01/16 3.121 Chapter 3.5 04/01/16 3.106 Chapter 3 28/12/15 3.134 Chapter 2 19/12/15 3.147 Chapter 1 11/12/15 3.348

WARC-Identified-Payload-Type: text/html

Cache-Control: no-cache

Pragma: no-cache

Content-Type: text/html; charset=utf-8


markdown:
Help
Ghi nhớ?
Diễn đàn
Bài gửi hôm nay
Hỏi - Đáp
Lịch
Cộng đồng
Danh sách Thành viên
Thiết lập nhanh
Đánh dấu diễn đàn là đã đọc
Liên kết Nhanh
Danh sách Ban Quản Trị
Có gì mới?
Tìm kiếm nâng cao
Danh sách Thành viên
W164
Tổ chức chụp ảnh áo dài, áo yếm ở làng cổ Phước Lộc Thọ ngày 08/06/2014
Hình ảnh chuyến đi chụp ảnh dã ngoại Đà Lạt tháng 12/2013
Hướng dẫn upload và chèn nhiều ảnh & text vào diễn đàn một cách nhanh chóng
If this is your first visit, be sure to check out the
FAQ
by clicking the link above. You may have to
register
before you can post: click the register link above to proceed. To start viewing messages, select the forum that you want to visit from the selection below.
W164
Member
Tìm tất cả bài viết
Tìm tất cả Bài bắt đầu
Mini Statistics
Ngày tham gia
11-11-2013
Hoạt động cuối
30-11-2015
02:43 PM
Khách ghé thăm gần đây
10 khách ghé thăm gần đây là:
amthuc2234kld
,
lamh590 vsvs32
,
minhl541vsvs34
,
tanthienchi@yahoo.com.vn
,
vaobong22
,
vaobong23
,
vaobong25
,
vaobong29
,
vaobong31
,
vaobong32
Trang này được xem 1,316 lần
Tab Content
Về tôi
Về W164
Thông tin cơ bản
Thống kê
Tổng số bài gửi
Tổng số bài gửi
74
Gửi mỗi ngày
0.04
Total Thanks
Total Thanks
155
Thanked 124 Times in 68 Posts
Thông tin chung
Hoạt động cuối
30-11-2015
02:43 PM
Ngày tham gia
11-11-2013
-- Default Style
-- Default Mobile Style
-- English (US)
-- VN
Nhiếp Ảnh Miễn Phí
Lưu trữ
Lên trên
Múi giờ GMT +8. Bây giờ là
06:29 AM
.
Diễn đàn sử dụng vBulletin® Phiên bản 3.1.1
Bản quyền của 2018 vBulletin Solutions, Inc. Tất cả quyền được bảo lưu.
Ban quản trị không chịu trách nhiệm về nội dung do thành viên đăng.


markdown:
# 福山教室の講座を探す
![ 他の教室の講座も検索する](/common_rn13b/images/btn_other_n.png)
![ 他の教室の講座も検索する](/common_rn13b/images/btn_other_n.png)
- 残席状況は変動しますので、申し込み時には異なる場合があります。
- ほとんどの講座は入会が必要です。
新しく始まる講座
18時以降に始まる講座
電話か窓口にてお問い合わせください。
2講座中 1〜2件目を表示しています。
[石鹸と自然派コスメ](/nhkcc-webapp/web/WKozaClick.do?rd=10177066&kd=1661687&td=986251)
〜季節の天然素材で作る〜
〜季節の天然素材で作る〜
- 講師
- ハンドメイド石けんシニアソーパー 岩本 尚美
[開運 九星気学](/nhkcc-webapp/web/WKozaClick.do?rd=10177066&kd=1661713&td=1132405)
〜人生を楽しく〜
〜人生を楽しく〜
- 講師
- 高野山真言宗正光寺住職 占い処「文殊坊」鑑定士 野田 泰洋

Content-Type: text/html; charset=Shift_JIS


markdown: Tr Thiết kế shop cửa hàng , tủ kệ trưng bày sản phẩm- Nội thất Shop Tiếng Việt English Tiếng Việt Hotline: 0911.133.855 Địa chỉ: 164 Nguyễn Xiển, Thanh Xuân, Hà Nội Trang chủ Giới thiệu Giới thiệu về Nội thất Shop Nhà xưởng Dịch vụ Thiết kế và thi công nội thất Thi công , sản xuất nội thất Gói thiết kế - thi công nội thất Tin liên quan Xả Kho Nội Thất Sản phẩm Kệ trưng bày Kệ tủ trưng bày điện thoại Kệ tủ trưng bày mỹ phẩm Kệ tủ trưng bày giày dép Kệ trưng bày trái cây Kệ tủ trưng bày túi xách Kệ trưng bày thiết bị vệ sinh Kệ tủ trưng bày trang sức Kệ trưng bày gạch Kệ trưng bày siêu thị Kệ tủ trưng bày quần áo Biển quảng cáo Kiosk Booth quảng cáo Kệ điện máy Phụ kiện cho shop Các mẫu kệ trưng bày khác Online Friday Giang hàng hội chợ Dự án Mẫu dự án bán chạy Cửa hàng thời trang Thiết kế shop thời trang Cửa hàng điện thoại, điện máy, máy tính Cửa hàng mỹ phẩm, làm đẹp, spa Cửa hàng thực phẩm chức năng, thực phẩm sạch Siêu thị mini Showroom rượu Cửa hàng vật liệu xây dựng và thiết bị vệ sinh Gian hàng hội chợ, TTTM Cửa hàng mẹ và bé Cửa hàng nữ trang, vàng bạc Cửa hàng dụng cụ thể thao Tư vấn khách hàng Phong thủy Hướng dẫn đặt hàng Hỏi đáp - FAQs Tin tức Không gian sống Ý kiến khách hàng Đối tác Tuyển dụng Liên hệ Chăm sóc khách hàng Khuyến mại Tìm kiếm: Lĩnh vực/Mảng Công sởDược phẩmMáy tínhMỹ phẩmQuà tặngSiêu thịThế thaoThời trangThực phẩmTruyền thôngXây dựngĐiện thoại Tỉnh/Thành phố An GiangBà Rịa - Vũng TàuBắc GiangBắc KạnBạc LiêuBắc NinhBến TreBình ĐịnhBình DươngBình PhướcBình ThuậnCà MauCần ThơCao BằngĐà NẵngĐắk LắkĐắk NôngĐiện BiênĐồng NaiĐồng ThápGia LaiHà GiangHà NamHà NộiHà TĩnhHải DươngHải PhòngHậu GiangHòa BìnhHưng YênKhánh HòaKiên GiangKon TumLai ChâuLâm ĐồngLạng SơnLào CaiLong AnNam ĐịnhNghệ AnNinh BìnhNinh ThuậnPhú ThọPhú YênQuảng BìnhQuảng NamQuảng NgãiQuảng NinhQuảng TrịSóc TrăngSơn LaTây NinhThái BìnhThái NguyênThanh HóaThừa Thiên HuếTiền GiangTP.Hồ Chí MinhTrà VinhTuyên QuangVĩnh LongVĩnh PhúcYên Bái Diện tích Dưới 100 m2100 m2 - 200 m2200 m2 - 500 m2500m2 - 1000m2Trên 1000 m2 Vật liệu/Phụ kiện MikaInoxKính màuBaambooMDF woodMelamine An Cuong WoodPine WoodRussian oak woodVerneer wood Tìm kiếm Tiếng Việt English Tiếng việt Tìm kiếm: Lĩnh vực/Mảng Công sởDược phẩmMáy tínhMỹ phẩmQuà tặngSiêu thịThế thaoThời trangThực phẩmTruyền thôngXây dựngĐiện thoại Tỉnh/Thành phố An GiangBà Rịa - Vũng TàuBắc GiangBắc KạnBạc LiêuBắc NinhBến TreBình ĐịnhBình DươngBình PhướcBình ThuậnCà MauCần ThơCao BằngĐà NẵngĐắk LắkĐắk NôngĐiện BiênĐồng NaiĐồng ThápGia LaiHà GiangHà NamHà NộiHà TĩnhHải DươngHải PhòngHậu GiangHòa BìnhHưng YênKhánh HòaKiên GiangKon TumLai ChâuLâm ĐồngLạng SơnLào CaiLong AnNam ĐịnhNghệ AnNinh BìnhNinh ThuậnPhú ThọPhú YênQuảng BìnhQuảng NamQuảng NgãiQuảng NinhQuảng TrịSóc TrăngSơn LaTây NinhThái BìnhThái NguyênThanh HóaThừa Thiên HuếTiền GiangTP.Hồ Chí MinhTrà VinhTuyên QuangVĩnh LongVĩnh PhúcYên Bái Diện tích Dưới 100 m2100 m2 - 200 m2200 m2 - 500 m2500m2 - 1000m2Trên 1000 m2 Vật liệu/Phụ kiện MikaInoxKính màuBaambooMDF woodMelamine An Cuong WoodPine WoodRussian oak woodVerneer wood Tìm kiếm Giới thiệu Nội thất Shop Nội thất Shop là doanh nghiệp cung cấp dịch vụ thiết kế và thi công nội thất cửa hàng. showroom chuyên nghiệp và uy tín. Thành lập từ năm 2010, đến nay Nội thất Shop đã có hơn 5 năm hoạt động trong lĩnh vực trang trí nội thất, công ty đã có những bước tiến vững chắc, trở thành một điểm đến đáng tin cậy. Xem thêm Dịch vụ Nội thất Shop Thiết kế và thi công nội thất Nội Thất Shop đơn vị thiết kế - thi công nội ... Thi công , sản xuất nội thất Đội thợ thi công tay nghề cao với nhiều năm kinh ... Kệ tủ trưng bày điện thoại Thiết kế tủ kệ trưng bày điện thoại anh Toàn - Hà Đông MSP: Giá SP:4,000,000 đ Giá SP:4,000,000 đ Kệ thiết kế tủ trưng bày điện thoại công ty Semco MSP: Giá SP:8,000,000 đ Giá SP:8,000,000 đ Thiết kế kệ tủ trưng bày điện thoại Anh Hùng- Long Biên SKK 301 MSP: SKK301 Giá SP:6,600,000 đ Giá SP:6,000,000 đ Thiết kế kệ trưng bày điện thoại SKP309 MSP: SKP309 Giá SP:7,700,000 đ Giá SP:7,000,000 đ Tủ Kệ Trưng Bày Điện Thoại SKP308 MSP: SKP308 Giá SP:0 đ Giá SP:0 đ Tủ trưng bày điện thoại mã SKP307 MSP: SKP307 Giá SP:0 đ Giá SP:0 đ Tủ trưng bày điện thoại mã SKP306 MSP: SKP306 Giá SP:0 đ Giá SP:0 đ Mẫu bàn trải nghiệm điện thoại mã SKP305 MSP: SKP305 Giá SP:0 đ Giá SP:0 đ Kệ tủ trưng bày mỹ phẩm Thiết kế kệ tủ trưng bày mỹ phẩm anh Huy Giải Phóng MSP: Nghiệm thu sản phẩm xe đẩy mỹ phẩm MSP: Kệ trưng bày mỹ phẩm SKM310 MSP: SKM310 Giá SP:0 đ Giá SP:0 đ Tủ trưng bày mỹ phẩm SKM309 MSP: SKM309 Giá SP:0 đ Giá SP:0 đ Tủ trưng bày mỹ phẩm SKM308 MSP: SKM308 Giá SP:0 đ Giá SP:0 đ Tủ trưng bày mỹ phẩm SKM307 MSP: SKM307 Giá SP:0 đ Giá SP:0 đ Kệ trưng bày mỹ phẩm SKM306 MSP: SKM306 Giá SP:0 đ Giá SP:0 đ Tủ trưng bày mỹ phẩm SKM305 MSP: SKM305 Giá SP:0 đ Giá SP:0 đ Kệ trưng bày trái cây Thiết kế kệ trưng bày hoa quả - chị Ánh MSP: Thiết kế tủ kệ gỗ thông chị Vân Anh - Hà Nội MSP: Kệ trưng bày trái cây, hoa quả SKT182 MSP: SKT182 Giá SP:0 đ Kệ trưng bày trái cây, hoa quả SKT181 MSP: SKT181 Giá SP:0 đ Kệ trưng bày trái cây, hoa quả SKT174 MSP: SKT174 Giá SP:0 đ Kệ Trưng Bày Trái Cây SKT173 MSP: SKT173 Giá SP:0 đ Kệ tủ trưng bày trái cây, hoa quả SKT293 MSP: SKT293 Giá SP:0 đ Kệ trưng bày trái cây SKT178 MSP: SKT178 Giá SP:0 đ Sản phẩm bán chạy Đơn hàng thiết kế bàn làm việc công ty U&I Trường Chinh MSP: Thi công quầy bán vé Chị Hương-Bắc Ninh MSP: Kệ trưng bày điện thoại Vertu SKP284 MSP: SKP284 Giá SP:0 đ Giá SP:0 đ Mẫu kệ giữa nhà trưng bày điện thoại SKP283 MSP: SKP283 Giá SP:0 đ Giá SP:0 đ Tủ trưng bày quần áo trẻ em SKQ282 MSP: SKQ282 Giá SP:0 đ Giá SP:0 đ Kệ trưng bày quần áo thời trang trẻ em SKQ281 MSP: SKQ281 Giá SP:0 đ Giá SP:0 đ Kệ trưng bày mỹ phẩm SKM280 MSP: SKM280 Giá SP:0 đ Giá SP:0 đ Quầy bar thanh toán SKK279 MSP: SKK279 Giá SP:0 đ Giá SP:0 đ Đối tác Video clip Dự án Thi công nội thất cửa hàng điện thoại anh Vịnh - Hoài Đức Hỗ trợ trực tuyến Hà Nội 0911.133.855 Nghệ An 0912.414.855 Email: phongbanhang@cosp.com.vn Hỗ trợ trực tuyến Hỗ trợ trực tuyến Nhà xưởng Chứng chỉ ISO9001-2015 Máy dán cạnh trong sản xuất đồ nội thất Nhà xưởng của nội thất Shop Máy cưa trượt trong sản xuất nội thất Quy trình sản xuất sản phẩm tại phân xưởng Máy cưa panel phục vụ sản xuất Máy cắt CNC và máy máy cắt Laze phục vụ sản xuất Đóng gói- vận chuyển và trách nhiệm bảo vệ môi trường tại nội thất Shop Giới thiệu xưởng sản xuất tại miền Nam Giới thiệu xưởng sản xuất tại Miền Bắc Tư vấn khách hàng Tổng hợp mẫu kệ điện thoại được yêu thích nhất tháng 6 Top 5 mẫu thiết kế nội thất shop thời trang hot nhất 2017 10m2 có thiết kế được shop quần áo thời trang không? Kệ trưng bày sát tường cho cửa hàng mỹ phẩm So sánh ưu điểm của gỗ MDF và gỗ MFC trong thi công nội thất shop Tin tức nổi bật Lựa chọn loại tủ trưng bày cửa hàng đẹp cho shop thời trang mẹ và bé Sắc màu tuyệt đẹp cho phòng ngủ Nội thất đẹp cho thanh niên chưa vợ Phòng khách ấn tượng cho căn hộ lớn “Bí kíp” làm chậu xương rồng cực xinh từ vải nỉ Ý kiến khách hàng Anh Lâm Văn Ninh Vụ Bản - Lạc Sơn - Hòa Binh "Tôi thấy rất hài lòng với đội ngũ nhân viên giúp tôi có được không gian trưng bày sản phẩm như ý. Thay mặt công ty, tôi xin cảm ơn." Đỗ Văn Tuấn CEO của Kid PLAZA "" Nội thất Nội thất Shop thực sự là một đơn vị uy tín trong tư vấn, thiết kế và thi công nội thất cửa hàng. Chúng tôi sẽ tiếp tục hợp tác với Nội thất Shop trong các dự án mở rộng hệ thống sau này". " Thiết bị vệ sinh Cotto Thiết bị vệ sinh Cotto "" Chúng tôi may mắn vì có cơ hội được làm việc với Nội thất Shop và chúng tôi thực sự hài lòng với chất lượng của công trình thiết kế và thi công. Rât mong sẽ có cơ hội được hợp tác với Shop nhiều hơn "" Ông Lý Xuân Hải Tổng giám đốc ACB " Qua một người bạn, tôi biết đến Công ty Quảng cáo và Thiết kế nội thất Hoàng Gia. Chính Hoàng Gia đã tư vấn và giải quyết tất cả những vấn đề của tôi ." Hỗ trợ khách hàng Hà Nội 0911.133.855 Nghệ An 0912.414.855 Email: phongbanhang@cosp.com.vn Hỗ trợ trực tuyến Hỗ trợ trực tuyến Dịch vụ Thiết kế và thi công nội thất Thi công , sản xuất nội thất Gói thiết kế - thi công nội thất Tin liên quan Xả Kho Nội Thất Sản phẩm Kệ trưng bày Kệ tủ trưng bày điện thoại Kệ tủ trưng bày mỹ phẩm Kệ tủ trưng bày giày dép Kệ trưng bày trái cây Kệ tủ trưng bày túi xách Kệ trưng bày thiết bị vệ sinh Kệ tủ trưng bày trang sức Kệ trưng bày gạch Kệ trưng bày siêu thị Kệ tủ trưng bày quần áo Biển quảng cáo Kiosk Booth quảng cáo Kệ điện máy Phụ kiện cho shop Các mẫu kệ trưng bày khác Online Friday Giang hàng hội chợ Dự án Mẫu dự án bán chạy Cửa hàng thời trang Cửa hàng điện thoại, điện máy, máy tính Cửa hàng mỹ phẩm, làm đẹp, spa Cửa hàng thực phẩm chức năng, thực phẩm sạch Siêu thị mini Showroom rượu Cửa hàng vật liệu xây dựng và thiết bị vệ sinh Gian hàng hội chợ, TTTM Cửa hàng mẹ và bé Cửa hàng nữ trang, vàng bạc Cửa hàng dụng cụ thể thao Fanpage Văn phòng Hà Nội * VPTK: 164 Nguyễn Xiển - Thanh Xuân - Hà Nội * Showroom: 164 Nguyễn Xiển, Thanh Xuân, Hà Nội * Tel : 024.8587.5748; * Hotline: 0946-545-299 Mr.Huỳnh * Email : noithatshop@gmail.com PHÂN XƯỞNG: *Tổ 37, Đông Anh, Hà Nội. *Tel: 024.3839.3366 HỖ TRỢ KHÁCH HÀNG Quy định và hình thức thanh toán Chính sách bảo hành Chính sách đổi trả hàng Chính sách gia hàng và vận chuyển Chính sách bảo mật Chi Nhánh Tp. Hồ Chí Minh * Lầu 3, Phòng 301, Số 158 Đường Điện Biên Phủ , Phường 6 ,Quận 3 ,TP HCM * Tel : 0283.9333.568; * Hotline: 0913.545.299 * Email: noithatshop@gmail.com PHÂN XƯỞNG: * Địa chỉ: 5/51, Ấp 3, Xã Xuân Thới Thượng, Huyện Hóc Môn, Hồ Chí Minh ------------------------------------ CHI NHÁNH NGHỆ AN VPGD: Tòa nhà Minh Phú Plaza Xã Diễn Ngọc, Huyện Diễn Châu, Tỉnh Nghệ An, Việt Nam Hotline: 0912.414.855 Xưởng sản xuất: Đường Bộ Quốc Phòng, sát đền Cuông, Nghệ An chi nhánh đà nẵng * Số 68, Nguyễn Đăng Tuyển, phường Thọ Quang, Q. Sơn Trà, TP.ĐN * Tel: 0511.651.8580 * Hotline: 0912.167.686 * Email: noithatshop@gmail.com PHÂN XƯỞNG: * Số 68, Nguyễn Đăng Tuyển, phường Thọ Quang, Q. Sơn Trà, TP.ĐN * Tel: 0511.651.8580 * Hotline: 0912.167.686 * Email: noithatshop@gmail.com Bản quyền 2015 © Công ty cổ phần Quảng cáo và Trang trí Nội thất Hoàng Gia Mst: 0104913886 Trụ sở: Châu Phong- Liên Hà- Đông Anh- Hà Nội FOLLOW US Facebook Google + Twitter Youtube Linkedin Trang chủ Giới thiệu -- Nhà xưởng -- Giới thiệu về Nội thất Shop Dịch vụ -- Xả Kho Nội Thất -- Tin liên quan -- Gói thiết kế - thi công nội thất -- Thi công , sản xuất nội thất -- Thiết kế và thi công nội thất Sản phẩm -- Giang hàng hội chợ -- Online Friday -- Các mẫu kệ trưng bày khác -- Phụ kiện cho shop -- Kệ điện máy -- Booth quảng cáo -- Kiosk -- Biển quảng cáo -- Kệ tủ trưng bày quần áo -- Kệ trưng bày siêu thị -- Kệ trưng bày gạch -- Kệ tủ trưng bày trang sức -- Kệ trưng bày thiết bị vệ sinh -- Kệ tủ trưng bày túi xách -- Kệ trưng bày trái cây -- Kệ tủ trưng bày giày dép -- Kệ tủ trưng bày mỹ phẩm -- Kệ tủ trưng bày điện thoại -- Kệ trưng bày Dự án -- Cửa hàng dụng cụ thể thao -- Cửa hàng nữ trang, vàng bạc -- Cửa hàng mẹ và bé -- Gian hàng hội chợ, TTTM -- Cửa hàng vật liệu xây dựng và thiết bị vệ sinh -- Showroom rượu -- Siêu thị mini -- Cửa hàng thực phẩm chức năng, thực phẩm sạch -- Cửa hàng mỹ phẩm, làm đẹp, spa -- Cửa hàng điện thoại, điện máy, máy tính -- Cửa hàng thời trang Thiết kế shop thời trang -- Mẫu dự án bán chạy Tư vấn khách hàng -- Hỏi đáp - FAQs -- Hướng dẫn đặt hàng -- Phong thủy Tin tức -- Đối tác -- Ý kiến khách hàng -- Không gian sống Tuyển dụng Liên hệ Chăm sóc khách hàng Khuyến mại Mẫu dự án bán chạy Ý kiến khách hàng Đối tác

Cache-Control: no-store, no-cache, must-revalidate, post-check=0, pre-check=0

P3P: CP="CAO DSP COR CURa ADMa DEVa OUR IND PHY ONL UNI COM NAV INT DEM PRE"


markdown:
PGS.TS Lê Quang Minh, Phó giám đốc ĐH Quốc gia TP.HCM cho biết, sự phát triển kinh tế quá nhanh đang diễn ra tại Việt Nam là một thách thức trong phát triển bền vững. Có không ít cơ hội trong việc sử dụng nguồn năng lượng mới như năng lượng sinh học, gió và năng lượng mặt trời. Vấn đề bài toán môi trường, phát triển bền vững cần phải có sự phối hợp của nhiều ngành, các cấp. Hiện, các bộ ngành Việt Nam chỉ mới hợp tác theo hàng dọc, chưa có sự phối hợp đồng bộ; cần phải có những quy trình hợp tác rõ ràng, học hỏi kinh nghiệm từ những nước phát triển để đảm bảo cho mục tiêu này.
TS. Minh nhấn mạnh, hệ thống pháp luật chưa hoàn chỉnh, còn thụ động đi sau những tình huống gây ô nhiễm môi trường đã và đang xảy ra. Trong số 16 lưu vực sông được điều tra thì 5 lưu vực sông bị ô nhiễm nghiêm trọng, 5 lưu vực có chất lượng khá, còn lại ở mức trung bình. Chúng ta chưa có phần dự báo để có những điều luật hoàn thiện hơn.
Bên cạnh đó, đến năm 2015, DuPont sẽ tăng doanh thu hàng năm lên ít nhất 2 tỷ đô la Mỹ từ những sản phẩm tối ưu nguồn năng lượng và giảm thiểu tối đa quá trình thải khí đốt do nhà kính. DuPont ước tính các sản phẩm này sẽ góp phần giảm ít nhất 40 triệu tấn chất khí carbon dioxide mà khách hàng, đối tác của DuPont có thể thải ra. Bên cạnh đó, với những sản phẩm và dịch vụ an toàn, bảo hộ lao động của mình, DuPont sẽ tăng cường tập trung vào việc bảo vệ con người. Dự kiến từ nay đến năm 2015, công ty sẽ giới thiệu ra thị trường ít nhất 1.000 sản phẩm hoặc dịch vụ mới, giúp con người an toàn hơn trên toàn cầu.
Hồng Vinh – Hải Thanh


markdown:
[Abertamy 17
](/register/soap-pn/abertamy-17?language=cs)
Archiv:
Plzeň
Místo:
Bludná (Irrgang), Hřebečná (Hengstererben), Pstruží (Salmthal), Rybná (Fischbach)
Datace:
1855-1884 +, +i
Index nově foliován; B - Z = fol. 118 - 129
Seznam dětí narozených mimo Hřebečnou = fol. 130
Seznam osob, které v letech 1855 - 1884 zemřely a dožily se nejvyššího věku = poslední fol. ( správně 131 - tímto číslem chybně foliována předsádka )
Nová foliace indexu a přepsání ev. listu 16. 1. 2003

Server: Apache/2.4.7 (Ubuntu)

X-Content-Type-Options: nosniff

X-Generator: Drupal 7 (http://drupal.org)

X-Frame-Options: SAMEORIGIN

Cache-Control: no-cache, must-revalidate

Content-Language: cs