Lecture 4: Exponential Smoothing for Trended and Seasonal Time Series
|
|
- Godwin Day
- 6 years ago
- Views:
Transcription
1 NATCOR: Forecasing & Predicive Anayics Lecure 4: Exponenia Smoohing for Trended and Seasona Time Series John Boyan Lancaser Cenre for Forecasing Deparmen of Managemen Science
2 Srucure of his Lecure Par I Par II Par III Opimising parameers and iniia vaues Roing origin evauaion (ou-of-sampe Predicion Inervas Par IV Par V Par VI Trend Forecasing Seasona Forecasing Mehod Seecion Side 2 NATCOR Exp Smooh: Trend and Seasonaiy
3 Par I: Opimising parameers and iniia vaues Side 3 NATCOR Exp Smooh: Trend and Seasonaiy
4 Parameer Esimaion for SES Saes (A Forecas (F SKU A - SES α = Saes (A Forecas (F SKU A - SES α = Saes Saes Week Week MSE(in-sampe = 27.3 < MSE(in-sampe = 35.6 This parameer is beer, resuing in ower in-sampe MSE Based on his idea we can opimise parameers (for any exponenia smoohing mehod or generay. Side 4 NATCOR Exp Smooh: Trend and Seasonaiy
5 MSE as a funcion of apha We can cacuae he in-sampe MSE for various vaues of apha and idenify he vaue ha gives he owes error MSE Minimum In-sampe MSE α = Apha For more compex mehods, invoving rend and seasonaiy, we wi need o opimise severa parameers simuaneousy. Side 5 NATCOR Exp Smooh: Trend and Seasonaiy
6 MSE: funcion of Apha & Iniia Leve The same principe can be appied o choose he Iniia Leve (= Iniia Forecas for SES. Now we vary boh iniia forecas and smoohing parameer and joiny opimise boh Saes (A Forecas (F SKU B - SES 4000 Saes Week Side 6 Join opimisaion can be done in R Leve = 1697 Apha = NATCOR Exp Smooh: Trend and Seasonaiy
7 Is opimisaion aways he bes? No! MSE is very reacive o exreme errors, which may disor he error surface ha we are searching on for he opima vaues. Sampe imiaions as he number of parameers o opimise increases (especiay for seasona modes. Opimisaion is done in-sampe. The correaion beween insampe error and ou-of-sampe may be ow. Opimisaion is (ypicay done using 1 MSE. In pracice we forecas for onger horizons. MSE may no reae we o a company s financia objecives. Opimisaion is very usefu for auomaion, bu human expers may override idenified parameers if hey vioae heory or objecives. Side 7 NATCOR - Exp Smooh: Trend and Seasonaiy
8 Par II: Roing Origin Evauaion (Ou-of-sampe Side 8 NATCOR Exp Smooh: Trend and Seasonaiy
9 The hisorica observaions can be spi ino wo subses: In-sampe: In- and Ou-of-Sampe used for mode buiding and parameerisaion. Ou-of-sampe: used for mode evauaion. No used in buiding he mode and no seen by he mode. Used o simuae rue forecass, insead of waiing for new unobserved vaues, o evauae he forecasing performance of aernaive forecasing mehods In-sampe Ou-of-sampe 550 Unis Side Monh NATCOR Exp Smooh: Trend and Seasonaiy
10 Remarks on Singe Origin Evauaion Use o fi he mode In-sampe Use o evauae he mode Ou-of-sampe (Hodou Origin Forecass Limiaions: 1. Forecas suscepibe o srange origins. 2. One forecas per horizon 3. Averaging over differen horizons no recommended. Side 10 NATCOR Exp Smooh: Trend and Seasonaiy
11 Roing Origin Evauaion. In-sampe Ou-of-sampe Ro origin In-sampe (increased Origin Ou-of-sampe Ro origin In-sampe (increased Origin Ou-of-sampe Origin Coninue roing uni ou-of-sampe shorer han forecas horizon.
12 Visuaising Roing Origin Forecass Apha = 0.20 Apha = Side 12 In-sampe Ou-of-sampe In-sampe Ou-of-sampe Back dos are forecas origins Visuaising he roing origin forecass makes i easier o appreciae he imporance of smooh forecass ha fier noise NATCOR Exp Smooh: Trend and Seasonaiy
13 Remarks on Roing Origin Evauaion Hodou: 5 Hodou: 10 Forecas Fixed Origin Roing Origin Fixed Origin Roing Origin Horizon Number of forecass Number of forecass Number of forecass Number of forecass Provides more forecass per origin 2. Overcomes imiaions of fixed origin evauaion Provides more hisory per horizon for same ou-of-sampe Does no need o average over horizons Can overcome srange origins or arges Side 13 NATCOR Exp Smooh: Trend and Seasonaiy
14 Par III: Predicion Inervas Side 14 NATCOR Exp Smooh: Trend and Seasonaiy
15 Predicion Inervas SKU C SKU A 80% and 90% predicion inervas UK Android Marke Share US Air Passengers Observe ha he predicion inervas vary wih he series, mehod and horizon. Side 15 We have more confidence o forecass wih igh PIs.
16 Require Error Disribuion Daa SES We can use he error disribuion o formuae predicion inervas of he forecasing mehods. Probabiiy Error PDF e =A -F σ-2σ -σ µ σ 2σ 3σ 99.7% 95.5% 68.2% x 10-3 Side 16 NATCOR Exp Smooh: Trend and Seasonaiy
17 Which Errors? Three Approaches 1. Mode Errors (Irreguar Componens Requires ransaion of esimaes of he variance of mode errors (epsion erms, using a modeing framework, o esimaes of Predicion Inervas. This approach is adoped in es (R package and is based on anaysis in Hyndman e a (2008. Or Forecas Errors Two possibiiies o cacuae Predicion Inerva (PI : 2. In-sampe MSE 3. Ou-of-sampe MSE Side 17 NATCOR Exp Smooh: Trend and Seasonaiy
18 PI based on Ou-of-Sampe MSE We can use ou-of-sampe MSE: h = F h ± za / 2se = F h ± PI z h 2 a / MSE h This can be obained by cacuaing he roing origin ou-of-sampe MSE of he reevan forecas horizon: Advanages over oher approaches Genuiney ou-of-sampe forecas errors. Aows for parameers or mode ou-of-sampe being differen from parameers or mode in-sampe. Disadvanage over oher approaches Can be inaccurae if sma ou-of-sampe (es se. Side 18 NATCOR Exp Smooh: Trend and Seasonaiy
19 Par IV: Trend Forecasing Side 19 NATCOR Exp Smooh: Trend and Seasonaiy
20 Principes of Ho s Mehod Smooh he eve AND Smooh he rend Aows for differen raes of smoohing for eve and rend. Originay proposed in 1957 (Research Memorandum reprined in Inernaiona Journa of Forecasing in Impemened in mos commercia forecasing sofware packages. Forms he basis for more sophisicaed smoohing mehods. Side 20 NATCOR Exp Smooh: Trend and Seasonaiy
21 Ho s Mehod b y = αy h = = β ( (1 α( 1 hb 1 b 1 (1 β b 1 Firs equaion updaes he eve, using he aes observaion and he previous eve, adjused by he previous rend (b erm. Second equaion updaes he rend, using he aes difference in eves (oca rend and he previous smoohed rend. Third equaion provide a forecas h-seps-ahead by adding h imes he updaed rend o he updaed eve. Side 21 NATCOR Exp Smooh: Trend and Seasonaiy
22 Iniiaising Leve and Trend The iniia eve and rend are usuay opimised, aong wih he smoohing parameers. (As for SES in Par I of his ecure. However, some simpe empirica cacuaions incude: Iniia eve Firs acua observaion Simpe average of k firs observaions Iniia rend Difference of he firs wo observaions (Y 2 - Y 1 (Y 3 - Y 1 /2 ec Side 22 NATCOR Exp Smooh: Trend and Seasonaiy
23 Ho s Mehod: Linear Exrapoaion Does a sraigh ine forecas makes sense? WhDo any facors indicae ha he forecas shoud be differen o a sraigh ine? If as sraigh ine does no make sense, hen a damped rend mehod may be used insead Unis Daa Trend EXSM True Trend Monh Side 23 NATCOR Exp Smooh: Trend and Seasonaiy
24 Par V: Seasona Forecasing Side 24 NATCOR Exp Smooh: Trend and Seasonaiy
25 Need o Fier Noise Daa January January PDF Unis Probabiiy Jan96 Jan97 Jan98 Jan99 Jan00 Jan01 Jan02 Jan03 Jan04 Jan05 Jan06 Jan07 Jan08 Jan09 Jan10 Monh 100 Observed True Simiar o eve or rend, he seasona componen is affeced by noise, which needs o be fiered The deviaion around he rue January is due o noise Side 25 NATCOR Exp Smooh: Trend and Seasonaiy
26 Winers Mehod (Addiive Seasonaiy Side 26 NATCOR Exp Smooh: Trend and Seasonaiy Firs equaion updaes he eve, using he difference beween he aes observaion and previous seasona effec, adjused by he previous eve. Second equaion updaes he season, using he difference in aes observaions and eves, and he previous seasona. Third equaion provide a forecas h-seps-ahead by adding mos recen seasona o he aes eve. h m h m m s y s y s s y = = = 1 (1 ( (1 ( γ γ α α
27 Winers Mehod (Muipicaive Seasonaiy Side 27 NATCOR Exp Smooh: Trend and Seasonaiy These equaions work on he same basis as for addiive seasonaiy. Again, he seasona componen is smoohed hrough he gamma parameer (smoohing consan. h m h m m s y s y s s y = = = 1 (1 (1 γ γ α α
28 Principes of Ho-Winers Mehod Smooh he eve AND Smooh he rend AND Smooh he seasonaiy Aows for differen raes of smoohing for eve (apha, rend (bea and seasonaiy (gamma. Impemened in mos commercia forecasing sofware packages. Side 28 NATCOR Exp Smooh: Trend and Seasonaiy
29 Ho- Winers Mehod (Addiive Seasonaiy Side 29 NATCOR Exp Smooh: Trend and Seasonaiy Combines he Ho s and Winers (addiive seasonaiy mehods. Aows for hree differen eves of smoohing. Common o observe boh rend and seasonaiy ogeher h m h m m s hb y s y s b b b s y = = = = (1 ( (1 ( ( (1 ( γ γ β β α α
30 Ho- Winers Mehod (Muipicaive Seasonaiy Side 30 NATCOR Exp Smooh: Trend and Seasonaiy Combines he Ho s and Winers (muipicaive seasonaiy mehods. Again, aows for hree differen eves of smoohing. h m h m m s hb y s y s b b b s y = = = = ( (1 (1 ( ( ( γ γ β β α α
31 Forecasing Compeiions Two famous Forecasing Compeiions: M1 Compeiion ( rea series. M3 Compeiion ( rea series. Up o as 18 periods wihhed from compeiors o ensure genuine ou of sampe resus. Concusions Simper mehods (eg exponenia smoohing ofen compeiive wih more sophisicaed mehods. Performance depends on horizon and choice of accuracy meric. Averages of differen mehods can perform we. Side 31 NATCOR Exp Smooh: Trend and Seasonaiy
32 Par VI: Mehod Seecion for Exponenia Smoohing Side 32 NATCOR Exp Smooh: Trend and Seasonaiy
33 Famiy of Exponenia Smoohing Modes There are 30 Exponenia Smoohing Modes 15 wih addiive errors 15 wih muipicaive errors Each se of 15 modes comprise combinaions of: 5 ypes of rend 3 ypes of seasonaiy Forecas mehods do no differ beween addiive and muipicaive error modes (provided he same parameer vaues are used. Formuae for Predicion Inervas do differ. Side 33 NATCOR Exp Smooh: Trend and Seasonaiy
34 Famiy of Exponenia Smoohing Mehods
35 Seecing Exponenia Smoohing Mehod Temping o simpy compare in-sampe MSEs for differen Exponenia Smoohing mehods. This fais o ake ino accoun he number of parameers being esimaed. There is a risk of over-fiing. Side 35 NATCOR Exp Smooh: Trend and Seasonaiy
36 Over-Fiing Suppose we fi on a saionary ime series a Trend-Seasona Exponenia Smoohing mode: Trend seasona inerpreed noise as rend and seasonaiy! Trend seasona exaggeraed noise in-sampe! Daa Trend-Seasona SES NATCOR Exp Smooh: Trend and Seasonaiy
37 Comparison of Mehods using Informaion Crieria Le p = number of parameers (incuding Iniia Vaues used in he mode. Informaion Crieria penaize according o size of p Le n = sampe size afer aowing for he observaions used o esabish saring vaues. Akaike s Informaion Crierion Sma-sampe correced Bayesian Informaion Crierion AIC = n( MSE 2 p / n 2 p AICc = n( MSE n p 1 BIC = n( MSE p n( n / n The mode wih he owes Informaion Crierion is seeced. es has AIC, AICc and BIC; AICc is he defau. Side 37 NATCOR Exp Smooh: Trend and Seasonaiy
38 Summary Can opimize over whoe se of parameers (incuding iniia vaues. Roing Origins make more efficien use of he ou-ofsampe (es or hod-ou se. Predicion Inervas can be cacuaed using ou-of-sampe MSE for fuure forecass. Exponenia smoohing can be exended nauray for rend and/or seasona daa. Exponenia Smoohing mehods can be compared using Informaion Crieria which penaise he number of parameers. Side 38 NATCOR Exp Smooh: Trend and Seasonaiy
39 References Exponenia Smoohing Mehods Ho CC (2004 Forecasing seasonas and rends by exponeniay weighed moving averages. In J Forecasing, 20, Goodwin P (2010 The Ho-Winers approach o exponenia smoohing: 50 years od and going srong. Foresigh, Fa 2010, Forecasing Compeiions Makridakis S e a (1982 The accuracy of exrapoaion (ime series mehods: resus of a forecasing compeiion. J Forecasing, 1, Makridakis S and Hibon M (2000 The M3 Compeiion: resus, concusions and impicaions, In J Forecasing, 16, Side 39 NATCOR Exp Smooh: Trend and Seasonaiy
Lecture 3: Exponential Smoothing
NATCOR: Forecasing & Predicive Analyics Lecure 3: Exponenial Smoohing John Boylan Lancaser Cenre for Forecasing Deparmen of Managemen Science Mehods and Models Forecasing Mehod A (numerical) procedure
More informationLinear Gaussian State Space Models
Linear Gaussian Sae Space Models Srucural Time Series Models Level and Trend Models Basic Srucural Model (BSM Dynamic Linear Models Sae Space Model Represenaion Level, Trend, and Seasonal Models Time Varying
More informationTypes of Exponential Smoothing Methods. Simple Exponential Smoothing. Simple Exponential Smoothing
M Business Forecasing Mehods Exponenial moohing Mehods ecurer : Dr Iris Yeung Room No : P79 Tel No : 788 8 Types of Exponenial moohing Mehods imple Exponenial moohing Double Exponenial moohing Brown s
More informationExponential Smoothing
Exponenial moohing Inroducion A simple mehod for forecasing. Does no require long series. Enables o decompose he series ino a rend and seasonal effecs. Paricularly useful mehod when here is a need o forecas
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
More informationSmoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T
Smoohing Consan process Separae signal & noise Smooh he daa: Backward smooher: A an give, replace he observaion b a combinaion of observaions a & before Simple smooher : replace he curren observaion wih
More informationTime series Decomposition method
Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,
More information3.1 More on model selection
3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of
More informationChapter 15. Time Series: Descriptive Analyses, Models, and Forecasting
Chaper 15 Time Series: Descripive Analyses, Models, and Forecasing Descripive Analysis: Index Numbers Index Number a number ha measures he change in a variable over ime relaive o he value of he variable
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationSummer Term Albert-Ludwigs-Universität Freiburg Empirische Forschung und Okonometrie. Time Series Analysis
Summer Term 2009 Alber-Ludwigs-Universiä Freiburg Empirische Forschung und Okonomerie Time Series Analysis Classical Time Series Models Time Series Analysis Dr. Sevap Kesel 2 Componens Hourly earnings:
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationReliability of Technical Systems
eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy
More informationPROC NLP Approach for Optimal Exponential Smoothing Srihari Jaganathan, Cognizant Technology Solutions, Newbury Park, CA.
PROC NLP Approach for Opimal Exponenial Smoohing Srihari Jaganahan, Cognizan Technology Soluions, Newbury Park, CA. ABSTRACT Esimaion of smoohing parameers and iniial values are some of he basic requiremens
More informationFinancial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS
Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If
More informationFrequency independent automatic input variable selection for neural networks for forecasting
Universiä Hamburg Insiu für Wirschafsinformaik Prof. Dr. D.B. Preßmar Frequency independen auomaic inpu variable selecion for neural neworks for forecasing Nikolaos Kourenzes Sven F. Crone LUMS Deparmen
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationWisconsin Unemployment Rate Forecast Revisited
Wisconsin Unemploymen Rae Forecas Revisied Forecas in Lecure Wisconsin unemploymen November 06 was 4.% Forecass Poin Forecas 50% Inerval 80% Inerval Forecas Forecas December 06 4.0% (4.0%, 4.0%) (3.95%,
More informationUncertainty in predictive modelling
Uncerainy in predicive modelling Nikolaos Kourenzes Lancaser Universiy OR59 14/09/2017 Agenda 1. Forecasing & uncerainy 2. Models, assumpions & uncerainy 3. Managing modelling uncerainy 4. Approach I:
More informationA Two Stage Mixture Model for Predicting EAD
A wo Sage Mixure Mode for Predicing EAD Credi Scoring & Credi Conro, Edinburgh 2013 Mindy Leow & Jonahan Crook Credi Research Cenre, Universiy of Edinburgh, UK 1 Moivaion Reguaory Capia ase II requires
More informationMethodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.
Mehodology. Uni Roo Tess A ime series is inegraed when i has a mean revering propery and a finie variance. I is only emporarily ou of equilibrium and is called saionary in I(0). However a ime series ha
More informationSimultaneous Localisation and Mapping. IAR Lecture 10 Barbara Webb
Simuaneous Locaisaion and Mapping IAR Lecure 0 Barbara Webb Wha is SLAM? Sar in an unknown ocaion and unknown environmen and incremenay buid a map of he environmen whie simuaneousy using his map o compue
More informationEstimation Uncertainty
Esimaion Uncerainy The sample mean is an esimae of β = E(y +h ) The esimaion error is = + = T h y T b ( ) = = + = + = = = T T h T h e T y T y T b β β β Esimaion Variance Under classical condiions, where
More informationIntroduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.
Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since
More informationAffine term structure models
Affine erm srucure models A. Inro o Gaussian affine erm srucure models B. Esimaion by minimum chi square (Hamilon and Wu) C. Esimaion by OLS (Adrian, Moench, and Crump) D. Dynamic Nelson-Siegel model (Chrisensen,
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationRecent Developments In Evolutionary Data Assimilation And Model Uncertainty Estimation For Hydrologic Forecasting Hamid Moradkhani
Feb 6-8, 208 Recen Developmens In Evoluionary Daa Assimilaion And Model Uncerainy Esimaion For Hydrologic Forecasing Hamid Moradkhani Cener for Complex Hydrosysems Research Deparmen of Civil, Consrucion
More informationForecasting optimally
I) ile: Forecas Evaluaion II) Conens: Evaluaing forecass, properies of opimal forecass, esing properies of opimal forecass, saisical comparison of forecas accuracy III) Documenaion: - Diebold, Francis
More informationDistribution of Least Squares
Disribuion of Leas Squares In classic regression, if he errors are iid normal, and independen of he regressors, hen he leas squares esimaes have an exac normal disribuion, no jus asympoic his is no rue
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More informationIntermittent Demand Forecast and Inventory Reduction Using Bayesian ARIMA Approach
Proceedings of he 00 Inernaional Conference on Indusrial Engineering and Operaions Managemen Dhaka, Bangladesh, January 9 0, 00 Inermien Demand orecas and Invenory Reducion Using Bayesian ARIMA Approach
More informationDistribution of Estimates
Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationLicenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A
Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationNature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.
Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike
More informationACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.
ACE 56 Fall 5 Lecure 8: The Simple Linear Regression Model: R, Reporing he Resuls and Predicion by Professor Sco H. Irwin Required Readings: Griffihs, Hill and Judge. "Explaining Variaion in he Dependen
More informationWEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x
WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile
More informationBox-Jenkins Modelling of Nigerian Stock Prices Data
Greener Journal of Science Engineering and Technological Research ISSN: 76-7835 Vol. (), pp. 03-038, Sepember 0. Research Aricle Box-Jenkins Modelling of Nigerian Sock Prices Daa Ee Harrison Euk*, Barholomew
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationDynamic models for largedimensional. Yields on U.S. Treasury securities (3 months to 10 years) y t
Dynamic models for largedimensional vecor sysems A. Principal componens analysis Suppose we have a large number of variables observed a dae Goal: can we summarize mos of he feaures of he daa using jus
More informationRL Lecture 7: Eligibility Traces. R. S. Sutton and A. G. Barto: Reinforcement Learning: An Introduction 1
RL Lecure 7: Eligibiliy Traces R. S. Suon and A. G. Baro: Reinforcemen Learning: An Inroducion 1 N-sep TD Predicion Idea: Look farher ino he fuure when you do TD backup (1, 2, 3,, n seps) R. S. Suon and
More informationUsing the Kalman filter Extended Kalman filter
Using he Kalman filer Eended Kalman filer Doz. G. Bleser Prof. Sricker Compuer Vision: Objec and People Tracking SA- Ouline Recap: Kalman filer algorihm Using Kalman filers Eended Kalman filer algorihm
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationForward guidance. Fed funds target during /15/2017
Forward guidance Fed funds arge during 2004 A. A wo-dimensional characerizaion of moneary shocks (Gürkynak, Sack, and Swanson, 2005) B. Odyssean versus Delphic foreign guidance (Campbell e al., 2012) C.
More informationQuarterly ice cream sales are high each summer, and the series tends to repeat itself each year, so that the seasonal period is 4.
Seasonal models Many business and economic ime series conain a seasonal componen ha repeas iself afer a regular period of ime. The smalles ime period for his repeiion is called he seasonal period, and
More informationUnit Root Time Series. Univariate random walk
Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More informationCH Sean Han QF, NTHU, Taiwan BFS2010. (Joint work with T.-Y. Chen and W.-H. Liu)
CH Sean Han QF, NTHU, Taiwan BFS2010 (Join work wih T.-Y. Chen and W.-H. Liu) Risk Managemen in Pracice: Value a Risk (VaR) / Condiional Value a Risk (CVaR) Volailiy Esimaion: Correced Fourier Transform
More informationExcel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand
Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338
More informationResearch Article Developing a Series Solution Method of q-difference Equations
Appied Mahemaics Voume 2013, Arice ID 743973, 4 pages hp://dx.doi.org/10.1155/2013/743973 Research Arice Deveoping a Series Souion Mehod of q-difference Equaions Hsuan-Ku Liu Deparmen of Mahemaics and
More informationAir Traffic Forecast Empirical Research Based on the MCMC Method
Compuer and Informaion Science; Vol. 5, No. 5; 0 ISSN 93-8989 E-ISSN 93-8997 Published by Canadian Cener of Science and Educaion Air Traffic Forecas Empirical Research Based on he MCMC Mehod Jian-bo Wang,
More informationEnsamble methods: Boosting
Lecure 21 Ensamble mehods: Boosing Milos Hauskrech milos@cs.pi.edu 5329 Senno Square Schedule Final exam: April 18: 1:00-2:15pm, in-class Term projecs April 23 & April 25: a 1:00-2:30pm in CS seminar room
More informationEnsamble methods: Bagging and Boosting
Lecure 21 Ensamble mehods: Bagging and Boosing Milos Hauskrech milos@cs.pi.edu 5329 Senno Square Ensemble mehods Mixure of expers Muliple base models (classifiers, regressors), each covers a differen par
More informationLecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance
Lecure 5 Time series: ECM Bernardina Algieri Deparmen Economics, Saisics and Finance Conens Time Series Modelling Coinegraion Error Correcion Model Two Seps, Engle-Granger procedure Error Correcion Model
More information20. Applications of the Genetic-Drift Model
0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0
More informationGDP Advance Estimate, 2016Q4
GDP Advance Esimae, 26Q4 Friday, Jan 27 Real gross domesic produc (GDP) increased a an annual rae of.9 percen in he fourh quarer of 26. The deceleraion in real GDP in he fourh quarer refleced a downurn
More informationProblem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims
Problem Se 5 Graduae Macro II, Spring 2017 The Universiy of Nore Dame Professor Sims Insrucions: You may consul wih oher members of he class, bu please make sure o urn in your own work. Where applicable,
More informationAugmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004
Augmened Realiy II Kalman Filers Gudrun Klinker May 25, 2004 Ouline Moivaion Discree Kalman Filer Modeled Process Compuing Model Parameers Algorihm Exended Kalman Filer Kalman Filer for Sensor Fusion Lieraure
More informationStationary Time Series
3-Jul-3 Time Series Analysis Assoc. Prof. Dr. Sevap Kesel July 03 Saionary Time Series Sricly saionary process: If he oin dis. of is he same as he oin dis. of ( X,... X n) ( X h,... X nh) Weakly Saionary
More informationModeling the Volatility of Shanghai Composite Index
Modeling he Volailiy of Shanghai Composie Index wih GARCH Family Models Auhor: Yuchen Du Supervisor: Changli He Essay in Saisics, Advanced Level Dalarna Universiy Sweden Modeling he volailiy of Shanghai
More informationTHE UNIVERSITY OF TEXAS AT AUSTIN McCombs School of Business
THE UNIVERITY OF TEXA AT AUTIN McCombs chool of Business TA 7.5 Tom hively CLAICAL EAONAL DECOMPOITION - MULTIPLICATIVE MODEL Examples of easonaliy 8000 Quarerly sales for Wal-Mar for quarers a l e s 6000
More informationSequential Importance Resampling (SIR) Particle Filter
Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle
More information12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =
1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of
More informationuser! 2009 July 8-10
user! 009 July 8-0 . A Parnership beween firm and universiy. The company STEF-TFE 3. Frozen goods supply chain 4. Forecasing goals 5. Forecasing is no realiy 6. Forecasing model 7. FSS Horizons Universiy
More informationProbabilistic Robotics
Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae
More informationIntroduction to Simulation - Lecture 16. Methods for Computing Periodic Steady-State - Part II Jacob White
Inroducion o Simuaion - ecure 6 Mehods for Compuing Periodic Seady-Sae - Par II Jacob Whie hanks o Deepak Ramaswamy, Micha Rewienski, and Karen Veroy Ouine hree Mehods so far ime inegraion uni seady-sae
More informationLinear Combinations of Volatility Forecasts for the WIG20 and Polish Exchange Rates
Eliza Buszkowska Universiy of Poznań, Poland Linear Combinaions of Volailiy Forecass for he WIG0 and Polish Exchange Raes Absrak. As is known forecas combinaions may be beer forecass hen forecass obained
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationIn this paper the innovations state space models (ETS) are used in series with:
Time series models for differen seasonal paerns Blaconá, M.T, Andreozzi*, L. and Magnano, L. Naional Universiy of Rosario, (*)CONICET - Argenina Absrac In his paper Innovaions Sae Space Models (ETS) are
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationExplaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015
Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become
More informationVectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1
Vecorauoregressive Model and Coinegraion Analysis Par V Time Series Analysis Dr. Sevap Kesel 1 Vecorauoregression Vecor auoregression (VAR) is an economeric model used o capure he evoluion and he inerdependencies
More information1 Evaluating Chromatograms
3 1 Evaluaing Chromaograms Hans-Joachim Kuss and Daniel Sauffer Chromaography is, in principle, a diluion process. In HPLC analysis, on dissolving he subsances o be analyzed in an eluen and hen injecing
More informationChapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws
Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species
More informationMonitoring and data filtering II. Dynamic Linear Models
Ouline Monioring and daa filering II. Dynamic Linear Models (Wes and Harrison, chaper 2 Updaing equaions: Kalman Filer Discoun facor as an aid o choose W Incorporae exernal informaion: Inervenion General
More informationComparing Means: t-tests for One Sample & Two Related Samples
Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion
More informationCointegration and Implications for Forecasting
Coinegraion and Implicaions for Forecasing Two examples (A) Y Y 1 1 1 2 (B) Y 0.3 0.9 1 1 2 Example B: Coinegraion Y and coinegraed wih coinegraing vecor [1, 0.9] because Y 0.9 0.3 is a saionary process
More informationA First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18
A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly
More informationAn recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes
WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,
More informationForecasting. Summary. Sample StatFolio: tsforecast.sgp. STATGRAPHICS Centurion Rev. 9/16/2013
STATGRAPHICS Cenurion Rev. 9/16/2013 Forecasing Summary... 1 Daa Inpu... 3 Analysis Opions... 5 Forecasing Models... 9 Analysis Summary... 21 Time Sequence Plo... 23 Forecas Table... 24 Forecas Plo...
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationTwo Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017
Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =
More information13.3 Term structure models
13.3 Term srucure models 13.3.1 Expecaions hypohesis model - Simples "model" a) shor rae b) expecaions o ge oher prices Resul: y () = 1 h +1 δ = φ( δ)+ε +1 f () = E (y +1) (1) =δ + φ( δ) f (3) = E (y +)
More informationForecast of Adult Literacy in Sudan
Journal for Sudies in Managemen and Planning Available a hp://inernaionaljournalofresearch.org/index.php/jsmap e-issn: 2395-463 Volume 1 Issue 2 March 215 Forecas of Adul Lieracy in Sudan Dr. Elfarazdag
More informationRepeat Sales Regression on Heterogeneous Properties. Liang Peng
Repea Saes Regression on Heerogeneous Properies Liang Peng Absrac This paper proposes a generaized repea saes regression (GRSR) ha uses repea saes from he enire marke, in which properies may have heerogeneous
More informationSTA 114: Statistics. Notes 2. Statistical Models and the Likelihood Function
STA 114: Saisics Noes 2. Saisical Models and he Likelihood Funcion Describing Daa & Saisical Models A physicis has a heory ha makes a precise predicion of wha s o be observed in daa. If he daa doesn mach
More informationPET467E-Analysis of Well Pressure Tests/2008 Spring Semester/İTÜ Midterm Examination (Duration 3:00 hours) Solutions
M. Onur 03.04.008 PET467E-Analysis of Well Pressure Tess/008 Spring Semeser/İTÜ Miderm Examinaion (Duraion 3:00 hours) Soluions Name of he Suden: Insrucions: Before saring he exam, wrie your name clearly
More informationLinear Dynamic Models
Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference
More informationLecture Notes 2. The Hilbert Space Approach to Time Series
Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More informationIB Physics Kinematics Worksheet
IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?
More informationReady for euro? Empirical study of the actual monetary policy independence in Poland VECM modelling
Macroeconomerics Handou 2 Ready for euro? Empirical sudy of he acual moneary policy independence in Poland VECM modelling 1. Inroducion This classes are based on: Łukasz Goczek & Dagmara Mycielska, 2013.
More informationVector autoregression VAR. Case 1
Vecor auoregression VAR So far we have focused mosl on models where deends onl on as. More generall we migh wan o consider oin models ha involve more han one variable. There are wo reasons: Firs, we migh
More informationUnobserved Component Model with Observed Cycle Use of BTS Data for Short-Term Forecasting of Industrial Production
Sławomir Dudek Dawid Pachucki Research Insiue for Economic Developmen (RIED) Warsaw School of Economics (WSE) Unobserved Componen Model wih Observed Cycle Use of BTS Daa for Shor-Term Forecasing of Indusrial
More informationSpeaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis
Speaker Adapaion Techniques For Coninuous Speech Using Medium and Small Adapaion Daa Ses Consaninos Boulis Ouline of he Presenaion Inroducion o he speaker adapaion problem Maximum Likelihood Sochasic Transformaions
More informationLecture 33: November 29
36-705: Inermediae Saisics Fall 2017 Lecurer: Siva Balakrishnan Lecure 33: November 29 Today we will coninue discussing he boosrap, and hen ry o undersand why i works in a simple case. In he las lecure
More information