Outoregressiewe Bewegende Gemiddelde Reeks

Outoregressiewe bewegende gemiddelde ARMA (p, q) Modelle vir Tydreeksanalise - Deel 3 Deur Michael Saal-Moore op 7 September 2015 Dit is die derde en laaste pos in die mini-reeks oor outoregressiewe bewegende gemiddelde (ARMA) modelle vir tydreeks ontleding. Weve bekendgestel outoregressiemodelle en bewegende gemiddelde modelle in die twee vorige artikels. Nou is dit tyd om hulle te kombineer om 'n meer gesofistikeerde model te produseer. Uiteindelik sal hierdie ons lei tot die ARIMA en GARCH modelle wat ons sal toelaat om bate opgawes en voorspelling wisselvalligheid voorspel. Hierdie modelle sal die basis vir handel seine en risikobestuur tegnieke vorm. As jy het gelees Deel 1 en Deel 2 sal jy gesien het dat ons geneig is om 'n patroon vir ons ontleding van 'n tydreeks model volg. Siek herhaal dit kortliks hier: Rasionaal - Hoekom is ons belangstel in hierdie spesifieke model Definisie - 'n wiskundige definisie vir dubbelsinnigheid te verminder. Correlogram - plot van 'n monster correlogram 'n modelle gedrag te visualiseer. Simulasie en Fitting - Pas die model om simulasies, ten einde weve verseker verstaan ​​die model korrek. Real finansiële inligting - Pas die model om werklike historiese batepryse. Voorspelling - Voorspelling daaropvolgende waardes te handel seine of filters te bou. Ten einde hierdie artikel volg, is dit raadsaam om 'n blik op die vorige artikels oor tydreeksanalise neem. Hulle kan al hier gevind word. Bayes inligting maatstaf in Deel 1 van hierdie artikel reeks het ons gekyk na die Akaike Inligting Criterion (AIC) as 'n manier om ons te help kies tussen afsonderlike beste tyd reeks modelle. A nou verwant instrument is die Bayes inligting Kriterium (BIC). In wese is dit het 'n soortgelyke gedrag by die AIC deurdat dit penaliseer modelle vir die feit dat te veel parameters. Dit kan lei tot overfitting. Die verskil tussen die BIC en AIC is dat die BIC is strenger met sy penalisering van addisionele parameters. Bayes inligting Criterion As ons die waarskynlikheid funksie vir 'n statistiese model wat k parameters het, en L maksimeer die waarskynlikheid. dan die Bayes inligting Criterion word gegee deur: waar n die aantal datapunte in die tyd reeks. Ons sal met behulp van die AIC en BIC hieronder by die keuse van geskikte ARMA (p, q) modelle. Ljung-Box toets in Deel 1 van hierdie artikel reeks Rajan genoem in die Disqus kommentaar dat die Ljung-Box toets was meer gepas as die gebruik van die Akaike Inligting Criterion van die Bayes inligting Kriterium om te besluit of 'n ARMA model was 'n goeie passing vir 'n tyd reeks. Die Ljung-Box toets is 'n klassieke hipotese toets wat ontwerp is om te toets of 'n stel van outokorrelasies van 'n toegeruste tydreeksmodel aansienlik verskil van nul. Die toets nie elke individu lag vir willekeur te toets nie, maar eerder toets die willekeur oor 'n groep van lags. Ljung-Box Toets Ons definieer die nulhipotese soos: Die tydreeksdata by elke lag is i. i.d .. dit is die korrelasies tussen die bevolking reeks waardes is nul. Ons definieer die alternatiewe hipotese as: Die tydreeksdata is nie i. i.d. en besit serial korrelasie. Ons bereken die volgende toetsstatistiek. V: Waar N is die lengte van die tyd reeks monster, hoed k is die monster outokorrelasie op lag k en h die aantal lags onder die toets. Die besluit reël om te bepaal of die nulhipotese verwerp is om vas te stel of Q GT Chi2, vir 'n chi-kwadraat verspreiding met h grade van vryheid aan die 100 (1-alfa) ste persentiel. Terwyl die besonderhede van die toets effens kompleks mag lyk, kan ons in werklikheid gebruik R tot die toets vir ons te bereken, vereenvoudig die prosedure ietwat. Autogressive bewegende gemiddelde (ARMA) Models van orde p, q Noudat weve die BIC en die Ljung-Box toets bespreek, was gereed om ons eerste gemengde model, naamlik die outoregressiewe bewegende gemiddelde van orde p, q, of ARMA (bl bespreek, Q). Rasionaal Tot op datum het ons outoregressiewe prosesse beskou en bewegende gemiddelde prosesse. Die voormalige model beskou sy eie verlede gedrag as insette vir die model en as sodanig pogings om die mark deelnemer effekte, soos momentum en gemiddelde-terugkeer in-beurs vang. Laasgenoemde model word gebruik om skok inligting kenmerk van 'n reeks, soos 'n verrassing verdienste aankondiging of onverwagte gebeurtenis (soos die BP Horizon Deep oliestorting). Dus, 'n ARMA model poog om beide hierdie aspekte te vang wanneer modellering finansiële tydreekse. Let daarop dat 'n ARMA model nie in ag neem wisselvalligheid groepering, 'n belangrike empiriese verskynsels van baie finansiële tydreekse. Dit is nie 'n voorwaardelik heteroscedastic model. Vir wat sal ons moet wag vir die boog en GARCH modelle. Definisie Die ARMA (p, q) model is 'n lineêre kombinasie van twee lineêre modelle en dus is self nog lineêre: outoregressiewe bewegende gemiddelde Model van orde p, q 'n tydreeksmodel, is 'n outoregressiewe bewegende gemiddelde model van orde p, q . ARMA (p, q), indien: begin xt alfa1 x alfa2 x ldots wt beta1 w beta2 w ldots betaq w end Waar is wit geraas met E (WT) 0 en variansie sigma2. As ons kyk na die agterste Shift-operateur. (Sien 'n vorige artikel) dan kan ons herskryf bogenoemde as 'n funksie theta en phi van: Ons kan reguit sien dat ons deur die oprigting van p neq 0 en Q0 herstel die AR (p) model. Net so as ons 'p 0 en Q neq 0 herstel ons die MA (Q) model. Een van die belangrikste kenmerke van die ARMA model is dat dit karig en oorbodig in sy parameters. Dit wil sê, 'n ARMA model sal dikwels minder parameters as 'n AR (p) of MA (Q) model alleen vereis. Daarbenewens, as ons herskryf die vergelyking in terme van die BSO, dan is die theta en phi polinome kan soms 'n gemeenskaplike faktor wat sal lei tot 'n eenvoudiger model. Simulasies en Correlograms Soos met die outoregressiewe en bewegende gemiddelde modelle sal ons nou simuleer verskeie ARMA reeks en dan probeer om ARMA modelle te pas by hierdie realisasies. Ons dra dit uit, want ons wil om te verseker dat ons verstaan ​​die gepaste prosedure, insluitend hoe om vertrouensintervalle bereken vir die modelle, asook te verseker dat die proses eintlik redelik skattings vir die oorspronklike ARMA parameters nie herstel. In Deel 1 en Deel 2 hand ons die AR en MA-reeks gebou deur 'N monsters van 'n normale verspreiding en dan knutselen die spesifieke tydreeksmodel behulp lags van hierdie monsters. Daar is egter 'n meer eenvoudige manier om AR, MA, ARMA en selfs ARIMA data, na te boots bloot deur die gebruik van die arima. sim metode in R. Kom ons begin met die eenvoudigste moontlike nie-triviale ARMA model, naamlik die ARMA (1,1 ) model. Dit wil sê, 'n outoregressiewe model van orde een gekombineer met 'n bewegende gemiddelde model van orde een. So 'n model het slegs twee koëffisiënte, Alpha en Beta, wat die eerste lags van die tydreeks self en die skok wit geraas terme verteenwoordig. So 'n model word gegee deur: Ons moet die koëffisiënte voor spesifiseer om simulasie. Kom ons neem 'n alfa 0,5 en beta -0,5: Die produksie is soos volg: Kom ook plot die correlogram: Ons kan sien dat daar geen beduidende outokorrelasie, wat verwag kan word van 'n ARMA (1,1) model. Ten slotte, kan probeer bepaal die koëffisiënte en hul standaard foute met behulp van die ARIMA funksie: Ons kan die vertrouensintervalle bereken vir elke parameter gebruik van die standaard foute: Die vertrouensintervalle doen bevat die ware parameter waardes vir beide gevalle egter moet ons daarop let dat die 95 vertrouensintervalle is baie breed ( 'n gevolg van die redelike groot standaard foute). Kom nou probeer om 'n ARMA (2,2) model. Dit wil sê, 'n AR (2) model gekombineer met 'n MA (2) model. Ons moet vier parameters vir hierdie model spesifiseer: alfa1, alfa2, beta1 en beta2. Kom ons neem alfa1 0.5, alpha2-0.25 beta10.5 en beta2-0.3: Die uitset van ons ARMA (2,2) model is soos volg: En die ooreenstemmende autocorelation: Ons kan nou probeer pas 'n ARMA (2,2) model te die data: Ons kan ook bereken die vertrouensintervalle vir elke parameter: Let daarop dat die vertrouensintervalle vir die koëffisiënte vir die bewegende gemiddelde komponent (beta1 en beta2) die oorspronklike parameter waarde nie eintlik bevat. Dit gee 'n uiteensetting van die gevaar van 'n poging om modelle te pas om data, selfs wanneer ons weet die ware parameterwaardes egter vir doeleindes van handeldryf ons net nodig het om 'n voorspellende krag wat kans oorskry en produseer genoeg wins bo transaksiekoste het, ten einde winsgewend in te wees die lang termyn. Nou dat weve 'n paar voorbeelde van gesimuleerde ARMA modelle gesien moet ons meganisme vir die keuse van die waardes van p en q toe pas om die modelle te real finansiële data. Die keuse van die beste ARMA (p, q) Model, Q van die ARMA model is geskik vir 'n reeks Ten einde vas te stel watter volgorde p, moet ons die AIC (of BIC) te gebruik in 'n subset van waardes vir p, q, en dan van toepassing die Ljung-Box toets om te bepaal of 'n goeie passing is bereik, vir bepaalde waardes van p, q. Om hierdie metode gaan ons eerstens na te boots 'n bepaalde ARMA (p, q) proses wys. Ons sal dan lus oor die hele paarsgewyse waardes van p in en Q in en bereken die AIC. Ons sal die model met die laagste AIC kies en dan hardloop 'n Ljung-Box toets op die residue om te bepaal of ons 'n goeie passing bereik het. Kom ons begin deur simuleer 'n ARMA (3,2) reeks: Ons sal nou 'n voorwerp finale om die beste model pas en laagste AIC waarde te stoor. Ons loop oor die verskillende p, q kombinasies en gebruik die huidige voorwerp om die pas van 'n ARMA (i, j) model stoor, vir die herhaling veranderlikes i en j. As die huidige AIC minder as 'n voorheen bereken AIC is het ons die finale AIC om hierdie huidige waarde en kies daardie volgorde. By beëindiging van die lus het ons die einde van die ARMA model gestoor in final. order en die ARIMA (p, d, q) pas self (met die Geïntegreerde d komponent ingestel op 0) gestoor as final. arma: Kom uitset die AIC , orde en ARIMA koëffisiënte: Ons kan sien dat die oorspronklike bevel van die gesimuleerde ARMA model verhaal, naamlik met P3 en q2. Ons kan die corelogram van die residue van die model plot om te sien of hulle lyk soos 'n verwesenliking van diskrete wit geraas (DWN): Die corelogram inderdaad lyk soos 'n verwesenliking van DWN. Ten slotte, ons voer die Ljung-Box toets vir 20 lags om dit te bevestig: Let daarop dat die p-waarde groter as 0.05, wat bepaal dat die residue is onafhanklik op die vlak 95 en dus 'n ARMA (3,2) model bied 'n goeie model pas. Dit is duidelik dat indien dit die geval wees, aangesien weve gesimuleerde die data onsself Dit is egter juis die proses sal ons gebruik wanneer ons kom ARMA (p, q) modelle om die SampP500 indeks in die volgende artikel te pas. Finansiële data nou dat weve beskryf die prosedure vir die keuse van die optimale tyd reeks model vir 'n gesimuleerde reeks, is dit eerder eenvoudig om dit toe te pas om finansiële data. Vir hierdie voorbeeld gaan ons weer kies die SampP500 VSA Equity Index. Kom ons laai die daaglikse sluitingspryse behulp quantmod en dan skep die log opbrengste stroom: Kom uit te voer dieselfde pas prosedure as vir die gesimuleerde ARMA (3,2) reeks bo op die puntelys opbrengste reeks van die SampP500 met behulp van die AIC: Die beste pas model het einde ARMA (3,3): Kom ons plot die residue van die toegeruste model om die SampP500 teken daaglikse opgawes stroom: Let daarop dat daar 'n paar beduidende hoogtepunte, veral by hoër lags. Dit is 'n aanduiding van 'n swak passing. Kom ons doen 'n Ljung-Box toets om te sien as ons statistiese bewyse vir hierdie: Terwyl ons vermoed, die p-waarde is minder as 0,05 en as sodanig kan ons nie sê dat die residue is 'n verwesenliking van diskrete wit geraas. Daar is dus bykomende outokorrelasie in die residue wat nie verklaar word deur die ingeboude ARMA (3,3) model. Volgende stappe soos weve al langs in hierdie artikel reeks het ons bewyse van voorwaardelike heteroskedastisiteit (wisselvalligheid groepering) gesien word in die SampP500 reeks bespreek, veral in die tydperke rondom 2007-2008. Wanneer ons later gebruik 'n GARCH model in die artikel reeks sal ons sien hoe hierdie outokorrelasies skakel. In die praktyk, ARMA modelle is nooit oor die algemeen goed pas vir log aandele opbrengste. Ons moet rekening hou met die voorwaardelike heteroskedastisiteit en gebruik 'n kombinasie van ARIMA en GARCH. Die volgende artikel sal oorweeg ARIMA en wys hoe die Geïntegreerde komponent verskil van die ARMA model het ons oorweeg in hierdie artikel. Michael Saal-Moore Mike is die stigter van QuantStart en is betrokke by die kwantitatiewe finansiële sektor vir die afgelope vyf jaar, in die eerste plek as 'n quant ontwikkelaar en later as 'n quant handelaar konsultasie vir verskansingsfondse. Verwante ArticlesAutoregressive bewegende gemiddelde ARMA (p, q) Modelle vir Tydreeksanalise - Deel 2 Deur Michael Saal-Moore op 24 Augustus 2015 in Deel 1 beskou ons die outoregressiewe model van orde p, ook bekend as die AR (p) model. Ons lei dit as 'n uitbreiding van die ewekansige loop model in 'n poging om bykomende reeks korrelasie in finansiële tydreekse verduidelik. Uiteindelik het ons besef dat dit was nie buigsaam genoeg om werklik al die outokorrelasie te vang in die laaste pryse van Amazon Inc. (AMZN) en die SampP500 VSA Equity Index. Die primêre rede hiervoor is dat beide van hierdie bates is voorwaardelik heteroskedastic. wat beteken dat hulle nie-stasionêre en het tydperke van wisselende variansie of wisselvalligheid groepering, wat nie in ag geneem word deur die AR (p) model geneem. In toekomstige artikels sal ons uiteindelik opbou tot die outoregressiewe geïntegreerde bewegende gemiddelde (ARIMA) modelle, asook die voorwaardelik heteroskedastic modelle van die boog en GARCH families. Hierdie modelle sal ons met ons eerste realistiese pogings tot vooruitskatting batepryse. In hierdie artikel, maar ons gaan die bewegende gemiddelde van orde q model, bekend as MA (Q) bekend te stel. Dit is 'n komponent van die meer algemene ARMA model en as sulks het ons nodig het om dit te verstaan ​​voordat verdere beweeg. Ek raai jy die vorige artikels gelees in die Tydreeksanalise versameling as jy dit nog nie gedoen. Hulle kan al hier gevind word. Bewegende gemiddelde (MA) Models van orde q Rasionaal 'n bewegende gemiddelde model is soortgelyk aan 'n outoregressiewe model, behalwe dat in plaas daarvan om 'n lineêre kombinasie van die verlede tyd reeks waardes, dit is 'n lineêre kombinasie van die afgelope wit geraas terme. Intuïtief, beteken dit dat die MA model sien soos ewekansige wit geraas skokke direk by elke huidige waarde van die model. Dit is in teenstelling met 'n AR (p) model, waar die wit geraas skokke slegs indirek gesien. via regressie op vorige terme van die reeks. 'N Belangrike verskil is dat die MA-model net ooit sal sien die laaste Q skokke vir 'n spesifieke MA (Q) model, terwyl die AR (p) model al voor skokke in ag sal neem, al is dit in 'n decreasingly swak wyse. Definisie Wiskundig die MA (Q) is 'n lineêre regressiemodel en is insgelyks gestruktureer om AR (p): Moving Gemiddelde Model van orde q 'n tydreeksmodel, is 'n bewegende gemiddelde model van orde q. MA (Q), indien: begin xt wt beta1 w ldots betaq w end Waar is wit geraas met E (WT) 0 en variansie sigma2. As ons kyk na die agterste Shift-operateur. (Sien 'n vorige artikel) dan kan ons herskryf bogenoemde as 'n funksie phi van: begin xt (1 beta1 beta2 2 ldots betaq Q) wt phiq () wt einde Ons sal gebruik maak van die phi-funksie in die latere artikels te maak. Tweede Orde Properties Soos met AR (p) die gemiddelde van 'n MA (Q) proses is nul. Dit is maklik om te sien as die gemiddelde is bloot 'n som van middel van wit geraas terme, wat al self nul is. begin teks enspace MUX E (xt) som E (Wi) 0 einde begin teks enspace sigma2w (1 beta21 ldots beta2q) einde teks enspace rhok links 1 teks enspace k 0 som betai beta / sumq beta2i teks enspace k 1, ldots, Q 0 teks enspace k GT Q einde reg. Waar beta0 1. Was nou gaan 'n paar gesimuleerde data te genereer en gebruik dit om correlograms skep. Dit sal die formule hierbo vir rhok ietwat meer beton. Simulasies en Correlograms MA (1) Kom ons begin met 'n MA (1) proses. As ons 'beta1 0.6 verkry ons die volgende model: Soos met die AR (p) modelle in die vorige artikel kan ons R te gebruik om so 'n reeks te boots en dan trek die correlogram. Sedert weve het 'n baie oefening in die vorige Tydreeksanalise artikel reeks van die uitvoering van erwe, sal ek die R-kode skryf ten volle, eerder as om te verdeel dit: Die produksie is soos volg: Soos ons hierbo gesien het in die formule vir rhok , vir k GT Q, al outokorrelasies moet nul wees. Sedert Q 1, moet ons 'n beduidende hoogtepunt op k1 en dan onbelangrik pieke na daardie sien. As gevolg van steekproefneming vooroordeel ons moet verwag om 5 (effens) beduidende pieke sien op 'n monster outokorrelasie plot. Dit is presies wat die correlogram wys vir ons in hierdie geval. Ons het 'n beduidende hoogtepunt op k1 en dan onbelangrik pieke vir k GT 1, behalwe by K4 waar ons 'n effens beduidende piek. Trouens, dit is 'n nuttige manier om te sien of 'n MA (Q) model toepaslik is. Deur die neem van 'n blik op die correlogram van 'n bepaalde reeks kan ons sien hoeveel opeenvolgende nie-nul lags bestaan. As Q so lags bestaan ​​dan kan ons tereg probeer om 'n MA (Q) model geskik is om 'n bepaalde reeks. Aangesien ons bewyse uit ons gesimuleerde data van 'n MA (1) proses, is nou van plan om te probeer en pas 'n MA (1) model vir ons gesimuleerde data. Ongelukkig is daar isnt 'n ekwivalente ma opdrag om die outoregressiewe model ar opdrag in R. In plaas daarvan, moet ons die meer algemene ARIMA opdrag gebruik en stel die outoregressiewe en geïntegreerde komponente aan nul. Ons doen dit deur die skep van 'n 3-vektor en die opstel van die eerste twee komponente (die autogressive en geïntegreerde parameters, onderskeidelik) na nul: Ons ontvang 'n paar nuttige uitset van die ARIMA opdrag. Eerstens, kan ons sien dat die parameter is beraam as hoed 0,602, wat baie naby aan die werklike waarde van beta1 0.6. In die tweede plek is die standaard foute reeds bereken vir ons, maak dit maklik om vertrouensintervalle bereken. Derdens, ontvang ons 'n geskatte variansie, log-waarskynlikheid en Akaike Inligting Criterion (wat nodig is vir model vergelyking). Die groot verskil tussen ARIMA en ar is dat ARIMA skat 'n onderskepdrie termyn omdat dit nie die gemiddelde waarde van die reeks af te trek. Vandaar ons nodig het om versigtig te wees wanneer die uitvoering van voorspellings met behulp van die ARIMA opdrag. Wel later terug te keer na hierdie punt. As 'n vinnige check op pad was om vertrouensintervalle vir hoed bereken: Ons kan sien dat die 95 vertrouensinterval bevat die ware parameter waarde van beta1 0.6 en daarom het ons die model kan oordeel 'n goeie passing. Dit is duidelik dat hierdie moet verwag word, aangesien ons die data nageboots in die eerste plek Hoe dinge verander as ons die teken van beta1 om -0,6 verander Kom dieselfde analise uit te voer: Die produksie is soos volg: Ons kan sien dat by k1 ons het 'n beduidende hoogtepunt in die correlogram, behalwe dat dit toon negatiewe korrelasie, as wed verwag van 'n MA (1) model met negatiewe eerste koëffisiënt. Weereens al pieke buite k1 is onbelangrik. Kom ons pas 'n MA (1) model en skat die parameter: hoed -0,730, wat is 'n klein onderskat van beta1 -0,6. Ten slotte, kan bereken die vertroue interval: Ons kan sien dat die ware parameter waarde van beta1-0.6 is vervat in die 95 vertrouensinterval, die verskaffing van ons met bewyse van 'n goeie model pas. MA (3) Kom ons loop deur dieselfde prosedure vir 'n MA (3) proses. Hierdie keer moet ons beduidende hoogtepunte op k verwag, en onbelangrik pieke vir k GT 3. Ons gaan die volgende koëffisiënte gebruik: beta1 0.6, beta2 0.4 en beta3 0.2. Kom na te boots 'n MA (3) proses van hierdie model. Ive het die aantal ewekansige monsters tot 1000 in hierdie simulasie, wat dit makliker maak om die ware outokorrelasie struktuur sien, ten koste van die maak van die oorspronklike reeks moeiliker om te interpreteer: Die produksie is soos volg: Soos verwag die eerste drie pieke is beduidende . Maar so is die vierde. Maar ons kan tereg daarop dui dat hierdie mag wees as gevolg van steekproefneming vooroordeel soos ons verwag om te sien 5 van die pieke wat beduidende buite KQ. Kom nou pas 'n MA (3) model om die data te probeer en skatting parameters: Die skattings hoed 0,544, hoed 0,345 en hoed 0,298 is naby aan die ware waardes van beta10.6, beta20.4 en beta30.3, onderskeidelik. Ons kan produseer ook vertrouensintervalle gebruik van die onderskeie standaard foute: In elk geval nie die 95 vertrouensintervalle die ware parameter waarde bevat en kan ons aflei dat ons 'n goeie passing met ons MA (3) model, soos verwag kan word. Finansiële inligting in Deel 1 beskou ons Amazon Inc. (AMZN) en die SampP500 VSA Equity Index. Ons toegerus die AR (p) model vir beide en gevind dat die model nie in staat was om effektief te vang die kompleksiteit van die reeks korrelasie, veral in die rolverdeling van die SampP500, waar langtermyn-geheue effekte blyk teenwoordig te wees. Ek sal nie stip die kaarte weer vir die pryse en outokorrelasie, in plaas Siek verwys u na die vorige post. Amazon Inc. (AMZN) Kom ons begin deur te probeer om 'n seleksie van MA (Q) pas modelle om AMZN, naamlik met Q in. Soos in Deel 1, goed gebruik quantmod om die daaglikse pryse vir AMZN aflaai en dan sit hulle in 'n log opbrengste stroom sluitingstyd pryse: Noudat ons die log opbrengste stroom kan ons die ARIMA opdrag gebruik om in te pas MA (1), MA (2) en MA (3) modelle en dan skat die parameters van elke. Vir MA (1) ons het: Ons kan die residue van die daaglikse log opbrengste en die toegeruste model plot: Let daarop dat ons 'n paar beduidende hoogtepunte op lags k2, K11, K16 en k18, wat aandui dat die MA (1) model is onwaarskynlik dat 'n goeie passing vir die gedrag van die AMZN log opbrengste wees, aangesien dit lyk nie soos 'n verwesenliking van wit geraas. Kom ons probeer 'n MA (2) model: Beide van die skattings vir die beta koëffisiënte is negatief. Kom ons plot die residue weer: Ons kan sien dat daar byna nul outokorrelasie in die eerste paar lags. Ons het egter vyf effens beduidende hoogtepunte op lags K12, K16, K19, k25 en K27. Dit is suggestief dat die MA (2) model is die opneem van 'n groot deel van die outokorrelasie, maar nie almal van die lang-geheue effekte. Hoe gaan dit met 'n MA (3) model Weereens, kan ons die residue Plot: Die MA (3) residue plot lyk byna identies aan dié van die MA (2) model. Dit is nie verbasend nie, as 'n nuwe parameter is die toevoeging van 'n model wat skynbaar weg verduidelik baie van die korrelasies met korter lags, maar dit sal nie veel van 'n uitwerking op die langer termyn loop. Al hierdie getuienis is suggestief van die feit dat 'n MA (Q) model is onwaarskynlik handig al die korrelasie in isolasie te wees. ten minste vir AMZN. SampP500 As jy onthou, in Deel 1 het ons gesien dat die eerste orde differenced daaglikse log opbrengste struktuur van die SampP500 besit baie beduidende pieke op verskillende lags, beide kort en lang. Dit verskaf bewyse van beide voorwaardelike heteroskedasticity (dit wil sê wisselvalligheid groepering) en langtermyn-geheue effekte. Dit lei ons tot die gevolgtrekking dat die AR (p) model onvoldoende is om al die outokorrelasie teenwoordig te vang was. Soos weve bo die MA (Q) model gesien was onvoldoende om bykomende reeks korrelasie in die residue van die toegeruste model om die eerste orde vang differenced daaglikse log prys reeks. Ons sal nou probeer om die MA (Q) model om die SampP500 pas. Mens kan vra waarom ons doen dit as ons weet dat dit is onwaarskynlik dat 'n goeie passing wees. Dit is 'n goeie vraag. Die antwoord is dat ons nodig het om te sien presies hoe dit is nie 'n goeie passing, want dit is die uiteindelike proses sal ons volgende wanneer ons teëkom baie meer gesofistikeerd modelle, wat potensieel moeiliker om te interpreteer. Kom ons begin deur die verkryging van die data en dit na 'n eerste orde differenced reeks logaritmies getransformeer daaglikse sluitingspryse soos in die vorige artikel: Ons gaan nou 'n MA (1), MA (2) en MA (3) model aan te pas die reeks, soos ons hierbo gedoen het vir AMZN. Kom ons begin met MA (1): Kom ons maak 'n plot van die residue van hierdie toegeruste model: Die eerste beduidende piek plaasvind op k2, maar daar is baie meer aan k in. Dit is duidelik nie 'n besef van wit geraas en so moet ons die MA (1) model as 'n potensiële goeie passing vir die SampP500 verwerp. Maak die situasie te verbeter met MA (2) Weereens, kan 'n plot van die residue van hierdie toegeruste MA (2) model: Terwyl die piek by K2 verdwyn (soos wed verwag), is ons nog steeds links met die beduidende hoogtepunte op baie meer lags in die residue. Weereens, vind ons die MA (2) model is nie 'n goeie passing. Ons moet verwag nie, want die MA (3) model, minder korrelasie by K3 as vir die MA (2) te sien, maar ons moet ook verwag weereens geen vermindering in verdere lags. Ten slotte, kan 'n plot van die residue van hierdie toegeruste MA (3) model: Dit is presies wat ons sien in die correlogram van die residue. Vandaar die MA (3), met die ander modelle hierbo, is nie 'n goeie passing vir die SampP500. Volgende stappe Weve nou ondersoek twee groot tydreeksmodelle in detail, naamlik die Autogressive model van orde p, AR (p) en dan bewegende gemiddelde van orde q, MA (Q). Weve gesien dat hulle is albei in staat te verduidelik weg van die outokorrelasie in die residue van eerste orde differenced daaglikse log pryse van aandele en indekse, maar wisselvalligheid groepering en langtermyn-geheue effek voortduur. Dit is uiteindelik tyd om ons aandag te draai na die kombinasie van hierdie twee modelle, naamlik die outoregressiewe bewegende gemiddelde van orde p, q, ARMA (p, q) om te sien of dit die situasie verder sal verbeter. Ons sal egter moet wag tot die volgende artikel vir 'n volledige bespreking Michael Saal-Moore Mike is die stigter van QuantStart en is betrokke by die kwantitatiewe finansiële sektor vir die afgelope vyf jaar, in die eerste plek as 'n quant ontwikkelaar en later as 'n quant handelaar konsultasie vir verskansing funds. A Rima staan ​​vir outoregressiewe geïntegreerde bewegende gemiddelde modelle. Eenveranderlike (enkele vektor) ARIMA is 'n vooruitskatting tegniek wat die toekomstige waardes van 'n reeks ten volle gebaseer op sy eie traagheid projekte. Die belangrikste aansoek is op die gebied van korttermyn voorspelling wat ten minste 40 historiese data punte. Dit werk die beste wanneer jou data toon 'n stabiele of konsekwent patroon met verloop van tyd met 'n minimum bedrag van uitskieters. Soms genoem word Posbus-Jenkins (ná die oorspronklike skrywers), ARIMA is gewoonlik beter as gladstrykingstegnieke eksponensiële wanneer die data is redelik lank en die korrelasie tussen die verlede waarnemings is stabiel. As die data is kort of baie volatiel, dan kan 'n paar smoothing metode beter te presteer. As jy nie ten minste 38 datapunte het, moet jy 'n ander metode as ARIMA oorweeg. Die eerste stap in die toepassing van ARIMA metode is om te kyk vir stasionariteit. Stasionariteit impliseer dat die reeks bly op 'n redelik konstante vlak met verloop van tyd. As 'n tendens bestaan, soos in die meeste ekonomiese of besigheid aansoeke, dan is jou data nie stilstaan. Die data moet ook 'n konstante stryd in sy skommelinge oor tyd te wys. Dit is maklik gesien met 'n reeks wat swaar seisoenale en groei teen 'n vinniger tempo. In so 'n geval, sal die wel en wee van die seisoen meer dramaties met verloop van tyd. Sonder hierdie stasionariteit voorwaardes voldoen word, baie van die berekeninge wat verband hou met die proses kan nie bereken word nie. As 'n grafiese plot van die data dui stationariteit, dan moet jy verskil die reeks. Breukmetodes is 'n uitstekende manier om die transformasie van 'n nie-stationaire reeks om 'n stilstaande een. Dit word gedoen deur die aftrekking van die waarneming in die huidige tydperk van die vorige een. As hierdie transformasie slegs een keer gedoen word om 'n reeks, sê jy dat die data het eers differenced. Hierdie proses elimineer wese die tendens as jou reeks groei teen 'n redelik konstante tempo. As dit groei teen 'n vinniger tempo, kan jy dieselfde prosedure en verskil die data weer aansoek doen. Jou data sal dan tweede differenced. Outokorrelasies is numeriese waardes wat aandui hoe 'n data-reeks is wat verband hou met self met verloop van tyd. Meer presies, dit meet hoe sterk datawaardes op 'n bepaalde aantal periodes uitmekaar gekorreleer met mekaar oor tyd. Die aantal periodes uitmekaar is gewoonlik bekend as die lag. Byvoorbeeld, 'n outokorrelasie op lag 1 maatreëls hoe waardes 1 tydperk uitmekaar gekorreleer met mekaar oor die hele reeks. 'N outokorrelasie op lag 2 maatreëls hoe die data twee periodes uitmekaar gekorreleer regdeur die reeks. Outokorrelasies kan wissel van 1 tot -1. 'N Waarde naby aan 1 dui op 'n hoë positiewe korrelasie, terwyl 'n waarde naby aan -1 impliseer 'n hoë negatiewe korrelasie. Hierdie maatreëls is meestal geëvalueer deur middel van grafiese plotte genoem correlagrams. A correlagram plotte die motor - korrelasie waardes vir 'n gegewe reeks by verskillende lags. Dit staan ​​bekend as die outokorrelasie funksie en is baie belangrik in die ARIMA metode. ARIMA metode poog om die bewegings in 'n stilstaande tyd reeks beskryf as 'n funksie van wat is outoregressiewe en bewegende gemiddelde parameters genoem. Dit is waarna verwys word as AR parameters (autoregessive) en MA parameters (bewegende gemiddeldes). 'N AR-model met slegs 1 parameter kan geskryf word as. X (t) 'n (1) X (t-1) E (t) waar x (t) tydreekse wat ondersoek word 'n (1) die outoregressiewe parameter van orde 1 X (t-1) die tydreeks uitgestel 1 periode E (t) die foutterm van die model beteken dit eenvoudig dat enige gegewe waarde X (t) kan verduidelik word deur 'n funksie van sy vorige waarde, X (t-1), plus 'n paar onverklaarbare ewekansige fout, E (t). As die beraamde waarde van A (1) was 0,30, dan is die huidige waarde van die reeks sal wees met betrekking tot 30 van sy waarde 1 periode gelede. Natuurlik, kan die reeks word wat verband hou met meer as net 'n verlede waarde. Byvoorbeeld, X (t) 'n (1) X (t-1) A (2) X (t-2) E (t) Dit dui daarop dat die huidige waarde van die reeks is 'n kombinasie van die twee onmiddellik voorafgaande waardes, X (t-1) en X (t-2), plus 'n paar random fout E (t). Ons model is nou 'n outoregressiewe model van orde 2. bewegende gemiddelde modelle: 'n Tweede tipe Box-Jenkins model is 'n bewegende gemiddelde model genoem. Hoewel hierdie modelle lyk baie soortgelyk aan die AR model, die konsep agter hulle is heel anders. Bewegende gemiddelde parameters verband wat gebeur in tydperk t net om die ewekansige foute wat plaasgevind het in die verlede tyd periodes, naamlik E (t-1), E (t-2), ens, eerder as om X (t-1), X ( t-2), (xt-3) as in die outoregressiewe benaderings. 'N bewegende gemiddelde model met 'n MA termyn kan soos volg geskryf word. X (t) - B (1) E (t-1) E (t) Die term B (1) genoem word 'n MA van orde 1. Die negatiewe teken voor die parameter is slegs vir konvensie en word gewoonlik gedruk uit motor - dateer deur die meeste rekenaarprogramme. Bogenoemde model eenvoudig sê dat enige gegewe waarde van X (t) direk verband hou net aan die ewekansige fout in die vorige tydperk, E (t-1), en die huidige foutterm, E (t). Soos in die geval van outoregressiemodelle, kan die bewegende gemiddelde modelle uitgebrei word na 'n hoër orde strukture wat verskillende kombinasies en bewegende gemiddelde lengtes. ARIMA metode kan ook modelle gebou word dat beide outoregressiewe en gemiddelde parameters saam beweeg inkorporeer. Hierdie modelle word dikwels na verwys as gemengde modelle. Hoewel dit maak vir 'n meer ingewikkelde voorspelling instrument, kan die struktuur inderdaad die reeks beter na te boots en produseer 'n meer akkurate skatting. Suiwer modelle impliseer dat die struktuur bestaan ​​slegs uit AR of MA parameters - nie beide. Die ontwikkel deur hierdie benadering modelle word gewoonlik genoem ARIMA modelle omdat hulle 'n kombinasie van outoregressiewe (AR) te gebruik, integrasie (I) - verwys na die omgekeerde proses van breukmetodes die voorspelling te produseer, en bewegende gemiddelde (MA) operasies. 'N ARIMA model word gewoonlik gestel as ARIMA (p, d, q). Dit verteenwoordig die orde van die outoregressiewe komponente (p), die aantal breukmetodes operateurs (d), en die hoogste orde van die bewegende gemiddelde termyn. Byvoorbeeld, ARIMA (2,1,1) beteken dat jy 'n tweede orde outoregressiewe model met 'n eerste orde bewegende gemiddelde komponent waarvan die reeks is differenced keer om stasionariteit veroorsaak. Pluk die reg spesifikasie: Die grootste probleem in die klassieke Box-Jenkins probeer om te besluit watter ARIMA spesifikasie gebruik - i. e. hoeveel AR en / of MA parameters in te sluit. Dit is wat die grootste deel van Box-Jenkings 1976 is gewy aan die identifikasieproses. Dit was afhanklik van grafiese en numeriese eval - uation van die monster outokorrelasie en gedeeltelike outokorrelasiefunksies. Wel, vir jou basiese modelle, die taak is nie te moeilik. Elk outokorrelasiefunksies dat 'n sekere manier te kyk. Maar wanneer jy optrek in kompleksiteit, die patrone is nie so maklik opgespoor. Om sake nog moeiliker maak, jou data verteenwoordig slegs 'n voorbeeld van die onderliggende proses. Dit beteken dat steekproeffoute (uitskieters, meting fout, ens) die teoretiese identifikasie proses kan verdraai. Dit is waarom tradisionele ARIMA modellering is 'n kuns eerder as 'n science.2.1 bewegende gemiddelde modelle (MA modelle) tydreeksmodelle bekend as ARIMA modelle kan die volgende insluit outoregressiewe terme en / of bewegende gemiddelde terme. In Week 1, het ons geleer 'n outoregressiewe term in 'n tydreeks model vir die veranderlike x t is 'n vertraagde waarde van x t. Byvoorbeeld, 'n lag 1 outoregressiewe termyn is x t-1 (vermenigvuldig met 'n koëffisiënt). Hierdie les definieer bewegende gemiddelde terme. 'N bewegende gemiddelde termyn in 'n tydreeks model is 'n verlede fout (vermenigvuldig met 'n koëffisiënt). Laat (WT omslaan N (0, sigma2w)), wat beteken dat die w t is identies, onafhanklik versprei, elk met 'n normaalverdeling met gemiddelde 0 en dieselfde afwyking. Die 1 ste orde bewegende gemiddelde model, aangedui deur MA (1) is (xt mu wt theta1w) Die 2de orde bewegende gemiddelde model, aangedui deur MA (2) is (xt mu wt theta1w theta2w) Die Q de orde bewegende gemiddelde model , aangedui deur MA (Q) is (xt mu wt theta1w theta2w kolle thetaqw) Nota. Baie handboeke en sagteware programme definieer die model met negatiewe tekens voor die terme. Dit nie die geval verander die algemene teoretiese eienskappe van die model, hoewel dit flip die algebraïese tekens van beraamde koëffisiënt waardes en (unsquared) terme in formules vir ACFs en afwykings. Jy moet jou sagteware kyk om te kontroleer of negatiewe of positiewe tekens is gebruik om korrek te skryf die beraamde model. R gebruik positiewe tekens in sy onderliggende model, soos ons hier doen. Teoretiese Eienskappe van 'n tydreeks met 'n MA (1) Model Let daarop dat die enigste nie-nul waarde in die teoretiese ACF is vir lag 1. Alle ander outokorrelasies is 0. So 'n monster ACF met 'n beduidende outokorrelasie net by lag 1 is 'n aanduiding van 'n moontlike MA (1) model. Vir belangstellende studente, bewyse van hierdie eienskappe is 'n bylae tot hierdie opdragstuk. Voorbeeld 1 Veronderstel dat 'n MA (1) model is x t 10 w t 0,7 w t-1. waar (WT omslaan N (0,1)). So het die koëffisiënt 1 0.7. Die teoretiese ACF gegee word deur 'n plot van hierdie volg ACF. Die plot net aangedui is die teoretiese ACF vir 'n MA (1) met 1 0.7. In die praktyk, 'n monster gewoond gewoonlik verskaf so 'n duidelike patroon. Die gebruik van R, gesimuleerde ons N 100 monster waardes gebruik te maak van die model x t 10 w t 0,7 w t-1 waar w t IID N (0,1). Vir hierdie simulasie, 'n tydreeks plot van die steekproefdata volg. Ons kan nie sê baie van hierdie plot. Die monster ACF vir die gesimuleerde data volg. Ons sien 'n skerp styging in lag 1 gevolg deur die algemeen nie-beduidende waardes vir lags afgelope 1. Let daarop dat die monster ACF kom nie ooreen met die teoretiese patroon van die onderliggende MA (1), en dit is dat al outokorrelasies vir lags afgelope 1 sal wees 0 . 'n ander voorbeeld sou 'n effens verskillende monster ACF hieronder getoon, maar sal waarskynlik dieselfde breë funksies. Theroretical Eienskappe van 'n tydreeks met 'n MA (2) model vir die MA (2) model, teoretiese eienskappe is soos volg: Let daarop dat die enigste nie-nul waardes in die teoretiese ACF is vir lags 1 en 2. outokorrelasies vir hoër lags is 0 . So, 'n monster ACF met 'n beduidende outokorrelasies by lags 1 en 2, maar nie-beduidende outokorrelasies vir hoër lags dui op 'n moontlike MA (2) model. IID N (0,1). Die koëffisiënte is 1 0.5 en 2 0.3. Want dit is 'n MA (2), sal die teoretiese ACF nul waardes het net by lags 1 en 2. Waardes van die twee nie-nul outokorrelasies is 'n plot van die teoretiese ACF volg. Soos byna altyd die geval is, monster data gewoond te tree heeltemal so perfek as teorie. Ons gesimuleerde N 150 monster waardes vir die model x t 10 w t 0,5 w t-1 0,3 w t-2. waar w t IID N (0,1). Die tydreekse plot van die data volg. Soos met die tydreeks plot vir die MA (1) voorbeeld van die data, kan nie vir jou sê baie daaruit. Die monster ACF vir die gesimuleerde data volg. Die patroon is tipies vir situasies waar 'n MA (2) model nuttig kan wees. Daar is twee statisties beduidende spykers by lags 1 en 2, gevolg deur nie-beduidende waardes vir ander lags. Let daarop dat as gevolg van steekproeffout, die monster ACF nie die teoretiese patroon presies ooreenstem. ACF vir Algemene MA (Q) Models n eiendom van MA (Q) modelle in die algemeen is dat daar nie-nul outokorrelasies vir die eerste Q lags en outokorrelasies 0 vir alle lags GT q. Nie-uniekheid van verband tussen waardes van 1 en (rho1) in MA (1) Model. In die MA (1) model, vir enige waarde van 1. die wedersydse 01/01 gee dieselfde waarde vir so 'n voorbeeld, gebruik 0,5 vir 1. en gebruik dan 1 / (0,5) 2 vir 1. Jy sal kry (rho1) 0.4 in beide gevalle. Om 'n teoretiese beperking genoem inverteerbaarheid bevredig. Ons beperk MA (1) modelle om waardes met absolute waarde minder as 1. In die voorbeeld net gegee, 1 0.5 sal 'n toelaatbare parameter waarde wees nie, terwyl 1 1 / 0.5 2 nie. Inverteerbaarheid van MA modelle 'n MA-model word gesê omkeerbare te wees indien dit algebraïes gelykstaande aan 'n konvergerende oneindige orde AR model. Bevestig deur die, bedoel ons dat die AR koëffisiënte daal tot 0 as ons terug beweeg in die tyd. Inverteerbaarheid is 'n beperking geprogrammeer in die tyd reeks sagteware wat gebruik word om die koëffisiënte van modelle te skat met MA terme. Dit is nie iets wat ons gaan vir die data-analise. Bykomende inligting oor die inverteerbaarheid beperking vir MA (1) modelle word in die bylaag. Gevorderde teorie Nota. Vir 'n MA (Q) model met 'n bepaalde ACF, daar is net een omkeerbare model. Die noodsaaklike voorwaarde vir inverteerbaarheid is dat die koëffisiënte waardes sodanig dat die vergelyking 1- 1 y. - Q y q 0 het oplossings vir y wat buite die eenheidsirkel val. R-kode vir die voorbeelde in Voorbeeld 1, ons geplot die teoretiese ACF van die model x t 10 w t. 7W t-1. en dan nageboots N 150 waardes van hierdie model en geplot die monster tydreekse en die monster ACF vir die gesimuleerde data. Die R bevele gebruik word om die teoretiese ACF plot was: acfma1ARMAacf (Mac (0,7), lag. max10) 10 lags van ACF vir MA (1) met theta1 0.7 lags0: 10 skep 'n veranderlike genaamd lags wat wissel van 0 tot 10. plot (lags, acfma1, xlimc (1,10), ylabr, typeh, hoof ACF vir MA (1) met theta1 0.7) abline (H0) voeg n horisontale as om die plot die eerste opdrag bepaal die ACF en slaan dit in 'n voorwerp vernoem acfma1 (ons keuse van naam). Die plot opdrag (die 3de gebod) erwe lags teenoor die ACF waardes vir lags 1 tot 10. Die ylab parameter etikette die y-as en die belangrikste parameter sit 'n titel op die plot. Om te sien die numeriese waardes van die ACF net gebruik die opdrag acfma1. Die simulasie en erwe is gedoen met die volgende opdragte. xcarima. sim (N150, lys (Mac (0,7))) Simuleer N 150 waardes van MA (1) xxc10 voeg 10 tot gemiddelde 10. Simulasie gebreke maak beteken 0. plot (x, typeb, mainSimulated MA (1) data) ACF (x, xlimc (1,10), mainACF vir gesimuleerde steekproefdata) In Voorbeeld 2, ons geplot die teoretiese ACF van die model xt 10 wt 0,5 w t-1 0,3 w t-2. en dan nageboots N 150 waardes van hierdie model en geplot die monster tydreekse en die monster ACF vir die gesimuleerde data. Die R bevele gebruik was acfma2ARMAacf (Mac (0.5,0.3), lag. max10) acfma2 lags0: 10 plot (lags, acfma2, xlimc (1,10), ylabr, typeh, hoof ACF vir MA (2) met theta1 0.5, theta20.3) abline (H0) xcarima. sim (N150, lys (Mac (0.5, 0.3))) xxc10 plot (x, typeb, hoof Gesimuleerde MA (2) Series) ACF (x, xlimc (1,10), mainACF vir gesimuleerde MA (2) Data) Bylae: Bewys van eiendomme van MA (1) vir belangstellende studente, hier is bewyse vir teoretiese eienskappe van die MA (1) model. Variansie: (teks (xt) teks (mu wt theta1 w) 0 teks (WT) teks (theta1w) sigma2w theta21sigma2w (1theta21) sigma2w) Wanneer h 1, die vorige uitdrukking 1 W 2. Vir enige h 2, die vorige uitdrukking 0 . die rede hiervoor is dat per definisie van onafhanklikheid van die WT. E (w k w j) 0 vir enige k j. Verder, omdat die w t het intussen 0, E (w j w j) E (w j 2) w 2. Vir 'n tydreeks, Pas hierdie resultaat aan die ACF hierbo kry. 'N omkeerbare MA model is die een wat geskryf kan word as 'n oneindige orde AR model wat konvergeer sodat die AR koëffisiënte konvergeer na 0 as ons oneindig terug in die tyd beweeg. Wel demonstreer inverteerbaarheid vir die MA (1) model. Ons het toe plaasvervanger verhouding (2) vir w t-1 in vergelyking (1) (3) (ZT wt theta1 (Z - theta1w) wt theta1z - theta2w) op tydstip t-2. vergelyking (2) word Ons het toe plaasvervanger verhouding (4) vir w t-2 in vergelyking (3) (ZT wt theta1 Z - theta21w wt theta1z - theta21 (Z - theta1w) wt theta1z - theta12z theta31w) As ons voortgaan ( oneindig), sou ons die oneindige orde AR model kry (ZT wt theta1 Z - theta21z theta31z - theta41z kolletjies) Nota egter dat as 1 1, die koëffisiënte die lags van Z vermenigvuldig sal toeneem (oneindig) in grootte as ons terug beweeg in tyd. Om dit te voorkom, moet ons 1 LT1. Dit is die voorwaarde vir 'n omkeerbare MA (1) model. Oneindige Bestel MA model In week 3, goed sien dat 'n AR (1) model kan omgeskakel word na 'n oneindige orde MA model: (xt - mu wt phi1w phi21w kolle phik1 w kolle som phij1w) Hierdie opsomming van verlede wit geraas terme is bekende as die oorsaaklike voorstelling van 'n AR (1). Met ander woorde, x t is 'n spesiale tipe MA met 'n oneindige aantal terme terug gaan in die tyd. Dit is 'n oneindige orde MA of MA () genoem. 'N Eindige orde MA is 'n oneindige orde AR en enige eindige orde AR is 'n oneindige orde MA. Onthou in Week 1, het ons opgemerk dat 'n vereiste vir 'n stilstaande AR (1) is dat 1 LT1. Kom ons bereken die Var (x t) met behulp van die oorsaaklike verteenwoordiging. Die laaste stap gebruik 'n basiese feit oor meetkundige reeks wat vereis (phi1lt1) anders sal die reeks divergeer. NavigationDocumentation is die onvoorwaardelike gemiddelde van die proses, en x03C8 (L) is 'n rasionele, oneindige-graad lag operateur polinoom, (1 x03C8 1 L x03C8 2 L 2 x2026). Let wel: Die konstante eienskap van 'n ARIMA model voorwerp ooreenstem met c. en nie die onvoorwaardelike gemiddelde 956. Deur Wolds ontbinding 1. Vergelyking 5-12 ooreenstem met 'n stilstaande stogastiese proses op voorwaarde dat die koëffisiënte x03C8 Ek is absoluut summable. Dit is die geval wanneer die AR polinoom, x03D5 (L). is stabiel. wat beteken dat al sy wortels lê buite die eenheidsirkel. Daarbenewens het die proses is kousale op voorwaarde dat die MA polinoom is omkeerbaar. wat beteken dat al sy wortels lê buite die eenheidsirkel. Ekonometrie Gereedskap dwing stabiliteit en inverteerbaarheid van ARMA prosesse. Wanneer jy 'n ARMA model spesifiseer met behulp van ARIMA. jy 'n fout as jy koëffisiënte wat nie ooreenstem met 'n stabiele AR polinoom of omkeerbare MA polinoom betree. Net so, skat lê stasionariteit en inverteerbaarheid beperkings tydens beraming. Verwysings 1 Wold, H. 'n studie in die ontleding van tydreekse. Uppsala, Swede: Almqvist amp Wiksell, 1938. Kies 'n land