MASARYKOVA UNIVERZITA PŘÍRODOVĚDECKÁ FAKULTA ÚSTAV MATEMATIKY A STATISTIKY
Diplomová práce
Brno 2015
Ekaterina Pushkareva
MASARYKOVA UNIVERZITA PŘÍRODOVĚDECKÁ FAKULTA ÚSTAV MATEMATIKY A STATISTIKY
R/S analýza kursů, finančních aktiv Diplomová práce
Ekaterina Pushkareva
Vedoucí práce: RNDr. Václav Studený, Ph.D.
Brno 2015
Bibliografický záznam Autor:
Ekaterina Pushkareva Přírodovědecká fakulta, Masarykova univerzita Ústav matematiky a statistiky
Název práce:
R/S analýza kursů, finančních aktiv
Studijní program:
Matematika
Studijní obor:
Statistika a analýza dat
Vedoucí práce:
RNDr. Václav Studený, Ph.D
Akademický rok:
2014/2015
Počet stran:
vii + 64
Klíčová slova:
R/S analýza; finance; dlouhodobé časové řady; finanční trh; Hurstův exponent; V-statistika; persistentní vlastnosti časový řad
Bibliographic Entry Author:
Ekaterina Pushkareva Faculty of Science, Masaryk University Department of Mathematics and Statistics
Title of Thesis:
R/S analysis of financial assets
Degree programme:
Mathematics
Field of Study:
Statistics and Data Analysis
Supervisor:
RNDr. Václav Studený, Ph.D
Academic Year:
2014/2015
Number of Pages:
vii + 64
Keywords:
R/S analysis; finances; long-term time series; financial market; Hurst exponent; V-statistics; persistence properties of time series
Abstrakt Táto diplomová práce se věnuje R/S analýze finančních aktiv, jimiž jsou v tomto případě ceny akcií. První kapitola je zaměřena na historii problematiky, především na stručný popis existujících druhů analýz finančních aktiv. V další části této práce se odvádí vztahy k výpočtu Hurstova exponentu a V-statistiky. Také jsou zmíněny důležité vlastnosti R/S-analýzy. Poslední kapitola je věnována R/S analýze dlouhodobých časových řad reprezentujících ceny akcií vybraných společností. Prozkoumané a shrnuté výsledky analýzy dávají údaje, jejž mohou být použité při práce na finančním trhu.
Abstract The thesis is dedicated to R/S analysis of financial assets, which are mainly stock prices. The first chapter is devoted to the historical development of other analyses of financial assets. The following chapter describes the development of formulas for Hurst exponent and V-statistics. Main properties of R/S analysis are also given in this chapter. The last chapter is focused on R/S analysis that studies long term time series, which represent stock prices of world well-known corporations. Studied and summarized results of the analysis are given as certain data that could be used in stock markets and sectors.
Poděkování Na tomto místě bych chtěla poděkovat vedoucímu mé diplomové práce RNDr. Václavu Studenému, Ph.D. za odborné rady, pomoc a trpělivost, za čas, který mi věnoval. V neposlední řadě bych také chtěla poděkovat svým rodičům, sestře a manželovi za podporu, pomoc a lásku.
Prohlášení Prohlašuji, že jsem svoji diplomovou práci vypracovala samostatně s využitím informačních zdrojů, které jsou v práci citovány.
Brno 5. ledna 2015
……………………………… Ekaterina Pushkateva
Obsah Úvod ................................................................................................................................... 1 Kapitola 1. Finanční trh, metody analýzy finančních aktiv ............................................... 2 1.1. Finanční trh a riziko ................................................................................................ 2 1.2. Statistická analýza a Teorie efektivních trhů .......................................................... 3 1.3. Technická analýza ................................................................................................... 8 Kapitola 2. R/S-analýza finančních aktiv ........................................................................ 18 2.1. Objevení a vývoj R/S-analýzy............................................................................... 18 2.1.1. Wienerův proces ......................................................................................... 19 2.1.2. Výpočet Hurstova exponentu ..................................................................... 21 2.2. Druhy časových řad podle hodnoty Hurstova exponentu ..................................... 25 2.3. Vyhledávání cyklů časové řady pomocí R/S-analýzy ........................................... 28 Kapitola 3. Použití R/S-analýzy ....................................................................................... 35 3.1. Analýza cen akcií společnosti Apple Inc. ............................................................. 36 3.2. Analýza cen akcií společnosti IBM Corp. ............................................................. 39 3.3. Analýza cen akcií společnosti BP plc. .................................................................. 42 3.4. Analýza cen akcií společnosti Wal-Mart Stores Inc. ............................................ 46 3.5. Analýza cen akcií společnosti McDonald’s Corporation ...................................... 48 3.6. Analýza cen akcií společnosti General Electric .................................................... 52 3.7. Analýza cen akcií společnosti Johnson & Johnson ............................................... 55 3.8. Zobecnění .............................................................................................................. 57 Závěr ................................................................................................................................ 60 Seznam použité literatury................................................................................................. 61 Příloha 1 ........................................................................................................................... 63
vii
Úvod Tato diplomová práce se zaměřuje na analýzu dlouhodobých časových řad reprezentujících ceny akcií mezinárodních korporací. Předmět výzkumu byl zvolen tak proto, že současná teorie analýzy finančních aktiv není úplná. Předpoklad, že změny cen jsou nezávislé a mají Gaussovo rozdělení, neodpovídá v plné míře skutečné situaci na finančních trzích. Rozvoj teorie fraktálů vnesl nové nástroje k analýze finančních aktiv. Fraktaln9 model změny cen a kurzů měn dal impulz k výzkumu financí z hlediska uspořádaného chaosu. Fraktální model dobře popisuje změny cen finančních aktiv. Hlavním představitelem nové teorie, a zároveň i její zakladatelem, se stal Benoît Mandelbrot. Jeho výzkum navazuje na publikace Edwina Hursta, Vilfreda Pareta, na Cantorovy, Peanovy, Hilbertovy, Hausdorfovy a Kochovy výsledky ze začátku století a vlastní teorii fraktálů. Mandelbrot jako první přišel s novým způsobem analýzy finančních aktiv, která se zabývala především dlouhodobými časovými řadami. Později se věnoval této teorii i Edgar Peters, jež ji rozšířil o nové poznatky. Fraktální analýza je založená na předpokladu, že změny hodnot finančních aktiv jsou chaotické, ale jsou podřízené zákonitostem. Jedním z nástrojů fraktální analýzy je R/S-analýza, která je popsána v dané diplomové práci. První kapitola této diplomové práce se zaměřuje na vysvětlení základních pojmů, jako jsou např. finanční trh a riziko. Dále se zabývá historií rozvoje analýzy finančních trhů se stručným popisem vybraných druhů analýz. Z hlediska technické analýzy, kterou používají tradeři ve svých analýzách, budou v této části práce analyzovány dlouhodobé časové řady. Další kapitola se zabývá R/S-analýzou, jejím objevením a rozvojem a popisem jedné její důležité vlastnosti, která pomáhá objevit neperiodicky se opakující kolísání cen v dlouhodobých časových řadách. Část kapitoly se věnuje dvěma základním přístupům výpočtu Hurstova exponentu a jeho hodnotě pro persistentní a antipersistentní procesy. Třetí kapitola se zabývá cenami akcií sedmi vybraných společností působících na mezinárodním trhu. K analýze a objevení cyklů ve zkoumaných časových řadách je využita hodnota Hurstova exponentu, která je popsána ve druhé kapitole. Práce je založená především na výzkumech Benoîta Mandelbrota ([7], [8], [9]) a Edgara Petersa ([15], [16]) s použitím některých definicí z Otevřené encyklopedie Wikipedie (https://cs.wikipedia.org)
1
Kapitola 1 Finanční trh, metody analýzy finančních aktiv Tato kapitola je věnována základním pojmům finančního trhu, metodám analýzy finančních aktiv a historii jejich rozvoje. Kapitola je zpracována na základě [7] a [16] s použitím obrázků z [16] nebo nakreslených pomocí softwaru Maple 18. Soubory zanalyzovaných dat jsou stáhnuty se serveru www.google.com/finance.
1.1. Finanční trh a riziko Definice 1.1. Finanční trh je systém institucí a instrumentů zabezpečující pohyb peněz a kapitálu (nabízeného ve formě cenných papírů) ve všech jeho formách mezi různými ekonomickými subjekty; a to na základě poptávky a nabídky. Poptávka na finančním trhu má pro pořizovaný kapitál tři kritéria, kterými jsou riziko, likvidita a výnosnost [11]. Definice 1.2. Riziko představuje pravděpodobnost, že návratnost investice bude jiná než návratnost očekávaná [11]. Definice 1.3. Trend je dominantní směr na trhu v rámci zobrazeného časového rámce. Rozlišuje se vzestupný trend, sestupný trend a pohyb do strany [22]. Během 20. století se finančníci a ekonomové snažili zanalyzovat a pochopit rizika finančních trhů, objevit jejich původ, kvantitativně je ohodnotit a získat výhodu. Pro tyto potřeby byla analytiky dlouhou dobu využívána tzv. fundamentální analýza. Definice 1.4. Fundamentální analýza je jedna z nejstarších analýz finančních aktiv. Říká, že hlavní důvody, proč kurz cenných papírů roste či klesá, je třeba hledat ve společnosti, která je vlastní, v oblasti, ve které společnost působí, či je kurz závisí na ekonomické situaci v dané oblasti. Důkladný průzkum pomůže nejen objevit důvod, ale i předpovědět další změnu ceny. Kurzy akcií, obligací, opcí a měn se mění proto, že trh je vystaven vnějším vlivům. Tím pádem je fundamentální analýza založena na předpokladu, že pokud je znán důvod změny ceny, lze předpovědět další změny a předcházet rizikům [7]. Bohužel skutečnost není tak jednoduchá. V současné době jsou totiž předmětem obchodu i informace. Proto mohou být některá fakta, data a události přístupné jen malému počtu subjektů. Důležité informace jsou často zatajovány, či převraceny. Přesto 2
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________3
stále působí na finanční trh a objekty obchodu. Zprávy o válkách, přírodních katastrofách, spojení společností, výrobě nových produktů se hodnotí každým účastníkem trhu subjektivně. Někdo předpokládá růst akčních cen nebo hodnoty měn, jiní naopak předpovídají pokles. Co se však děje ve skutečnosti? Fundamentální analýza se toto zkouší předpovídat, avšak takový odhad není zcela přesný. Je pochopitelné, že účastníci finančního trhu vždy potřebovali co největší přesnost finančních předpovědí. A tedy i nové nástroje k analýze trhu. Nová generace analytiků začala používat metody z teorie pravděpodobnosti a statistiky. Základní koncepce nové školy prohlašovala, že ceny jsou náhodné a nepředvídatelné, avšak jejich kolísání lze popsat matematickými zákony. Následující paragraf se zabývá podrobným popisem tohoto přístupu, jeho historií, výhodami a nevýhodami a pro lepší pochopení je uvedeno několik příkladů.
1.2. Statistická analýza a Teorie efektivních trhů Statistická analýza finančních aktiv je založená na výpočtu statistických měr, kterými jsou střední hodnota, směrodatná odchylka, kvantily atd. K tomu se používá především normální rozdělení. Nejprve si připomeňme základní definici. 2 Definice 1.5. Normální (Gaussovo) rozdělení N ( , ) s parametry a 0 je rozdělení určené hustotou
1 f x e 2
( x )2 2 2
, x (, )
(1.1)
Rozdělení N (0,1) s parametry 0 a 1 se nazývá normované normální rozdělení. Tedy hustota je určená vztahem 2
1 x2 x e , x (, ) 2
(1.2)
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________4
Obrázek 1.1. Hustota normálního rozdělení pravděpodobnosti pro různá ( , ) .
Za zakladatele nové analýzy finančního trhu se považuje francouzský vědec Louis Bachelier, který použil model Brownova pohybu pro finanční aktiva ve své disertační práci Teorie spekulací (francouzsky Théorie de la spéculation). Bachelier psal, že se ceny mohou zvýšit nebo klesnout se stejnou pravděpodobností, a porovnal tento proces s házením mince, kdy orel nebo panna mohou padnout se stejnou pravděpodobností. Z toho vyplývá, že změny cen na finančním trhu jsou náhodným procesem s nulovou střední hodnotou a přírůstky cen tak mohou být odhadnuty [7], [16]. Za tohoto předpokladu zkusil Bachelier použít normální rozdělení k popsání změny cen. Je 68 % změn menších než jedna směrodatná odchylka od střední hodnoty, 95 % změn patří do intervalu ( 2 , 2 ) a 98 % do intervalu (3 ,3 ) .
Obrázek 1.2. Rozdělení dat do
-intervalů
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________5
Ze zkušenosti víme, že relativní četnost změn větších než 3 je vetší, než by měla být, kdyby se tyto změny řídily normálním rozdělením. Je známo, že výnosy či změny cen mají více velkých výkyvů a proto neodpovídají normálnímu rozdělení přesně. Ale při analýze aktiv se takové údaje zanedbávají, aby bylo možné pokračovat s použitím normálního rozdělení. Ačkoliv Bachelier navrhl statistickou analýzu finančního trhu již na začátku 20. století, do praxe se dostala až o několik desítek let později. Dalším krokem rozvoje statistické analýzy se stala Teorie efektivních trhů (anglicky Efficient Market Theory – EMH), navržena americkým ekonomem Eugenem Famou. Definice 1.6. Uvažujeme prostor (, , P) , kde je množina elementárních tržních situací, je σ-algebra podmnožin prostoru , P je pravděpodobnostní míra na . Nechť {St }t 0 je prostor cen indexů cenných papírů. Kapitálový trh (, , P) se nazývá efektivním, jestliže E St St k 0 , kde St k St k St k 1 [5]. Teorie efektivních trhů je jednou z teorií, která se pokouší popsat chování kurzů cenných papírů se zaměřením na akcie. Tato teorie předpokládá, že kurzy cenných papírů jsou ovlivňovány pouze objektivními informacemi, očekávanými zisky, dividendami, možnými riziky a dalšími kurzotvornými informacemi. Tržní cena akcií na trhu pak představuje objektivní hodnotu, akcie jsou v každém okamžiku správně oceněny a na trhu nelze najít podhodnocené nebo nadhodnocené tituly. Z toho vyplývá, že úspěšnost obchodování není možno zvýšit fundamentální či technickou analýzou, ani studiem historických údajů. Trh reaguje jen na nové informace a je tak zcela nepředvídatelný [14]. Je hodně stoupenců a odpůrců Teorie efektivních trhů. Hlavním argumentem odpůrců je fakt, že na rozdíl od teorie, skutečnost není ideální. V podstatě jde o to, že ceny nejsou objektivní, neboť se na trhu nachází velké množství nadhodnocených a podhodnocených cenných papírů. Změny cen nejsou náhodné veličiny. Ceny mají tzv. „paměť“, tj. dnešní cena má vliv na budoucí. Jestliže se dnes cena významně zvýší, pravděpodobnost, že bude růst pokračovat i zítra, je velmi vysoká [7]. Nejlépe je prozkoumána krátkodobá závislost. Pod tímto pojmem rozumíme vliv kolísání ceny v daném okamžiku na kolísání ceny v určitém blízkém časovém období. Blízkost budoucího okamžiku záleží na časovém horizontu výzkumu. Nejčastěji se využívá časový horizont několika hodin, dní, týdnů či let. Zkoumání krátkodobé závislosti změny cen ukazuje na „efekt setrvačnosti“ [7]. Ten spočívá v tendenci zachování trendu: jestli cena akcií roste/klesá v daném časovém období, patrně bude růst/pokles pokračovat ještě nějakou dobu. Přičemž čím rychlejší
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________6
byl růst/pokles v prvním období, tím bude pravděpodobnost zachování trendu v následujícím období vyšší. Avšak ve střednědobé perspektivě (tři až osm let) panuje opačná situace. Akcie, jejichž ceny během několika předchozích let rostly, budou s největší pravděpodobností v následujících letech klesat, a naopak. Tuto zákonitost objevili Eugen Fama a Kenneth French v roku 1988 [2]. Dalším důvodem ke kritice klasické statistické analýzy je rozdíl mezi rozdělením změn ceny a grafem hustoty normálního rozdělení. Riziko, že dojde k velkým změnám hodnot, je větší, než ukazuje normální rozdělení. Podle normálního rozdělení, pravděpodobnost, že nastane událost větší než 3 , činí 0,5 %. Pravděpodobnost, že nastane událost větší než 4 by měla být 0,01 %. Nicméně, skutečné změny neodpovídají daným předpokladům. Edgar Peters [15] prozkoumal velké soubory dat a zjistil, že pravděpodobnost, že nastane událost větší než 3 , je 2,4 %, 4 – 1 %. Jako příklad lze uvést změnu hodnot Dow Jonesova indexu za období od 2. ledna 1888 do 31. prosince 1991. Do grafu vyneseme hodnoty hustoty rozdělení 5-denních a 90-denních výnosů a porovnáme je s hustotou normálního rozdělení. Výsledek lze vidět na obrázku 1.3.
Obrázek 1.3. Hustota rozdělení výnosů, Dow Jonesův index 1888-1991 [16]
Obě dvě rozdělení mají ostrý vrchol ve střední hodnotě a tlusté konce. Navíc 5denní a 90-denní výnosy mají skoro stejný tvar rozdělení s ostrou špičkou a relativně velkými pravděpodobnostmi výrazných změn. Podle Petersa mají shodný tvar rozdělení výnosy i pro jiné časové horizonty. Obrázky 1.4-1.6 ukazují rozdíl mezi rozdělením hustoty výnosů a hustoty normálního rozdělení.
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________7
Obrázek 1.4. Rozdíl hustoty rozdělení jednodenních výnosů a hustoty normálního rozdělení, Dow Jonesův index 1888-1991 [16]
Obrázek 1.5. Rozdíl hustoty rozdělení 30-denních výnosů a hustoty normálního rozdělení, Dow Jonesův index 1888-1991 [16]
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________8
Obrázek 1.6. Rozdíl hustoty rozdělení 90-denních výnosů a hustoty normálního rozdělení, Dow Jonesův index 1888-1991 [16]
Na každém grafu lze vidět, že pravděpodobnost změny ceny v hodnotě 4σ je dost vysoká. Nezáleží na tom, se kterým časovým horizontem pracuje trader, všichni mají skoro stejné riziko narazit na velký výkyv cen. Většina traderů nepoužívá v praxi statistickou analýzu, ale dává přednost technické analýze trhů.
1.3. Technická analýza Definice 1.7. Technická analýza je systematick7m zkoumáním, analyzováním a vyhodnocováním starších i současných dat. Používá se na předpovídání budoucích cenových pohybů. Tato metoda je používána u všech finančních produktů, včetně cenných papírů, futures kontraktů a úrokových produktů. Na rozdíl od fundamentální analýzy využívá pouze údaje tvořené trhem, jako je např. cena, objem, množství otevřených kontraktů na trhu, popřípadě mezitržní korelace. Technická analýza se proto nezabývá takovými jevy a skutečnostmi, jako jsou zveřejněná ekonomická data, politická situace, daňová politika státu nebo ekonomické prostředí [12]. Ve skutečnosti se technická analýza stává analýzou grafů. Existuje velký počet různých druhů grafů, které ukazují změny cen, jejich přírůstků, objem transakcí, jejich střední hodnoty, volatilitu, atd. Analytici se řídí pravidly která říkají „Ceny
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________9
se nepohybují náhodně, ale v trendech„ a „Historie má tendenci se opakovat“. Proto hledají trendy a na základě svých odhadů se snaží předpovědět vývoj ceny. V následující části práce budou použity metody technické analýzy a prozkoumány vybrané časové řady. Vzhledem k tomu, že další výzkum se bude zabývat časovými řadami s velkým počtem pozorování, bude také provedena grafická analýza těchto časových řad. Pro účely této práce byly vybrány záznamy o cenách akcií sedmi mezinárodních společností, jež působí na trhu minimálně třicet let. Těmi jsou
Apple Inc., IBM Corp., BP plc. (dřív British Petroleum), Wal-Mart Stores Inc., McDonald’s Corporation, General Electric, Johnson & Johnson.
Společnosti byly zvoleny tak, aby se jejich činnost týkala různých oblastí. Apple Inc. a IBM Corp. působí na trhu digitálních technologií, BP plc. je energetická společnost, Wal-Mart Stores Inc. je řetězec obchodních domů, McDonald’s Corporation je světově proslulý řetězec restaurací, General Electric vyniká v oblasti technologie a Johnson & Johnson je globální americká farmaceutická firma. Čtyři z sedmi výše uvedených společností patří do struktury Dow Jonesova indexu. Byly analyzovány ceny akcií společnosti Apple Inc. za období od 2. ledna 1981 do 31. října 2014. Jedná se o soubor 8 534 nejvyšších denních cen, jejichž hodnoty jsou zobrazeny na obrázku 1.7.
Obrázek 1.7. Ceny akcií Apple Inc., 1981-2014
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________10
Z vývoje cen akcií je vidět, že kolísání cen není na začátku sledovaného období výrazné a v porovnání s dnešními hodnotami, je tehdejší cena akcií velmi nízká, v podstatě zanedbatelná. Pro tuto situaci je vhodnější zobrazit změny na logaritmické stupnici. Výhodou používání logaritmické stupnice je, že umožňuje názorně zobrazovat veličiny v rozpětí mnoha řádů. Díky tomu diagramy uvádí skutečnou tržní situaci během daného období. Na obrázku 1.8 logaritmická funkce nabývá nejen kladné, ale i záporné hodnoty, což odpovídá skutečnosti, kdy akcie Apple Inc. stály méně než jeden americký dolar. Kolísání cen je tak výraznější než na obrázku 1.7, přičemž můžeme pozorovat značné skoky ceny, kdy např. v roce 1987 došlo k výraznému nárůstu a naopak v roce 2001 k poklesu.
Obrázek 1.8. Ceny akcií Apple Inc. na logaritmické stupnici, 1981-2014
Obrázek 1.9. Přírůstky cen akcií Apple Inc., 1981-2014
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________11
Graf na obrázku 1.9 znázorňuje turbulentní vlastnost finančního trhu: čím dále se pohybujeme po časové ose, tím můžeme pozorovat větší přírůstky v cenách. Lze vidět, že během zkoumaného období nastávaly výrazné poklesy a růsty. Podle níže uvedeného vzorce byl poté vytvořen graf relativních přírůstků: Yt
X t 1 X t Xt
(1.3)
Obrázek 1.10. Relativní přírůstky cen akcií Apple Inc., 1981-2014
Na obrázku 1.10 lze zřetelně pozorovat, že ostré změny ceny akcií nastávaly nepředvídatelně po obdobích relativního klidu. Pro další společnosti budou vykresleny pouze grafy cen na logaritmické stupnici a grafy relativních přírůstků. Analýza vývoje cen akcií společnosti IBM Corp. byla provedena za období od 3. ledna 1978 do 31. října 2014. Celkem bylo použito 9 293 hodnot.
Obrázek 1.11. Ceny akcií IBM Corp. na logaritmické stupnici, 1978-2014
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________12
Na obrázku 1.11 můžeme vidět docela výrazné změny ceny. Přičemž je těžké předpovědět, jakým směrem bude vývoj ceny pokračovat v následujících letech.
Obrázek 1.12. Relativní přírůstky cen akcií IBM Corp., 1978-2014
Obrázek 1.12 názorně ukazuje, že ostré změny nastávaly docela často a občas i nečekaně. Stejné časové období jako v případě společnosti IBM Corp. je předmětem analýzy i u následující společnosti, kterou je BP plc. Celkem bylo využito 9 294 záznamů o výši denních cen akcií této společnosti. Zase začneme z pohybu cen za dané období na logaritmické stupnici.
Obrázek 1.13. Ceny akcií BP plc. na logaritmické stupnici, 1978-2014
Situace je stejná: existují zde náhlé změny ceny a to jak směrem dolů, tak i nahoru. Menších změn je poté evidováno poměrně značné množství. To potvrzuje i graf relativních přírůstků.
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________13
Obrázek 1.14. Relativní přírůstky cen akcií BP plc., 1978-2014
Společnost Wal-Mart Stores Inc. působí na trhu od roku 1962. Pro analýzu cen akcií bylo vybráno období od 3. ledna 1984 do 31. října 2014, což činí celkem 7 776 záznamů. Toto období je zvoleno tak proto, aby odpovídalo našemu požadavku délky (nejméně 30 let) a ceny při tom byly vyšší než 1 dolar. Situace, kdy jsou ceny akcií nižší než jeden dolar, byla znázorněna na příkladu společnosti Apple Inc. Na logaritmické stupnici je níže znázorněn pohyb cen akcií společnosti Wal-Mart Stores Inc.
Obrázek 1.15. Ceny akcií Wal-Mart Stores Inc. na logaritmické stupnici, 1984-2014
V porovnání s předchozími grafy (obrázky 1.8, 1.11, 1.13) nemá graf na obrázku 1.15 natolik výrazné poklesy nebo růsty. V posledním desetiletí je kolísání cen poměrně stabilní. Podobnou situaci lze sledovat také v období 1993 až 1998. Od roku 1998 do roku 2000 pak období relativního klidu střídá období, kdy dochází k poměrně výraznému nárůstu cen akcií. Můžeme očekávat náhlý růst i v současné době? Nebo dojde naopak k poklesu? Jsou akcie Wal-Mart Stores Inc. vhodné k dlouhodobému investování? Nebo nemá vůbec smysl přidávat dané aktivity do portfolia cenných papírů? Bohužel ani technická analýza nám nedokáže odpovědět na tyto otázky. Bez ohledu na růst cen na konci 20.století jsou relativní přírůstky docela nízké.
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________14
Obrázek 1.16. Relativní přírůstky cen akcií Wal-Mart Stores Inc., 1984-2014
Také pro společnost McDonald’s Corporation bylo zvoleno období cen akcií od 3. ledna 1984 do 31. října 2014, jedná se celkem o 7 776 záznamů. Kreslí se pohyb cen za dané období na logaritmické stupnici.
Obrázek 1.17. Ceny akcií McDonald’s Corporation na logaritmické stupnici, 1984-2014
Vývoj ceny akcií ukazuje jak na její růst, tak i na její pokles. Na obrázku 1.18 jsou pak vidět velké skoky v relativních přírůstcích cen.
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________15
Obrázek 1.18: Relativní přírůstky cen akcií McDonald’s Corporation, 1984-2014
Pro ceny akcií společnosti General Electric bylo analyzováno období od 3. ledna 1984 do 31. října 2014, počet záznamů činí 7 777. Na rozdíl od předchozích grafů pohybu cen akcií, je situace na grafu 1.19 odlišná. V letech 2008 a 2009 lze sledovat výrazný pokles cen akcií, který nebyl, v takové míře, u předchozích společností zaznamenán.
Obrázek 1.19. Ceny akcií General Electric na logaritmické stupnici, 1984-2014
Obrázek 1.20. Relativní přírůstky cen akcií General Electric, 1984-2014
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________16
Na obrázku 1.20 je vidět, že většinou byly malé změny cen, ale došlo i k docela výrazným skokům směrem nahoru a dolů. Tedy je zřejmá turbulence kolísání a lze vidět velké výkyvy ceny. Poslední analýza se týká společnosti Johnson & Johnson, pro kterou je k dispozici celkem 9 293 záznamů z období od 3. ledna 1978 do 31. října 2014.
Obrázek 1.21. Ceny Johnson & Johnson na logaritmické stupnici, 1978-2014
Obrázek 1.22. Relativní přírůstky cen akcií Johnson & Johnson, 1978-2014
Na obrázcích 1.21-1.22 je patrná jistá turbulence v kolísání cen akcií. Ačkoliv má graf na obrázku 1.21 rostoucí trend, neznamená to však, že bude tento trend zachován i v následujícím období. Na výše uvedených grafech je vidět trendy, ale na základě historických údajů není možné přesně předpovědět, kdy začne korekce ceny a směr trendu se změní. Dobrým příkladem je společnost General Electric. V období, kdy akcie ostatních firem rostou, u GE dochází k výraznému poklesu. Kdybychom se však na tyto grafy dívali
Kapitola 1. Finanční trh, metody analýzy finančních aktiv____________________________17
v roce 2000, mohli bychom předpokládat, že růst akcií General Electric bude pokračovat. Z výše uvedených grafů lze vypozorovat, že k výraznému růstu a kolísání cen akcií docházelo ve všech případech v 90. letech 20. století. Jak by to bylo možné vysvětlit? Jedním z předpokladů je fakt, že v 80. letech zavedli některé burzy systémy prodeje malého množství aktiv. V 90. letech pak, s rozmachem informačních technologií, získalo velké množství lidí snadný přístup k internetu a pomocí on-line brokerů mohlo začít obchodovat na finančních burzách. Tím pádem se počet subjektů trhu mnohonásobně zvýšil a kontrola cen akcií se stala komplikovanější. Odhad rizika je komplikovaný a investor nemůže vědět, jak velké budou výnosy a ztráty. Nesmíme zapomenout, že oceňování rizika je nejdůležitější věc pro investora. Statistická a technická analýzy mohou pomoci zkušenému triedru odhadnout trend. Ale nemohou dát odpověď na otázku, kolik investor ztratí v případě špatného odhadu. Odpověď na tuto otázku nám dává R/S-analýza. V této kapitole jsme se seznámili s historií rozvoje analýzy finančního trhu a uvedli několik příkladů ukazujících, že se klasická technická analýza ne vždy hodí k pochopení chování cen akcií. Také byla provedena základní analýza časových řad odpovídajících cenám akcií sedmi korporací, které působí na trhu již delší dobu. V následujících kapitolách se seznámíme s docela novou metodou analýzy velkých časových řad – R/S-analýzou a zkusíme prozkoumat dané časové řady tímto způsobem.
Kapitola 2 R/S-analýza finančních aktiv V této kapitole je popsána historie objevení a rozvoje R/S-analýzy, metody výpočtu hlavního parametru R/S-analýzy, tzv. Hurstova exponentu (H), třídění časových řad podle hodnoty H. Dále jsou zde uvedeny příklady použití R/S-analýzy a zmíněny jsou také její důležité vlastnosti. Kapitola je vypracována především na základě [6], [7] a [16] s použitím obrázků z těchto zdrojů.
2.1. Objevení a vývoj R/S-analýzy Historie R/S-analýzy a Hurstova exponentu sahá již do začátku 20. století, avšak tehdy měla s finanční matematikou jen pramálo společného. V roce 1906 přijel anglický vědec Harold Hurst do Egypta. Účelem jeho návštěvy bylo postavit na řece Nil velkou Asuánskou přehradu. Hlavním úkolem Hursta bylo stanovit výšku přehrady tak, aby v případě záplav nepřekročila hladina výšku hráze a také aby bylo v přehradě dostatečné množství vody i v obdobích sucha. Náročnost úkolu spočívala v tom, že nebylo možné rozsah povodně předpovědět. Nil má velmi složitý vodní režim, průtok se nepředvídatelně mění každý rok. Deštivé roky se seskupují bez zřejmých period, což platí i v případě roků suchých [13]. Hurst měl k dispozici roční data minimálních hladin řeky od roku 622 do roku 1469 (celkem 847 údajů) a všiml si, že existují neperiodické cykly ve změně průtoku. Hurst tedy došel k závěru, že hodnoty hladiny řek nejsou nezávislé a záplavy z minulých let mají vliv na současnou situaci [16]. Standardní analýza však neprokázala vzájemný vztah mezi pozorováními. Hurst tak vynalezl novou metodu analýzy, která pak byla modifikovaná Benoîtem Mandelbrotem. Své výklady Hurst začal ze standardního předpokladu: hodnoty časové řady jsou náhodné a nezávislé. To znamená, že daná časová řada je náhodnou procházkou. Víme, že jednorozměrný případ náhodné procházky můžeme simulovat házením mince. Pravděpodobnost, že padne panna nebo orel je stejná a rovná se 0.5. Vždy, kdy padne panna, zvětšíme svůj výsledek o jeden bod, kdy padne orel, zmenšíme výsledek o jeden bod. Při dost velkém počtů házení vidíme, že existují periody, kdy pořád padá panna, a také dochází k periodám, kdy několikrát po sobě padne orel. Jak můžeme spočítat
18
Kapitola 2. R/S-analýza finančních aktiv__________________________________________19
výsledek? Použijme Wienerův proces, který je matematickým modelem Brownova pohybu. Ale předem si uveďme několik základních definicí. Definice 2.1. Rozsah je rozdíl mezi nejmenší a největší hodnotou veličiny, která se v daném případě může objevit. Definice 2.2. Směrodatná odchylka je statistickou mírou, která ukazuje průměr druhých mocnin odchylek hodnot znaku od jejich aritmetického průměru a je určená vztahem
E ( X 2 ) ( E ( X ))2
(2.1)
Definice 2.3. R/S-analýza je jedna z metod analýzy finančních trhů. Je založena na výpočtu poměru mezi rozsahem časové řady na daném intervalu a odhadem směrodatné odchylky, který je úměrný počtu pozorování na daném intervalu umocněnému na nějakou konstantu H: R c NH , S
(2.2)
kde R je rozsah časové řady, S je směrodatná odchylka, N je počet pozorování, c je konstanta. Definice 2.4. Parametr H ze vztahu (2.2) se nazývá Hurstův exponent, nabývá hodnoty 0 H 1 , je charakteristickou mírou pro dlouhou paměť v časových řadách.
2.1.1. Wienerův proces Definice 2.5. Wienerův proces Wt je stochastický proces spojitého času s přírůstky nezávislými na poloze, který splňuje: 1. W0 0 . 2. Wt je téměř skoro jistě spojitý. 3. Wt
má
na
poloze
nezávislé
přírůstky
s rozdělením
Wt Ws ~ N (0, 2 (t s)) pro 0 s t . Wienerův proces Wt lze chápat jako limitu náhodné procházky při zmenšování časového a prostorového kroku
x a t (tj. x 0 , t 0 ). Ukažme, že
t ( x) 2
(2.3)
Kapitola 2. R/S-analýza finančních aktiv__________________________________________20
Nechť P( X i 1) P( X i 1)
1 , kde 2
X i ,.. X n ,... jsou stejně rozdělené
nezávislé náhodné veličiny s E ( X i ) 0 a Var ( X i ) 1 , S 0 je standardní symetrická náhodná procházka. Nechť
Sn S0 X1 X 2 .. X n Zvolíme délku časového kroku t a velikost prostorového kroku t t n t (tedy n ) definujeme proces: t
St Sn t ( X 1 X 2 .. X n ) x
(2.4)
x . Pro
(2.5)
Z nezávislosti přírůstků X j plyne, že E ( St ) 0 a Var ( St ) ( x) 2 n ( x) 2
t t
(2.6)
Zjistíme, jak se chová tento proces při x 0 a t 0 . Uvažujeme mocninnou závislost mezi x a t . Stanovme pro p 0
t ( x) p
(2.7)
Tudíž pro t 0 máme možnosti:
0, pro p 2 2 ( x) Var ( St ) t t , pro p 2 t , pro p 2
(2.8)
Tedy konečný nenulový rozptyl dostaneme pro p 2 , tj. platí (2.3) [4]. Z toho plyne, že
x ( t)
1 2
(2.9)
Zobecněním (2.9) dostaneme vzorec k výpočtu vzdálenosti mezi nejlepším a nejhorším výsledkem: 1
R cT 2
kde R je rozsah, T je čas, c je kladná konstanta.
(2.10)
Kapitola 2. R/S-analýza finančních aktiv__________________________________________21
Hurst použil tento vzorec na hodnoty výšky hladiny Nilu a zjistil, že vztah neodpovídá skutečnosti. V podstatě se rozsah zvětšoval rychleji než druhá odmocnina z počtu pozorování. Tento fakt podnítil Angličana k vymyšlení vlastního vztahu závislosti mezi rozsahem a časem. V následující části je popsán výpočetní postup.
2.1.2. Výpočet Hurstova exponentu Uvažujeme časovou řadu X X i , i 1,.., N . Střední hodnota se počítá podle vzorce:
EX
X 1 .. X N 1 N N
N
X
(2.11)
i
i 1
Směrodatná odchylka je: SN
1 2
N
( X
i
EX )2
(2.12)
i 1
Normovanou řadu dostaneme pomocí vzorce: Zi X i EX , i 1,.., N
(2.13)
Pomocí (2.11) ukažme, že získaná řada Z Zi , i 1,.., N má nulovou střední hodnotu: 1 N
EZ
N
Zi i 1
1 N
N
X i i 1
1 N
N
1
EX EX N N EX 0
(2.14)
i 1
V dalším kroku vytvoříme kumulovanou časovou řadu: t
t
i 1
i 1
Yi Zi ( X i EX ) X 1 X 2 .. X t t EX (2.15)
t
X i t EX , t 1,.., N i 1
Všimneme si, že poslední člen takto definované časové řady vždy bude nulový: N
N
i 1
i 1
YN X i EX X i N EX N N EX N EX 0
1 N X i N EX N i 1
(2.16)
Kapitola 2. R/S-analýza finančních aktiv__________________________________________22
Definujeme upravený rozsah RN vztahem:
RN max Y1 ,.., YN min(Y1 ,.., YN ) t
t
max(X i t EX ) min(X i t EX ) 1t N
i 1
1t N
(2.17)
i 1
Kumulovaná časová řada Y je normovaná k nulové střední hodnotě. Proto maximální hodnota Y je vždy větší nebo se rovná nule a minimální hodnota Y je vždy menší nebo se rovná nule. Z toho vyplývá, že upravený rozsah RN může nabývat pouze kladné nebo nulové hodnoty. Upravený rozsah RN je délkou cesty, kterou prochází systém během času N. Dále pomocí (2.10) objevil Hurst nový vztah mezí rozsahem a časem: R/S
RN cNH , S
(2.18)
kde c je konstanta. Definice 2.6. Veličina R / S ze vztahu (2.18) se nazývá normovaným rozsahem časové řady [16]. Hurstův exponent, který tak později nazval Benoît Mandelbrot, lze přibližně odhadnout prostřednictvím grafu s hodnotami log( R / S ) nakreslenými proti hodnotám
log( N ) . Používáme metodu nejmenších čtverců na základě rovnice: R log C H log N , S
(2.19)
kde C log c je konstanta. Kdyby zkoumaná časová řada byla náhodná a její členy nezávislé výsledkem výpočtů by byla hodnota H 0.5 . Avšak se normovaný vztah zvětšuje rychleji než druhá odmocnina z času. To znamená, že systém prochází větší distancí, než za stejný čas prochází náhodný proces. Toto je možné jen v případě, kdy změny hodnot řady mají vzájemný vliv, tj. hodnoty jsou korelovány [16]. Připomeňme si, že cílem práce Hursta byla stavba vhodné vodní přehrady a jeho vzorec je stručným popisem matematického modelu pro výpočet optimální výšky přehrady a objemu nádrže. Z tradičního hlediska by měl začít z výpočtu rozsahu jako rozdílu mezi maximálními a minimálními hodnotami. Ale Hurst standardní postup vylepšil. Nejdříve odstranil trend pro časové řady s rozdílnými počátečními okamžiky a
Kapitola 2. R/S-analýza finančních aktiv__________________________________________23
různou délkou a následně spočítal normovaný rozdíl. Tím zjistil, že takový rozsah se zvětšuje rychleji než druhá odmocnina z počtu pozorování (totiž z délky času) [7]. Normování upraveného rozsahu dělením na směrodatnou odchylku odstraňuje problém trendu v dlouhodobých časových řadách. Z hlediska finančních aktiv to znamená, že bez ohledu na inflaci můžeme porovnávat mezi sebou periody rozdělené dlouhou dobou [7].
Obrázek 2.1. Model cenového diagramu [7]
Pro lepší a přesnější výsledky při výpočtu Hurstova exponentu se zkoumá samotná časová řada a všechny její podmnožiny. Náhodně se volí počáteční hodnota t a délka intervalu , kde 0 N t . Pro každý pár (t , ) se provádí výpočty podle vzorců (2.11)-(2.18). Tímto způsobem se vypočítá pro každou hodnotu soubor hodnot R/S, odpovídající každému t. Dále, pomocí bodového diagramu a regresní analýzy se zjišťuje hodnota parametru H. Zobecníme výpočet pro dlouhodobou časovou řadu
X X i , i 1,.., N .
Náhodně zvolíme počáteční bod t, 0 t N 1 , a délku podmnožiny , 1 N t . Normovaná časová řada Z se definuje vztahem: Zi X t i EX , i 1,..,
Zkonstruujme kumulovanou časovou řadu:
(2.20)
Kapitola 2. R/S-analýza finančních aktiv__________________________________________24 u
u
i 1
i 1
Yi Zi ( X t i EX ) X t 1 X t 2 .. X t u u EX (2.21)
u
X t i u EX , u 1,.., i 1
t u
u
Sumu
X
t i
X . Tím pádem pro částečnou řadu
můžeme zapsat ve tvaru
i
i t 1
i 1
X X i , i t 1,.., t , vztah (2.17) můžeme zapsat ve tvaru:
R(t , ) max Y1 ,.., Y min(Y1 ,.., Y ) t u
t u
(2.22)
max( X i u EX ) min( X i t EX ) 1u
1u
i t 1
i t 1
Střední hodnota celé časové řady může být nicméně neznámá, kdy pracujeme jen s podmnožinou této řady. Nahradíme střední hodnotu řady X střední hodnotou částečné řady X X i , i t 1,.., t :
( EX )
X t 1 .. X t
1
t
X
(2.23)
i
i t 1
Přeměníme (2.22) s úvahou (2.23): t u
R(t , ) max( X i 1 u
i t 1
u
t u
t u
X ) min( X i
1u
i t 1
i
i t 1
u
t u
X )
(2.24)
i
i t 1
Podle (2.12) bychom měli normovat rozsah R t , směrodatnou odchylkou řady X . Nahradíme ji přepočítanou odchylkou S t , ze stejného důvodu, proč jsme nahradili EX střední hodnotou částečné řady X X i , i t 1,.., t .
S (t , )
1
t
X
i t 1
i
1
2
t
X
j t 1
j
1 2
t 1 Xi i t 1
1
2 X j j t 1 t
2
(2.25)
Dostáváme:
R/S
R(t , ) S (t , )
(2.26)
Kapitola 2. R/S-analýza finančních aktiv__________________________________________25
Pro 1 máme R(t , ) S t , 0 a poměr R / S je neurčitý. Pro 2 máme R t , 2 / S t , 2 2 bez ohledu na proces. Tím pádem nás zajímají hodnoty 3 [8].
2.2. Druhy časových řad podle hodnoty Hurstova exponentu Podle hodnoty parametru H můžeme rozlišovat tři druhy procesů:
0 H 0.5 ukazuje anti-persistentní proces. To znamená, že řada mění znaménko častěji, než bychom očekávali pro náhodný proces, a prochází menší cestu než nezávislý náhodný proces. Pokud v daném časovém okamžiku systém roste, můžeme s určitou jistotou předpokládat, že v následujícím časovém okamžiku bude systém klesat a naopak; H 0.5 znamená, že časová řada je nezávislý náhodný proces; 0.5 H 1 znamená, že řada je persistentní a následuje lokální trendy. Takové řadě se říká řada s dlouhodobou pamětí, protože změny řady mají od začátku do současnosti vliv na budoucí hodnoty.
Na obrázku 2.2 jsou uvedené diagramy chování časových řad pro tří odlišné hodnoty Hurstova exponentu. Každý krok nahoru nebo dolů ukazuje změnu časové řady od jednoho časového okamžiku k dalšímu.
Obrázek 2.2. Kolísání hodnot časových řad pro různá H [7]
Kapitola 2. R/S-analýza finančních aktiv__________________________________________26
Ukažme si výsledné grafy R/S-analýzy pro různé druhy časových řad. Použijme příklady z Mandelbrot & Wallis [8]. Modely časových řad reprezentují tři druhy procesů: persistentní, anti-persistentní a nezávislý náhodný proces. Obrázek 2.3 ukazuje závislost log( R / S ) na log( ) . Daná časová řada má N 10000 .
Hodnoty
log( R / S )
jsou
vypočtené
pro
{3, 4, 5, 6, 7,10 , 20, 40, 70,100, 200, 400, 700,1000, 2000, 4000, 7000, 9000} . Hodnoty jsou zvoleny tak, aby byly hodnoty log( ) rozdělené rovnoměrně na ose a byly více méně stejně vzdálené. Na grafu je vidět několik hodnot log( R / S ) pro každé . Veličina logaritmu je vypočtená pro stejně vzdálená t. Tento výběr je způsoben časovou náročností vyčíslení log( R / S ) pro všechna t, 1 t N , pro dané . Ovšem zvětšením počtu hodnot t zlepšujeme přesnost výpočtů a stabilitu výsledků, protože používáme střední hodnotu všech hodnot R / S . Střední hodnota každé posloupnosti hodnot log( R / S ) je zdůrazněná na grafu malým čtvercem. Pro malé hodnoty má poměr R / S komplikované chování. Pro 20 lze vidět trend. Metoda nejmenších čtverců nám dává výsledek H 0.9 .
Obrázek 2.3. Výsledný graf R/S-analýzy pro persistentní proces, H 0.9 [8] Dále uvažujeme stejně dlouhou časovou řadu, která je nezávislým náhodným procesem. Výsledný graf a trend jsou znázorněné na obrázku 2.4.
Kapitola 2. R/S-analýza finančních aktiv__________________________________________27
Obrázek 2.4. Výsledný graf R/S-analýzy nezávislého náhodného procesu, H 0.5 [8] Výrazný trend je vidět od 20 , pro menší hodnoty je chování nejisté. Na obrázku 2.5 je vidět výsledný graf R/S-analýzy časové řady s parametrem H 0.3 .
Obrázek 2.5. Výsledný graf R/S-analýzy pro anti-persistentní proces, H 0.3 [8] V tomto případě není trend výrazný ani pro hodnoty 20 . Jen zhruba můžeme připustit trend pro 70 a odhadnout H.
Kapitola 2. R/S-analýza finančních aktiv__________________________________________28
2.3. Vyhledávání cyklů časové řady pomocí R/S-analýzy V dané části je popsána důležitá vlastnost R/S-analýzy, kterou je objevení cyklů v časové řadě. Definice 2.7. Cyklus změny cen akcií označuje pravidelné střídání relativního růstu a relativního poklesu, které trvají určitou dobu. Periodický cyklus se opakuje pravidelně s určitou frekvencí. Neperiodický cyklus nemá žádnou pevnou frekvenci. Doposud bylo provedeno mnoho pokusů, ve kterých se vědci snažili najít cykly ve změnách cen na finančním trhu. K tomu se používaly veškeré známé druhy analýz. Pomocí fundamentální analýzy, statistické analýzy, technické analýzy, Fourierovy analýzy, spektrální analýzy a ani dalších metod se nepodařilo odhalit periodické nebo neperiodické cykly. R/S-analýza je odhaduje. Důkaz této vlastnosti je založen na použití fraktální dimenzi. Lze najít ve práci B. Mandelbrota Statistical Methodology For Nonperiodic Cycles: From The Covariance To R/S-Analysis [6]. Výchozí podmínkou je předpoklad, že rozsah nemůže nikdy překročit meze amplitudy. Proto veličina R / S dosáhne maximální hodnoty po skončení cyklu [16]. Nechť je časová řada definovaná vztahem Yt sin t .
(2.27)
Na každém intervalu [2 k , 2 2 k ], k 0,1, 2,.. uvažujeme 100 pozorování. Výsledný graf R/S-analýzy časové řady (2.27) je zobrazen na obrázku 2.6. Při t 100 je změna trendu zřejmá.
Kapitola 2. R/S-analýza finančních aktiv__________________________________________29
Obrázek 2.6. Výsledný graf R/S-analýzy pro časovou řadu Yt sin t , trend se mění po 100 pozorování [16] Toto byl příklad časové řady s periodickým cyklem. Následující příklad ukazuje časovou řadu s neperiodickým cyklem. Zkoumá se denní hodnotu Dow Jonesova indexu za období od 2. ledna 1888 do 31. prosince 1991. Uvažují se 20-denní přírůstky, celkem je 1 320 záznamů. 20-denní údaje tvoří přibližně jeden kalendářní měsíc. Na výsledném grafu je vidět změnu trendu. Ostrý pohyb nahoru nastává přibližně při t 52 , (log(52) 1.7) . Je možné předpokládat, že daná časová řada má cyklus rovný přibližně padesáti dvěma provozním měsícům, což je 1 040 provozních dnů.
Kapitola 2. R/S-analýza finančních aktiv__________________________________________30
Obrázek 2.7. Výsledný graf R/S-analýzy 20-denních výnosů Dow Jonesova indexu [16] Ovšem taková mezera v trendu není výrazná a potřebuje se ověřit, zda předpoklad o existenci cyklu je správný nebo ne. K tomu se používá jiný nástroj R/Sanalýzy, který se nazývá V-statistika. Definice 2.8. V-statistika je statistická míra časové řady, která je definována vztahem R V (t , ) ( ) / , S
kde t je čas,
(2.28)
R je normovaný rozsah časové řady na intervalu [t , t ] , je S
počet pozorování. Na graf se vynáší hodnoty V (t , ) proti hodnotám log( ) . Dále se bodový graf aproximuje vhodnou přímkou a zjišťuje se hodnota H Hurstova exponentu. Pokud časová řada prokazuje persistentní vlastnost ( H 0.5 ), graf V-statistiky roste. V okamžiku, kdy časová řada změní svou persistentní vlastnost a přemění se v náhodnou procházku ( H 0.5 ) či antipersistentní proces ( H 0.5 ), směr grafu se změní k vodorovné přímce nebo bude klesat [16].
Kapitola 2. R/S-analýza finančních aktiv__________________________________________31
Teď nakreslíme graf V-statistiky pro 20-denní výnosy Dow Jonesova indexu. Na obrázku 2.8 je vidět očekávanou změnu směru grafu pro t 52 .
Obrázek 2.8. V-statistika 20-denních výnosů Dow Jonesova indexu [16] Tím pádem je možné pokračovat v tvrzení, že výnosy Dow Jonesova indexu mají cyklus přibližně 1 040 provozních dnů, což tvoří téměř 4 roky. Čtyřletý cyklus může souviset s ekonomickými cykly, například s Kitchinovým cyklem [20]. Avšak, jestli je tento čtyřletý cyklus opravdovým cyklem a neobjevil se kvůli výpočtové chybě, měl by být nezávislý na zvolené délce časové řady. Zkoumá se stejná časová řada denních záznamu za období od 2. ledna 1888 do 31. prosince 1991. Pro 5-denní přírůstky je 5 280 hodnot. Na obrazcích 2.9 a 2.10 jsou R/S-analýza a V-statistika 5-denních výnosů. Je vidět mezeru v trendu Hurstova exponentu a změnu směru grafu V-statistiky pro hodnotu log(t ) 2.32 , tj. pro t 209 . Přepočítáním na provozní dny se dojde k výsledku 1 045 dnů.
Kapitola 2. R/S-analýza finančních aktiv__________________________________________32
Obrázek 2.9. Výsledný graf R/S-analýzy 5-denních výnosů Dow Jonesova indexu [16]
Obrázek 2.10. V-statistika 5-denních výnosů Dow Jonesova indexu [16] Tudíž se čtyřletý cyklus objevuje nezávisle na zvoleném přírůstku času. Proto lze říct, že objevení tohoto cyklu není výpočtovou chybou nebo důsledkem šumu.
Kapitola 2. R/S-analýza finančních aktiv__________________________________________33
Zajímavé by bylo zjistit, zda-li existují cykly s periodou menší než 4 roky. Prozkoumají se jednodenní výnosy Dow Jonesova indexu. Je 24 900 pozorování a výsledný graf R/S-analýzy je na obrázku 2.11.
Obrázek 2.11. Výsledný graf R/S-analýzy denních výnosů Dow Jonesova indexu [16]
Obrázek 2.12. V-statistika denních výnosů Dow Jonesova indexu [16]
Kapitola 2. R/S-analýza finančních aktiv__________________________________________34
Na obrázku 2.11 nejsou mezery pro t 40 , t 250 a t 1250 výrazné, přesto je tam můžeme nalézt. Na obrázku 2.12, který ukazuje výsledný graf V-statistiky, jsou zlomy v trendu pro t 40 a t 1250 zřejmé. Tudíž můžeme vyčlenit cyklus s periodou 1 250 dnů, což těsně odpovídá již dříve zmíněnému cyklu s periodou 4 roky. Zároveň můžeme předpokládat, že existuje cyklus s délkou přibližně 40 dnů [16]. Tento předpoklad bychom mohli ověřit, kdybychom měli například hodinové záznamy Dow Jonesova indexu a prozkoumali je. Zatím však tento cyklus ověřit nemůžeme, avšak nemáme ani důvod k jeho odmítnutí. Bod t 250 je označen na grafu z toho důvodu, že se trend V-statistiky také mění při procházení tohoto bodu. Pro 40 t 250 je náklon grafu shodný s náhodnou procházkou. Pro 250 t 1250 graf ostře roste a pak začíná klesat [16]. Provedla se R/S-analýza Dow Jonesova indexu za období od 2. ledna 1888 do 31. prosince 1991. Objevily se dva cykly s periodou 40 provozních dnů a 1 250 dnů. Tyto cykly je možné použít k technické analýze nebo ke konstrukci vhodného modelu pro testování historických údajů [16]. R/S-analýza je součástí fraktální analýzy, která se používá ke zkoumání chaotických struktur. Tady jsme ukázali jen malou část možností použití R/S-analýzy a jen některé výsledky. Celkový popis všech možností tohoto sice důležitého, ale docela málo známého nástroje analýzy časových řad, je předmětem rozsáhlejšího výzkumu. Edgar Peters provedl ve výzkumu R/S-analýzy obrovský kus práce, jejímž výsledkem jsou publikace Chaos and order in the capital markets : a new view of cycles, prices, and market volatility (1996) [15] a Fractal Market Analysis: Applying Chaos Theory to Investment and Economics (2003) [16]. Avšak nejsou prozkoumané všechny oblasti aplikace R/S-analýzy a různé aspekty problematiky stále čekají na své výzkumníky. V následující kapitole bude provedena R/S-analýza cen akcií sedmi společností zmíněných v první kapitole a pokusí se najít dlouhou paměť daných časových řad, jejich trend a cykly.
Kapitola 3 Použití R/S-analýzy V dané kapitole se provádí R/S-analýza cen akcií sedmi společností. K výpočtu Hurstova exponentu byl napsán program (víz příloha 1) v softwaru Maple 18. Soubory dat mají různé počty prvků, proto se k analýze časových řad používají rozdílné soubory hodnot . Pro společnosti Wal-Mart Stores Inc., McDonald’s Corporation a General Electric jsou hodnoty z souboru (3.1), pro Apple Inc. z (3.2) a pro IBM Corp., BP plc. a Johnson & Johnson z (3.3).
{3,5,7,10,20,25,40,60,80,140,200,400,700,1000,2000,4000,7000}
(3.1)
{3,5,7,10,20,25,40,60,80,140,200,400,700,1000,2000,4000,7000,8500}
(3.2)
{3,5,7,10,20,25,40,60,80,140,200,400,700,1000,2000,4000,7000,9000}
(3.3)
Jak bylo zmíněno v kapitole 2, hodnota poměru R / S se vypočítá jako střední hodnota všech hodnot R / S (t , ) pro daný pár (t , ) . Proto čím víc je počet hodnot t pro dané , tím je přesnější výpočty. Avšak je provádět výpočty pro všechna t časově náročné. Máme zvolit vhodný počet hodnot t tak, aby přesnost byla docela velká a zároveň výpočty trvaly rozumnou dobu. Podle Petersa [16], Mandelbrota & Wallisa [8] a vlastních pokusů byly zvolené následující kroky pro t:
pro časové řady s počtem prvků N 7000,..,9000 (původní časové řady)
N se hodnota R / S počítá pro ti t0 300i , kde t0 1 , i 1,.., ; 300 pro časové řady s počtem prvků N 1500,..,1800 (5-denní přírůstky) N se hodnota R / S počítá pro ti t0 50i , kde t0 1 , i 1,.., , přičemž 50 uvažujeme dvacet čtyři hodnot :
{3,5,7,10,15,20,25,30,40,50,60,70,80,100,140,200,300,400,500,700,750,1000,1200,1500}
pro časové řady s počtem prvků N 380,.., 465 (20-denní přírůstky)
N se hodnota R / S počítá pro ti t0 10i , t0 1 , i 1,.., , přičemž 10 uvažujeme devatenáct hodnot :
{3,5,7,10,12,15,20,25,30,40,50,60,70,80,100,140,200,300,350} 35
Kapitola 3. Použití RS-analýzy__________________________________________________36
K znázornění výsledku taky byly použité funkce softwaru Maple. Kapitola je rozdělena na osm částí, z nichž je sedm věnováno R/S-analýze cen akcií, v osmé pokusíme najít zákonitostí a shrneme výsledky analýzy.
3.1. Analýza cen akcií společnosti Apple Inc. Začneme analýzou cen akcií společnosti Apple Inc. za období od 2. ledna 1981 do 31. října 2014. Celkem máme 8 534 záznamů. Vypočteme hodnotu Hurstova exponentu pro celou časovou řadu a aproximujeme vhodnou přímkou.
Obrázek 3.1. Výsledný graf analýzy cen akcií Apple Inc.
Na obrázku 3.1 jsou nakreslené dvě přímky. Červená přímka aproximuje krabicový graf původní časové řady. Dále časová řada byla náhodně promíchána (k tomu se používá funkce Shuffle z softwaru Maple). Modrá přímka ukazuje trend pro promíchanou časovou řadu. Tím ukazujeme, že pořadí hodnot cen akcií je důležité a dlouhá paměť existuje jen pro toto pořadí. V případě míšení řady dlouhá paměť neexistuje a řada si chová skoro jako náhodná procházka. Lze vidět, že hodnota Hurstova exponentu je hodně velká a rovná se 0.97. Sestavíme novou časovou řadu:
Zi X i 1 X i , i 1,.., N 1 ,
(3.4)
kde N je počet prvků původní časové řady. Časová řada (3.4) reprezentuje denní přírůstky cen akcií. Z obrázku 3.2 je vidět, že se hodnota Hurstova exponentu blíží k hodnotě H pro nezávislý náhodný proces.
Kapitola 3. Použití RS-analýzy__________________________________________________37
Obrázek 3.2.Výsledný graf analýzy denních změn cen akcií Apple Inc.
Zkusíme najít cykly v dané časové řadě. Na obrázku (3.2) změny trendu jsou, ale nejsou výrazné. Prozkoumáme V-statistiku denních, pětidenních a dvacetidenních přírůstků cen akcií Apple Inc. Hodnoty V-statistiky jsou vypočtený podle vzorce (2.28) a zapsány do tabulky 3.1. Tučným písmem jsou zdůrazněny hodnoty lokálních maxim a odpovídajících jím počtů dnů. Na obrázcích 3.3-3.5 jsou znázorněny grafy V-statistiky.
Obrázek 3.3. V-statistika denních změn cen akcií Apple Inc.
Kapitola 3. Použití RS-analýzy__________________________________________________38
Obrázek 3.4. V-statistika 5-denních změn cen akcií Apple Inc.
Obrázek 3.5. V-statistika 20-denních změn cen akcií Apple Inc.
Na obrázku 3.3 vidíme ostrou špičku pro 7 a v tabulce 3.1 a) je zaznamenána lokální maximální hodnota V-statistiky. Analýza časových řad s větším přírůstkem času zřejmě neukáže tento cyklus. Tedy se můžeme spolehnout jen na denní přírůstky v objevení cyklu s takovou malou délkou. 1-denní přírůstky
5-denní přírůstky
20-denní přírůstky
Kapitola 3. Použití RS-analýzy__________________________________________________39 δ 3 5 7 10 20 25 40 60 80 140 200 400 700 1,000 2,000 4,000 7,000 8,500
Dny 3 5 7 10 20 25 40 60 80 140 200 400 700 1,000 2,000 4,000 7,000 8,500
V 0.754 0.895 0.978 0.967 1.053 1.074 1.127 1.168 1.163 1.233 1.317 1.306 1.447 1.499 1.516 1.635 2.223 2.313
δ 3 5 7 10 15 20 25 30 40 50 60 70 80 100 140 200 300 400 500 700 750 1,000 1,200 1,500 1,700
Dny 15 25 35 50 75 100 125 150 200 250 300 350 400 500 700 1,000 1,500 2,000 2,500 3,500 3,750 5,000 6,000 7,500 8,500
V 0.784 0.844 0.881 0.978 0.986 1.035 1.050 1.114 1.165 1.182 1.211 1.217 1.217 1.247 1.282 1.343 1.333 1.318 1.303 1.406 1.437 1.600 1.808 2.135 1.938
δ 3 5 7 10 12 15 20 25 30 40 50 60 70 80 100 140 200 300 350 400 420
Dny 60 100 140 200 240 300 400 500 600 800 1,000 1,200 1,400 1,600 2,000 2,800 4,000 6,000 7,000 8,000 8,400
V 0.779 0.867 0.927 1.006 1.059 1.117 1.089 1.156 1.147 1.184 1.220 1.235 1.216 1.185 1.148 1.196 1.303 1.618 1.807 2.215 1.969
a) b) c) Tabulka 3.1. Hodnoty V-statistik přírůstků cen akcií Apple Inc.
Srovnáním hodnot V-statistiky pro časové řady z denních, pětidenních a dvacetidenních přírůstku najdeme další možné cykly. Jestli cyklus existuje, měl by se objevit ve všech třech případech. Proto lze mluvit o cyklech s délkou přibližně 350 dnů a 8 000 dnů s ohledem, že se délka prvního určí z analýzy pěti a dvaceti denních výnosů. V případě denních přírůstků je vidět, že změna trendu nastává na intervalu [200, 400] , přesnější hodnotu by bylo možné určit zvětšením přesnosti výpočtů. Cyklus s délkou 8 000 dnů bereme v úvahu s ohledem na to, že máme ceny akcií jen za období 34 let. Kdyby naším cílem bylo sestavení modelu chování cen v budoucím období, měli bychom používat údaje o tomto cyklu opatrně. Apple Inc. existuje jen 39 let a bylo by domýšlivě prohlašovat, že třicetiletý cyklus je nesporný fakt.
3.2. Analýza cen akcií společnosti IBM Corp.
Kapitola 3. Použití RS-analýzy__________________________________________________40
Uvažujeme 9 293 záznamů za období od 3. ledna 1978 do 31. října 2014.
Obrázek 3.6. Výsledný graf analýzy cen akcií IBM Corp.
Z grafu na obrázku 3.6 je vidět, že daná časová řada má vysokou hodnotu parametru H. Opět, promíchaná časová řada se blíží k náhodné procházce, což odpovídá myšlence, že pořadí prvků má vliv na dlouhou paměť. Prozkoumáme V-statistiku časových řad sestavených z přírůstků cen za jeden den, pět dnů a dvacet dnů.
Obrázek 3.7. V-statistika denních změn cen akcií IBM Corp.
Kapitola 3. Použití RS-analýzy__________________________________________________41
Obrázek 3.8. V-statistika 5-denních změn cen akcií IBM Corp.
Obrázek 3.9. V-statistika 20-denních změn cen akcií IBM Corp.
Z analýzy V-statistiky denních přírůstků plyne cyklus s délkou 20 dnů. Tento cyklus nemůžeme potvrdit na dalších grafech, proto bereme v úvahu teď. Je možné předpokládat cyklus s délkou 100 až 200 dnů. V-statistika všech tří časových řad mění svou hodnotu na tomto intervalu. Lze se domnívat, že se délka cyklu blíží k hodnotě 150 dnů.
Kapitola 3. Použití RS-analýzy__________________________________________________42 1-denní přírůstky δ Dny V 3 3 0.784 5 5 0.863 7 7 0.908 10 10 0.941 20 20 1.075 25 25 1.049 40 40 1.112 60 60 1.176 80 80 1.248 140 140 1.271 200 200 1.340 400 400 1.323 700 700 1.355 1,000 1,000 1.396 2,000 2,000 1.370 4,000 4,000 1.522 7,000 7,000 1.175 9,000 9,000 1.042
5-denní přírůstky δ Dny V 3 15 0.785 5 25 0.816 7 35 0.882 10 50 0.962 15 75 0.999 20 100 1.087 25 125 1.093 30 150 1.129 40 200 1.109 50 250 1.113 60 300 1.178 70 350 1.183 80 400 1.188 100 500 1.148 140 700 1.158 200 1,000 1.206 300 1,500 1.223 400 2,000 1.229 500 2,500 1.286 700 3,500 1.382 750 3,750 1.391 1,000 5,000 1.299 1,200 6,000 1.158 1,500 7,500 1.026
20-denní přírůstky δ Dny V 3 60 0.777 5 100 0.869 7 140 0.866 10 200 0.937 12 240 0.976 15 300 1.006 20 400 1.026 25 500 1.018 30 600 1.034 40 800 1.065 50 1,000 1.085 60 1,200 1.120 70 1,400 1.121 80 1,600 1.106 100 2,000 1.085 140 2,800 1.199 200 4,000 1.296 300 6,000 1.119 350 7,000 1.043 400 8,000 0.969 450 9,000 0.971
a) b) c) Tabulka 3.2. Hodnoty V-statistik přírůstků cen akcií IBM Corp.
Zřejmě existuje cyklus s délkou přibližně 4 000 provozních dnů. Cykly s délkami 400 a 1 000 dnů plynou jen ze dvou analýz V-statistiky. Ověřit jejich existenci je možné zvětšením přesnosti výpočtů, buď analýzou časové řady sestavenou z desetidenních nebo jiných přírůstků. Tím pádem můžeme předpokládat, že daná časová řada má tři cykly s délkami přibližně 20, 150 a 4 000 dnů.
3.3. Analýza cen akcií společnosti BP plc. Máme ceny akcií BP plc. od 3. ledna 1978 do 31. října 2014, tedy je 9 294 záznamů. Daná časová řada ukazuje hodnotu Hurstova exponentu H 0.97 , tj. je
Kapitola 3. Použití RS-analýzy__________________________________________________43
persistentním procesem. Modrá přímka na obrázku 3.10 aproximuje krabicový diagram R/S-analýzy promíchané časové řady. Tedy vidíme, že pořadí prvků má vliv na persistentní vlastnosti procesu.
Obrázek 3.10. Výsledný graf analýzy cen akcií BP plc.
Prozkoumají se časové řady denních, pětidenních a dvacetidenních přírůstků.
Obrázek 3.11. V-statistika denních změn cen akcií BP plc.
Kapitola 3. Použití RS-analýzy__________________________________________________44
Obrázek 3.12. V-statistika 5-denních změn cen akcií BP plc.
Obrázek 3.13. V-statistika 20-denních změn cen akcií BP plc.
Z analýzy denních přírůstků je vidět cyklus s délkou 20 dnů. Také je možné mluvit o cyklu s délkou 50 až 60 dnů, který plyne z analýzy denních a pětidenních změn cen. V-statistika 20-denních přírůstků tento cyklus zřejmě neodhadne.
Kapitola 3. Použití RS-analýzy__________________________________________________45 1-denní přírůstky δ Dny V 3 3 0.796 5 5 0.862 7 7 0.899 10 10 0.980 20 20 1.133 25 25 1.105 40 40 1.172 60 60 1.185 80 80 1.162 140 140 1.209 200 200 1.210 400 400 1.194 700 700 1.184 1,000 1,000 1.282 2,000 2,000 1.321 4,000 4,000 1.125 7,000 7,000 0.956 9,000 9,000 1.003
5-denní přírůstky δ Dny V 3 15 0.775 5 25 0.882 7 35 0.924 10 50 0.941 15 75 0.920 20 100 0.967 25 125 0.965 30 150 0.976 40 200 1.020 50 250 0.988 60 300 0.998 70 350 1.013 80 400 0.982 100 500 1.009 140 700 1.038 200 1,000 1.102 300 1,500 1.175 400 2,000 1.210 500 2,500 1.201 700 3,500 1.109 750 3,750 1.090 1,000 5,000 1.106 1,200 6,000 1.018 1,500 7,500 0.987
20-denní přírůstky δ Dny V 3 60 0.780 5 100 0.865 7 140 0.901 10 200 0.953 12 240 0.978 15 300 0.971 20 400 1.015 25 500 0.993 30 600 1.025 40 800 1.069 50 1,000 1.137 60 1,200 1.161 70 1,400 1.190 80 1,600 1.204 100 2,000 1.232 140 2,800 1.253 200 4,000 1.149 300 6,000 1.054 350 7,000 1.000 400 8,000 1.020 450 9,000 0.987
a) b) c) Tabulka 3.3. Hodnoty V-statistik přírůstků cen akcií BP plc.
Je možné předpokládat cykly s délkami 200 až 240 provozních dnů a přibližně 2 000-2 500 dnů. Ověřit cyklus a zjistit přesnou hodnotu je možné zvýšením přesností výpočtu V-statistiky, což se zatím neprovádí z důvodu velké časové náročnosti. Tím pádem pro danou časovou řadu předpokládáme cykly s periodami 20, 60, 200 a 2 500 dnů. Cykly s větší délkou jsme zatím neobjevili, ale to neznamená, že neexistují. Je možné, že jsou cykly s rozsahem víc než máme počtů pozorování v daném výzkumu.
Kapitola 3. Použití RS-analýzy__________________________________________________46
3.4. Analýza cen akcií společnosti Wal-Mart Stores Inc. K analýze akcií společnosti Wal-Mart Stores Inc. používáme ceny akcií od 3. ledna 1984 do 31. října 2014. Celkem je 7 776 záznamů. Proveďme R/S-analýzu celé časové řady.
Obrázek 3.14. Výsledný graf analýzy cen akcií Wal-Mart Stores Inc.
Na obrázku 3.14 je červená přímka aproximací krabicového diagramu původní časové řady, je modrá přímka aproximací krabicového diagramu promíchané časové řady. Hodnota Hurstova exponentu se blíží k jedničce a rovna se 0.96. Stejně jak v předchozích případech se analýzují časové řady sestavené z denních, pětidenních a dvacetidenních přírůstků.
Kapitola 3. Použití RS-analýzy__________________________________________________47 Obrázek 3.15. V-statistika denních změn cen akcií Wal-Mart Stores Inc.
Obrázek 3.16. V-statistika 5-denních změn cen akcií Wal-Mart Stores Inc.
Obrázek 3.17. V-statistika 20-denních změn cen akcií Wal-Mart Stores Inc.
Z analýzy V-statistiky denních přírůstků opět je vidět cyklus s délkou 20 provozních dnů. Ve všech třech případech se opakují cykly s délkami přibližně 200 až 240 dnů a 4 000 až 5 000 dnů.
Kapitola 3. Použití RS-analýzy__________________________________________________48 1-denní přírůstky δ Dny V 3 3 0.787 5 5 0.851 7 7 0.908 10 10 0.991 20 20 1.100 25 25 1.092 40 40 1.116 60 60 1.170 80 80 1.213 140 140 1.198 200 200 1.235 400 400 1.148 700 700 1.151 1,000 1,000 1.176 2,000 2,000 1.183 4,000 4,000 1.316 7,000 7,000 1.087
5-denní přírůstky δ Dny V 3 15 0.772 5 25 0.843 7 35 0.865 10 50 0.946 15 75 0.945 20 100 0.974 25 125 1.022 30 150 1.060 40 200 1.121 50 250 1.093 60 300 1.065 70 350 1.057 80 400 1.081 100 500 1.060 140 700 1.098 200 1,000 1.116 300 1,500 1.130 400 2,000 1.198 500 2,500 1.252 700 3,500 1.425 750 3,750 1.435 1,000 5,000 1.448 1,200 6,000 1.388 1,500 7,500 1.511
20-denní přírůstky δ Dny V 3 60 0.782 5 100 0.861 7 140 0.881 10 200 0.918 12 240 0.944 15 300 0.933 20 400 0.918 25 500 0.946 30 600 0.950 40 800 0.980 50 1,000 1.018 60 1,200 1.044 70 1,400 1.068 80 1,600 1.067 100 2,000 1.128 140 2,800 1.191 200 4,000 1.309 300 6,000 1.188 350 7,000 1.119
a) b) c) Tabulka 3.4. Hodnoty V-statistik přírůstků cen akcií Wal-Mart Stores Inc.
Tím pádem můžeme závažně tvrdit, že daná časová řada má minimálně dva cykly s periodami 200 provozních dnů a 4 000 provozních dnů, což je téměř deset měsíců a šestnáct let. Taky můžeme předpokládat cyklus délky 20 dnů. Ověřit existenci tohoto cyklu je možné zvětšením přesnosti výpočtů a analýzou časové řady třeba čtyřhodinových přírůstků.
3.5. Analýza cen akcií společnosti McDonald’s Corporation K analýze ceny akcií McDonald’s Corporation máme 7 776 záznamů za období od 3. ledna 1984 do 31. října 2014. Opět začneme R/S-analýzou dané časové řady a
Kapitola 3. Použití RS-analýzy__________________________________________________49
taktéž ukažme trend pro promíchanou časovou řadu (modrá přímka). Zase máme velkou hodnotu parametru H 0.97 , což znamená, že máme proces s dlouhou pamětí.
Obrázek 3.18. Výsledný graf analýzy cen akcií McDonald’s Corporation
V-statistika časových řad je znázorněna na obrázcích 3.19-3.21 a její hodnoty jsou zapsány do tabulky 3.5.
Obrázek 3.19. V-statistika denních změn cen akcií McDonald’s Corporation
Kapitola 3. Použití RS-analýzy__________________________________________________50
Obrázek 3.20. V-statistika 5-denních změn cen akcií McDonald’s Corporation
Obrázek 3.21. V-statistika 20-denních změn cen akcií McDonald’s Corporation
Z analýzy V-statistiky denních přírůstků je možné předpokládat cykly s délkami 7 a 20 provozních dnů. Ověřit je by bylo možné v případě analýzy časové řady sestavenou z menších přírůstku, například hodinových nebo čtyřhodinových.
Kapitola 3. Použití RS-analýzy__________________________________________________51 1-denní přírůstky δ Dny V 3 3 0.775 5 5 0.903 7 7 1.012 10 10 1.009 20 20 1.123 25 25 1.077 40 40 1.072 60 60 1.144 80 80 1.118 140 140 1.092 200 200 1.132 400 400 1.098 700 700 1.135 1,000 1,000 1.150 2,000 2,000 1.232 4,000 4,000 1.475 7,000 7,000 1.500
5-denní přírůstky δ Dny V 3 15 0.774 5 25 0.844 7 35 0.866 10 50 0.948 15 75 0.947 20 100 0.975 25 125 1.022 30 150 1.060 40 200 1.121 50 250 1.093 60 300 1.065 70 350 1.057 80 400 1.082 100 500 1.060 140 700 1.098 200 1,000 1.115 300 1,500 1.130 400 2,000 1.197 500 2,500 1.253 700 3,500 1.425 750 3,750 1.435 1,000 5,000 1.449 1,200 6,000 1.388 1,500 7,500 1.511
20-denní přírůstky δ Dny V 3 60 0.786 5 100 0.874 7 140 0.933 10 200 0.966 12 240 0.961 15 300 0.988 20 400 0.999 25 500 1.028 30 600 1.046 40 800 1.095 50 1,000 1.105 60 1,200 1.102 70 1,400 1.108 80 1,600 1.123 100 2,000 1.202 140 2,800 1.312 200 4,000 1.452 300 6,000 1.410 350 7,000 1.515
a) b) c) Tabulka 3.5. Hodnoty V-statistik přírůstků cen akcií McDonald’s Corporation
Analýza V-statistik denních a pětidenních přírůstků ukazuje cyklus délky 50 až 60 dnů. Není možné ověřit jej pomocí V-statistiky 20 denních přírůstků. Proto se tento cyklus bere v úvahu, ale s ohledem, že by bylo možné ověřit stejným způsobem jak 7denní cyklus. Je zřejmá existence cyklu s délkou 200 provozních dnů. Cyklus délky 4 000 až 5 000 provozních dnů plyne z analýzy pětidenních a dvacetidenních přírůstků. Ačkoliv se neobjevil při analýze denních přírůstku, ale to může souviset s přesností výpočtů. Všimněme si, že na obrázcích 3.19-3.21 graf V-statistiky roste na posledním úseku. Můžeme předpokládat, že daná časová řada má ještě vetší cyklus, který zatím nemůžeme objevit z důvodů nedostatečného počtu pozorování.
Kapitola 3. Použití RS-analýzy__________________________________________________52
3.6. Analýza cen akcií společnosti General Electric Uvažujeme ceny akcií General Electric od 3. ledna 1984 do 31. října 2014, celkem je 7 777 cen. Prozkoumáme danou časovou řadu stejně jak předchozí a prověříme, jaké vlastnosti ukáže R/S-analýza. Za prvé vidíme vysokou persistenci procesu, hodnota Hurstova exponentu je H 0.96 . Ověříme, že pořadí prvků je důležité a po promíchání řady se vyskytují vlastností náhodné procházky.
Obrázek 3.22. Výsledný graf analýzy cen akcií General Electric
Obrázek 3.23. V-statistika denních změn cen akcií General Electric
Kapitola 3. Použití RS-analýzy__________________________________________________53
Obrázek 3.24. V-statistika 5-denních změn cen akcií General Electric
Obrázek 3.25. V-statistika 20-denních změn cen akcií General Electric
Stejně jako v předchozích případech z analýzy denních přírůstků plyne existence cyklu délky 20 dnů. Je možné ten ověřit jen analýzou časové řady sestavenou z menších přírůstků cen.
Kapitola 3. Použití RS-analýzy__________________________________________________54
Cyklus délky 60 dnů, jenž je vidět v tabulce 3.6 a), neobjevuje se dále, proto se ten neuvažuje. Avšak cyklus s délkou 140 až 200 provozních dnů je zřejmý. 1-denní přírůstky δ Dny V 3 3 0.790 5 5 0.889 7 7 0.927 10 10 1.050 20 20 1.091 25 25 1.063 40 40 1.051 60 60 1.109 80 80 1.102 140 140 1.162 200 200 1.159 400 400 1.130 700 700 1.212 1,000 1,000 1.318 2,000 2,000 1.489 4,000 4,000 1.670 7,000 7,000 1.759
5-denní přírůstky δ Dny V 3 15 0.793 5 25 0.837 7 35 0.914 10 50 0.963 15 75 1.001 20 100 1.023 25 125 1.030 30 150 1.061 40 200 1.048 50 250 1.066 60 300 1.082 70 350 1.090 80 400 1.062 100 500 1.078 140 700 1.117 200 1,000 1.196 300 1,500 1.334 400 2,000 1.403 500 2,500 1.410 700 3,500 1.500 750 3,750 1.526 1,000 5,000 1.580 1,200 6,000 1.639 1,500 7,500 1.702
20-denní přírůstky δ Dny V 3 60 0.769 5 100 0.859 7 140 0.928 10 200 0.942 12 240 0.941 15 300 0.991 20 400 0.953 25 500 0.984 30 600 1.017 40 800 1.090 50 1,000 1.141 60 1,200 1.174 70 1,400 1.231 80 1,600 1.288 100 2,000 1.355 140 2,800 1.407 200 4,000 1.582 300 6,000 1.625 350 7,000 1.649
a) b) c) Tabulka 3.6. Hodnoty V-statistik přírůstků cen akcií General Electric
Je možné se domnívat, že taky existuje cyklus s délkou 300 až 350 dnů. Objevil se při analýze pětidenních a dvacetidenních přírůstků cen. Bere se v úvahu s ohledem na kontrolu zvětšením přesnost R/S-analýzy denních změn cen akcií. Není vidět dlouhodobé cykly. Avšak grafy V-statistik ve všech případech rostou pro velké hodnoty a je možné předpokládat existenci cyklu délky vetší, než je počet pozorování v dané časové řadě. Tím pádem se předpokládají cykly s periodami přibližně 20, 150 a 350 dnů.
Kapitola 3. Použití RS-analýzy__________________________________________________55
3.7. Analýza cen akcií společnosti Johnson & Johnson Uvažujeme 9 293 záznamů za období od 3. ledna 1978 do 31. října 2014. Na obrázku 3.26 jsou nakreslené přímky aproximací R/S-analýzy původní časové řady a promíchané. Znovu vidíme, že pořadí prvků má vliv na dlouhou paměť procesu. Hodnota Hurstova exponentu se blíží k jedničce a je H 0.97 . Časová řada ukazuje vysokou úroveň persistence
Obrázek 3.26. Výsledný graf analýzy cen akcií Johnson&Johnson
Obrázek 3.27. V-statistika denních změn cen akcií Johnson&Johnson
Kapitola 3. Použití RS-analýzy__________________________________________________56
Obrázek 3.28. V-statistika 5-denních změn cen akcií Johnson&Johnson
Obrázek 3.29. V-statistika 20-denních změn cen akcií Johnson&Johnson
Cyklus délky 40 dnů se objevuje při analýze denních přírůstků. Je jasné, že se neobjeví při analýze 20-denních změn, ale mohl by plynout z analýzy pětidenních přírůstků. Avšak se tam neobjevuje. Proto se tento cyklus nebere v úvahu za současné podmínky. Nicméně by bylo možné ten ověřit zvětšením přesnosti výpočtů nebo analýzou časové řady sestavenou z menších časových přírůstků.
Kapitola 3. Použití RS-analýzy__________________________________________________57 1-denní přírůstky δ Dny V 3 3 0.781 5 5 0.852 7 7 0.906 10 10 0.988 20 20 1.106 25 25 1.128 40 40 1.170 60 60 1.121 80 80 1.147 140 140 1.175 200 200 1.182 400 400 1.146 700 700 1.205 1,000 1,000 1.211 2,000 2,000 1.169 4,000 4,000 1.230 7,000 7,000 1.111 9,000 9,000 1.000
5-denní přírůstky δ Dny V 3 15 0.790 5 25 0.868 7 35 0.889 10 50 0.943 15 75 0.952 20 100 0.992 25 125 1.039 30 150 1.055 40 200 1.092 50 250 1.080 60 300 1.073 70 350 1.069 80 400 1.074 100 500 1.090 140 700 1.134 200 1,000 1.166 300 1,500 1.156 400 2,000 1.119 500 2,500 1.110 700 3,500 1.143 750 3,750 1.161 1,000 5,000 1.198 1,200 6,000 1.097 1,500 7,500 1.049
20-denní přírůstky δ Dny V 3 60 0.771 5 100 0.864 7 140 0.917 10 200 0.979 12 240 0.962 15 300 1.003 20 400 1.037 25 500 1.040 30 600 1.062 40 800 1.097 50 1,000 1.103 60 1,200 1.114 70 1,400 1.108 80 1,600 1.093 100 2,000 1.087 140 2,800 1.104 200 4,000 1.199 300 6,000 1.159 350 7,000 1.104 400 8,000 1.055 450 9,000 0.959
a) b) c) Tabulka 3.7. Hodnoty V-statistik přírůstků cen akcií Johnson&Johnson
Na rozdíl od předchozích případů ostatní cykly je vidět při analýze všech tří časových řad. Jsou cykly s délkami přibližně 200, 1 200 a 4 000 provozních dnů.
3.8. Zobecnění V dané kapitole byly prozkoumány časové řady reprezentující ceny akcií sedmi společností působících na mezinárodním trhu. Každá z časových řad má vysokou hodnotu parametru H. Hurstův exponent ve všech případech se blíží k jedničce, což ukazuje vysokou persistenci procesů. V praxe velká hodnota parametru H ukazuje, že řada má trend a dlouhou paměť. Sice není možné předpovědět cenu akcií za deset nebo
Kapitola 3. Použití RS-analýzy__________________________________________________58
dvacet let, ale je možné zmodelovat chování řady a odhadnout možné výnosy nebo ztráty. Takže byly prozkoumány i podmnožiny daných časových řad délky 5-6 let. Sice během různých period Hurstův exponent nabývá odlišné hodnoty, ale jsou téměř rovny jedničce.
Obrázek 3.30. Hurstův exponent pro pětiletá období
Ve všech prozkoumaných časových řadách se objevily cykly, jež je možné rozdělit na tři třídy: krátkodobé (míň než jeden kalendářní rok), střednědobé (1 až 10 kalendářních roků) a dlouhodobé (víc než 10 let). Znázornění výsledků je v následující tabulce: Délka cyklů (v provozních dnech) Název společnosti
krátkodobé
střednědobé
Dlouhodobé
Hurstův exponent
Apple Inc. IBM Corp. BP plc. Wal-Mart Stores Inc.
7 20, 150 20, 60, 200 20, 200 7, 20, 60, 200 20, 150 200
350 2 500 -
8 000 4 000 4 000
0.97 0.96 0.97 0.96
-
5 000
0.97
350 1 200
4 000
0.96 0.97
McDonald’s Corporation General Electric Johnson&Johnson
Tabulka 3.8. Délky cyklů časových řad a hodnoty Hurstova exponentu
Kapitola 3. Použití RS-analýzy__________________________________________________59
Z tabulky 3.8 je vidět, že krátkodobé cykly mají vše časové řady, dlouhodobý cyklus nemají jen akcie BP plc. a General Electric. Nelze závažně tvrdit, že dané časové řady zcela nemají dlouhodobé cykly. Jsou možné cykly s délkou víc než 30 prozkoumaných let. Je vidět, že se střednědobé cykly objevily jen u čtyř korporací ze sedmi. Přičemž tyhle společnosti působí v různých oblastech. Tudíž cykly a jejich délky nezáleží na oboru firmy. V části 3.1 bylo zmíněno, že sice se objevil u společnosti Apple Inc. cyklus s periodou 30 let, ale je nutné pamatovat, že se zkoumají údaje jen za období 34 let. Avšak ostatní procesy ukazují cykly s délkou 4 000 provozních dnů, což je téměř 16 let. Tato doba odpovídá ekonomickému cyklu, tzv. Kuznetsovu cyklu, který má rozsah délky 15 až 25 let. Můžeme předpoklad, že změna cen akcií v určité míře odráží současnou ekonomickou situaci ve světě.
Závěr Cílem dané diplomové práce bylo seznámení s R/S-analýzou dlouhodobých časových řad. Tato metoda je založená na hledání poměru mezi normovaným rozsahem časové řady a její směrodatnou odchylkou. Hodnota poměru je úměrná délce časové řady umocněné na parametr H, jenž se nazývá Hurstovým exponentem a ukazuje persistentní vlastnost časové řady. Jinými slovy, pomocí Hurstova exponentu lze usoudit, jestli časová řada má dlouhou paměť. Existence dlouhé paměti pomáhá pochopit chování časové řady, což je důležité pro investory, jestli zkoumaná časová řada je pozorováním cen akcií nebo kurzů měn. Ve třetí kapitole jsme zkusili použít teoretickou informaci z druhé kapitoly a vypočetli hodnotu Hurstova exponentu pro sedm dlouhodobých časových řad. Taky jsme se pokusili najít cykly v změnách cen a interpretovat je. Sekundárním cílem diplomové práce se bylo naučit programovat v softwaru Maple 18. Tento cíl je dosažen, potvrzením je programovací kód pro výpočet hodnoty parametru H (víz příloha 1) a různé druhy obrázků, nakreslených pro danou práci. Ovšem existují další neprostudované možnosti tohoto softwaru. Nicméně ovládání Maple 18 bylo natolik zajímavým, že lze předpokládat další použití mozností softwaru ve studiu a práci.
60
Seznam použité literatury [1]
Butakov V., Grakovsky A. Оценка уровня стохастичности временных рядов произвольного происхождения при помощи показателя Херста. Computer Modelling and New Technologies, 2005, Vol.9, No.2, 27-32.
[2]
Fama E., French K. Permanent and Temporary Components of Stock Prices. The Journal of Political Economu, Vol. 96, No. 2 (Apr., 1988), 246-273.
[3]
Chernushkin V. Использование среды Maple для решения задач квантовой механики. – [on-line], Воронеж, 2005. Dostupné [13.09.2014] z http://tchernouchkine.narod.ru/maple/docums/maple_tchern_web.pdf
[4]
Kolář M. Stochastické procesy ve financní matematice. Učební text.
[5]
Krištoufek L., Vošvrda M. Efektivita kapitálových trhů: fraktální dimenze, Hurstův exponent a entropie. Politická ekonomie №2, 2012, 208-221.
[6]
Mandelbrot B. Statistical Methodology For Nonperiodic Cycles: From The Covariance To R/S-Analysis. Annals of Economic and Social Measurement, Volume 1, number 3, 1972, 259-290
[7]
Mandelbrot B., Hudson R. (Не)послушные рынки. Фрактальная революция в финансах. – М.: Издательский дом «Вильямс», 2006, 35-86, 137-138, 196197, 217-250.
[8]
Mandelbrot B., Wallis J. Computer experiments with fractional Gaussian noises. Water Resources Research: 5, 1969, 228-241.
[9]
Mandelbrot B., Wallis J. Computer experiments with fractional Gaussian noises. Water Resources Research: 3, 1969, 260-267.
[10]
Naiman E. Расчет показателя Херста с целью выявления трендовости (персистентности) финансовых рынков и макроэкономических показателей. Економіст № 10, 2009. 18-28.
[11]
Otevřená encyklopedie Wikipedie [on-line]. z http://cs.wikipedia.org/wiki/Finanční_trh
Dostupné
[10.10.2014]
[12]
Otevřená encyklopedie Wikipedie [on-line]. z https://cs.wikipedia.org/wiki/Technická_analýza
Dostupné
[15.11.2014]
[13]
Otevřená encyklopedie Wikipedie z http://cs.wikipedia.org/wiki/Nil
Dostupné
[24.03.2014]
[on-line].
61
Seznam použité literatury______________________________________________________62
[14]
Otevřená encyklopedie Wikipedie [on-line]. Dostupné z https://cs.wikipedia.org/wiki/Teorie_efektivních_trhů
[15]
Peters E. Chaos and order in the capital markets : a new view of cycles, prices, and market volatility. New York: Wiley, 1996.
[16]
Peters E. Фрактальный анализ финансовых рынков: применение теории Хаоса в инвестициях и экономике. – М.: Интернет-трейдинг, 2004, 22-124.
[17]
Prochorov G., Kolbeev V., Zhelnov K., Ledenev M. Математический пакет Maple V Release 4: Руководство пользователя. – Калуга: Облиздат, 1998.
[18]
Savotchenko S., Kuzmicheva T. Методы решения математических задач в Maple: Учебное пособие – Белгород: Изд. Белаудит, 2001.
[19]
Schroeder M. Фракталы, хаос, степенные законы. Миниатюры из бесконечного рая. – Ижевск: НИЦ «Регулярная и хаотическая динамика», 2001.
[20]
Strygin A. Анализ фрактальных свойств финансово-экономических процессов в экономике РФ. – Санкт-Петербург, 2004. Diplomová práce. Státní univerzita, Výsoká Škola Ekonomická.
[21]
Teplov S. Применение R/S анализа на фондовых рынках. Финансы и Бизнес №1, 2008, 130-137.
[22]
Zhdanov I. Выбор акций на основе показателя Херста. Фрактальный анализ Российского фондового рынка [on-line], 2010. Dostupné [12.04.2014] z http://www.beintrend.ru/2010-09-01-13-03-40
[24.03.2014]
Příloha 1
>
>
> > >
63
Příloha 1____________________________________________________________________64