Zagadnienia

15. Metody stochastyczne w finansach – cd

15.1. Value at Risk – inne spojrzenie na ryzyko finansowe

15.1.1. Wprowadzenie i definicja

Value at Risk (wartość zagrożona ryzykiem, w skrócie VaR) jest przykładem miary ryzyka opartej na nieco innej filozofii niż miary omówione w wykładzie 14. VaR, to miara, za pomocą której nie tyle określa się, na ile niedokładne są prognozy końcowego wyniku inwestycji, ale uzyskuje informacje jak duże środki należy przeznaczyć na zabezpieczenie danej inwestycji, a nawet całej instytucji finansowej. Krótko mówiąc VaR uprzedza nas, ile możemy stracić.

Miary takie jak VaR wykorzystuje się m.in.:

  1. jako element kontroli, np.:
    – poszczególnych departamentów i oddziałów przez zarząd banku,
    – poszczególnych banków przez nadzór bankowy,
    – zarządu spółki akcyjnej przez radę nadzorczą i akcjonariuszy;

  2. przy ustalaniu limitów dotyczących adekwatności kapitałowej instytucji finansowych, czyli określaniu ,,jak duże powinny być środki własne firmy, aby zabezpieczona była jej wypłacalność”;

  3. przy ustalaniu wysokości depozytów wymaganych przez izby rozrachunkowe działające na rynku instrumentów pochodnych;

  4. jako benchmark przy podejmowaniu decyzji inwestycyjnych, np.:
    ,,Ustala się pewną wartość graniczną d i zarządzający portfelem inwestycyjnym może tylko tak kształtować skład portfela, aby dana miara ryzyka jej nie przekroczyła”.

Ogólnie, Value at Risk to odpowiedź probabilisty na pytanie ,,ile można stracić?”. Jest to największa strata, jaką można ponieść przy zadanym poziomie ufności c.

Ustalmy pewien horyzont czasowy T. Niech X będzie zmienną losową, za pomocą której modelujemy przyszły stan należności netto danej instytucji finansowej (pasywa minus aktywa) albo stratę z danej inwestycji (aktualna wartość inwestycji minus jej wartość po upływie czasu T). Jeśli X ma rozkład ciągły i jej dystybuanta jest ściśle rosnąca (np. X ma rozkład normalny), to VaRcX określa się w następujący sposób

PXVaRcX=c.

W przypadku ogólnym (np. gdy zmienna losowa X jest dyskretna) powyższe równanie może nie mieć rozwiązań lub może mieć ich więcej niż jedno. Dlatego też zdefiniujemy VaR jako tzw. dolny kwantyl.

Definicja 15.1

Dla dowolnej zmiennej losowej X i dowolnego poziomu ufności c0,1

VaRcX=supV:PXV<c,

Oczywiście można też definiować VaR jako górny kwantyl

VaRcuX=supV:PX<Vc,

lub jako średnią ważoną obu kwantyli.

Kwestią otwartą pozostaje wybór właściwego poziomu ufności c. Analitycy banku J.P.Morgan, twórcy pojęcia Value at Risk, przyjmowali c=95% ([37]). Natomiast Komitet Bazylejski postuluje przyjęcie c=99% ([5]). Ta różnica wynika między innymi z faktu, że według analityków z banku J.P.Morgan Value at Risk miało być narzędziem kontrolnym używanym przez zarząd banku, natomiast Komitet Bazylejski zaproponował użycie Value at Risk w ramach nadzoru bankowego. Oczywiste jest, że nadzór zewnętrzny jest bardziej wymagający (większe c oznacza większą wartość VaR).

Omówimy pokrótce podstawowe własności VaR. Zauważmy, że oba zbiory V:PX<Vc i V:PXV<c są półprostymi, wobec tego kres górny każdego z nich jest równy kresowy dolnemu odpowiednio każdego z dopełnień. Zatem VaRcuX i VaRcX są równe odpowiednio

VaRcuX=infV:PX<V>c,VaRcX=infV:PXVc.

Korzystając z powyższego, udowodnimy następujący lemat.

Lemat 15.1

Dla dowolnej zmiennej losowej X i dowolnego poziomu ufności c0,1:
1. PX<VaRcuXcPXVaRcX;
2. Jeśli x jest większy od VaRcuX, to PX<x>c;
3. Jeśli x jest mniejszy od VaRcX, to PXx<c.
Ponadto, gdy VaRcuX nie jest równy VaRcX, to:
4. PX<VaRcuX=PXVaRcX=c;
5. Dla dowolnego x z przedziału otwartego VaRcX,VaRcuX,
PX<x=PXx=c;
6. PXVaRcX,VaRcuX=0.

Dowód.
Punkt 1 wynika z ciągłości prawdopodobieństwa ([21] §Twierdzenie 7).

ω:Xω<VaRcuX=n=1ω:Xω<VaRcuX-1n,

zatem, ponieważ PX<VaRcuX-1nc, to

PX<VaRcuX=limnPX<VaRcuX-1nc.

Podobnie

ω:XωVaRcX=n=1ω:XωVaRcX+1n.

Ponieważ PXVaRcuX+1nc, to

PXVaRcX=limnPXVaRcX+1nc.

Punkt 2 wynika z przedstawienia VaRcuX jako kres dolny, a punkt 3 z przedstawienia VaRcX jako kres górny.

Punkt 4 wynika z nierówności z punktu 1 i z faktu, że VaRcX jest mniejszy od VaRcuX.

cPXVaRcXPX<VaRcuXc,

wobec tego wszystkie nierówności ,,” w tym wzorze są równościami.

Podobnie pokazuje się punkt 5:

c=PXVaRcXPX<xPXxPX<VaRcuX=c.

Punkt 6 wynika z punktu 4.

PXVaRcX,VaRcuX=
=P(X<VaRcu(X))-P(XVaRc(X))=c-c=0.

15.1.2. Alokacja kapitału

Na zakończenie przeanalizujmy możliwość zastosowania Value at Risk do wyznaczania optymalnej alokacji posiadanych środków. Rozważymy następujący prosty przykład.

Udziałowcy Towarzystwa Ubezpieczeniowego TU SA zastanawiają się, czy zachodzi potrzeba dokapitalizowania spółki. TU SA sprzedało szereg polis ubezpieczeniowych co może skutkowac wypłatą wysokich odszkodowań. Jeśli kwota odszkodowań przekroczy posiadane fundusze, TU SA będzie zagrożone bankructwem. Aby go uniknąć udziałowcy będą zmuszeni do poniesienia dodatkowych kosztów. Z drugiej strony dokapitalizowanie też wiąże się z pewnymi kosztami. W tej sytuacji udziałowcy muszą rozstrzygnąć, jaki jest optymalny poziom funduszy posiadanych przez TU SA.

Przyjmiemy następujące oznaczenia:
Y – wielkość odszkodowań pomniejszoną o środki posiadane przez TU SA (modelujemy ją jako zmienną losową);
x – kwota, którą udziałowcy chcą dokapitalizować spółkę;
Kx,Y – koszty w zależności od x i Y.

Kx,Y=ax+bY-x+,

gdzie: a – koszt kapitału (np. stopa procentowa), b – koszty, które poniosą udziałowcy na pokrycie roszczeń przekraczających fundusze TU SA (w przeliczeniu na 1 jednostkę monetarną). Oczywiście 0<a<b.

Racjonalni akcjonariusze będą się starać zminimalizować oczekiwane koszty. A więc wybiorą x, które zminimalizuje wartość oczekiwaną relatywnej straty, czyli różnicy Kx,Y-K0,Y,

EKx,Y-K0,Ymin.

Pokażemy, że optymalny poziom dokapitalizowania jest równy Value at Risk Y, wyznaczonemu dla poziomu ufności c=1-ab.

Zauważmy, że w odróżnieniu od straty Kx,Y, relatywna strata

Vx=Kx,Y-K0,Y=ax+bY-x+-bY+,

jest ograniczona i niezależnie od tego, jaki jest rozkład Y i ile wynosi x, posiada skończoną wartość oczekiwaną, która jest funkcją ciągłą zmiennej x. Ponadto, gdy Kx,Y mają skończoną wartość oczekiwaną, to EVx i EKx,Y przyjmują minimum w tych samych punktach.

Twierdzenie 15.1

VaRcY minimalizuje wartość oczekiwaną relatywnej straty EVx dla c=1-ab. Ponadto, jeśli dystrybuanta Y jest ściśle rosnąca w pewnym otoczeniu punktu x=VaRcY, to jest to jedyne minimum. W przeciwnym wypadku zbiór minimów EVx jest przedziałem domkniętym

VaRcY,VaRcuY.

Dowód.
Niech x1 i x2, x1<x2, będą dowolnymi liczbami rzeczywistymi. Zapiszemy różnicę relatywnych strat za pomocą funkcji charakterystycznych I

Vx2-Vx1=ax2+bY-x2+-Y+-ax1-bY-x1++Y+=
=a(x2-x1)-b(x2-x1)IYx2-b(Y-x1)IYx1,x2=
=ax2-x1-bx2-x1IY>x1+bx2-YIYx1,x2.

Zatem przyrost wartości oczekiwanej relatywnej straty wynosi

EVx2-EVx1=EVx2-Vx1=
=(x2-x1)(a-bP(Yx2))-bE((Y-x1)IYx1,x2)=
=x2-x1a-bPY>x1+bEx2-YIYx1,x2.

Rozważymy trzy przypadki:
A. Oba punkty leżą na prawo od VaRcu,

VaRcuY<x1<x2;

B. Oba punkty leżą na lewo od VaRc,

x1<x2<VaRc;

C. Oba punkty leżą pomiędzy VaRc i VaRcu,

VaRcY<x1<x2<VaRcuY.

Przypadek A.
Na mocy lematu 5.4.1

PY>x1=1-PYx1<1-c=ab.

Wobec tego

a-bPY>x1>a-bab=0.

Zatem, ponieważ wartość oczekiwana nieujemnej zmiennej losowej jest nieujemna, to

E(V(x2))-E(V(x1))=(x2-x1)(a-bP(Y>x1))+bE((x2-Y)IYx1,x2)0.

W ten sposób pokazaliśmy, że funkcja EVx jest ściśle rosnąca na półprostej
VaRcuY,+. Ale wiemy, że jest ona ciągła, zatem jest ściśle rosnąca na półprostej domkniętej VaRcuY,+.

Przypadek B.
Na mocy lematu 5.4.1

PYx2=1-PY<x2>1-c=ab.

Wobec tego

a-bPYx2<a-bab=0.

Zatem, ponieważ wartość oczekiwana nieujemnej zmiennej losowej jest nieujemna, to

E(V(x2))-E(V(x1))=(x2-x1)(a-bP(Yx2))-bE((Y-x1)IYx1,x2)0.

W ten sposób pokazaliśmy, że funkcja EVx jest ściśle malejąca na półprostej -,VaRcY. Z ciągłości wynika, że jest ona ściśle malejąca na półprostej domkniętej -,VaRcY.

Przypadek C.
Gdy VaRcY jest mniejszy od VaRcuY, to na mocy lematu 5.4.1

PYx2=1-PY<x2=1-c=ab.

Wobec tego

a-bPYx2=a-bab=0.

Ponadto

PYx1,x2=0,

czyli

EY-x1IYx1,x2=0.

Z powyższego wynika, że

E(V(x2))-E(V(x1))=(x2-x1)(a-bP(Yx2))-bE((Y-x1)IYx1,x2)=0.

W ten sposób pokazaliśmy, że funkcja EVx jest stała na odcinku
VaRcY,VaRcuY. Z ciągłości wynika, że jest ona stała na odcinku domkniętym VaRcY,VaRcuY.

Podsumowując, gdy dystrybuanta Y jest ściśle monotoniczna w pewnym otoczeniu VaRcY, to VaRcuY=VaRcY i EVx przyjmuje minimum dokładnie w jednym punkcie. W przeciwnym przypadku każdy punkt z przedziału VaRcY,VaRcuY minimalizuje EVx.

15.2. Teoretyczne podstawy pomiaru ryzyka dla inwestycji finansowych.

Każdą inwestycję A charakteryzujemy za pomocą dwóch wielkości:
wartości początkowej (nakładów) ozn. W0A
i wartości końcowej (wypłaty, bogactwa końcowego) ozn. W1A.
W0A>0 i jest znane w momencie rozpoczęcia inwestycji.
W1A0 modelujemy jako zmienną losową określoną na pewnej przestrzeni probabilistycznej Ω,F,P.
W tej konwencji
W1A-W0A oznacza zysk,
W0A-W1A stratę,
a stopa zwrotu to RA=W1A-W0AW0A.

Ponadto ze względów technicznych zakładamy, że zbiór dopuszczalnych inwestycji G jest stożkiem wypukłym

i.AGλ>0CGWiC=λWiA,i=0,1;
ii.A,BGCGWiC=WiA+WiB,i=0,1;

oraz zawiera inwestycję o zerowym zysku

iii.CGW1C=W0C.

Inwestycje C z punktów i. oraz ii. będziemy oznaczać odpowiednio C=λA i C=A+B.

Z punktu widzenia metodologicznego wyróżniamy dwie grupy miar ryzyka: miary dyspersji i miary monotoniczne. Osobną grupę stanowią miary efektywności, ściśle związane z ryzykiem nieosiągnięcia zakładanego poziomu zysku albo stopy zwrotu.

15.2.1. Miary dyspersji.

Miary dyspersji ostrzegają, o ile stopa zwrotu RA albo wypłata W1A różnią się od prognozy (np. wartości oczekiwanej lub mediany). Ostrzeżenie to może dotyczyć zarówno:
1. Możliwości wystąpienia efektu niezgodnego z przewidywaniami. Nieważne, czy jest to przykra niespodzianka, czy przyjemne zaskoczenie; ważne, że prognoza była niedokładna.
2. Możliwości poniesienia straty. Uwzględniamy tylko przykre niespodzianki.

Najbardziej popularne miary ryzyka-dyspersji to omówione na poprzednim wykładzie (patrz 14.3.1):

1. Odchylenie standardowe i wariancja.

VR=ER-ER2,σR=VR

Wyznaczamy prognozę błędu prognozy.

2. Ujemne semiodchylenie standardowe i ujemna semiwariancja.

SV-R=ER-ER-2,σ-R=SV-R

Są to odpowiedniki odchylenia standardowego i wariancji w przypadku drugiej interpretacji ryzyka.

3. Odchylenie przeciętne i ujemne semiodchylenie przeciętne.

dR=ER-ER,sd-R=ER-ER-.

Więcej przykładów i informacji o miarach dyspersji (dewiacji) czytelnik może znaleźć w pracy [40].

15.2.2. Miary monotoniczne – zgodne z dominacją stochastyczną.

Miary monotoniczne mają za zadanie oceniać, która inwestycja może przynieść większą stratę. Funkcja

ρ:GR+

jest monotoniczą miarą ryzyka jeżeli spełnia następujące warunki

A1.W0AW0BW1AW1BρAρB;
A2.W0A=W1AρA=0.

Poza tym ,,mile widziane” są następujące własności:

A3.λ1,λ20,λ1+λ2=1ρλ1A1+λ2A2λ1ρA1+λ2ρA2;
A4.cRW1B-W0B=W1A-W0A+cρB=ρA-c;
A5.λ>0ρλA=λρA.

Miary spełniające aksjomaty A1, A2, A3 i A4 nazywa się wypukłymi (patrz [14]). Natomiast miary spełniające wszystkie pięć aksjomatów to tzw. koherentne miary ryzyka (patrz [2, 3]). Podstawowym zastosowaniem takich miar jest wyznaczanie wielkości ”kapitału ekonomicznego” potrzebnego do zabezpieczenia danej pozycji.

Przykłady miar monotonicznych.

1. Prawdopodobieństwo nieosiągnięcia wyznaczonego poziomu T (poziomu aspiracji)

SFTA=PW1A-W0A<T.

Gdy T=0 jest to prawdopodobieństwo poniesienia straty.
Gdy T0, to SFT spełnia tylko warunki A1 i A2.

2. Value at Risk - największa strata przy zadanym poziomie istotności α (α=0.05,0.01,).

VaRαA=infV:PW0A-W1A>V<α.

VaRα spełnia warunki A1, A2, A4 i A5, ale na ogół nie spełnia A3 – patrz lemat 15.1.

3. Conditional VaR (zwany też Probable Maximum Loss lub Expected Shortfall) - warunkowa oczekiwana wartość straty pod warunkiem, że osiągnie ona lub przekroczy VaR.

CVaRα(A)=E(W0(A)-W1(A)|W0(A)-W1(A)VaRα(A)).

Dla zmiennych losowych W1() o ciągłym rozkładzie CVaR jest koherentną miarą ryzyka i spełnia wszystkie pięć aksjomatów (patrz [32, 39, 1]).

4. Spadek oczekiwanej użyteczności.

ρφA=φW0A-EφW1A,

gdzie φ wklęsła i rosnąca funkcja na R+ (funkcja użyteczności) np. lnt.
ρφ spełnia warunki A1, A2 i A3.

5. Miary ryzyka indukowane przez zmianę miary probabilistycznej.
Niech Yi, i=1,k, nieujemne zmienne losowe o skończonej wartości oczekiwanej (np. kursy waluty lokalnej w i-tej walucie lub koszyku walut).

ρ(A)=sup{EYiW0A-W1AEYi},i=1,,k}.

Miary indukowane są koherentnymi miarami ryzyka i spełniają wszystkie pięć aksjomatów.

6. Miary wyznaczone przez poziom krytyczny r taki, że wartość oczekiwana jego przekroczenia przez stratę jest stała i wynosi γ, γ>0.

ργA=r+γ, gdzie EW0A-W1A-r+=γ.

Miary ργ są wypukłe i spełniają warunki A1, A2, A3 i A4 ([22]).

7. Miary wyznaczone przez poziom krytyczny r taki, że wartość oczekiwana przekroczenia r jest proporcjonalna do r i wynosi γr, γ>0.

ργA=0,gdyPW1A<W0A=0,r,gdyPW1A<W0A>0,

gdzie

r>0 i EW0A-W1Ar-1+=γ.

Miary ργ spełniają warunki A1, A2, A3 i A5 ([22]).

8. Miary wyznaczone przez poziom krytyczny r taki, że dla ustalonego parametru, γ>0.

rEW0A-W1A i EW0A-W1A-r+=γr-EW0A-W1A.

Miary ργ są koherentne i spełniają warunki A1, A2, A3, A4 i A5 ([22]).

9. Składka holenderska.

ρcA=EW0A-W1A+cEW1A-EW1A-,c0,1.

Miary ργ są koherentne i spełniają warunki A1, A2, A3, A4 i A5.

15.2.3. Miary efektywności

Miary efektywności (performance measures) mają za zadanie ocenić w jakim stopniu dana strategia inwestycyjna gwarantuje uzyskanie (i przekroczenie) zadanego poziomu aspiracji T. Zwykle benchmark T jest wielkością deterministyczną, ale może być też zależny od czynników losowych.

Miary efektywności najczęściej konstruowane są jako iloraz dwóch miar:
miary nadwyżki stopy zwrotu nad benchmarkiem T
i miary dyspersji stopy zwrotu względem benchmarku T.
Można tu zauważyć pewną analogię z testami statystycznymi dotyczącymi wartości oczekiwanej.

Przykłady miar efektywności ([31]).

Przyjmiemy następujące oznaczenia:
R – stopa zwrotu z inwestycji, a T – benchmark czyli wymagany poziom stopy zwrotu (poziom aspiracji).

1. Uogólniony współczynnik Sharpe'a.

λR=ER-TσR-T.

2. Współczynnik Sortino.

SoRR=ER-TσR-T-.

3. Wskaźnik potencjału nadwyżki (Upside Potential Ratio).

UPRR=ER-T+σR-T-.

4. Stosunek zysku do straty (Gain-Loss Ratio), zwany też funkcją omega.

ΩR=ER-T+ER-T-.

Miary 2, 3, i 4, w odróżnieniu od ,,symetrycznego” współczynnika Sharpe'a, określane są w literaturze jako ”Downside Risk Measures”.

15.3. Ćwiczenia

Ćwiczenie 15.1

Wyznacz VaRc dla straty X o rozkładzie normalnym, XNμ,σ2, dla c=0,95 i dla c=0,99.

Rozwiązanie.

VaRcX=FX-1c=μ+σF-1c,

gdzie F jest dystrybuantą standardowego rozkładu normalnego. Zatem VaR0.95=μ+1,65σ a VaR0.99=μ+2,33σ.

Odpowiedź.VaR0.95=μ+1,65σ a VaR0.99=μ+2,33σ.

Ćwiczenie 15.2

Wyznacz VaRc dla straty X o rozkładzie lognormalnym, X=expY YNμ,σ2, dla c0,1.

Rozwiązanie.

VaRcX=expVaRcY=expμ+σF-1c,

gdzie F jest dystrybuantą standardowego rozkładu normalnego.

Odpowiedź.VaRcX=expμ+σF-1c, gdzie F jest dystrybuantą standardowego rozkładu normalnego.

Ćwiczenie 15.3

Wyznacz VaRc dla straty X o rozkładzie Pareto, FXx=1-1xk, x1, k>0 dla c=0,95 i dla c=0,99.

Rozwiązanie.

FX-1y=1-y-1/k.

Zatem

VaR0,95X=FX-10,95=0,05-1/k=201/k,
VaR0,99X=FX-10,99=0,01-1/k=1001/k.

Odpowiedź.VaR0,95X=0,05-1/k=201/k, a VaR0,99X=0,01-1/k=1001/k.

Ćwiczenie 15.4

Przybliżamy łączny rozkład dziennych stóp zwrotu kursów USD i EUR w PLN rozkładem normalnym o parametrach: m1=m2=0, σ1=0,01, σ2=0,02 i ρ=0,5. Aktualny kurs USD wynosi 3,5 PLN, a EUR – 4 PLN. Wyznaczyć dzienny Value at Risk dla portfela walutowego złożonego z 20 tys. USD i 10 tys. EUR, dla c=0,95.

Rozwiązanie. Niech r1 i r2 oznaczają dzienne stopy zwrotu (relatywne przyrosty) kursów USD i EUR. Oznaczmy przez K wyrażoną w PLN, możliwą stratę, tzn. różnicę wartości aktualnej portfela i w dniu następnym.

K=-3,520000r1+410000r2=-70000r1+40000r2.

K ma rozkład normalny o zerowej wartości oczekiwanej i wariancji

σ2K=700002σ12+400002σ22+27000040000ρσ1σ2=1690000.

Odchylenie standardowe K wynosi 1300. Po przybliżeniu kwantyla rozkładu normalnego standardowego przez 1,65 otrzymujemy, że VaRK wynosi 13001,65=2145.

Odpowiedź. Dla poziomu ufności 0,95 Value at Risk wynosi 2145 PLN.

Treść automatycznie generowana z plików źródłowych LaTeXa za pomocą oprogramowania wykorzystującego LaTeXML.

Projekt współfinansowany przez Unię Europejską w ramach Europejskiego Funduszu Społecznego.

Projekt współfinansowany przez Ministerstwo Nauki i Szkolnictwa Wyższego i przez Uniwersytet Warszawski.