1
Z2069 Statistické metody a zpracování dat II
Analýza časových řad
• vývoj cen akcií
• objem obchodování na burze
• průměrný roční odtok vody z povodí
• vývoj počtu obyvatelstva určité lokality
• maximální denní srážkové úhrny na určité stanici
Příklady časových řad a jejich použití
Obchodní den
0
50
100
150
200
250
300
2 4 6 8 10 12 14 16 18 20
Objem obchodu
(intervalová řada)
Kurz akcie
(okamžiková řada)
http://www.cru.uea.ac.uk/cru/info/warming/
Základní pojmy
nyyy ,,, 21 L
ty , kde t=1, 2, …, n
y = ukazatel
t = časová proměnná
n = počet členů řady
( )tfyt =
Časová řada je chronologicky uspořádaná posloupnost
hodnot určitého statistického ukazatele.
Pomocí časových řad můžeme zkoumat dynamiku
jevů v čase.
Mají základní význam pro analýzu příčin, které na
tyto jevy působily a ovlivňovaly jejich chování v
minulosti, tak pro předvídání jejich budoucího vývoje.
…
• Problém volby časových bodů pozorování
• Problémy s délkou časové řady
• Problémy s kalendářem
• Problémy s nesrovnatelností jednotlivých měření
Problémy při sestavování časových řad
Uvedené problémy
mohou vést k narušení
homogenity časové řady
b
0
10
20
30
40
50
60
70
80
90
1921 1931 1941 1951 1961 1971 1981 1991
a
10
15
20
25
30
35
40
45
1921 1931 1941 1951 1961 1971 1981 1991
Maximální denní nárazy větru a počty dnů s
nárazy větru na stanici Praha, Karlov v období
1921-1990
Transformace časové řady
Jedná se o úpravu původní časové řady, tak aby
1. splňovala podmínky pro následnou analýzu (např.
stacionarita atd.)
2. zvýrazňovala dále analyzovanou složku
• přidání konstanty y = y + C
• linearizace řady y = ln(y)
• odečtení průměru
• standardizace
• odečtení hodnot trendové funkce (viz.
stacionarita)





 −
=
ds
yy
y
yyy −=
Běžné druhy transformací:
Základní typy časových řad
Časové řady deterministické - neobsahují prvek náhody
(sin(x)) a stochastické (realizace náhodného procesu)
Časové řady absolutních veličin (přímo zjišťovaných)
• okamžikové (počet obyvatel – k datu sčítání)
• intervalové (denní úhrn srážek)
Časové řady odvozené
• průměrných veličin (řada klouzavých průměrů)
• poměrných – relativních veličin (řada hektarových výnosů)
Časové řady ekvidistantní a neekvidistantní
2
Okamžikové časové řady
Jsou spojité v čase, záleží u nich na rozhodném okamžiku šetření.
Hodnota nezávisí na délce intervalu, za který je znak zjišťován.
Okamžikové ukazatele za několik intervalů nesčítáme. Je však pro
ně typické počítání průměrů v čase. Průměr okamžikové veličiny za
určité období označujeme jako tzv. chronologický průměr.
Nejprve spočteme průměr za časové okamžiky ti-1 a ti, pro i=2 až n.
Z těchto hodnot určíme průměr pro celou řadu:
Uvedený vztah platí v případě, že délka všech intervalů je
konstantní. Pokud ne, je nutné jednotlivé dílčí průměry vážit
délkami intervalů a vypočítat vážený chronologický průměr.
1
2
1
...
2
1
121
−
++++
=
−
n
yyyy
y
nn
Intervalové časové řady
Jednotlivé hodnoty se vztahují k časovým úsekům a přímo závisí na
jejich délce.
Hodnotu intervalového ukazatele zjištěnou za časový interval (ti-1, ti)
označme qi a přiřazujeme-ji ke středu časového intervalu.
Časovou řadu hodnot qi označujeme intervalovou řadou běžných
hodnot.
Za delší časové období lze intervalové ukazatele shrnovat a vytvářet
součtové (kumulativní) řady. Součtová řada vznikne postupným
sčítáním hodnot za sebou jdoucích časových intervalů. Podle průběhu
součtové řady můžeme posoudit rovnoměrnost vývoje hodnot znaku.
Požadavkem sestavování intervalových časových řad je konstantnost
délky časového intervalu. V řadě případů tento požadavek není splněn
(např. počet dnů v měsíci).
Dalším typem součtových časových řad jsou řady klouzavých úhrnů.
Jsou vhodné ke srovnání úrovně řady ve sledovaném období s úrovní
řady období předešlého.
Odvozené ukazatele časové řady
Při práci s časovými řadami je typické, že často
pracujeme ne přímo s původní časovou řadou, ale
s nějakou její transformací.
Absolutní přírůstek (první diference)
Jsou-li členy v řadě absolutních přírůstků prakticky
konstantní, potom hodnoty řady lineárně rostou (klesají).
Relativní přírůstek
Informuje nás o rychlosti (tempu) růstu
1−−=∆ ttt yyy
1
11
1
1
−=
−
=
∆
=
−−
−
− i
i
i
ii
i
i
i
y
y
y
yy
y
y
δ
Odvozené ukazatele časové řady
Koeficient růstu (řetězový index): vyjadřuje, o kolik
procent vzrostla hodnota časové řady v okamžiku ti
ve srovnání s hodnotou řady v čase ti-1.
Průměrný koeficient růstu: pro celou řadu se
vypočte jako geometrický průměr jednotlivých
hodnot koeficientů růstu.
Uvedený výpočet je vhodný pouze v případě stálého a
přibližně stejného růstu hodnot řady.
(%)1001
1
⋅=+=
−i
i
ii
y
y
k δ
1
1
1
12
3
1
21
121 ...... −−
−
−
− =⋅⋅=⋅= n
n
n
n
nn
n
y
y
y
y
y
y
y
y
kkkk
Odvozené ukazatele časové řady
Pro účely srovnání různých časových řad se jejich
hodnoty převádějí na tzv. bazické indexy (indexy
se stálým základem):
Hodnota yz je obvykle prvním nebo posledním
členem časové řady (základ).
(%)100'
⋅=
z
i
i
y
y
k
Odvozené ukazatele časové řady
• Absolutní přírůstek =D8-D7
• Relativní přírůstek =E8/D7*100
• Řetězový index =D8/D7*100
• Bazický index =D8/$D$7*100
3
Odvozené ukazatele časové řady Základy analýzy časových řad
Hlavní cíle analýzy časových řad
1. odhalení zákonitostí a příčin
dosavadního vývoje
2. prognóza chování časových
řad
Každá řada může obsahovat čtyři základní složky:
• trend (Tt)
• periodická (sezónní) složka (St)
• cyklická složka (Ct)
• náhodná složka (εt)
První tři složky tvoří systematickou část řady.
Trendová složka časové řady
• Trend je obecná tendence vývoje zkoumaného jevu za
dlouhé období.
• Je výsledkem dlouhodobých a stálých procesů (v
měřítku posuzované délky časové řady).
• Trend může být lineární či nelineární.
• Trend může být rostoucí, klesající nebo může
existovat řada bez trendu (s nulovým trendem).
• Časové řady bez trendu se označují jako stacionární.
Periodická složka časové řady
• Periodická složka je pravidelně se opakující odchylka od
trendové složky s pevnou délkou periody T.
• Perioda této složky je menší než celková velikost sledovaného
období.
• Typickým případem jsou sezónní kolísání a nebo řady denních,
měsíčních, čtvrtletních ukazatelů.
• Příčiny sezónnosti jsou různé, většinou však dobře definovatelné.
• Sezónnost je typická pro časové řady ekonomických ukazatelů.
)()( Ttftf ii +=
Cyklická složka
• Cyklická složka udává kolísání okolo trendu v důsledku
dlouhodobého cyklického vývoje.
• Cyklická složka může vykazovat změny v délce a amplitudě cyklu.
• Délka cyklu je tedy většinou neznámá. (př. demografický trend,
kolísání teploty vzduchu).
• Délka cyklu je tedy delší než 1 rok. V některých případech se
označuje jako „střednědobý trend“.
• Bývá typickou součástí časových řad meteorologických prvků (př.
problém globálního oteplování) či hydrologických jevů.
6,0
7,0
8,0
9,0
10,0
11,0
12,0
1961 1966 1971 1976 1981 1986
)()( Ttftf ii +≈
Náhodná složka časové řady
• Náhodná (stochastická) složka se nedá popsat
žádnou funkcí času.
• "Zbývá" po vyloučení trendu, sezónní a cyklické
složky.
• Jejím zdrojem jsou v jednotlivostech
nepostižitelné jevy.
• Lze ji však popsat pravděpodobnostně.
4
• Prvotní analýza spočívá v grafickém znázornění průběhu řady.
• Graf slouží k prvotnímu posouzení tendence změn či k hledání
opakujících se jevů („patterns“).
• I tyto jednoduché metody umožňují velmi krátkodobou
předpověď.
• Graf však velmi dobře může znázorňovat nehomogenity,
porovnávat dvě či více řad mezi sebou, …
• Slouží k výběru vhodné metody analýzy.
Grafické metody analýzy časových řad
Index severoatlantské cirkulace (NAOI), XII-II
-4
-3
-2
-1
0
1
2
3
1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000
Vývoj kurzu akcií – příklad výskytu jednoduchých obrazců
(patterns) v časové řadě
Grafické metody analýzy časových řad
Klasický (formální) model analýzy časových řad
• Klasický model je pouze popisem jednotlivých složek
časové řady jako forem pohybu, ne poznáním příčin.
• Jedná se o dekompozici na jednotlivé složky a jejich
formální popis např. tzv. aditivním modelem:
• Základem je popis systematické složky (trendu,
cyklických a periodických kolísání).
ttttttt CSTYy εε +++=+=
Model časové řady s aditivní sezónní složkou
Analýza trendu
A. Klasický přístup založený na matematickostatistickém
modelování. Modelované parametry jsou
KONSTANTNÍ v čase. Neadaptivní metody – např.
regresní modely. Umožňují snadnou předpověď.
B. Adaptivní přístup – parametry se v čase VYVÍJEJÍ.
Například charakter lineárního trendu se mění (mění
se směrnice trendu). Za jednoduchou adaptivní
metodu lze považovat i metodu klouzavých průměrů
(viz. dále).
Analýza trendu – základní metody
vyrovnávání:
• analytické (popis časové řady funkcí)
• mechanické (klouzavé průměry)
• …
analytické vyrovnání trendu
Analytické vyrovnávání trendu
matematickou křivkou
• Patří mezi neadaptivní metody. Vychází z
předpokladu, že se trend po celou sledovanou dobu
nemění a že je možné ho popsat některým typem
matematické křivky.
• Identifikace trendu se redukuje na výběr správného
typu matematické křivky a odhad jejích parametrů.
• Na problém analýzy trendu lze pohlížet jako na
speciální případ regresní závislosti, kdy nezávisle
proměnnou je čas.
• Časovou řadu vyrovnáváme křivkou, která nejlépe
vystihuje její vývojový trend. Výpočet parametrů křivky
se děje metodou nejmenších čtverců.
ttt ETry +=
5
Lineární trend
Parametr b1 představuje přírůstek hodnoty y
připadající na jednotkovou změnu časové proměnné.
Řada se vyznačuje konstantními absolutními přírůstky
(první diference).
tbbyt 10 +=
2
1
2
1 1
1
tnt
ytty
b n
t
n
t
n
t
tt
−
−
=
∑
∑ ∑
=
= =
Lineární trend
Hodnoty parametrů b0 a b1 získáme metodou
nejmenších čtverců obdobně jako v případě jednoduché
lineární regrese, tedy:
tbyb 10 −=
Předpověď budoucí hodnoty (bodová předpověď)
má tvar:
TbbyT 10
ˆ +=
Exponenciální trend
Parametr b1 představuje průměrný přírůstek hodnot yt.
Ty se chovají jako členy geometrické posloupnosti.
Protože se již nejedná o funkci lineární v parametrech,
lze k odhadu exponenciálního trendu využít metody
nejmenších čtverců pouze po její logaritmické
transformaci:
t
t bby 10 ⋅=
10 logloglog btbyt ⋅+=
Polynomický trend
Při volbě stupně polynomu je třeba postupovat
opatrně. Vyšší stupeň zajišťuje těsnější proložení
empirických hodnot křivkou, vede ale k nestabilitě
trendu.
Vyšší polynomy se většinou vůbec nehodí k
extrapolacím.
K odhadu parametrů lze využít MNČ.
k
kt tbtbtbby ++++= ...2
210
Polynomický trend
Rozdíly v proložení časové řady polynomem 2. (nahoře) a 6. (dole) stupně
Polynomický trend
Proložení časové řady polynomem 5 stupně v programu Statistica
Grafy – 2 D bodové grafy
6
Logistická křivka
Křivka má tři úseky, první je charakterizován
pozvolným vzestupem, druhá v okolí inflexního bodu
prudkým růstem a třetí určitou vrcholovou stagnací.
(patří mezi tzv. S-křivky).
tt
bbk
y
10
1
⋅+
=
Verifikace modelu
Je zapotřebí zhodnotit statistickou významnost
odhadnutých parametrů modelu i modelu jako celku.
MNČ – podstatou je, že model vždy vysvětlí pouze část
variability (proměnlivosti) pozorovaných dat.
Je nutné zjistit (testovat), zda model jako celek dává
lepší vysvětlení, než je možné očekávat jako důsledek
náhody a to na jisté hladině významnosti.
Koeficient determinance R2 – základní ukazatel
vhodnosti použitého modelu (vzorec a interpretace viz.
korelační počet)
Analýza rozptylu
A. Rozptyl empirických hodnot
(celkový)
B. Rozptyl vyrovnaných hodnot
(modelový)
C. Rozptyl reziduální ∑=
− −=
n
t
tiyy yy
n
s
1
22
ˆ )ˆ(
1
∑=
−=
n
t
ty yy
n
s
1
22
)(
1
∑=
−=
n
t
ty yy
n
s
1
22
ˆ )ˆ(
1
2
ˆ
2
ˆ
2
yyyy sss −+=
Analýza rozptylu
AB
C
y
)( yy −
Interpretace výsledků analýzy rozptylu
Interpretace: p < 0,05 - existuje statisticky významný
rozdíl mezi rozptylem vysvětleným regresní přímkou (tedy
modelem trendu) a zbytkovým (reziduálním) rozptylem –
zvolený model trendu je vhodný
p-hodnota
Model vysvětluje více než 63 % proměnlivosti studované
charakteristiky v čase
Kritéria pro volbu vhodného
modelu trendové funkce I.
A. Volba vhodné trendové funkce by v prvé řadě měla
vycházet z věcné analýzy zkoumaného jevu.
Ta nám umožní zaměřit se na určité typy (skupiny)
funkcí či některé jiné předem vyloučit – jde o
funkci rostoucí či klesající, má inflexní bod či je
nekonečně rostoucí.
Pro použitou trendovou funkci je důležité, zda má
(logistický trend) či nemá (lineární trend – růst
řady není ničím omezen) asymptotu. Je to důležité
pro předpovídání chování časové řady.
Kritéria pro volbu vhodného modelu
trendové funkce II.
B. Analýza grafu časové řady a analýza reziduí.
-
yt
– empirické hodnoty
– teoretické hodnoty –
vyrovnané trendovou
funkcí
tyˆ
tyˆ
7
Kritéria pro volbu vhodného modelu
trendové funkce III (objektivní)
Spočívají v minimalizaci předem zvoleného kritéria
(jako v případě regresní analýzy). Za toto kritérium
se nejčastěji bere součet čtverců odchylek
empirických hodnot yt od hodnot vyrovnaných
(součet čtvercových chyb):
Z uvažovaných funkcí se vybírá ta s nejmenší
hodnotou reziduálního součtu čtverců.
POZOR – jde o formální kritérium. Např. použijeme-li
polynom vysokého stupně, může být reziduální součet
čtverců i nulový, avšak zcela nepoužitelný.
( )∑=
−=
n
t
tt yySSE
1
2
ˆ
tyˆ
Objektivní kritéria pro volbu vhodného
modelu trendové funkce III
Počítačové programy obvykle nabízejí následující
míry vhodnosti zvolené trendové funkce:
Střední chyba odhadu (M.E. – Mean Error)
( )
n
yy
EM
n
t
tt∑=
−
= 1
ˆ
..
( )
n
yy
ESM
n
t
tt∑=
−
= 1
2
ˆ
...
Střední čtvercová chyba odhadu (M.S.E. – Mean
Square Error)
Je to nejpoužívanější kritérium.
Informativní testy pro volbu vhodné
trendové křivky:
Trend Informativní test
lineární První diference (yt+1 - yt) jsou přibližně konstantní
kvadratický Druhé diference (yt+2 - 2yt+1 + yt) jsou přibližně
konstantní
exponenciální Podíly sousedních hodnot (yt+1/yt) resp. První
diference logaritmů tvaru (log yt+1 - log yt) jsou
přibližně konstantní
logistický Křivka prvních diferencí (yt+1-yt) se podobá
křivce normální hustoty, podíly (1/yt+2 -
1/yt+1)/(1/yt+1 - 1/yt) jsou přibližně konstantní
Gompertzova
křivka
Podíly (log yt+2 – log yt+1)/(log yt+1 – log yt) jsou
přibližně konstantní
Mechanické vyrovnávání trendu
metodami klouzavých průměrů
Používá se v případě, že se trend mění a nelze ho
vyrovnat „globálně“ jednou matematickou křivkou.
Metoda je vhodná pro neperiodické řady, neumožňuje
extrapolaci hodnot.
Vlastní průměry se používají jako prosté či vážené.
V některých případech lze použít klouzavých mediánů.
Klouzavé průměry mohou být necentrované a centrované
Metody klouzavých průměrů
Jako klouzavé průměry obecně označujeme lineární
kombinace členů původní řady.
)222(
8
1
2112 ++−− ++++ ttttt yyyyy
Patří mezi tzv. adaptivní přístupy k trendové složce
časové řady.
Tzv. polynomické klouzavé
průměry umožňují vyrovnání
hodnot na počátku a konci
časové řady
)(
5
1
2112 ++−− ++++ ttttt yyyyyprosté
vážené
Volba řádu klouzavých průměrů
• Subjektivní posouzení charakteru dat
• Délka klouzavých průměrů by měla odpovídat
periodě sezónních či cyklických fluktuací
• Vzorce pro výpočet optimální délky
Obsahuje-li řada sezónní složku, je vhodné volit řád klouzavých
průměrů tak, aby zahrnoval celou délku periody sezónní složky.
8
Centrované klouzavé průměry
Ve většině případů se používají klouzavé průměry liché
délky, u sudé délky je problém s přiřazením hodnot
časovému okamžiku.
V ekonomických časových řadách, které často obsahují
sezónní složku délky 4 (řady čtvrtletních hodnot) či 12
(řady měsíčních hodnot), se tento problém řeší tzv.
centrováním.
Výsledné klouzavé průměry pro sudou délku klouzavé
části vypočteme jako průměry dvou sousedních
klouzavých průměrů liché délky.
Centrované klouzavé průměry
Příklad: Abychom vystihli roční chod určitého
ukazatele, chceme pro řadu měsíčních hodnot použít
klouzavých průměrů délky 12.
Shlazená hodnota však spadá doprostřed mezi
„červen“ a „červenec“. Další shlazená hodnota pak
mezi „červenec“ a „srpen“.
Tyto dva jednoduché klouzavé průměry vezmeme a
zprůměrňujeme. Výsledek pak už můžeme přiřadit k
„červencové“ hodnotě. Tedy vytváříme klouzavé
průměry o délce 13:
)2...22(
24
1
))...(
12
1
)...(
12
1
(
2
1
ˆ
65456
645556
++−−−
+−−+−−
+++++
=+++++++=
ttttt
ttttttt
yyyyy
yyyyyyy
Centrované klouzavé průměry
Obecně místo jednoduchých klouzavých průměrů délky
2m vytváříme centrované klouzavé průměry délky 2m+1
podle tohoto obecného vzorce:
)2...2(
4
1
ˆ 11 mtmtmtmtt yyyy
m
y +−++−− ++++=
m - polovina řádu klouzavých průměrů (shlazovacího okna)
Vážené klouzavé průměry
• Jednotlivé členy úseku řady přiřazeny váhy.
• Tyto váhy většinou lineárně klesají směrem od
středního (vyrovnávaného) členu.
• Váhy mohou mít také např. podobu tzv. gaussova filtru.
Člen
řady
váha
yt-4 0,014
yt-3 0,048
yt-2 0,117
yt-1 0,201
yt 0,241
yt+1 0,201
yt+2 0,117
yt+3 0,048
yt+4 0,014
Gaussův filtr pro m=4
Statistiky-Pokročilé lineární/ nelineární modely-Časové řady/Predikce
Analýza časových řad v programu Statistica
Výpočet klouzavých průměrů
Pozor na vykreslování klouzavých
průměrů v programu EXCEL !!!
… a EXCEL
1. klasický přístup k sezónní dekompozici
2. úvod do autokorelační analýzy
Sezónní složka St je typická pro časové řady, jejichž interval pozorování
je kratší než jeden rok (sezóna může mít délku týden, měsíc, roční
období).
Objevuje se v řadách ekonomických (tržby, produkce, …), ale i v řadách
meteorologických prvků (roční chod teploty vzduchu).
Řada obsahující sezónní složku se vyznačuje pravidelným opakováním
hodnot kolem trendu a toto opakování může mít délku např. 7 dnů (do
týdne), 12 měsíců či 4 roční období (do roku).
Sezónní složka může mít aditivní resp. multiplikativní charakter
Analýza sezónní složky časových řad
(sezónní očišťování)
9
Obecný model řady při sezónním
očišťování
Trendovou a cyklickou složku považujeme za jeden
celek:
aditivní model: tttt STCY ε++=
Yt je pozorovaná hodnota časové řady v čase t.
Statistiky-Pokročilé lineární/
nelineární modely-Časové
řady/Predikce
Jednotlivé kroky analýzy sezónní složky
1. Z originální řady obsahující
sezónní složku je vypočtena
řada klouzavých průměrů
s délkou klouzavých průměrů
rovnou délce sezónní složky.
2. Vytvoříme novou řadu jako
rozdíl řady původní a řady
shlazené.
Jednotlivé kroky analýzy sezónní složky
3. Tzv. sezónní komponenty
jsou vypočteny jako průměr pro
každý člen v rámci sezóny.
Výsledné hodnoty představují
průměrnou sezónní složku
v časové řadě.
4. Sezónně očištěná řada (tedy
řada obsahující vedle náhodné
složky ještě složku TCt) se potom
vyjádří jako rozdíl řady
originální a sezónní
komponenty.
Jednotlivé kroky analýzy sezónní složky
5. Složka TCt se většinou
aproximuje řadou shlazenou
váženým klouzavým
průměrem délky 5 se
symetrickými vahami (1, 2,
3, 2, 1).
6. Obdobně lze izolovat
náhodnou složku jako rozdíl
řady sezónně očištěné a řady
se zvýrazněnou složkou TCt (
viz. bod 5).
Autokorelace časových řad
Autokorelační analýza - metoda, kterou lze zkoumat
vzájemné vztahy mezi hodnotami jedné časové řady.
Může sloužit jako metoda k definování sezónní a
cyklické složky časových řad.
Jejím základem je výpočet autokorelačního
koeficientu, resp. autokorelační funkce.
Autokorelační koeficient
Autokorelační koeficient rk je relativní míra
proměnlivosti členů časové řady posunutých o určitou
hodnotu k. Definuje vztah mezi členy časové řady yt a
yt+k.
Posun k se z angličtiny označuje jako lag. Je to tedy
korelační koeficient vypočtený mezi jednotlivými členy
časové řady, mezi kterými je k-1 jiných pozorování tedy
lag = k a označujeme ho jako autokorelační koeficient ktého
řádu.
Pro k = 0 je hodnota r0 = 1 - je to vlastně hodnota
korelačního koeficientu.
10
Základní pojmy
Rozptyl (variance) – míra variability
(proměnlivosti) statistického znaku x 1
)(
1
2
2
−
−
=
∑=
n
xx
s
n
i
i
x
Kovariance – absolutní míra
vzájemné variability dvou
statistických znaků x; y 1
))((
1
−
−−
=
∑=
n
yyxx
s
n
i
ii
xy
Korelace - relativní míra vzájemné
variability dvou statistických znaků
x; y yx
xy
xy
ss
s
r
⋅
=
Základní vztahy
Autokorelační funkce –
hodnoty ry(k) pro k=1,2,…M,
kde M < N/2, N – délka řady
Autokorelace – relativní míra
proměnlivosti členů časové řady y
posunutých o určitou hodnotu k.
Autokovariance – absolutní míra
proměnlivosti členů časové řady y
posunutých o určitou hodnotu k. 1
))((
1
−−
−−
=
∑
−
=
+
kn
yyyy
c
kn
i
kii
k
2
0
)(
y
kk
y
s
c
c
c
kr ==
Autokorelační funkce
Autokorelační funkce (ACF) je potom závislost mezi
hodnotami autokorelačního koeficientu a hodnotami
posunu k.
Vyjadřuje se formou grafu – tzv. korelogramu (viz.
obrázek). Na ose x jsou hodnoty lag (k), na ose y
hodnoty autokorelačního koeficientu.
Hodnoty autokorelační funkce se pohybují v intervalu –
1,1.
ACF je vhodným nástrojem k posouzení, zda časová
řada obsahuje cyklickou či periodickou složku a také
zda je či není řadou náhodných čísel – tedy do jaké míry
je možné ji extrapolovat (předpovídat).
Princip autokorelace časové řady
Interpretace ACF I
Korelogram bývá doplňován intervaly spolehlivosti,
kterými lze hodnotit statistickou významnost
autokorelačních koeficientů.
95 % interval
spolehlivosti ACF lze z
dostatečnou přesností
zkonstruovat ze
vztahu:
N
2±
N – délka časové řady
Časová řada náhodných čísel (bílý šum) a její
autokorelační funkce
11
Časová řada bez periodické složky se silnou
autokorelací a její autokorelační funkce
Časová řada obsahující výraznou sezónní
složku a její autokorelační funkce