Пластичность волатильности:Приложение:Меры волатильности — различия между версиями
WikiSysop (обсуждение | вклад) |
WikiSysop (обсуждение | вклад) |
||
(не показаны 4 промежуточные версии этого же участника) | |||
Строка 8: | Строка 8: | ||
Для положительной случайной величины <math>\textstyle z>0</math> ширину её распределения можно характеризовать относительной ошибкой <math>\textstyle \sigma_z/\bar{z}</math>, где <math>\textstyle \sigma_z</math>, как обычно, обозначает среднеквадратичное отклонение <math>\textstyle \sigma^2_z=\overline{(z-\bar{z})^2}</math>. | Для положительной случайной величины <math>\textstyle z>0</math> ширину её распределения можно характеризовать относительной ошибкой <math>\textstyle \sigma_z/\bar{z}</math>, где <math>\textstyle \sigma_z</math>, как обычно, обозначает среднеквадратичное отклонение <math>\textstyle \sigma^2_z=\overline{(z-\bar{z})^2}</math>. | ||
− | Заметим, что относительные ширины распределений для <math>\textstyle z</math> и <math>\textstyle z^2</math> различны, поэтому, вообще говоря, существуют различные критерии оптимальности мер измерения волатильности. Например, для вычисления стационарной волатильности обычно используется, усреднение либо квадратов доходностей, либо квадратов амплитуд размаха <ref name="Parkinson1980"/>: | + | Заметим, что относительные ширины распределений для <math>\textstyle z</math> и <math>\textstyle z^2</math> различны, поэтому, вообще говоря, существуют различные критерии оптимальности мер измерения волатильности. Например, для вычисления стационарной волатильности обычно используется, усреднение либо квадратов доходностей, либо квадратов амплитуд размаха |
+ | <ref name="Parkinson1980"> | ||
+ | M. Parkinson, 1980 ''The extreme value method for estimating the variance of the rate of return'', | ||
+ | The Journal of Business}, Vol.53, No.1. | ||
+ | </ref>: | ||
{| 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(34)'''</div> |
|} | |} | ||
Так как в этой статье мы изучаем нестационарный характер волатильности, и для её сглаживания используем нелинейный HP-фильтр, то удобнее усреднять непосредственно волатильности <math>\textstyle \sigma</math>, а не их квадраты, которые, как мы увидим ниже, при малых <math>\textstyle n</math> дают, к тому же, смещённое значение <math>\textstyle \sigma</math>. Тем не менее, говоря о различных мерах волатильности, мы будем вычислять относительную ширину как самой величины, так и её квадрата. | Так как в этой статье мы изучаем нестационарный характер волатильности, и для её сглаживания используем нелинейный HP-фильтр, то удобнее усреднять непосредственно волатильности <math>\textstyle \sigma</math>, а не их квадраты, которые, как мы увидим ниже, при малых <math>\textstyle n</math> дают, к тому же, смещённое значение <math>\textstyle \sigma</math>. Тем не менее, говоря о различных мерах волатильности, мы будем вычислять относительную ширину как самой величины, так и её квадрата. | ||
− | Приведём некоторые известные меры. Базовой будет мера Parkinson (1980) <ref name="Parkinson1980"/>, равная амплитуде размаха <math>\textstyle v_P=a</math>. Garman and Klass (1980) <ref name="Garman1980"/>, в классе аналитических функций по <math>\textstyle h</math>, <math>\textstyle l</math>, <math>\textstyle r</math>, предложили следующую оптимальную комбинацию, которая лучше меры Parkinson: | + | Приведём некоторые известные меры. Базовой будет мера Parkinson (1980) <ref name="Parkinson1980"/>, равная амплитуде размаха <math>\textstyle v_P=a</math>. Garman and Klass (1980) |
+ | <ref name="Garman1980"> | ||
+ | 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. | ||
+ | </ref>, в классе аналитических функций по <math>\textstyle h</math>, <math>\textstyle l</math>, <math>\textstyle r</math>, предложили следующую оптимальную комбинацию, которая лучше меры Parkinson: | ||
{| 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(35)'''</div> |
|} | |} | ||
Строка 30: | Строка 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(36)'''</div> |
|} | |} | ||
Строка 37: | Строка 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(37)'''</div> |
|} | |} | ||
Строка 44: | Строка 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">'''( | + | | <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">'''( | + | | <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> | ||
Строка 64: | Строка 71: | ||
Видно, что коэффициент при <math>\textstyle \mu^2</math> в случае <math>\textstyle \beta=1/2</math> в четыре раза меньше, чем в случае <math>\textstyle \beta=0</math> (<math>\textstyle v_P=a</math>). Соответственно, в четыре раза меньше и зависимость от <math>\textstyle \mu</math>. При <math>\textstyle \beta=2/3</math> (далее <math>\textstyle v_{2/3}</math>) коэффициент при <math>\textstyle \mu^2</math> становится равным нулю и зависимость от <math>\textstyle \mu</math> ослабевает ещё сильнее, хотя полностью она исчезает только для меры Rogers and Satchell. | Видно, что коэффициент при <math>\textstyle \mu^2</math> в случае <math>\textstyle \beta=1/2</math> в четыре раза меньше, чем в случае <math>\textstyle \beta=0</math> (<math>\textstyle v_P=a</math>). Соответственно, в четыре раза меньше и зависимость от <math>\textstyle \mu</math>. При <math>\textstyle \beta=2/3</math> (далее <math>\textstyle v_{2/3}</math>) коэффициент при <math>\textstyle \mu^2</math> становится равным нулю и зависимость от <math>\textstyle \mu</math> ослабевает ещё сильнее, хотя полностью она исчезает только для меры Rogers and Satchell. | ||
− | <math>\textstyle \bullet</math> Сравним теперь статистические параметры различных мер волатильности: : | + | <math>\textstyle \bullet</math> Сравним теперь статистические параметры различных мер волатильности: |
+ | |||
+ | <center>[[File:volat_tbl1h.png]]</center> | ||
+ | |||
+ | Курсивным шрифтом приведены значения Монте-Карло моделирования по 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.
Сравним теперь статистические параметры различных мер волатильности:

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

Таким образом, мера имеет достаточно узкое распределение, и, следовательно, ошибку измерения волатильности. При этом её очевидным преимуществом является простота, по сравнению с мерами и . Кроме этого она не смещена при малых размерах выборки, что существенно при исследовании эффектов нестационарности.
Примчания
- ↑ Перейти к: 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,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.
- Перейти ↑ 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.
Приложение: Броуновское блуждание << | Оглавление | >> Приложение: Моделирование блуждания |
---|