Работа над выработкой индексов активов, разрешенных для инвестирования, несомнено, будет продолжена. Она законодательно вменена уполномоченным на это органам управления фондовым рынком (в рамках Закона РФ «Об инвестировании...» [3]. Продолжится работа и над формированием индексов активов, не разрешенных для инвестирования в них накопительной составляющей трудовых пенсий.
№ пп | Тикер индекса модельного класса | Краткое описание модельного класса |
1 | SBS Rus Govt | Государственные обязательства России |
2 | SBS Rus Muni | Обязательства субъектов РФ |
3 | SBS Rus Corp | Корпоративные облишгации российских эмитентов |
4 | RTSI RUB | Акции российских эмитентов (1-й эшелон) |
5 | AK&M-2 | Акции российских эмитентов (2-й эшелон) |
6 | CBR Rus CD | Банковские депозиты в российских рублях |
7 | TYX RUB Cum | Государственные долгосрочные обязательства США |
8 | Moody AAA RUB Cum | Облигации высоконадежных корпораций США |
9 | S&P500 RUB | Акции крупнейших корпораций США |
10 | FED US CD RUB Cum | Банковские депозиты в долларах США |
11 | USD_RUB | Доллары США на банковских счетах |
12 | BE Gilts RUB Cum | Государственные обязательства европейских стран |
13 | BE CD RUB Cum | Банковские депозиты в европейских странах |
14 | DAX RUB | Акции крупнейших эмитентов Еврозоны |
15 | EURO_RUB | Евро на банковских счетах |
16 | BOJ Japan Govt RUB Cum | Государственные обязательства Японии |
17 | BOJ Japan CD RUB Cum | Банковские депозиты Японии |
18 | Japan Nikkei Equity RUB | Акции крупнейших корпораций Японии |
19 | JPY_RUB | Японские иены на счетах в банках |
20 | GBP_RUB | Английские функты стерлингов на счетах в банках |
21 | MSCI Emerging RUB | Фондовые активы развивающихся стран |
3.2. Нечетко-множественная оценка доходности и риска индексов
Традиционной вероятностной моделью поведения индекса является модель винеровского случайного процесса c постоянными параметрами m (коэффициент сноса, по смыслу – предельная курсовая доходность) и s (коэффикциент диффузии, по смыслу – стандартное уклонение от среднего значения предельной доходности). Аналитическое описание винеровского процесса [115]:
(3.1)где z(t) – стандартный винеровский процесс (броуновское движение, случайное блуждание) с коэффициентом сноса 0 и коэффициентом диффузии 1.
В приращениях запись (3.1) приобретает вид
(3.2)Из (3.1) – (3.2) следует, что доходность, как ее понимает модель винеровского процесса, имеет нормальное распределение с матожиданием m и среднеквадратическим отклонением s. Обозначим плотность этого распределения j(r,m,s), где r – расчетное значение доходности.
Однако, если пронаблюдать фактическое ценовое поведение индексов, то мы увидим, что текущая доходность индексов не колеблется вокруг постоянной случайной величины, но образует динамический тренд. Очень характерным для анализа в этом смысле является интервал 1998-2002 г.г., когда тренд доходности поменял знак, и винеровская модель оказалась абсолютно неадекватной.
Чтобы повысить достоверность оценки доходности и риска индексов, необходимо отказаться от винеровской модели и перейти к нечеткой модели финальной (конечной) доходности следующего вида:
S(t) = S(t0) ´ (1+r(t)´(t-t0)), (3.3)
где t – текущее время, t0 – начальный отсчет времени, S(t) - прогнозный уровень индекса – треугольная нечеткая функция, r(t) – расчетный коридор доходности индекса - треугольная нечеткая функция. В каждый момент t случайная величина r(t) имеет нормальное распределение j(r,m,s) с треугольно-нечеткими параметрами m,s. Подробно такое нормальное распределение описано в Приложении 1 к настоящей монографии.
Оценим треугольные параметры m,s по принципу максимума правдобия. Пусть у нас есть квазистатистика доходностей (r1, …rN) мощности N и соответствующая ей гистограмма (n1,...,nM) мощности M. Для этой квазистатистики мы подбираем двупараметрическое нормальное распределение, руководствуясь критерием правдоподобия
, (3.4)где ri – отвечающее i-му столбцу гистограммы расчетное значение доходности, Dr – уровень дискретизации гистограммы.
Задача (3.4) – это задача нелинейной оптимизации, которое имеет решение
, (3.5)причем m0, s0 – аргументы максимума F(m,s), представляющие собой контрольную точку.
Выберем уровень отсечения F1 < F0 и признаем все вероятностные гипотезы правдоподобными, если соответствующий критерий правдоподобия лежит в диапазоне от F1 до F0. Тогда всем правдоподобным вероятностным гипотезам отвечает множество векторов À’, которое в двумерном фазовом пространстве представляет собой выпуклую область с нелинейными границами.
Впишем в эту область прямоугольник максимальной площади, грани которого сориентированы параллельно фазовым осям. Тогда этот прямоугольник представляет собой усечение À’ и может быть описан набором интервальных диапазонов по каждой компоненте
À’’ = (mmin, mmax; smin, smax) Î À’. (3.6)
Назовем À’’зоной предельного правдоподобия. Разумеется, контрольная точка попадает в эту зону , то есть выполняется
mmin< m0 <mmax, smin < s0 < smax , (3.7)
что вытекает из унимодальности и гладкости функции правдоподобия. Тогда мы можем рассматривать числа m = (mmin, m0, mmax), s = (smin, s0, smax) как треугольные нечеткие параметры плотности распределения j(·), которая и сама в этом случае имеет вид нечеткой функции.
Рассмотрим пример. Пусть по результатам наблюдений за индексом сформирована квазистатистика мощностью N=100 отсчетов, представленная в диапазоне –5 ¸ +15 процентов годовых следующей гистограммой c уровнем дискретизации 2% годовых мощностью M=10 интервалов (таблица 3.2):
Таблица 3.2. Гистограмма квазистатистики
Расчетная доходность ri, % годовых (середина интервала) | Число попавших в интервал отсчетов квазистатистики ni | Частота ni = ni/N |
-4 | 5 | 0.05 |
-2 | 2 | 0.02 |
0 | 3 | 0.03 |
2 | 8 | 0.08 |
4 | 10 | 0.1 |
6 | 20 | 0.2 |
8 | 28 | 0.28 |
10 | 19 | 0.19 |
12 | 5 | 0.05 |
14 | 0 | 0 |
Оценить параметры нормального распределения доходности.
Таблица 3.3. Гистограмма квазистатистики
m | F(m,s) ´ 10000 при s = | ||||
2 | 2.5 | 3 | 3.5 | 4 | |
6 | -214 | -120 | -79 | -66 | -67 |
6.5 | -151 | -76 | -49 | -45 | -52 |
7 | -104 | -46 | -29 | -32 | -44 |
7.5 | -77 | -31 | -22 | -29 | -43 |
8 | -76 | -34 | -28 | -36 | -49 |
8.5 | -100 | -56 | -47 | -52 | -62 |
Видно, что при данном уровне дискретизации параметров можно построить зону предельного правдоподобия двумя путями:
À’’1 = (7.5,8.0; 2.5,3.5), À’’2 = (7.0,8.0; 3.0,3.5), (3.8)
причем контрольная точка попадает в оба эти прямоугольника. Точное же решение этой задачи, разумеется, единственное:
À’’ = (6.8,8.3; 2.3,3.8), (3.9)
и m = (6.8, 7.55, 8.3), s = (2.3, 2.95, 3.8) – искомая нечеткая оценка параметров распределения.
Теперь, когда мы научились получать достоверные оценки доходности и риска фондовых индексов, можно переходить к решению задачи оптимизации портфеля на модельных активах.
3.3. Нечетко-множественная оптимизация модельного портфеля
Исторически первым методом оптимизации фондового портфеля был метод, предложенный Гарри Марковицем в [134]. Суть его в следующем.
Пусть портфель содержит N типов ценных бумаг (ЦБ), каждая из которых характеризуется пятью параметрами:
- начальной ценой Wi0 одной бумаги перед помещением ее в портфель;
- числом бумаг ni в портфеле;
- начальными инвестициями Si0 в данный портфельный сегмент, причем
Si0 = Wi0 ´ ni; (3.10)
- среднеожидаемой доходностью бумаги ri;
- ее стандартным отклонением si от значения ri.
Из перечисленных условий ясно, что случайная величина доходности бумаги имеет нормальное распределение с первым начальным моментом ri и вторым центральным моментом si. Это распределение не обязательно должно быть нормальным, но из условий винеровского случайного процесса нормальность вытекает автоматически.