Пластичность волатильности:Приложение:Меры волатильности — различия между версиями

Материал из synset
Перейти к: навигация, поиск
 
(не показана 1 промежуточная версия этого же участника)
Строка 16: Строка 16:
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> \sigma^2_R = \frac{1}{n-1}\sum^n_{t=1} (r_t-\bar{r})^2,\;\;\;\;\;\;\;\;\;\;\;\sigma^2_P = \frac{1}{n}\sum^n_{t=1} \frac{a_t^2}{4\ln 2}. </math>
 
  | width="90%" align="center"|<math> \sigma^2_R = \frac{1}{n-1}\sum^n_{t=1} (r_t-\bar{r})^2,\;\;\;\;\;\;\;\;\;\;\;\sigma^2_P = \frac{1}{n}\sum^n_{t=1} \frac{a_t^2}{4\ln 2}. </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(34)'''</div>
 
  |}
 
  |}
  
Строка 28: Строка 28:
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> v^2_{GK} = 0.511 \cdot a^2 - 0.019 \bigl(r\cdot (h-l)+ 2 h\cdot l\bigr) - 0.383 \cdot r^2. </math>
 
  | width="90%" align="center"|<math> v^2_{GK} = 0.511 \cdot a^2 - 0.019 \bigl(r\cdot (h-l)+ 2 h\cdot l\bigr) - 0.383 \cdot r^2. </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(35)'''</div>
 
  |}
 
  |}
  
Строка 37: Строка 37:
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> v^2_{RS} = h\cdot (h-r) + l\cdot (l+r). </math>
 
  | width="90%" align="center"|<math> v^2_{RS} = h\cdot (h-r) + l\cdot (l+r). </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(36)'''</div>
 
  |}
 
  |}
  
Строка 44: Строка 44:
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> v_\beta=a - \beta\cdot |r| </math>
 
  | width="90%" align="center"|<math> v_\beta=a - \beta\cdot |r| </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(37)'''</div>
 
  |}
 
  |}
  
Строка 51: Строка 51:
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> \frac{\overline{(a-\beta \cdot |r|)^2}}{\bigl(\;\bar{a}-\beta \cdot \overline{|r|}\;\bigr)^2} = min \;\;\;\;\;\;\;=>\;\;\;\;\beta = 6-8 \ln 2 \approx 0.455. </math>
 
  | width="90%" align="center"|<math> \frac{\overline{(a-\beta \cdot |r|)^2}}{\bigl(\;\bar{a}-\beta \cdot \overline{|r|}\;\bigr)^2} = min \;\;\;\;\;\;\;=>\;\;\;\;\beta = 6-8 \ln 2 \approx 0.455. </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(38)'''</div>
 
  |}
 
  |}
  
 
Однако, так как критерий <math>\textstyle \sigma_v/\bar{v}</math> не является единственным, и в силу слабой чувствительности относительной волатильности от <math>\textstyle \beta</math>, мы в статье использовали значение <math>\textstyle \beta=1/2</math> и обозначение <math>\textstyle v=a-|r|/2</math>. Далее <math>\textstyle v_\beta=a-0.455 \cdot|r|</math>.
 
Однако, так как критерий <math>\textstyle \sigma_v/\bar{v}</math> не является единственным, и в силу слабой чувствительности относительной волатильности от <math>\textstyle \beta</math>, мы в статье использовали значение <math>\textstyle \beta=1/2</math> и обозначение <math>\textstyle v=a-|r|/2</math>. Далее <math>\textstyle v_\beta=a-0.455 \cdot|r|</math>.
  
Заметим, что существует ещё одна простая мера волатильности, сравнимая по эффективности с (), следующего вида:
+
Заметим, что существует ещё одна простая мера волатильности, сравнимая по эффективности с (37), следующего вида:
  
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> v_F = \frac{a}{1+r^2/a^2}. </math>
 
  | width="90%" align="center"|<math> v_F = \frac{a}{1+r^2/a^2}. </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(39)'''</div>
 
  |}
 
  |}
  
Хотя вероятность нулевого значения <math>\textstyle a</math> для конечной длительности лага <math>\textstyle T</math> исчезающе мала, необходимо, тем не менее, доопределить <math>\textstyle v_F=0</math> при <math>\textstyle a=0</math>. Вообще, () и () не являются аналитическими функциями по <math>\textstyle a</math> и <math>\textstyle r</math> и выпадают из действия леммы приложения B. работы <ref name="Garman1980"/>.
+
Хотя вероятность нулевого значения <math>\textstyle a</math> для конечной длительности лага <math>\textstyle T</math> исчезающе мала, необходимо, тем не менее, доопределить <math>\textstyle v_F=0</math> при <math>\textstyle a=0</math>. Вообще, (37) и (39) не являются аналитическими функциями по <math>\textstyle a</math> и <math>\textstyle r</math> и выпадают из действия леммы приложения B. работы <ref name="Garman1980"/>.
  
<math>\textstyle \bullet</math> Кроме ширины распределения, в качестве критерия иногда используется отсутствие или слабая зависимость от сноса <math>\textstyle \mu</math>. Заметим, что для дневных или более коротких лагов <math>\textstyle \mu \ll \sigma</math>. Поэтому этот критерий не является столь значимым. Предложенная выше мера модифицированной амплитуды размаха, как и сама амплитуда, зависит от <math>\textstyle \mu</math>. Однако эта зависимость существенно слабее, чем у амплитуды. Если воспользоваться разложениями (), (), для <math>\textstyle v_\beta</math> можно записать:
+
<math>\textstyle \bullet</math> Кроме ширины распределения, в качестве критерия иногда используется отсутствие или слабая зависимость от сноса <math>\textstyle \mu</math>. Заметим, что для дневных или более коротких лагов <math>\textstyle \mu \ll \sigma</math>. Поэтому этот критерий не является столь значимым. Предложенная выше мера модифицированной амплитуды размаха, как и сама амплитуда, зависит от <math>\textstyle \mu</math>. Однако эта зависимость существенно слабее, чем у амплитуды. Если воспользоваться разложениями (32), (33), для <math>\textstyle v_\beta</math> можно записать:
  
 
:<center><math>\overline{v}_\beta = (2-\beta)\cdot \sqrt{\frac{2}{\pi}} + \frac{(2-3\beta)\mu^2}{3\sqrt{2\pi}} - \frac{(2-5\beta)\mu^4}{60\sqrt{2\pi}} +...\;\;\;\;\;</math></center>
 
:<center><math>\overline{v}_\beta = (2-\beta)\cdot \sqrt{\frac{2}{\pi}} + \frac{(2-3\beta)\mu^2}{3\sqrt{2\pi}} - \frac{(2-5\beta)\mu^4}{60\sqrt{2\pi}} +...\;\;\;\;\;</math></center>
Строка 77: Строка 77:
 
Курсивным шрифтом приведены значения Монте-Карло моделирования по 3.5 миллионам лагов, каждый из которых являлся блужданием из миллиона тиков. В этом случае для средних и волатильности в последнем знаке возможна ошибка порядка <math>\textstyle \pm 0.002</math>. Для определения остальных значений (не курсивных) использовались аналитические выражения.
 
Курсивным шрифтом приведены значения Монте-Карло моделирования по 3.5 миллионам лагов, каждый из которых являлся блужданием из миллиона тиков. В этом случае для средних и волатильности в последнем знаке возможна ошибка порядка <math>\textstyle \pm 0.002</math>. Для определения остальных значений (не курсивных) использовались аналитические выражения.
  
<math>\textstyle \bullet</math> В условиях нестационарности данных часто необходимо проводить усреднение по достаточно малому числу наблюдений <math>\textstyle n</math>. В этом случае начинает проявляться смещённость квадратичных мер для оценки волатильности <math>\textstyle \sigma</math>. Даже, если мы вычисляем классический квадрат волатильности <math>\textstyle \sigma^2_R</math> по несмещённой формуле (), значение <math>\textstyle \sigma_R</math> будет смещено, так как при усреднении по большому числу выборок, каждая из которых имеет размер <math>\textstyle n</math>, мы имеем: <math>\textstyle <\sigma^2_R>=\sigma^2</math>, но <math>\textstyle <\sqrt{\sigma^2_R}> \neq\sigma</math>. Если нас интересует значение именно волатильности, а не её квадрата, лучше для нестационарных данных использовать линейные, а не квадратичные меры.
+
<math>\textstyle \bullet</math> В условиях нестационарности данных часто необходимо проводить усреднение по достаточно малому числу наблюдений <math>\textstyle n</math>. В этом случае начинает проявляться смещённость квадратичных мер для оценки волатильности <math>\textstyle \sigma</math>. Даже, если мы вычисляем классический квадрат волатильности <math>\textstyle \sigma^2_R</math> по несмещённой формуле (34), значение <math>\textstyle \sigma_R</math> будет смещено, так как при усреднении по большому числу выборок, каждая из которых имеет размер <math>\textstyle n</math>, мы имеем: <math>\textstyle <\sigma^2_R>=\sigma^2</math>, но <math>\textstyle <\sqrt{\sigma^2_R}> \neq\sigma</math>. Если нас интересует значение именно волатильности, а не её квадрата, лучше для нестационарных данных использовать линейные, а не квадратичные меры.
  
Для иллюстрации эффекта смещения справа приведены графики средних значений волатильности, полученных усреднением по большому числу выборок по <math>\textstyle n</math> значений в каждой для стандартного определения <math>\textstyle \sigma_R</math> и меры <math>\textstyle \sigma_{RG}=\sqrt{v^2_{RG}}</math> () по сравнению с линейной мерой <math>\textstyle \sigma=(a-|r|/2)\sqrt{2\pi}/3</math>.  
+
Для иллюстрации эффекта смещения справа приведены графики средних значений волатильности, полученных усреднением по большому числу выборок по <math>\textstyle n</math> значений в каждой для стандартного определения <math>\textstyle \sigma_R</math> и меры <math>\textstyle \sigma_{RG}=\sqrt{v^2_{RG}}</math> (35) по сравнению с линейной мерой <math>\textstyle \sigma=(a-|r|/2)\sqrt{2\pi}/3</math>.  
  
 
<center>[[File:volat_tbl.png]]</center>
 
<center>[[File:volat_tbl.png]]</center>

Текущая версия на 20:42, 6 марта 2010

Приложение: Броуновское блуждание << Оглавление >> Приложение: Моделирование блуждания

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

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

(34)

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

Приведём некоторые известные меры. Базовой будет мера Parkinson (1980) [1], равная амплитуде размаха . Garman and Klass (1980) [2], в классе аналитических функций по , , , предложили следующую оптимальную комбинацию, которая лучше меры Parkinson:

(35)

Отметим также более простую и не зависящую от сноса меру Rogers and Satchell (1991) [3]:

(36)

Покажем, что простейшая линейная модификация меры Parkinson

(37)

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

(38)

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

Заметим, что существует ещё одна простая мера волатильности, сравнимая по эффективности с (37), следующего вида:

(39)

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

Кроме ширины распределения, в качестве критерия иногда используется отсутствие или слабая зависимость от сноса . Заметим, что для дневных или более коротких лагов . Поэтому этот критерий не является столь значимым. Предложенная выше мера модифицированной амплитуды размаха, как и сама амплитуда, зависит от . Однако эта зависимость существенно слабее, чем у амплитуды. Если воспользоваться разложениями (32), (33), для можно записать:

Видно, что коэффициент при в случае в четыре раза меньше, чем в случае (). Соответственно, в четыре раза меньше и зависимость от . При (далее ) коэффициент при становится равным нулю и зависимость от ослабевает ещё сильнее, хотя полностью она исчезает только для меры Rogers and Satchell.

Сравним теперь статистические параметры различных мер волатильности:

Volat tbl1h.png

Курсивным шрифтом приведены значения Монте-Карло моделирования по 3.5 миллионам лагов, каждый из которых являлся блужданием из миллиона тиков. В этом случае для средних и волатильности в последнем знаке возможна ошибка порядка . Для определения остальных значений (не курсивных) использовались аналитические выражения.

В условиях нестационарности данных часто необходимо проводить усреднение по достаточно малому числу наблюдений . В этом случае начинает проявляться смещённость квадратичных мер для оценки волатильности . Даже, если мы вычисляем классический квадрат волатильности по несмещённой формуле (34), значение будет смещено, так как при усреднении по большому числу выборок, каждая из которых имеет размер , мы имеем: , но . Если нас интересует значение именно волатильности, а не её квадрата, лучше для нестационарных данных использовать линейные, а не квадратичные меры.

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

Volat tbl.png

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

Примчания

  1. Перейти к: 1,0 1,1 M. Parkinson, 1980 The extreme value method for estimating the variance of the rate of return, The Journal of Business}, Vol.53, No.1.
  2. Перейти к: 2,0 2,1 M.B. Garman, M.J.Klass, 1980, On the estimation of security price volatilities from historical data, The Journal of Business, Vol.53, No.1.
  3. Перейти L.C.G. Rogers, S.E.Satchell, 1991, {Estimating variance from high, low and closing prices}, The Annals of Applied Probability, Vol.1. No.4, pp.504-512.

Приложение: Броуновское блуждание << Оглавление >> Приложение: Моделирование блуждания