Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
|
|
- Jonas McDowell
- 5 years ago
- Views:
Transcription
1 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number of pages = 21 Number of supplemenal figures = 10 Number of supplemenal ables = 2 S1
2 16 SUPPLEMENTAL MATERIALS, TABLES AND FIGURES Kriging Deails Consider he following model which decomposes an observaion (on he log scale) ino addiive signal and noise componens: Z S ; T, (1) where Z represens an observaion a ime wihin some coninuous ime inerval T, and S and represen he underlying signal and noise componens, respecively. In pracice, T is ypically discreized ino a uniform grid of, say, daily unis. The signal and noise componens are assumed o be uncorrelaed, and he noise componen is 2 assumed o be a zero-mean uncorrelaed ime-series process wih consan variance =. Suppose ha he signal componen can be furher decomposed ino a fixed regression componen and a random componen represening a zero-mean, second-order saionary ime-series process possessing a emporal correlaion srucure. Specifically, consider he following model: S x β ; T, (2) where x represens a (fixed) vecor of covariae values a ime and β represens he vecor of associaed regression coefficiens; and represens he zero-mean residual a ime. Assume ha he residual ime-series process has he following second-order saionary variance- covariance srucure: cov(, u) C( u ); for all, u T, (3) 33 where Ch ( ) represens he covariance beween residuals a disance of h apar; in paricular, 34 le C 2 (0). Second-order saionariy assumes a consan mean and variance, and a emporal S2
3 correlaion srucure ha depends only on he lag disance beween wo observaions, and his can be checked graphically. Therefore, we can express he observed process as: 38 where Z x β ; T, (4) and are assumed o be uncorrelaed, and r represens he observed residual a ime, and hence has he following second-order saionary variance-covariance srucure: 2 2 ; u cov( r, ru ) CZ ( u ), for all, u T. C( u ); u (5) Noe ha cov(, ) and cov( r, r ) are idenical excep in he case where = u (i.e., where he u u sampled and predicion locaions coincide). Suppose ha a sample of ime poins has been drawn, represened by A { : i 1,..., n} T {1,..., N}, where T represens he complee se of days in he ime i inerval of ineres (i.e., he N = 122 days of he usage season). Following Cressie (1993; Secion 3.2.1), we ypically wan o predic (some funcion of) he signal process S ( S1,..., SN ) over all days in T (i.e., we wan o filer ou he noise componen), using informaion abou 48 Z ( Z,..., Z ) obained from he sampled daa and informaion from he covariae(s) from he 1 n sampled and non-sampled daa. In he absence of noise (i.e, if 0), predicions a he sampled poins will coincide wih he Z-values a hose poins. The (heoreical) universal kriging predicor of he signal a ime 0 can be expressed as: ˆ ˆ 1 S ˆ x β c Σ Z Xβ ( ), (6) S3
4 where βˆ ( X Σ X) X Σ Z is he generalized leas squares (GLS) esimaor of β, 53 X ( x,..., x ), Σ var( Z), and 1 n c ( C( ),..., C( )). From linear-model heory n (Goldberger, 1962), we have ha ˆβ is he bes linear unbiased esimaor (BLUE) of β, and ha S is he bes linear unbiased predicor (BLUP) of ˆ0 S 0, where "bes" means ha he mean squared (predicion) error is minimized. The mean squared predicion error (MSPE) of S ˆ can 0 be expressed as: ˆ ˆ MSPE( S ) E( S S ) c Σ c ( x c Σ X)( X Σ X) ( x X Σ c ) (7) 58 hence he (1 α) 100% predicion inerval of S ˆ can be expressed as: 0 (S 0 z α 2 MSPE(S 0 ), S 0 + z α 2 MSPE(S 0 )) (8) where z α 2 is he upper α 2 criical value of he sandard normal disribuion. Typically, he elemens of Σ and c 0 are unknown, so we need o inser esimaed values in heir place; he esimaed values are obained from he following procedures involving semivariograms. The (heoreical) semivariogram of he observable (noisy) r residuals is defined as: 0, u 1 r (, u) 2 var( r ru ), for all, u T, 2 2 C( u ), u (9) where is ofen referred o as he "parial sill," and is ofen referred o as he "sill." 2 Also, (0 ) lim ( u, ) is ofen referred o as he "nugge," and he "range" ofen refers r u 0 r o he smalles lag disance for which C( u ) 0. A special case of universal kriging (called ordinary kriging) applies when x β (i.e., he mean funcion is consan). S4
5 In pracice, ( u, ) is firs esimaed a he u lags in he sampled daa. For example, a 7-day r sampling frequency would allow us o esimae he semivariogram only a lags 7, 14, 21, ec. The nex sep is o fi a semivariogram model o hose esimaes (see Cressie, 1993, for furher 71 deails). The semivariogram model will provide esimaes of 2, 2, and C( u ), which can hen be used as esimaes of he elemens of Σ and c. 0 Because of he limied sample and lack of daa a shor lags, only he piecewise linear semivariogram model (wih range = 14) was used for 14-day sampling. For 7-day sampling, we compared wo approaches: he firs was o successively fi exponenial, hen spherical models unil a sable fi was achieved, wih each semivariogram model fied o semivariogram esimaes successively derived from hree mehods described in Cressie (1993; Secion 2.4) (e.g., mehod-of-momens, hen robus, hen median approach). If no fi was sable hen a piecewise linear semivariogram model (wih range = 7) was used. The second was o fi only he linear semivariogram model. Also in pracice, universal kriging predicion ypically proceeds according o an ieraive process involving he following seps: 1) a regression model is fied o he daa by means of ordinary leas squares (OLS); 2) he elemens of Σ and c 0 are hen esimaed from he residuals of he resuling model; 3) he regression model is refied using he updaed elemens of Σ and he elemens of Σ and c 0 : 4) c 0 are hen re-esimaed from he residuals of he updaed model; and 5) his ieraive process is coninued unil he esimaes converge (e.g., Cressie, 1993; Secion 3.4.3). However, for his paper, he ieraive process was runcaed a he second sep, and hence he MSPE of S ˆ described above can be viewed as an approximaion. The jusificaions for 0 runcaing he ieraive process include logisical reasons (i.e., we needed a general simplified predicion approach ha could be applied across many sie years), and ha since our goal is S5
6 predicion (wih mehod performance evaluaed by comparison agains near-complee monioring daases), quesions abou he approximaions assumed in deriving he models are of less ineres. Ordinary kriging is hen applied o predic he residual ime series, and finally hose prediced values are added o he regression model o form he universal kriging predicions. The reason for OLS esimaion in he firs sep insead of GLS esimaion as described in he heoreical universal kriging predicor (6), is ha he GLS esimaor requires simulaneous knowledge of Σ and c 0, bu hese quaniies can only be esimaed in he second sep of he ieraive process. OLS esimaion is sill unbiased and is only slighly less precise han GLS esimaion, so he pracical applicaion in he wo seps described is ypically followed. Oher pracical consideraions included he use of rimming prediced values. For example, predicions ha fell below he LOD were rimmed o he LOD value, and ordinary kriging predicions of he residuals a non-sampled locaions were rimmed so ha hey could no exceed wo sandard deviaions from zero. Finally, exponeniaion of he kriging predicions and predicion limis was applied o back- ransform hese values o he original arazine scale. Exponeniaion may resul in a small amoun of bias, bu no correcion facor was applied (Aldworh and Cressie, 2003). SUPPLEMENTAL FIGURE CAPTIONS Supplemenal Figure S1. Sudy monioring sies, as idenified by he Arazine Ecological Monioring Program (AEMP) sie ID code Supplemenal Figure S2. Semivariogram models fied o empirical semivariograms for wo example sie-years wih daily daa. S6
7 Supplemenal Figure S3. RRMSPE (i.e., RMSPE/observed maximum rolling average) versus observed maximum rolling average for ordinary kriging (solid circle) and linear inerpolaion (open circle) for 7-day sampling and maximum 7-day rolling average (a), and 14-day sampling and maximum 7-day rolling average (b). Connecing lines are blue if kriging performs beer and red if i performs worse han linear inerpolaion Supplemenal Figure S4. RRMSPE (i.e., RMSPE/observed maximum rolling average) versus observed maximum rolling average for ordinary kriging (solid circle) and linear inerpolaion (open circle) for 7-day sampling for yearly maximum, maximum 7-day rolling average, maximum 14-day rolling average, and maximum 30-day rolling average. Connecing lines are blue if kriging performs beer and red if i performs worse han linear inerpolaion Supplemenal Figure S5. RRMSPE (i.e., RMSPE/observed maximum rolling average) versus observed maximum rolling average for ordinary kriging (solid circle) and linear inerpolaion (open circle) for 14-day sampling for yearly maximum, maximum 7-day rolling average, maximum 14-day rolling average, and maximum 30-day rolling average. Connecing lines are blue if kriging performs beer and red if i performs worse han linear inerpolaion Supplemenal Figure S6. RRMSPE (i.e., RMSPE/observed maximum rolling average) versus observed maximum rolling average for universal kriging wih PRZM-hybrid covariae (solid circle) and linear inerpolaion (open circle) for 7-day sampling for yearly maximum, maximum 7-day rolling average, maximum 14-day rolling average, and maximum 30-day rolling average. Connecing lines are blue if kriging performs beer and red if i performs worse han linear inerpolaion. S7
8 Supplemenal Figure S7. RRMSPE (i.e., RMSPE/observed maximum rolling average) versus observed maximum rolling average for universal kriging wih PRZM-hybrid covariae (solid circle) and linear inerpolaion (open circle) for 14-day sampling for yearly maximum, maximum 7-day rolling average, maximum 14-day rolling average, and maximum 30-day rolling average. Connecing lines are blue if kriging performs beer and red if i performs worse han linear inerpolaion Supplemenal Figure S8. RRMSPE (i.e., RMSPE/observed maximum rolling average) versus observed maximum rolling average for log-linear inerpolaion (solid circle) and linear inerpolaion (open circle) for 7-day sampling for yearly maximum, maximum 7-day rolling average, maximum 14-day rolling average, and maximum 30-day rolling average. Connecing lines are blue if log-linear inerpolaion performs beer and red if i performs worse han linear inerpolaion Supplemenal Figure S9. RRMSPE (i.e., RMSPE/observed maximum rolling average) versus observed maximum rolling average for log-linear inerpolaion (solid circle) and linear inerpolaion (open circle) for 14-day sampling for yearly maximum, maximum 7-day rolling average, maximum 14-day rolling average, and maximum 30-day rolling average. Connecing lines are blue if log-linear inerpolaion performs beer and red if i performs worse han linear inerpolaion Supplemenal Figure S10. Ordinary kriging predicions (red) and corresponding predicion inervals a he 90% level (orange) for each possible sample of every 14-day sampling for sie KS-02 in S8
9 154 S9
10 TABLES Supplemenal Table S1. Analysis daa se of sie-years, usage season sample sizes and percen days sampled in he Arazine Ecological Monioring Program (AEMP). 158 No. Days No. Days Sie Year Sampled Sampled Percen Percen (Max=122 Sie Year (Max=122 Complee Complee in April 1- in April 1- July 31) July 31) IA MO-01D IA MO IA MO IA MO-04A IA MO-04A IL MO IL MO-05B IL MO IL MO IL MO IL MO IN MO KS NE KS NE KS NE KS NE KS NE-08B KS NE MO NE MO-01C OH MO-01D OH S10
11 Supplemenal Table S2. Comparison of linear versus nonlinear semivariogram models using mean RRMSPE (i.e., RMSPE/observed maximum rolling average) for maximum rolling averages of duraion m=1,7,14,30, and for 7-day sampling. 164 Duraion of Rolling Average 1 Day 7 Days 14 Days 30 Days Linear semivariogram model Nonlinear semivariogram model Raio beween linear and nonlinear models S11
12 166 SUPPLEMENTAL FIGURES 167 Supplemenal Figure S S12
13 170 Supplemenal Figure S S13
14 174 Supplemenal Figure S S14
15 178 Supplemenal Figure S S15
16 183 Supplemenal Figure S5 184 S16
17 185 Supplemenal Figure S6 186 S17
18 187 Supplemenal Figure S7 188 S18
19 189 Supplemenal Figure S8 190 S19
20 191 Supplemenal Figure S9 192 S20
21 193 Supplemenal Figure S S21
State-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 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 informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
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 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 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 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 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 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 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 informationGeneralized Least Squares
Generalized Leas Squares Augus 006 1 Modified Model Original assumpions: 1 Specificaion: y = Xβ + ε (1) Eε =0 3 EX 0 ε =0 4 Eεε 0 = σ I In his secion, we consider relaxing assumpion (4) Insead, assume
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 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 informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
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 information(10) (a) Derive and plot the spectrum of y. Discuss how the seasonality in the process is evident in spectrum.
January 01 Final Exam Quesions: Mark W. Wason (Poins/Minues are given in Parenheses) (15) 1. Suppose ha y follows he saionary AR(1) process y = y 1 +, where = 0.5 and ~ iid(0,1). Le x = (y + y 1 )/. (11)
More informationReferences are appeared in the last slide. Last update: (1393/08/19)
SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be
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 informationHow to Deal with Structural Breaks in Practical Cointegration Analysis
How o Deal wih Srucural Breaks in Pracical Coinegraion Analysis Roselyne Joyeux * School of Economic and Financial Sudies Macquarie Universiy December 00 ABSTRACT In his noe we consider he reamen of srucural
More informationA Specification Test for Linear Dynamic Stochastic General Equilibrium Models
Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models
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 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 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 informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
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 informationInstitute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler
MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion
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 information4.1 Other Interpretations of Ridge Regression
CHAPTER 4 FURTHER RIDGE THEORY 4. Oher Inerpreaions of Ridge Regression In his secion we will presen hree inerpreaions for he use of ridge regression. The firs one is analogous o Hoerl and Kennard reasoning
More informationBias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé
Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070
More informationWednesday, November 7 Handout: Heteroskedasticity
Amhers College Deparmen of Economics Economics 360 Fall 202 Wednesday, November 7 Handou: Heeroskedasiciy Preview Review o Regression Model o Sandard Ordinary Leas Squares (OLS) Premises o Esimaion Procedures
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 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 informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More informationSolutions to Odd Number Exercises in Chapter 6
1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b
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 informationThe Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form
Chaper 6 The Simple Linear Regression Model: Reporing he Resuls and Choosing he Funcional Form To complee he analysis of he simple linear regression model, in his chaper we will consider how o measure
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
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 informationRegression with Time Series Data
Regression wih Time Series Daa y = β 0 + β 1 x 1 +...+ β k x k + u Serial Correlaion and Heeroskedasiciy Time Series - Serial Correlaion and Heeroskedasiciy 1 Serially Correlaed Errors: Consequences Wih
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 informationThe average rate of change between two points on a function is d t
SM Dae: Secion: Objecive: The average rae of change beween wo poins on a funcion is d. For example, if he funcion ( ) represens he disance in miles ha a car has raveled afer hours, hen finding he slope
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 informationACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin
ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model
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 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 informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
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 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 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 informationDynamic Models, Autocorrelation and Forecasting
ECON 4551 Economerics II Memorial Universiy of Newfoundland Dynamic Models, Auocorrelaion and Forecasing Adaped from Vera Tabakova s noes 9.1 Inroducion 9.2 Lags in he Error Term: Auocorrelaion 9.3 Esimaing
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 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 informationSTAD57 Time Series Analysis. Lecture 14
STAD57 Time Series Analysis Lecure 14 1 Maximum Likelihood AR(p) Esimaion Insead of Yule-Walker (MM) for AR(p) model, can use Maximum Likelihood (ML) esimaion Likelihood is join densiy of daa {x 1,,x n
More informationChristos Papadimitriou & Luca Trevisan November 22, 2016
U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationLesson 2, page 1. Outline of lesson 2
Lesson 2, page Ouline of lesson 2 Inroduce he Auocorrelaion Coefficien Undersand and define saionariy Discuss ransformaion Discuss rend and rend removal C:\Kyrre\sudier\drgrad\Kurs\series\lecure 02 03022.doc,
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 information0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED
0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable
More informationProperties of Autocorrelated Processes Economics 30331
Properies of Auocorrelaed Processes Economics 3033 Bill Evans Fall 05 Suppose we have ime series daa series labeled as where =,,3, T (he final period) Some examples are he dail closing price of he S&500,
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationSensors, Signals and Noise
Sensors, Signals and Noise COURSE OUTLINE Inroducion Signals and Noise: 1) Descripion Filering Sensors and associaed elecronics rv 2017/02/08 1 Noise Descripion Noise Waveforms and Samples Saisics of Noise
More informationExponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits
DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,
More informationOutline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests
Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange
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 information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationTesting the Random Walk Model. i.i.d. ( ) r
he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp
More informationMath 115 Final Exam December 14, 2017
On my honor, as a suden, I have neiher given nor received unauhorized aid on his academic work. Your Iniials Only: Iniials: Do no wrie in his area Mah 5 Final Exam December, 07 Your U-M ID # (no uniqname):
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 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 informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationEcon Autocorrelation. Sanjaya DeSilva
Econ 39 - Auocorrelaion Sanjaya DeSilva Ocober 3, 008 1 Definiion Auocorrelaion (or serial correlaion) occurs when he error erm of one observaion is correlaed wih he error erm of any oher observaion. This
More informationACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.
ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models
More informationm = 41 members n = 27 (nonfounders), f = 14 (founders) 8 markers from chromosome 19
Sequenial Imporance Sampling (SIS) AKA Paricle Filering, Sequenial Impuaion (Kong, Liu, Wong, 994) For many problems, sampling direcly from he arge disribuion is difficul or impossible. One reason possible
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 informationESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING
Inernaional Journal of Social Science and Economic Research Volume:02 Issue:0 ESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING Chung-ki Min Professor
More information6. COMPUTATION OF CENTILES AND Z-SCORES FOR VELOCITIES BASED ON WEIGHT, LENGTH AND HEAD CIRCUMFERENCE
6. COMPUTATION OF CENTILES AND Z-SCORES FOR VELOCITIES BASED ON WEIGHT, LENGTH AND HEAD CIRCUMFERENCE The same mehod used o calculae ceniles and -scores for he aained growh sandards based on weigh is used
More informationEcological Archives E A1. Meghan A. Duffy, Spencer R. Hall, Carla E. Cáceres, and Anthony R. Ives.
Ecological Archives E9-95-A1 Meghan A. Duffy, pencer R. Hall, Carla E. Cáceres, and Anhony R. ves. 29. Rapid evoluion, seasonaliy, and he erminaion of parasie epidemics. Ecology 9:1441 1448. Appendix A.
More informationCONFIDENCE LIMITS AND THEIR ROBUSTNESS
CONFIDENCE LIMITS AND THEIR ROBUSTNESS Rajendran Raja Fermi Naional Acceleraor laboraory Baavia, IL 60510 Absrac Confidence limis are common place in physics analysis. Grea care mus be aken in heir calculaion
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationChapter 11. Heteroskedasticity The Nature of Heteroskedasticity. In Chapter 3 we introduced the linear model (11.1.1)
Chaper 11 Heeroskedasiciy 11.1 The Naure of Heeroskedasiciy In Chaper 3 we inroduced he linear model y = β+β x (11.1.1) 1 o explain household expendiure on food (y) as a funcion of household income (x).
More informationOrdinary dierential equations
Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36
More informationNon-parametric techniques. Instance Based Learning. NN Decision Boundaries. Nearest Neighbor Algorithm. Distance metric important
on-parameric echniques Insance Based Learning AKA: neares neighbor mehods, non-parameric, lazy, memorybased, or case-based learning Copyrigh 2005 by David Helmbold 1 Do no fi a model (as do LDA, logisic
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationMathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013
Mahemaical Theory and Modeling ISSN -580 (Paper) ISSN 5-05 (Online) Vol, No., 0 www.iise.org The ffec of Inverse Transformaion on he Uni Mean and Consan Variance Assumpions of a Muliplicaive rror Model
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 informationFITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA
FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA N. Okendro Singh Associae Professor (Ag. Sa.), College of Agriculure, Cenral Agriculural Universiy, Iroisemba 795 004, Imphal, Manipur
More informationBlock Diagram of a DCS in 411
Informaion source Forma A/D From oher sources Pulse modu. Muliplex Bandpass modu. X M h: channel impulse response m i g i s i Digial inpu Digial oupu iming and synchronizaion Digial baseband/ bandpass
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
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 informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More informationA Bayesian Approach to Spectral Analysis
Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2
More informationMathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation
Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion
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 informationOrdinary Differential Equations
Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described
More informationLecture 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 informationRobot Motion Model EKF based Localization EKF SLAM Graph SLAM
Robo Moion Model EKF based Localizaion EKF SLAM Graph SLAM General Robo Moion Model Robo sae v r Conrol a ime Sae updae model Noise model of robo conrol Noise model of conrol Robo moion model
More information