Tidsserieanalys Exempel Vilka srdrag har tidsseriedata Varfr behvs
Tidsserieanalys Exempel: Vilka särdrag har tidsseriedata? Varför behövs nya metoder? 1
Vilka särdrag har tidsseriedata? Varför behövs nya metoder? • Observationerna är inte oberoende, eftersom det finns ett beroende i tiden • Observationerna ger ett mönster över tiden – en trend: fallande eller stigande värden med tiden – en periodisk variation över en tidsperiod av bestämd längd (säsongseffekter eller liknande) 2
Exempel på tidsseriedata Olika typer av ekonomiska data: – Arbetslöshetssiffror – Försäljningsvärden – Konsumentprisindex och andra index – Export- och importmängder Miljömätdata: – Fosforhalt i havsvattenbassänger – Ozonhalt i luftrummet över en storstad Medicinska data: – Antal fall av viss sjukdom (influensa, påssjuka. . . ) 3
4
5
Tidsserieregression: 6
sold time month x 1 x 2 x 3 x 4 x 5 x 6 x 7 x 8 x 9 x 10 x 11 2 1 1 1 0 0 0 0 0 6 2 2 0 1 0 0 0 0 0 5 3 3 0 0 1 0 0 0 0 5 4 4 0 0 0 1 0 0 0 0 10 5 5 0 0 1 0 0 0 8 6 6 0 0 0 1 0 0 0 . . . . 7 46 10 0 0 0 0 1 0 6 47 11 0 0 0 0 0 1 3 48 12 0 0 0 7
Tidsserieregression med enbart en linjär trend. 8
Regression Analysis The regression equation is Sold = 5, 78 + 0, 0430 time Predictor St. Dev T P 5, 7761 0, 9429 6, 13 0, 000 0, 04302 0, 03420 1, 26 0, 215 Constant Coef time S = 3, 181 R-Sq = 3, 4% R-Sq(adj) = 1, 2% Analysis of Variance Source DF SS MS F P 1 16, 00 1, 58 0, 215 Residual Error 45 455, 27 10, 12 Total 46 471, 28 Regression 9 Tidsserieregression fungerar statistiskt som vanlig regression.
10
Regression Analysis The regression equation is Sold = 3, 65 + 0, 0285 time - 1, 69 x 1 - 0, 47 x 2 + 2, 75 x 3 + 1, 22 x 4 + 6, 20 x 5 + 2, 42 x 6 + 8, 14 x 7 + 6, 36 x 8 + 0, 58 x 9 + 2, 55 x 10 + 1, 02 x 11 Predictor Constant time x 1 x 2 x 3 x 4 x 5 x 6 x 7 x 8 x 9 x 10 x 11 S = 1, 342 Coef 3, 6491 0, 02851 -1, 691 -0, 469 2, 752 1, 224 6, 195 2, 417 8, 138 6, 360 0, 581 2, 553 1, 024 St. Dev 0, 8526 0, 01481 1, 028 1, 027 1, 026 1, 025 1, 026 1, 027 1, 028 R-Sq = 87, 0% T 4, 28 1, 92 -1, 65 -0, 46 2, 68 1, 19 6, 04 2, 36 7, 94 6, 20 0, 57 2, 49 1, 00 P 0, 000 0, 063 0, 109 0, 651 0, 011 0, 241 0, 000 0, 024 0, 000 0, 575 0, 018 0, 326 R-Sq(adj) = 82, 4% Analysis of Variance Source DF SS MS F P Regression 12 410, 031 34, 169 18, 97 0, 000 Residual Error 34 61, 246 1, 801 Total 46 471, 277 11
Predictor Constant time x 1 x 2 x 3 x 4 x 5 x 6 x 7 x 8 x 9 x 10 x 11 Coef 3, 6491 0, 02851 -1, 691 -0, 469 2, 752 1, 224 6, 195 2, 417 8, 138 6, 360 0, 581 2, 553 1, 024 St. Dev 0, 8526 0, 01481 1, 028 1, 027 1, 026 1, 025 1, 026 1, 027 1, 028 T 4, 28 1, 92 -1, 65 -0, 46 2, 68 1, 19 6, 04 2, 36 7, 94 6, 20 0, 57 2, 49 1, 00 P 0, 000 0, 063 0, 109 0, 651 0, 011 0, 241 0, 000 0, 024 0, 000 0, 575 0, 018 0, 326 Tolkning av parametrar: – Antal sålda hus ökar i genomsnitt med 0, 0285 enheter per tidsenhet (månad) – I januari säljs det färre hus (-1. 69 hus) jämfört med december, i mars säljs det fler hus (+ 2. 75). . (observera att december är basperioden, eftersom dummyvariabeln för december inte finns med – decembernivån är 12 alltså inbakad i konstanten)
Inferens, som konfidensintervall, prognosintervall, t-test, F-test, partiellt F-test. . . kan göras på samma sätt som i vanlig regressionsanalys. Residualanalys bör göras för att kontrollera om villkoren för regression är uppfyllt: – Oberoende residualer – Normalfördelade residualer (för att kunna lita på testen) – Residualer med konstant varians (inga strutmönster) 13
Antagandet om oberoende residualer är ofta inte uppfyllt när det gäller tidsseriedata. Det kan också vara svårt att kolla detta antagande visuellt. 14
Enklare att se om observationerna är sammanbundna. Här ser man tydligt att en negativ residual vanligtvis följs av en positiv residual och tvärtom. Detta är ett tecken på autokorrelation. 15
Statistiskt test för att kontrollera om residualerna är oberoende: Durbin-Watson-test – Durbin-Watson-testet bedömer om autokorrelation (eller seriell korrelation) förekommer bland residualerna: Corr(et, et-1) – Vi skiljer mellan positiv autokorrelation och negativ autokorrelation 16
Negativ autokorrelation Positiv autokorrelation 17
Durbin-Watson-testvariabeln: I vårt exempel: Durbin-Watson statistic = 2. 66 18
Vi kan testa nollhypotesen: H 0: Det finns ingen autokorrelation i residualerna Om d > d. U, /2 eller (4 – d ) > d. U, /2 Ingen signifikant autokorrelation, H 0 kan ej förkastas Om d < d. L, /2 Signifikant positiv autokorrelation Om (4 – d ) < d. L, /2 Signifikant negativ autokorrelation Om d. L, /2 d d. U, /2 och d. L, /2 (4 – d ) d. U, /2 Inget uttalande kan göras 19
Om det inte finns någon autokorrelation i residualerna så kommer d att ligga nära 2. En approximativ kontroll kan göras genom att se om d är lägre än 1 eller högre än 3 Då finns autokorrelation i residualerna. 20
Vissa tidsserier har en så kallad exponentiell trend: Modell: Modellen kallas i boken ’growth curve model’ och jag har gått igenom de här modellerna i förra föreläsningen (avsnittet om 21 exponentiella modeller).
Klassisk komponentuppdelning: – Multiplikativ modell: – Additiv modell: där • • TRt=Trendkomponenten SNt=Säsongkomponenten CLt=Cykliska komponenten IRt=Slumpkomponenten 22
Trendkomponenten TR står för en (ofta) linjär funktion av tiden t Säsongkomponenten SN består av ett värde per säsong, som uttrycker skillnaden mellan denna säsong och årsgenomsnittet (jämför säsongsdummies) Cykliska komponenten CL står för en oregelbunden funktion som avspeglar konjunktursvängningar, alltså i form av en cykel. Slumpkomponenten är resten av variationen som är helt oregelbunden och som inte kan förklaras. 23
Multiplikativ eller additiv modell? Multiplikativ modell: Modellen används om säsongssvängningarna ökar med ökat nivå i serien. För ekonomiska data brukar denna modell ofta vara bäst. Additiv modell: Fungerar vid mer stabila tidsserier där säsongssvängningarna ej beror av nivån. 24
Skattning av komponenter, termer Säsongrensning: Borttagandet av säsongsvariation yt - SNt i den additiva modellen yt / SNt i den multiplikativa modellen Säsongsvariation överskuggar ofta andra relevanta komponenter. Genom säsongrensningen kan man alltså enklare se trender och andra komponenter. ’Detrending’: Borttagandet av trenden yt - TRt yt / TRt 25
Skattning av komponenter, steg-för-steg 1. Säsongrensning: Serien rensas från säsongkomponenten genom beräkning av centrerade och viktade glidande medelvärden (centered moving averages, CMA): där L=Antal säsonger i serien (L=2 för halvårsdata, 4 för kvartalsdata och 12 för månadsdata) 26
Exempel (sales data från tidigare) tid månad 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 1 14 2 15 3 antal 2 6 5 5 10 8 10 11 4 7 3 3 3 2 6 CMA * * * 6. 21 6. 08 5. 95. . 27
Trend och cyklisk komponent skattas grovt av CMAt. En första skattning av säsongkomponenterna erhålls genom att beräkna • yt/CMAt i en multiplikativ modell • yt – CMAt i en additiv modell och sen beräkna medelvärden för alla värden som avser samma säsong. (t. ex. alla januari-värden av yt/CMAt, etc. ) Totalt L medelvärden. 28
Medelvärdena måste dessutom justeras så att de vid multiplikativ modell får medelvärde 1, (dvs summan av alla justerade säsongmedelvärden ska bli L) vid additiv modell får medelvärde 0, (dvs summan av alla justerade säsongmedelvärden ska bli 0. ) De justerade värdena kallas för säsongskomponenter sn 1, . . . , sn. L 29
Exempel, forts 30
Medelvärden av grova säsongskomponenter: Juli: (1. 61074+2. 14013+1. 64571)/3 1. 7989 Aug: (1. 80822+1. 36709+1. 64571)/3 1. 6070 Sep: (0. 67133+0. 58896+0. 83237)/3 0. 6976 Okt: (1. 15862+1. 01818+0. 97110)/3 1. 0493 Nov: (0. 49655+1. 01205+0. 82286)/3 0. 7772 Dec: (0. 50350+0. 71006+0. 69767)/3 0. 6371 Jan: (0. 49315+0. 28571+0. 42353)/3 0. 4008 Feb: (0. 32432+0. 56805+0. 42105)/3 0. 4378 Mar: (0. 98630+1. 24138+1. 12941)/3 1. 1190 Apr: (0. 98630+0. 68182+0. 84706)/3 0. 8384 Maj: (1. 44000+1. 50857)/2 1. 4743 Juni: (1. 07692+1. 10345)/2 1. 0902 Obs! Bara två värden här! …och här! 31
Summan av de beräknade medelvärdena: 1. 7989 +1. 6070 + 0. 6976 + 1. 0493 + 0. 7772 + 0. 6371 + 0. 4008 + 0. 4378 + 1. 1190 + 0. 8384 + 1. 4743 + 1. 0902) 11. 9276 Summan skall bli L=12 För att få den till 12 multipliceras samtliga medelvärden med 12/11. 9276 1. 00607 32
Slutligt skattade säsongkomponenter: Jan: sn 1 = 0. 4008 · 1. 00607 0. 403 Feb: sn 2 = 0. 4378 · 1. 00607 0. 440 Mar: sn 3 = 1. 1190 · 1. 00607 1. 126 Apr: sn 4 = 0. 8384 · 1. 00607 0. 843 Maj: sn 5 = 1. 4743 · 1. 00607 1. 483 Juni: sn 6 = 1. 0902 · 1. 00607 1. 097 Juli: sn 7 = 1. 7989 · 1. 00607 1. 809 Aug: sn 8 = 1. 6070 · 1. 00607 1. 617 Sep: sn 9 = 0. 6976 · 1. 00607 0. 702 Okt: sn 10 = 1. 0493 · 1. 00607 1. 056 Nov: sn 11 = 0. 7772 · 1. 00607 0. 782 Dec: sn 12 = 0. 6371 · 1. 00607 0. 641 33
Tidsserien säsongrensas genom • vid multiplikativ modell • vid additiv modell där är något av värdena beroende på vilken av säsongerna som t motsvarar. 34
Exempel, forts 35
2. De säsongrensade värdena används för att skatta trendkomponenten. Skatta en linjär (eller kvadratisk) trend TRt med hjälp av regressionsanalys 36
3. Cyklisk och oregelbunden komponent: Om cyklisk komponent ej finns med: Residualerna från regressionsanalysen utgör skattning av termen IRt i den klassiska modellen. Om cyklisk komponent finns med: Skatta cyklisk och oregelbunden komponent som en komponent (CLIRt) 37
Den cykliska komponenten skattas nu genom ett centrerat oviktat glidande medelvärde: och den oregelbundna komponenten skattas slutligen som 38
Vilka glidande medelvärden ska användas? 2 m+1 väljs i regel till något av värdena 3, 5, 7, 9, 11, 13 Hur m skall väljas bestäms genom att titta på den slutliga skattningen av IRt m väljs så att autokorrelationen och variansen för dessa värden blir så låg som möjligt. 2 m+1 kallas antal punkter i det glidande medelvärdet 39
Minitab kan användas för komponentuppdelning med Stat Time series Decomposition Multiplikativ modell är dock något annorlunda: yt = TRt·SNt+IRt Val av modelltyp Möjlighet att välja komponenter, men dock begränsat 40
Säsongrensade data 41
Time Series Decomposition Data Sold Length 47, 0000 NMissing 0 Trend Line Equation Yt = 5, 77613 + 4, 30 E-02*t Seasonal Indices Period Index Dessa blir något annorlunda jämfört med handräkningen tidigare p g a att modellen är annorlunda 1 0, 425997 2 0, 425278 3 1, 14238 4 0, 856404 5 1, 52471 6 1, 10138 MAPE: 7 1, 65646 MAD: 0, 9057 8 1, 65053 MSD: 1, 6388 9 0, 670985 10 1, 02048 11 0, 825072 12 0, 700325 Accuracy of Model 16, 8643 42
43
44
45
Skattade trend- och säsongkomponenter har lagrats i kolumnerna TREN 1 resp. SEAS 1 Beräkning av kan göras genom att dividera originaldata med produkten av dessa två CLIR 1=Sold/(TREN 1· SEAS 1) Den cykliska komponenten skall nu skattas genom beräkning av glidande medelvärden på CLIR 1 46
Stat Time Series Moving Average… Antal punkter i det glidande medelvärdet 47
Sparar de glidande medelvärdena, dvs den skattade cykliska komponenten i en ny kolumn, som får namnet AVER 1 48
49
Den oregelbundna komponenten (IR) skattas slutligen genom att dividera CLIR 1 med AVER 1 De resulterade värdena studeras sedan med avseende på spridning, s och seriell korrelation, Corr ( irt , irt-1) 2 m+1 s Corr(irt, irt-1) 3 0. 219 -0. 686 5 0. 197 -0. 292 7 0. 173 -0. 343 9 0. 171 -0. 345 11 0. 181 -0. 276 13 0. 165 -0. 200 50
51
Seriella korrelationer kan enkelt beräknas med Stat Time series Lag och sedan Stat Basic statistics Correlation eller manuellt i Session window: MTB > lag ’RESI 4’ c 50 MTB > corr ’RESI 4’ c 50 52
Analys med additiv modell: 53
Time Series Decomposition Data Sold Length 47, 0000 NMissing 0 Trend Line Equation Yt = 5, 77613 + 4, 30 E-02*t Seasonal Indices Period Index 1 -4, 09028 2 -4, 13194 3 0, 909722 4 -1, 09028 5 3, 70139 6 0, 618056 7 4, 70139 MAPE: 8 4, 70139 MAD: 0, 9025 9 -1, 96528 MSD: 1, 6902 10 0, 118056 11 -1, 29861 12 -2, 17361 Accuracy of Model 16, 4122 54
Vad står måtten MAPE, MAD och MSD för? Alla tre är mått på hur bra anpassningen är och kan användas för att jämföra olika modeller. Den modell som har lägst MAPE, MAD och MSD har bäst anpassning. Oftast visar alla 3 måtten åt samma håll. Men i vissa fall kan man vara tvungen att välja en av dem. Vid val mellan t ex additiv modell och multiplikativ modell kan det hända att något av måtten är högre för den ena modellen medan ett annat mått är lägre. Det gäller alltså att tolka måtten med visst förnuft. 55
MSD kan också jämföras med MSE i den multipla regressionen: Mean Square Deviation Mean Square Error Formlerna är väldigt lika. Notera dock att vi dividerar med n och inte med n-k-1. Orsaken är att vi här inte har någon regressionsmodell med parametrar som måste skattas väntevärdesriktigt. Storleksmässigt kan dock MSD jämföras med MSE från tidsserieregressionen och är skillnaden markant kan vi också 56 se vilken av modellerna som ger bäst anpassning.
Mean Absolute Deviation Skillnaden mellan MAD och MSD är att MAD använder absolutavvikelser istället för kvadratiska avvikelser. MAD är mindre känslig för avvikande värden och blir mer användbar när vi har något enstaka värde som uppträder konstigt. Ytterligare en fördel med MAD är att dess värde är i samma skala som yt - observationerna själva, vilket gör det lättare att tolka. 57
Mean Absolute Percentage Error Måttet använder också absoluta avvikelser, men mäter dem relativt nivån hos y. Vi får alltså relativa (procentuella) avvikelser. Måttet är praktiskt för multiplikativa modeller där den oregelbundna komponenten (IRt ) är ganska betydande, eftersom avvikelserna då blir stora när vi har stora värden på y. 58
Multiplikativ Additiv 59
multiplikativ 60
Trend Line Equation Yt = 5. 77613 + 4. 30 E-02*t Seasonal Indices Period 1 2 3 4 5 6 7 8 9 10 11 12 Period Index 0. 425997 0. 425278 1. 14238 0. 856404 1. 52471 1. 10138 1. 65646 1. 65053 0. 670985 1. 02048 0. 825072 0. 700325 multiplikativ additiv 1 2 3 4 5 6 7 8 9 10 11 12 Index -4. 09028 -4. 13194 0. 909722 -1. 09028 3. 70139 0. 618056 4. 70139 -1. 96528 0. 118056 -1. 29861 -2. 17361 61
- Slides: 61