A Time-Varying Threshold STAR Model of Unemployment
|
|
- Jemimah Barton
- 6 years ago
- Views:
Transcription
1 A Time-Varying Threshold STAR Model of Unemloyment michael dueker a michael owyang b martin sola c,d a Russell Investments b Federal Reserve Bank of St. Louis c Deartamento de Economia, Universidad Torcuato di Tella d School of Economics, Mathematics Statistics, Birkbeck College, University of London January 2009 Abstract Smooth-transition autoregressive (STAR) models have roven to be worthy cometitors of Markov switching models of regime shifts, but the assumtion of a time-invariant threshold level does not seem realistic it holds back this class of models from reaching their otential usefulness. Indeed, an estimate of a time-varying threshold level of unemloyment, for examle, might serve as a meaningful estimate of the natural rate of unemloyment. More recisely, within a STAR framework, one might call the time-varying threshold the tiing level rate of unemloyment, at which the mean dynamics of the unemloyment rate shift. In addition, once the threshold level is allowed to be time-varying, one can add an error-correction term between the lagged level of unemloyment the lagged threshold level to the autoregressive terms in the STAR model. In this way, the time-varying latent threshold level serves dual roles: as a demarcation between regimes as art of an error-correction term. Keywords: Regime switching, smooth-transition autoregressive model, unemloyment, nonlinear models. JEL Classification: C22; E31; G12. 1
2 1 Introduction Our starting oint in modeling the rate of unemloyment is the contemoraneous smooth-transition autoregressive model from Dueker, Sola Sagnolo (2007), which weights the two regimes by the ex ante robability that the unemloyment rate will be below/above the contemoraneous value of the threshold: P (y t < y t I t 1, Θ 0 ) P (y t y t I t 1, Θ 1 ), where Θ 0 Θ 1 are the arameters that ertain to the two regimes the resective regime weights are P 0t = P 1t = P (y t < yt I t 1, Θ 0 ) P (y t < yt I t 1, Θ 0 ) + P (y t yt I t 1, Θ 1 ) P (y t yt I t 1, Θ 1 ) P (y t < yt I t 1, Θ 0 ) + P (y t yt I t 1, Θ 1 ) (1) (2) Note that P 0t + P 1t = 1 this way. With the time-varying threshold comes the ossibility of adding an errorcorrection term to a STAR model: y t = P 0t µ 0 + P 1 tµ 1 + i=1 i=1 Φ (i) 0 y t i + Γ 0 (y t 1 y t 1) + Φ (i) 1 y t i + Γ 1 (y t 1 y t 1) + ɛ t We model the time-varying threshold level of the unemloyment rate as an autoregressive rocess, although the coefficient θ could be close to one to accommodate a high level of ersistence in the threshold level: (4) y t = λ + θy t 1 + w t (5) 2
3 The threshold is endogenous in that its innovation is correlated with the shocks to the observable data: Cov(ɛ t, w t ) 0. the regime-secific covariance matrices are denoted Ω 0 Ω 1. Estimation algorithm for the endogenous, timevarying threshold model Because the vector of latent thresholds enters the regime weights, P 0t P 1t, the model cannot be estimated by maximum likelihood. Nor is it ossible to derive exact conditional distributions for Gibbs samling. Instead, a Metroolis-Hastings algorithm is needed for Bayesian estimation of this model. The key is to have a good, yet tractable, roosal draw for the latent threshold series, {y t } T t=1. The arameter grouings for Markov Chain Monte Carlo estimation of the model are: {yt } T t=1 latent threshold series M-H draw, unscented Kalman filter roosal φ, ρ, µ, γ, λ, θ regression coeffs. M-H draw Normal roosal Ω 0, Ω 1 cov. matrices inverted Wishart Estimation results for U.S. unemloyment We alied the model to the U.S. unemloyment rate since Figure 1 shows the relationshi we found between the unemloyment rate the estimated time-varying threshold level. We also include out-of-samle forecasts of both the unemloyment rate the threshold level. Starting from November 2008, when the unemloyment rate was 6.7 ercent, the STAR model redicts that the unemloyment rate would reach 7.5 ercent at its eak in November This rojected eak seems low, but it is imortant to bear in mind that this forecast is coming from a univariate model of the unemloyment rate. With additional information about the current state of the business cycle, the model should be able to rovide accurate forecasts 3
4 of unemloyment kee the interesting self-referential unemloyment dynamics catured by the model. Aendix: analysis using the unscented Kalman filter We can rewrite the model in the receding section in a state sace reresentation y t = Hz t + ε t, z t = α + g(z t 1 ) + u t (6) where the nonlinear function g(.) contains the cumulative density weighting function P at. Given this state-sace reresentation, inferred values for the latent variable z t can be obtained from the unscented Kalman filter (UKF). The UKF is a nonlinear filter that serves as an alternative to the extended Kalman filter, which uses first-order Taylor-series aroximations to any nonlinear functions in the measurement transition equations. The UKF tracks the state variable by comuting its distribution across a set of deterministic oints called sigma oints. We begin with a set of initial values. We then augment the mean with the exectation of the transition noise z a t 1 t 1 = zt 1 t 1 E u t augment the state covariance P a t 1 t 1 = Pt 1 t 1 0 N,N 0 N,N I N, where s t 1 t 1 P t 1 t 1 are the estimates of the state its covariance matrix at time t 1. Our task is to construct a set of 2L + 1 sigma oints, 4
5 where L is the dimension of s a t 1 t 1. t 1 t 1 = Let z a t 1 t 1, for = 0 z a t 1 t 1 + ( ) (L + λ) Pt 1 t 1 a, for = 1,..., L i z a t 1 t 1 ( ) (L + λ) Pt 1 t 1 a, for = L + 1,..., 2L define the initial sigma oints, where λ = a 2 (L + κ) L, a κ are user-chosen ( ) arameters that govern the sread scale, resectively. Here, X is the ith column of the lower triangular Cholesky factorization of i the square marix X. Given the set of initial sigma oints, we can then roagate { t 1 t 1} through the transition function g (.) to recover i L t t 1 = g ( t 1 t 1), for = 0,..., 2L. The redicted states covariances can then be extraolated from a weighted sum of the roagated sigma oints P t t 1 = wc =0 The weights are defined as w c = w s = { λ { ẑ t t 1 = =0 w s t t 1 t t 1 ẑ t t 1 t t 1 ẑ t t 1. λ, for = 0 L+λ 1 2(L+λ), for = 1,..., 2L + (1 L+λ a2 + b), for = 0 1, for = 1,..., 2L, 2(L+λ) where a is defined as above b is a arameter that incororates rior information of the distribution. 5
6 Filter recursions The udating ste roceeds in a similar manner. We can augment the redicted state covariance with the exectation covariance of the measurement noise z a t t 1 = ẑt t 1 v t P a t t 1 = Pt t 1 0 n,n 0 n,n Ω, where Ω is the covariance of the measurement noise. We again comute 2L + 1 sigma oints from t t 1 = z a t t 1, for = 0 z a t t 1 + ( ) (L + λ) Pt t 1 a, for = 1,..., L i z a t t 1 ( ). (L + λ) Pt t 1 a, for = L + 1,..., 2L i L We then roagate the udate sigma oints through the measurement equation: γ t t = H t t 1, = 0,..., 2L. The udated roagation oints, γ t t, are used to construct the redicted measurement error ŷ t = =0 w sγ t t, where the weights are defined as before. We then form the udated state similar to the stard Kalman filter ẑ t t = ẑ t t 1 + K t (y t ŷ t ), where K t is the Kalman gain defined by K t = P yz P 1 yy. 6
7 Here, P yz defines the cross-covariance P yz = wc =0 P yy is the redicted covariance P yy = wc =0 The udated covariance is defined by Smoothing t t 1 ẑ t t 1 γ t t ŷ t γ t t ŷ t γ t t ŷ t. P t t = P t t 1 K t P yy K t. Multi-move samling of the latent attractor, {s } T t=1, requires backwards samling from one-eriod smoothed inferences of the state vector. Here the unscented Rauch-Tung-Striebel smoother comes into lay. The URTS smoother (Sarkka, 2007) begins by augmenting the unscented Kalman filter with a smoothing ste that recomutes the state estimate. The smoothed state is with covariance matrix P s t ẑ s t = ẑ t + D t ẑ s t+1 ẑ t+1 t = P t + D t P s t+1 P t+1 t D t, where D t is the smoother gain defined by where with P zt,z t+1 = D t = P zt,z t+1 P 1 t+1 t, wc =0 t+1 t ẑ t+1 t t t ẑ t t, t t = t t 1 + K t (y t ŷ t t 1 ). 7
8 References 1 Julier, Simon J. Jeffery K. Uhlmann. A New Extension of the Kalman Filter to nonlinear Systems, in The Proceedings of AeroSense: The 11th International Symosium on Aerosace/Defense Sensing, Simulation Controls, Multi Sensor Fusion, Tracking Resource Management II, SPIE, Särkkä, Simo. Unscented Rauch-Tung-Striebel Smoothing, forthcoming IEEE, Wan, Eric van der Merwe, Rudolh. The Unscented Kalman Filter, in S. Hakin, Kalman Filtering Neural Networks, Wiley, Setember
9 Unemloyment rate latent threshold level with out-ofsamle forecasts /1/1968 2/1/1972 2/1/1976 2/1/1980 2/1/1984 2/1/1988 2/1/1992 2/1/1996 2/1/2000 2/1/2004 2/1/2008 ercent Latent threshold Actual Forecasts
Research Division Federal Reserve Bank of St. Louis Working Paper Series
Research Division Federal Reserve Bank of St. Louis Working Paper Series A Time-Varying Threshold STAR Model of Unemployment and the Natural Rate Michael Dueker Michael T. Owyang and Martin Sola Working
More informationResearch Division Federal Reserve Bank of St. Louis Working Paper Series
Research Division Federal Reserve Bank of St Louis Working Paper Series Kalman Filtering with Truncated Normal State Variables for Bayesian Estimation of Macroeconomic Models Michael Dueker Working Paper
More informationNonlinear and/or Non-normal Filtering. Jesús Fernández-Villaverde University of Pennsylvania
Nonlinear and/or Non-normal Filtering Jesús Fernández-Villaverde University of Pennsylvania 1 Motivation Nonlinear and/or non-gaussian filtering, smoothing, and forecasting (NLGF) problems are pervasive
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationThe Kalman Filter ImPr Talk
The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationNon-Markovian Regime Switching with Endogenous States and Time-Varying State Strengths
Non-Markovian Regime Switching with Endogenous States and Time-Varying State Strengths January 2004 Siddhartha Chib Olin School of Business Washington University chib@olin.wustl.edu Michael Dueker Federal
More informationDistributed Rule-Based Inference in the Presence of Redundant Information
istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced
More informationA new unscented Kalman filter with higher order moment-matching
A new unscented Kalman filter with higher order moment-matching KSENIA PONOMAREVA, PARESH DATE AND ZIDONG WANG Department of Mathematical Sciences, Brunel University, Uxbridge, UB8 3PH, UK. Abstract This
More informationPrediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother
Prediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother Simo Särkkä, Aki Vehtari and Jouko Lampinen Helsinki University of Technology Department of Electrical and Communications
More informationSelf-Driving Car ND - Sensor Fusion - Extended Kalman Filters
Self-Driving Car ND - Sensor Fusion - Extended Kalman Filters Udacity and Mercedes February 7, 07 Introduction Lesson Ma 3 Estimation Problem Refresh 4 Measurement Udate Quiz 5 Kalman Filter Equations
More informationin a Rao-Blackwellised Unscented Kalman Filter
A Rao-Blacwellised Unscented Kalman Filter Mar Briers QinetiQ Ltd. Malvern Technology Centre Malvern, UK. m.briers@signal.qinetiq.com Simon R. Masell QinetiQ Ltd. Malvern Technology Centre Malvern, UK.
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationNotes on Instrumental Variables Methods
Notes on Instrumental Variables Methods Michele Pellizzari IGIER-Bocconi, IZA and frdb 1 The Instrumental Variable Estimator Instrumental variable estimation is the classical solution to the roblem of
More informationNonlinear State Estimation! Particle, Sigma-Points Filters!
Nonlinear State Estimation! Particle, Sigma-Points Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Particle filter!! Sigma-Points Unscented Kalman ) filter!!
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Extended Kalman Filter Dr. Kostas Alexis (CSE) These slides relied on the lectures from C. Stachniss, J. Sturm and the book Probabilistic Robotics from Thurn et al.
More informationLearning Sequence Motif Models Using Gibbs Sampling
Learning Sequence Motif Models Using Gibbs Samling BMI/CS 776 www.biostat.wisc.edu/bmi776/ Sring 2018 Anthony Gitter gitter@biostat.wisc.edu These slides excluding third-arty material are licensed under
More informationShort course A vademecum of statistical pattern recognition techniques with applications to image and video analysis. Agenda
Short course A vademecum of statistical attern recognition techniques with alications to image and video analysis Lecture 6 The Kalman filter. Particle filters Massimo Piccardi University of Technology,
More informationCollaborative Place Models Supplement 1
Collaborative Place Models Sulement Ber Kaicioglu Foursquare Labs ber.aicioglu@gmail.com Robert E. Schaire Princeton University schaire@cs.rinceton.edu David S. Rosenberg P Mobile Labs david.davidr@gmail.com
More informationLecture 7: Optimal Smoothing
Department of Biomedical Engineering and Computational Science Aalto University March 17, 2011 Contents 1 What is Optimal Smoothing? 2 Bayesian Optimal Smoothing Equations 3 Rauch-Tung-Striebel Smoother
More informationFiltering and Likelihood Inference
Filtering and Likelihood Inference Jesús Fernández-Villaverde University of Pennsylvania July 10, 2011 Jesús Fernández-Villaverde (PENN) Filtering and Likelihood July 10, 2011 1 / 79 Motivation Introduction
More informationIntroduction to Unscented Kalman Filter
Introduction to Unscented Kalman Filter 1 Introdution In many scientific fields, we use certain models to describe the dynamics of system, such as mobile robot, vision tracking and so on. The word dynamics
More informationPerformance of lag length selection criteria in three different situations
MPRA Munich Personal RePEc Archive Performance of lag length selection criteria in three different situations Zahid Asghar and Irum Abid Quaid-i-Azam University, Islamabad Aril 2007 Online at htts://mra.ub.uni-muenchen.de/40042/
More informationLecture 6. 2 Recurrence/transience, harmonic functions and martingales
Lecture 6 Classification of states We have shown that all states of an irreducible countable state Markov chain must of the same tye. This gives rise to the following classification. Definition. [Classification
More informationDual Estimation and the Unscented Transformation
Dual Estimation and the Unscented Transformation Eric A. Wan ericwan@ece.ogi.edu Rudolph van der Merwe rudmerwe@ece.ogi.edu Alex T. Nelson atnelson@ece.ogi.edu Oregon Graduate Institute of Science & Technology
More informationDynamic Factor Models and Factor Augmented Vector Autoregressions. Lawrence J. Christiano
Dynamic Factor Models and Factor Augmented Vector Autoregressions Lawrence J Christiano Dynamic Factor Models and Factor Augmented Vector Autoregressions Problem: the time series dimension of data is relatively
More informationGeneral Linear Model Introduction, Classes of Linear models and Estimation
Stat 740 General Linear Model Introduction, Classes of Linear models and Estimation An aim of scientific enquiry: To describe or to discover relationshis among events (variables) in the controlled (laboratory)
More information8 STOCHASTIC PROCESSES
8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular
More informationFinancial Econometrics
Financial Econometrics Nonlinear time series analysis Gerald P. Dwyer Trinity College, Dublin January 2016 Outline 1 Nonlinearity Does nonlinearity matter? Nonlinear models Tests for nonlinearity Forecasting
More informationLower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data
Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment
More informationBayesian Dynamic Linear Modelling for. Complex Computer Models
Bayesian Dynamic Linear Modelling for Complex Computer Models Fei Liu, Liang Zhang, Mike West Abstract Computer models may have functional outputs. With no loss of generality, we assume that a single computer
More informationUniform Law on the Unit Sphere of a Banach Space
Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a
More informationON THE L P -CONVERGENCE OF A GIRSANOV THEOREM BASED PARTICLE FILTER. Eric Moulines. Center of Applied Mathematics École Polytechnique, France
O TH L P -COVRGC OF A GIRSAOV THORM BASD PARTICL FILTR Simo Särä Det. of lectrical ngineering and Automation Aalto University Finland ric Moulines Center of Alied Mathematics École Polytechnique France
More informationNon-Markovian Regime Switching with Endogenous States and Time-Varying State Strengths
WORKING PAPER SERIES Non-Markovian Regime Switching with Endogenous States and Time-Varying State Strengths Siddhartha Chib and Michael Dueker Working Paper 2004-030A http://research.stlouisfed.org/wp/2004/2004-030.pdf
More informationLabor-Supply Shifts and Economic Fluctuations. Technical Appendix
Labor-Supply Shifts and Economic Fluctuations Technical Appendix Yongsung Chang Department of Economics University of Pennsylvania Frank Schorfheide Department of Economics University of Pennsylvania January
More informationVAR models with non-gaussian shocks
VAR models with non-gaussian shocks Ching-Wai (Jeremy) Chiu Haroon Mumtaz Gabor Pinter September 27, 2016 Motivation and aims Density of the residuals from a recursive VAR(13) (1960m1-2015m6) Motivation
More informationConstrained State Estimation Using the Unscented Kalman Filter
16th Mediterranean Conference on Control and Automation Congress Centre, Ajaccio, France June 25-27, 28 Constrained State Estimation Using the Unscented Kalman Filter Rambabu Kandepu, Lars Imsland and
More informationScaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling
Scaling Multile Point Statistics or Non-Stationary Geostatistical Modeling Julián M. Ortiz, Steven Lyster and Clayton V. Deutsch Centre or Comutational Geostatistics Deartment o Civil & Environmental Engineering
More informationBayesian inference & Markov chain Monte Carlo. Note 1: Many slides for this lecture were kindly provided by Paul Lewis and Mark Holder
Bayesian inference & Markov chain Monte Carlo Note 1: Many slides for this lecture were kindly rovided by Paul Lewis and Mark Holder Note 2: Paul Lewis has written nice software for demonstrating Markov
More informationSAS for Bayesian Mediation Analysis
Paer 1569-2014 SAS for Bayesian Mediation Analysis Miočević Milica, Arizona State University; David P. MacKinnon, Arizona State University ABSTRACT Recent statistical mediation analysis research focuses
More informationASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE
J Jaan Statist Soc Vol 34 No 2004 9 26 ASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE Yasunori Fujikoshi*, Tetsuto Himeno
More informationChapter 7: Special Distributions
This chater first resents some imortant distributions, and then develos the largesamle distribution theory which is crucial in estimation and statistical inference Discrete distributions The Bernoulli
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationMODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS
MODEL-BASED MULIPLE FAUL DEECION AND ISOLAION FOR NONLINEAR SYSEMS Ivan Castillo, and homas F. Edgar he University of exas at Austin Austin, X 78712 David Hill Chemstations Houston, X 77009 Abstract A
More informationOpen Economy Macroeconomics: Theory, methods and applications
Open Economy Macroeconomics: Theory, methods and applications Lecture 4: The state space representation and the Kalman Filter Hernán D. Seoane UC3M January, 2016 Today s lecture State space representation
More informationFinite Mixture EFA in Mplus
Finite Mixture EFA in Mlus November 16, 2007 In this document we describe the Mixture EFA model estimated in Mlus. Four tyes of deendent variables are ossible in this model: normally distributed, ordered
More informationOnline appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US
Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US Gerdie Everaert 1, Lorenzo Pozzi 2, and Ruben Schoonackers 3 1 Ghent University & SHERPPA 2 Erasmus
More informationEinführung in Stochastische Prozesse und Zeitreihenanalyse Vorlesung, 2013S, 2.0h March 2015 Hubalek/Scherrer
Name: Mat.Nr.: Studium: Bitte keinen Rotstift verwenden! 15.593 Einführung in Stochastische Prozesse und Zeitreihenanalyse Vorlesung, 213S, 2.h March 215 Hubalek/Scherrer (Dauer 9 Minutes, Permissible
More informationUse of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek
Use of Transformations and the Reeated Statement in PROC GLM in SAS Ed Stanek Introduction We describe how the Reeated Statement in PROC GLM in SAS transforms the data to rovide tests of hyotheses of interest.
More informationOn Wrapping of Exponentiated Inverted Weibull Distribution
IJIRST International Journal for Innovative Research in Science & Technology Volume 3 Issue 11 Aril 217 ISSN (online): 2349-61 On Wraing of Exonentiated Inverted Weibull Distribution P.Srinivasa Subrahmanyam
More informationKalman filtering and friends: Inference in time series models. Herke van Hoof slides mostly by Michael Rubinstein
Kalman filtering and friends: Inference in time series models Herke van Hoof slides mostly by Michael Rubinstein Problem overview Goal Estimate most probable state at time k using measurement up to time
More informationA Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition
A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino
More informationApplied Fitting Theory VI. Formulas for Kinematic Fitting
Alied Fitting heory VI Paul Avery CBX 98 37 June 9, 1998 Ar. 17, 1999 (rev.) I Introduction Formulas for Kinematic Fitting I intend for this note and the one following it to serve as mathematical references,
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationThe Unscented Particle Filter
The Unscented Particle Filter Rudolph van der Merwe (OGI) Nando de Freitas (UC Bereley) Arnaud Doucet (Cambridge University) Eric Wan (OGI) Outline Optimal Estimation & Filtering Optimal Recursive Bayesian
More informationSTA 250: Statistics. Notes 7. Bayesian Approach to Statistics. Book chapters: 7.2
STA 25: Statistics Notes 7. Bayesian Aroach to Statistics Book chaters: 7.2 1 From calibrating a rocedure to quantifying uncertainty We saw that the central idea of classical testing is to rovide a rigorous
More informationNON-LINEAR NOISE ADAPTIVE KALMAN FILTERING VIA VARIATIONAL BAYES
2013 IEEE INTERNATIONAL WORKSHOP ON MACHINE LEARNING FOR SIGNAL PROCESSING NON-LINEAR NOISE ADAPTIVE KALMAN FILTERING VIA VARIATIONAL BAYES Simo Särä Aalto University, 02150 Espoo, Finland Jouni Hartiainen
More informationECE 534 Information Theory - Midterm 2
ECE 534 Information Theory - Midterm Nov.4, 009. 3:30-4:45 in LH03. You will be given the full class time: 75 minutes. Use it wisely! Many of the roblems have short answers; try to find shortcuts. You
More informationPlotting the Wilson distribution
, Survey of English Usage, University College London Setember 018 1 1. Introduction We have discussed the Wilson score interval at length elsewhere (Wallis 013a, b). Given an observed Binomial roortion
More informationPArtially observable Markov decision processes
Solving Continuous-State POMDPs via Density Projection Enlu Zhou, Member, IEEE, Michael C. Fu, Fellow, IEEE, and Steven I. Marcus, Fellow, IEEE Abstract Research on numerical solution methods for artially
More informationRECURSIVE OUTLIER-ROBUST FILTERING AND SMOOTHING FOR NONLINEAR SYSTEMS USING THE MULTIVARIATE STUDENT-T DISTRIBUTION
1 IEEE INTERNATIONAL WORKSHOP ON MACHINE LEARNING FOR SIGNAL PROCESSING, SEPT. 3 6, 1, SANTANDER, SPAIN RECURSIVE OUTLIER-ROBUST FILTERING AND SMOOTHING FOR NONLINEAR SYSTEMS USING THE MULTIVARIATE STUDENT-T
More informationOptimal filtering with Kalman filters and smoothers a Manual for Matlab toolbox EKF/UKF
Optimal filtering with Kalman filters and smoothers a Manual for Matlab toolbox EKF/UKF Jouni Hartikainen and Simo Särkkä Laboratory of Computational Engineering, Helsinki University of Technology, P.O.Box
More informationThe Scaled Unscented Transformation
The Scaled Unscented Transformation Simon J. Julier, IDAK Industries, 91 Missouri Blvd., #179 Jefferson City, MO 6519 E-mail:sjulier@idak.com Abstract This paper describes a generalisation of the unscented
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY PREFACE xiii 1 Difference Equations 1.1. First-Order Difference Equations 1 1.2. pth-order Difference Equations 7
More informationVolume 30, Issue 3. A note on Kalman filter approach to solution of rational expectations models
Volume 30, Issue 3 A note on Kalman filter approach to solution of rational expectations models Marco Maria Sorge BGSE, University of Bonn Abstract In this note, a class of nonlinear dynamic models under
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationX t = a t + r t, (7.1)
Chapter 7 State Space Models 71 Introduction State Space models, developed over the past 10 20 years, are alternative models for time series They include both the ARIMA models of Chapters 3 6 and the Classical
More informationVector Auto-Regressive Models
Vector Auto-Regressive Models Laurent Ferrara 1 1 University of Paris Nanterre M2 Oct. 2018 Overview of the presentation 1. Vector Auto-Regressions Definition Estimation Testing 2. Impulse responses functions
More informationState Estimation of Linear and Nonlinear Dynamic Systems
State Estimation of Linear and Nonlinear Dynamic Systems Part III: Nonlinear Systems: Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF) James B. Rawlings and Fernando V. Lima Department of
More informationVAR Models and Applications
VAR Models and Applications Laurent Ferrara 1 1 University of Paris West M2 EIPMC Oct. 2016 Overview of the presentation 1. Vector Auto-Regressions Definition Estimation Testing 2. Impulse responses functions
More informationPATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA
PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA Contents in latter part Linear Dynamical Systems What is different from HMM? Kalman filter Its strength and limitation Particle Filter
More informationSolutions to exercises on delays. P (x = 0 θ = 1)P (θ = 1) P (x = 0) We can replace z in the first equation by its value in the second equation.
Ec 517 Christohe Chamley Solutions to exercises on delays Ex 1: P (θ = 1 x = 0) = P (x = 0 θ = 1)P (θ = 1) P (x = 0) = 1 z)µ (1 z)µ + 1 µ. The value of z is solution of µ c = δµz(1 c). We can relace z
More informationBayesian Networks Practice
Bayesian Networks Practice Part 2 2016-03-17 Byoung-Hee Kim, Seong-Ho Son Biointelligence Lab, CSE, Seoul National University Agenda Probabilistic Inference in Bayesian networks Probability basics D-searation
More informationCausality Testing using Higher Order Statistics
Causality Testing using Higher Order Statistics Dr Sanya Dudukovic International Management Deartment Franklin College, Switzerland Fax: 41 91 994 41 17 E-mail : Sdudukov@fc.edu Abstract : A new causality
More informationIdentifying the Monetary Policy Shock Christiano et al. (1999)
Identifying the Monetary Policy Shock Christiano et al. (1999) The question we are asking is: What are the consequences of a monetary policy shock a shock which is purely related to monetary conditions
More informationState-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Fin. Econometrics / 53
State-space Model Eduardo Rossi University of Pavia November 2014 Rossi State-space Model Fin. Econometrics - 2014 1 / 53 Outline 1 Motivation 2 Introduction 3 The Kalman filter 4 Forecast errors 5 State
More informationEM-algorithm for Training of State-space Models with Application to Time Series Prediction
EM-algorithm for Training of State-space Models with Application to Time Series Prediction Elia Liitiäinen, Nima Reyhani and Amaury Lendasse Helsinki University of Technology - Neural Networks Research
More informationCOMMUNICATION BETWEEN SHAREHOLDERS 1
COMMUNICATION BTWN SHARHOLDRS 1 A B. O A : A D Lemma B.1. U to µ Z r 2 σ2 Z + σ2 X 2r ω 2 an additive constant that does not deend on a or θ, the agents ayoffs can be written as: 2r rθa ω2 + θ µ Y rcov
More informationBayesian Inference for DSGE Models. Lawrence J. Christiano
Bayesian Inference for DSGE Models Lawrence J. Christiano Outline State space-observer form. convenient for model estimation and many other things. Bayesian inference Bayes rule. Monte Carlo integation.
More informationA Comparison of the EKF, SPKF, and the Bayes Filter for Landmark-Based Localization
A Comparison of the EKF, SPKF, and the Bayes Filter for Landmark-Based Localization and Timothy D. Barfoot CRV 2 Outline Background Objective Experimental Setup Results Discussion Conclusion 2 Outline
More informationSmoothers: Types and Benchmarks
Smoothers: Types and Benchmarks Patrick N. Raanes Oxford University, NERSC 8th International EnKF Workshop May 27, 2013 Chris Farmer, Irene Moroz Laurent Bertino NERSC Geir Evensen Abstract Talk builds
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationBrownian Motion and Random Prime Factorization
Brownian Motion and Random Prime Factorization Kendrick Tang June 4, 202 Contents Introduction 2 2 Brownian Motion 2 2. Develoing Brownian Motion.................... 2 2.. Measure Saces and Borel Sigma-Algebras.........
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
More informationPDF hosted at the Radboud Repository of the Radboud University Nijmegen
PDF hosted at the Radboud Repository of the Radboud University Nijmegen The following full text is a preprint version which may differ from the publisher's version. For additional information about this
More informationChapter 7 Sampling and Sampling Distributions. Introduction. Selecting a Sample. Introduction. Sampling from a Finite Population
Chater 7 and s Selecting a Samle Point Estimation Introduction to s of Proerties of Point Estimators Other Methods Introduction An element is the entity on which data are collected. A oulation is a collection
More informationTAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω
ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.
More informationECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering
ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering Lecturer: Nikolay Atanasov: natanasov@ucsd.edu Teaching Assistants: Siwei Guo: s9guo@eng.ucsd.edu Anwesan Pal:
More informationIterative Filtering and Smoothing In Non-Linear and Non-Gaussian Systems Using Conditional Moments
Iterative Filtering and Smoothing In Non-Linear and Non-Gaussian Systems Using Conditional Moments Filip Tronarp, Ángel F. García-Fernández, and Simo Särkkä This is a post-print of a paper published in
More informationThe Recursive Fitting of Multivariate. Complex Subset ARX Models
lied Mathematical Sciences, Vol. 1, 2007, no. 23, 1129-1143 The Recursive Fitting of Multivariate Comlex Subset RX Models Jack Penm School of Finance and lied Statistics NU College of Business & conomics
More informationLecture 16: State Space Model and Kalman Filter Bus 41910, Time Series Analysis, Mr. R. Tsay
Lecture 6: State Space Model and Kalman Filter Bus 490, Time Series Analysis, Mr R Tsay A state space model consists of two equations: S t+ F S t + Ge t+, () Z t HS t + ɛ t (2) where S t is a state vector
More informationMultivariate Markov Switching With Weighted Regime Determination: Giving France More Weight than Finland
Multivariate Markov Switching With Weighted Regime Determination: Giving France More Weight than Finland November 2006 Michael Dueker Federal Reserve Bank of St. Louis P.O. Box 442, St. Louis, MO 63166
More informationPosterior Cramer-Rao Lower Bound for Mobile Tracking in Mixed Line-of-Sight/Non Line-of-Sight Conditions
Posterior Cramer-Rao Lower Bound for Mobile Tracking in Mixed Line-of-Sight/Non Line-of-Sight Conditions Chen Liang 1,2, Wu Lenan 2, Robert Piché 1 1 Tampere University of Technology, Finland 2 Southeast
More informationIntroduction to Machine Learning
Introduction to Machine Learning Brown University CSCI 1950-F, Spring 2012 Prof. Erik Sudderth Lecture 25: Markov Chain Monte Carlo (MCMC) Course Review and Advanced Topics Many figures courtesy Kevin
More informationElements of Asymptotic Theory. James L. Powell Department of Economics University of California, Berkeley
Elements of Asymtotic Theory James L. Powell Deartment of Economics University of California, Berkeley Objectives of Asymtotic Theory While exact results are available for, say, the distribution of the
More informationF. Augustin, P. Rentrop Technische Universität München, Centre for Mathematical Sciences
NUMERICS OF THE VAN-DER-POL EQUATION WITH RANDOM PARAMETER F. Augustin, P. Rentro Technische Universität München, Centre for Mathematical Sciences Abstract. In this article, the roblem of long-term integration
More informationKatsuhiro Sugita Faculty of Law and Letters, University of the Ryukyus. Abstract
Bayesian analysis of a vector autoregressive model with multiple structural breaks Katsuhiro Sugita Faculty of Law and Letters, University of the Ryukyus Abstract This paper develops a Bayesian approach
More informationClassical gas (molecules) Phonon gas Number fixed Population depends on frequency of mode and temperature: 1. For each particle. For an N-particle gas
Lecture 14: Thermal conductivity Review: honons as articles In chater 5, we have been considering quantized waves in solids to be articles and this becomes very imortant when we discuss thermal conductivity.
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More information