MODELY s ROZLOŽENÝMI ZPOŽDĚNÍMI V následující části popíšeme základní charakteristiky ekonomickomatematických modelů , jejichž základní vlastností je, že operují se zpožděnými hodnotami ekonomických veličin. Takové modely se vyskytují velmi frekventovaně v mnoha oblastech formalizovatelné ekonomické teorie, zpožděné veličiny se objevují zejména u spotřebních a investičních funkcí, dále u vztahů popisujících např. vývoj nezaměstnanosti, stavu zásob apod. Je tedy nesporně užitečné vyložit alespoň základní poznatky o tom, čím se přítomnost zpožděných proměnných vyznačuje a jaké specifické rysy s sebou přináší. Omezíme se přitom jen na nejcharakterističtější specifikace vztahů obsahujících zpožděné proměnné, aniž si budeme všímat dopadů výskytu zpožděných proměnných na statistické vlastnosti (zejména při odhadu modelových parametrů). Ty jsou ostatně předmětem často velmi podrobných analýz prováděných ekonometrickými prostředky, kde, jak známo, přítomnost zejména zpožděných endogenních proměnných vyžaduje aplikaci specifických postupů (u vícerovnicových modelů např. nasazení metod simultánního odhadu modelových parametrů, jako 2SLS, 3SLS, LIML apod. ) V následujícím výkladu budeme předpokládat obecný lineárně aditivní vztah mezi závisle proměnnou skalární veličinou Y a nezávisle proměnnou (rovněž skalární veličinou) X, tj. vztah (1) Y = a + bX. K rozlišení toho, zda z hlediska vývoje ekonomických proměnných budeme považovat čas přiřazený průběhu obou veličin (shodně) za diskrétní či spojitý, budeme v prvém případě psát : (1a) Y[t] = a + bX[t ] [ ], ve druhém pak (1b) Y( t ) = a + bX ( t ) [ ]. Pro oba tyto případy vyložíme níže základní poznatky vztahující se k přítomnosti zpožděných ekonomických veličin v lineárním kauzálním modelu. V dalším se ukáže se účelným zavést pojem potenciální hodnota závisle proměnné. K tomu účelu zavedeme vztah (2) Z[t] = a + bX[t ]jako veličinu, která popisuje lineární (nezpožděný) vztah mezi závisle proměnnou Y a nezávisle proměnnou X. Na vztah tento lze pohlížet jako na určitou aproximaci vysvětlované veličiny Z pomocí lineární transformace založené na vysvětlující veličině X.. 1. Situace s diskrétním časem 1A - Případ odkladu o pevnou dobu r Nejprve se zabývejme situací, kdy vliv působení nezávisle proměnné X na závisle proměnnou Y je jednorázový se zpožděním rovným pevné hodnotě r. Tento případ zapíšeme jako (3) Y[t] = a + bX [t-r , ] kde r je pevná délka zpoždění (odkladu), s nímž veličina X ovlivňuje závisle proměnnou Y. Nejčastějším případem (v ekonomických aplikacích vycházejících rovněž z obvyklých hodnot zpoždění) je zpoždění o 1 období odpovídající vztahu Y[t] = a + b.X [t-1], tedy pro r = 1. Ve vztahu k potenciální hodnotě závisle proměnné získáváme pro případ odkladu o pevnou dobu okamžitě vztah (3a) Y[t] = Z [t -r ] = a + bX[t -r ] , 1B - Případ obecného rozloženého zpoždění Případ uvedený v 1A rozšíříme o možnost, že hodnota veličina Y v čase t nemusí záviset toliko na jediné zpožděné hodnotě veličiny X, nýbrž obecně na celé posloupnosti (třeba nekonečné) takových hodnot. Znamená to tedy připustit platnost vztahu (4) Y[t] = a + b[1]X [t -1] +b[2]X [t -2] + b[3] X [t-3] +.......+ b[r]X [t-r] . . . . . . atd. s tím, že pro případ nekonečné posloupnosti musíme doplnit předpoklad (aby závisle proměnná byla konečná veličina) , že součet aa b[1] +b[2 ] +b[3 ] + ...... [ ]je konečný. Pokud bychom tímto vztahem popsali např. závislost mezi spotřebou C a důchodem Y, tzn. (5) C[t ] = a + b[1] Y[t -1] + b[2] Y [t -2] + b[3] Y [t -3 ] + . . . . . ( zde důchod Y vystupuje v roli nezávisle proměnné ) mohli bychom přirozeně příslušným parametrům přisoudit ekonomicky interpretovatelný význam. Tak bychom parametry b[i] mohli nazvat parciálními mezními sklony ke spotřebě, zatímco součtová hodnota b by vyjadřovala celkový sklon ke spotřebě. V obecném případě se součinitele b[k ] nazývají váhovými koeficienty (bez ohledu na to, zda je jich konečný nebo nekonečný počet). Dosti často (ne však nutně) se předpokládá, že součet koeficientů u všech zpožděných nezávisle proměnných je jedničkový neboli, že tyto koeficienty jsou normované, tzn. aa b[1] + b[2 ] + b[3 ] + ...... = 1, jindy pouze, že součet je konečný, tzn. aa b[1] + b[2 ] + b[3 ] + ...... = b. Normování na lze ostatně snadno provést. Váhové koeficienty však nemusí být nutně všechny kladné. Učiňme nyní několik poznámek : a) je zřejmé, že případ zpoždění o pevnou dobu je speciálním případem modelu rozloženého zpoždění, pro který platí b[k ] = 1 a současně b [j] = 0 pro všechna j ^1 k. b) pokud by všechny veličiny X[k-r ] = c (reálná konstanta), potom též Y = c. Pro obecné rozložené zpoždění pak dostaneme vztah (6) Y[t] = l[1] Z [t -1] + l[2] Z [t -2] + l [3] Z [t - 3 ] + . . . . . . . . , kde součet koeficientů l [i] dává dohromady 1, tj. l[1] + l[2] + l [3] + . . . . = 1. Ověření: Abychom to ukázali, dosaďme z transformace ( 2 ) Z[t] = a + bX[t ] zpětně za X[t ] = (Z[t] - a)/b . Po dosazení za X[1] , X[2] , X [3] v (4) atd. obdržíme : (7) Y[t] = a + b[1](Z[t-1] - a)/b +b[2](Z[t-2] - a)/b + b[3](Z[t-3] - a)/b +. . . . . atd. neboli po přeskupení členů (7*) Y[t] = b[1](Z[t-1])/b +b[2](Z[t-2])/b + b[3](Z[t-3])/b +....+ a - b[1.]a/b - b[2].a/b- b[3].a/b nebo též jinak zapsáno (8) Y[t] = l[1]Z[t-1] +l[2]Z[t-2] + l[3] Z[t-3] +. . . . . + a.( 1- b[1]/b - b[2]/b- b[3]/b + ....... ) ., kde - poněvadž b[1] + b[2] + b[3] + .... = b, je obsah závorky na pravé straně roven nule. OS . 1C - Případ geometricky rozloženého zpoždění Zvláštním ( a patrně i nejdůležitějším) typem rozloženého zpoždění je případ, kdy koeficienty b[1], b[2] , b[3 ] představují jednotlivé, po sobě jdoucí členy geometrické posloupnosti (tj. podíl b[k+1] / b[k ] je roven konstantní hodnotě-kvocientu) Potom lze závisle proměnnou zapsat jako (9) Y[t ] = a + d (1- c) ( X [t-1] +c[.]X [t-2] + c^2 X [t-3 ] + . . . . ) , kde váhové koeficienty příslušné jednotlivým zpožděným proměnným s rostoucí hodnotou zpoždění klesají. Kvocient je představován hodnotou c s podmínkou c<1 nutnou pro konvergenci geometrické řady. Máme tedy zaručeno, že platí aa b[1]+b[2]+b[3] +... = d, v případě d =1 pak jedničkový. Konstanta a (pouze) ovlivňuje střední hodnotu závisle proměnné Y[t]. Součet koeficientů u zpožděných proměnných je totiž 1 d (1-c) ( 1 + c + c^2 + c^3 + .,.. ) = d (1-c ) ----- = d. 1-c Konečně geometricky rozložené zpoždění můžeme charakterizovat podobnou závislostí mezi Y a Z : (10) Y[t] = (1-c) ( Z[t -1] +c Z [t -2] + c^2 Z [t - 3 ] + . . . . ) [ ]atd., což opět můžeme ověřit takto : Vyjdeme ze zápisu (9), v němž dosadíme za X [t] = ( Z [t] - a ) / d a obdobně za zpožděné hodnoty X [t-1] = ( Z [t-1]- a ) / d , X [t-2] = ( Z [t-2] - a ) / d atd. Dostaneme Y[t ] = a + d.(1- c) [ ( Z [t-1] - a ) / d + c ( Z [t-2] - a ) / d + c^2 ( Z [t-3] - a ) / d + .... ] , neboli po vykrácení všech výrazu v hranaté závorce konstantou d (11) Y[t ] = a + (1- c) (Z [t-1] - a + c[.]Z [t-2] - c. a + c. Z [t-3 ] - c^2. a + .... ) , až konečně po přeskupení členů Y[t ] = (1- c).( Z [t-1] + c[.]Z [t-2] + c^2[.]Z [t-3] + ... ) + a [ 1 - (1-c) - c.(1-c) - c^2.(1-c) - .... ] , Součet prvků v hranaté závorce je však nulový, neboť součet všech záporných členů uvnitř závorky je roven 1. Odtud tedy obdržíme hledaný vztah (10) Y[t ] = (1- c) ( Z [t-1] + c[.]Z [t-2] + c^2[.]Z [t-3] + ... ) Tento poslední případ můžeme dále upravit do vztahu popisujícího přírůstkovou relaci pro Y[t ] : Zaveďme nejprve změnu veličiny Y jako první diferenci DY [t] = Y[t +1] - Y[t ] a uplatněme výraz (10) DY [t ](11) ------ = ( Z[t] +c Z [t -1] + c^2 Z [t - 2 ] + . . . ) - ( Z[t -1] +c Z [t -2] + c^2 Z[t - 3 ] + ........ ) = 1-c = Z[t ] - (1-c) ( Z[t -1] +c Z [t -2] + c^2 Z [t - 3 ] + ........ ) = Z[t ] - Y[t ] . Položme nyní l = 1-c , kde 0 < l < 1 . Potom zřejmě platí (12) DY [t ] = - l (Y [t ] - Z [t] ) , kde Z [t] = a + b X [t ] Jestliže si dále uvědomíme, že jsme tento vztah odvodili pro pevnou délku období D [t] = 1, nebrání nám nic v tom, abychom uvažovali libovolnou hodnotu D [t ] místo jedničky. Vyplývá to z toho, že bychom mohli definovat jako DY[t ] = Y[t+Dt ] - Y[t ] a následně uvažovat [ ]průměrnou změnu Y[t] mezi časovými okamžiky t a t + Dt . Potom platí DY [t ](13) ------ = - l (Y [t ] - Z [t] ) D [t ] s následující věcnou interpretací : Marginální rychlost růstu (nebo průměrná za období D [t] ) je úměrná záporně vzatému rozdílu skutečné a potenciální hodnoty závisle proměnné. Koeficient této úměrnosti l se nazývá rychlost reakce (a je přirozeně ovlivněn hodnotou parametru c) Uvedený vztah tedy popisuje adaptaci /přizpůsobování závisle proměnné k její potenciální hodnotě, přičemž (při pevné rychlosti reakce) je přizpůsobení úměrné rozdílu mezi skutečnou a potenciální hodnotou této veličiny. 2. Situace se spojitým časem Zde předpokládáme, že veličina X, která je funkcí času t, ovlivňuje jinou veličinu Y (jež je také funkcí času) prostřednictvím vztahu (14) Y(t) = f (X(t) ) , kde f je funkce vhodných vlastností (např. spojitá). V dalším budeme předpokládat, že funkce f je lineární , takže vztah mezi oběma veličinami X,Y lze zapsat jako Y(t) = a + b.X (t), kde a, b jsou koeficienty lineární funkce (a úrovňová konstanta, b parametr sklonu). Pokud bude v tomto vztahu existovat zpoždění, označíme "aproximovanou" hodnotu veličiny Y pomocí veličiny X jako Z(t) = a + b.X(t) , kde Z opět hraje roli potenciální hodnoty veličiny Y. Skutečná hodnota Y se tedy bude opožďovat za svou potenciální hodnotou Z , hodnota zpoždění je r. 2A - Případ odkladu o pevnou dobu r Obdobně výše uvedenému diskrétnímu případu zapíšeme pro lineární vztah veličin X,Y a pro hodnotu zpoždění r vztah (15) Y(t) = Z ( t- r ) = a + b.X( t - r ) [ ]pro pevnou délku odkladu r, s nímž veličina X ovlivňuje závisle proměnnou Y. Potenciální hodnotu závisle proměnné Z(t) zavedeme analogicky jako v diskrétním případě (16) Z(t) = a + b.X(t) 2B - Případ obecného spojitého rozloženého zpoždění Místo nekonečné posloupnosti koeficientů b[1], b[2], ..., b[3], se kterou jsme pracovali v nespojitém případě, budeme uvažovat spojitou množinu hodnot nějaké funkce f(t) spojitě se měnícího argumentu t. Součet (konečný nebo nekonečný) se podobně nahradí integrálem . Dostaneme tedy YEN (17) Y(t) = a + b o f(t) X (t-[,]t) dt , přičemž funkce f musí splňovat 0 podmínku (obdobnou jako u hustoty spojitého pravděpodobnostního rozdělení definovaného na celém oboru nezáporných hodnot) YEN o f(t) dt = 1. 0 Omezení se na nezáporný definiční obor vyplývá z toho, že zpoždění uvažujeme jako kladné (v případě záporné hodnoty by šlo o předstih). Konkrétní tvar funkce f(t) udává (spojitě vyjádřené) časové rozložení zpoždění, jinými slovy závislost hodnoty proměnné Y v libovolném čase t na hodnotě proměnné X v čase (t-[,]t ) , tj. o dobu t dříve. Přejdeme-li od proměnné Y k její potenciální hodnotě Z, nabude výraz ( 17 ) tvaru YEN (18) Y(t) = a + b o f(t) Z(t-[,]t) dt se stejnými vlastnostmi funkce f jako dříve. 0 Poznámka : Odklad o pevnou dobu r by byl jistým speciálním případem vztahů ( 17 ), resp. ( 18 ), pokud bychom za f( t ) vzali "degenerovanou" váhovou funkci tvaru f( t ) = 1 pro t = r a f( t ) = 0 v jiných případech. ( Jde o obdobu situace, kdy "náhodná veličina" s nulovým rozptylem přechází s jednotkovou pravděpodobností v nenáhodnou-konstantní veličinu a kdy je její hustota degenerována do jediného bodu s jedničkovou hodnotou). Další, matematický problém ovšem vyvstává v tom, že příslušná "distribuční" funkce F(t) jako primitivní funkce k f(t) není v bodě t = r spojitá a že tedy by bylo místo f(t) dt nutno psát d F(t) a užít (místo Newtonova) Stieltjesova integrálu. 2C - Případ exponenciálně rozloženého zpoždění Jde o významný případ spojitého zpoždění analogický geometricky rozděleného zpoždění u nespojité situace. V tomto případě je váhová funkce představována hustotou exponenciálního rozdělení, která má tvar (19) f(t) = l.exp (-lt) , kde l je kladná konstanta. Lze snadno ověřit, že jde skutečně o hustotu , neboť integrál přes kladný obor hodnot je roven 1. Ověření: Zavedeme-li do vzorce (18 ) , v němž je Z = a + bX , hustotu exponenciálně rozděleného zpoždění, dostaneme YEN (20) Y(t) = l o exp (-lt) Z ( t-[,]t) dt 0 Při použití jednoduché substituce integrační proměnné x = t -[,]t , obdržíme t t (21) Y(t) = l o exp (-l(t-x)) Z (x) dx = l exp (-lt) o exp (lx).Z(x) dx - YEN - YEN Platí tedy t (22) Y(t) exp(lt) = l o exp (lx) Z(x) dx - YEN a derivujeme-li obě strany podle času t (mj. horní meze integrálu), dostaneme dY(t) t l.Y(t).exp(lt) + ------ exp(lt) = l . { o exp (-lx) Z(x) dx } = l .exp (lt). Z(t) dt - YEN Zkrátíme-li nyní obě strany výrazem exp(lt), dostaneme dY(t) (23) -------- = - l.(Y(t) - Z(t)) dt což je tvar akcelerátoru při exponenciálně rozděleném zpoždění. OS . Výsledný výraz je analogický vztahu ( 12 ) pro případ exponenciálně rozděleného zpoždění. Rychlost změny veličiny Z je vždy úměrná hodnotě -l.[ Y(t) - Z(t) ] , tj. záporně vzatému rozdílu mezi skutečnou a potenciální hodnotou proměnné Y. Koeficient úměrnosti l je rychlost reakce (akcelerátoru), jeho převrácená hodnota T = 1/l^ se nazývá časová konstanta zpoždění. Pojmenování lze vztáhnout jak k exponenciálně, tak k geometricky rozloženým zpožděním.