Timevarying VARs. Wouter J. Den Haan London School of Economics. c Wouter J. Den Haan
|
|
- Edwina Hunt
- 6 years ago
- Views:
Transcription
1 Timevarying VARs Wouter J. Den Haan London School of Economics c Wouter J. Den Haan
2 Time-Varying VARs Gibbs-Sampler general idea probit regression application (Inverted Wishart distribution Drawing from a multi-variate Normal in Matlab Time-varying VAR model specification Gibbs sampler 2 / 38
3 Gibbs Sampler Suppose x, y, and z are distributed according to f (x, y, z Suppose that drawing x, y, and z from f (x, y, z is diffi cult but you can draw x from f (x y, z and you can draw y from f (y x, z and you can draw y from f (z x, y 3 / 38
4 Gibbs Sampler - how it works Start with y 0, z 0 Draw x 1 from f (x y 0, z 0, Draw y 1 from f (y x 1, z 0, Draw z 1 from f (z x 1, y 1, Draw x 2 from f (x y 1, z 1 (x i, y i, z i is one draw from the joint density f (x, y, z Although series are constructed recursively, they are not time series 4 / 38
5 Gibbs Sampler - convergence The idea is that this sequence converges to a sequence drawn from f (x, y, z. Since convergence is not immediate, you have to discard beginning of sequence (burn-in period. See Casella and George (1992 for a discussion on why and when this works. 5 / 38
6 Gibbs Sampler - probit regression This example is from Lancaster (2004 y i is the i th observation of a binary variable, i.e., y i {0, 1} y i is an unobservable and given by y i = x i β + ε i, ε i N(0, 1 { 1 if y y = i 0 0 o.w. 6 / 38
7 Probit regression Parameters: β and y = [y 1, y 2,, y n] Data: X = [x 1, x 2,, x n ], Y = {y i, x i } n i=1 Objective: get p ( β Y, i.e., the distribution of β given Y. With the Gibbs sample we can get a sequence of obervations for ( β, y distributed according to p ( β, ŷ Y, from which we can get p ( β Y 7 / 38
8 Probit - Gibb sampler step 1 We need to draw from p ( β y, Y Given y and X β N ( (X X 1 X y, ( X X 1, since the standard deviation of ε i is known and equal to 1. 8 / 38
9 Probit - Gibb sampler step 2 We need to draw from p (y β, Y Since the y i s are independent, we can do this separately for each i y i N >0 (x i β, 1 if y i = 1 y i N <0 (x i β, 1 if y i = 0, where N >0 ( is a Normal distribution truncated on the left at 0 N <0 ( is a Normal distribution truncated on the right at 0 9 / 38
10 Wishart distribution generalization of Chi-square distribution to more variables X : n p matrix; each row drawn from N p (0, Σ, where Σ is the p p variance-covariance matrix W = X X W p (Σ, n, i.e., the p-dimensional Wishart with scale matrix Σ and degrees of freedom n You get the Chi-square if p = 1 and Σ = 1 10 / 38
11 Inverse Wishart distribution If W has a Wishart distribution with parameters Σ and n, then W 1 has an inverse Wishart with scale matrix Σ 1 and degrees of freedom n!!! In the assignment, the input of the Matlab Inverse Wishart function is Σ not Σ / 38
12 Inverse Wishart in Bayesian statistics Data: x t is a p 1 vector with i.i.d. random observations with distribution N (0, V prior of V : ( p (V = IW V 1, n posterior of V : p ( V X T ( = IW W 1, n + T W = V + V T V T = T x tx t t=1 Note that V T is like a sum of squares 12 / 38
13 Multivariate normal in Matlab x t is a p 1 vector and we want x t N(0, Σ C=chol(Σ Thus C is an upper-triangular matrix and C C = Σ e t is a p 1 vector with draws from N(0, I p E [ C e t e tc ] = Σ Thus, C e t is a p 1 vector with draws from N(0, Σ 13 / 38
14 Time-varying VARs - intro Idea: capture changes in model specification in a flexible way The analysis here is based on Cogley and Sargent (2002, CS 14 / 38
15 Model specification y t = X tθ t + ε t X t = [ 1, y t 1, y t 2,, y t p ] θ t = θ t 1 + v t ε t N(0, R v t N(0, Q 15 / 38
16 Model specification E t [ εt v t ] [ εt v t ] = V = ( R C C Q θ t : "parameters" R, C, and Q are the "hyperparameters" 16 / 38
17 Model specification - details Simplifying assumptions: CS impose that θ t is such that y t would be stationary if θ t+τ = θ t for all τ 0. This stationarity requirement is left out for transparency. C = / 38
18 Notation Y T = [ y 1,, ] y T θ T = [ θ 0, θ 1,, ] θ T 18 / 38
19 Priors Prior for initial condition: θ 0 N ( θ, P Prior for hyperparameters: p(v = IW (V 1, T 0 θ, P, V, T 0 are taken as given 19 / 38
20 Posterior The posterior is given by p (θ T, V Y T We can use the Gibbs sampler if we can draw from ( P θ T Y T, V and from P (V Y T, θ T 20 / 38
21 Stationarity CS exclude draws of θ t for which the dgp of y t is nonstationary: p (θ t is density without imposing stationarity and f (θ t is density with imposing stationarity This restriction is ignored in these slides 21 / 38
22 Gibbs part I: Posterior of theta given V Since f (A, B = f (A B f (B, we have ( ( p θ T Y T, V = f θ T Y T, V ( ( = f θ T 1 θ T, Y T, V f θ T Y T, V ( ( = f θ T 2 θ T, θ T 1, Y T, V f θ T 1 θ T, Y T, V ( f θ T Y T, V ( ( = f θ T 3 θ T, θ T 1, θ T 2, Y T, V f θ t 2 θ T, θ T 1, Y T, V ( ( f θ T 1 θ T, Y T, V f θ T Y T, V 22 / 38
23 Posterior of theta given V Since θ t = θ t 1 + v t, θ t+τ has no predictive power for θ t 1 for all τ 1 given Y T and θ t, Thus ( f θ T 2 θ T, θ T 1, Y T, V ( f θ T 3 θ T, θ T 1, θ T 2, Y T, V ( = f θ T 2 θ T 1, Y T, V ( = f θ T 3 θ T 2, Y T, V etc. 23 / 38
24 Posterior of theta given V Combining gives ( ( p θ T Y T, V = f θ T Y T T 1 (, V f θ t θ t+1, Y T, V t=1 All the densities are Gaussian = if we( know the means and the variances, then we can draw from p θ T Y T, V 24 / 38
25 Posterior of theta given V We need to find the means and variances of ( ( f θ T Y T, V & f θ t θ t+1, Y T, V Notation θ t t = E ( θ t Y t, V ( P t t 1 = VAR θ t Y t 1, V P t t = VAR ( θ t Y t, V θ t t+1 = E ( θ t θ t+1, Y t, V ( = E θ t θ t+1, Y T, V P t t+1 = VAR ( θ t θ t+1, Y t, V ( = VAR θ t θ t+1, Y T, V 25 / 38
26 Posterior of theta given V First, use Kalman filter to go forward start with θ 0 and P 0 0 Next, go backwards to get draws for θ t given θ t+1 26 / 38
27 Posterior of theta given V Kalman filter part: y t = X tθ t + ε t X t = [ 1, y t 1, y t 2,, y t p ] θ t = θ t 1 + v t ε t N(0, R v t N(0, Q Here: the p + 1 elements of X t are the known (time-varying coeffi cients of the state-space represenation the elements of θ t are the unobserved underlying state variables 27 / 38
28 Posterior of theta given V Kalman filter in the first period: and then iterate P 1 0 = P Q 1 K 1 = P 1 0 X 1 (X 1 P 1 0X 1 + R θ 1 1 = θ K 1 (y 1 X 1 θ 0 0 P t t 1 = P t 1 t 1 + Q 1 K t = P t t 1 X t (X tp t t 1 X t + R ( θ t t = θ t 1 t 1 + K t y t X tθ t 1 t 1 P t t = P t t 1 K t X tp t t 1 28 / 38
29 Posterior of theta given V In the Kalman filter part of the assignment: and we go up to TH(:,1 = θ 0 TH(:,t+1 = θ t t Pe(:,:,t = P t 1 t 1 Po(:,:,t = P t t 1 TH(:,T+1 = θ T T Pe(:,:,T+1 = P T T Po(:,:,T = P T T 1 29 / 38
30 Posterior of theta given V Distribution terminal state: ( f θ T Y T, V = N (θ T T, P T T From this we get a draw θ T 30 / 38
31 Posterior of theta given V Draws for θ T 1, θ T 2,, θ 1 are obtained recursively from ( f θ t θ t+1, Y T, V = N (θ t t+1, P t t+1 θ t t+1 = θ t t + P t t P (θ 1 t+1 t t+1 θ t t P t t+1 = P t t P t t P 1 t+1 t P t t The terms needed to calculate θ t t+1 and P t t+1 are generated by the Kalman filter (that is, from going forward and the standard projection formulas (and note that the covariance of θ t+1 and θ t is the variance of θ t 31 / 38
32 Details for previous slide E [y x] = µ y + Σ yx Σxx 1 (x µ x = E [θ t θ t+1 ; ] = E [θ t ] + Σ θt θ t+1 Σ 1 θ t+1 θ t+1 (θ t+1 E [θ t+1 ] = E [θ t ] + Σ θt θ t Σ 1 θ t+1 θ t+1 (θ t+1 E [θ t+1 ] = θ t t+1 = θ t t + P t t P 1 t+1 t (θ t+1 E [θ t + v t ] = θ t t + P t t P 1 t+1 t (θ t+1 E [θ t ] = θ t t + P t t P (θ 1 t+1 t t+1 θ t t Suppressing the dependence on Y t and V to simplify notation 32 / 38
33 Posterior of theta given V In the backward part of the assignment: Draw from TH(t-1 t In the for loop below t goes from high to low. At a particular t: 1 TH(:,t+1 it is a random draw from a normal that has already been determined (either in this loop or for T above 2 TH(:,t on the RHS of the mean equation is equal to theta_(t-1 (t-1 3 TH(:,t what we end up with is a random draw for theta(t-1 conditional on knowning theta in the next period 33 / 38
34 Why go forward & backward? The Kalman filter gives us E ( θ t Y t, V and VAR(θ t Y t, V With this information, we can also obtain draws for θ t ( However, we need draws from f θ T Y T, V not from ( f θ T Y t, V. The analysis above showed how to get draws ( from f θ T Y T, V recursively by going backward. 34 / 38
35 Relation to Kalman Smoother The Kalman smoother also goes backwards and resembles the procedure here. However, there is a difference. The Kalman smoother computes the mean and variance for f ( θ t Y T, V We need the mean and variance for f ( θ t θ t+1, Y T, V Since ( f θ t θ t+1, Y T, V = f ( θ t θ t+1, Y t, V, we can calculate these from Kalman filter without using the Kalman smoother 35 / 38
36 Gibbs part II: Posterior of V given theta Next step is to draw from the posterior given Y T and θ T, that is get a draw from p (V Y T, θ T The posterior combines the prior and information from the data = in each Gibbs iteration the prior is the same but the data set (i.e., θ T is different 36 / 38
37 Gibbs part II: Posterior of V given theta Given Y T and θ T, we can calcluate ε t and ν t. Both have mean zero and a Normal distribution Thus ( p V Y T, θ T = IW(V1 1, T 1 T 1 = T 0 + T V 1 = V + V T V T = T ( εt v t=1 t (ε t v t!!! Note that V, V T, &V 1 are like a sum of squares, whereas V (and R&Q are like a sum of squares divided by number of observations (same notation as in CS 37 / 38
38 References Casella, G., and E.I. George, 1992, Explaining the Gibbs Sampler, The American Statistician, Cogley, T. and T.J. Sargent, 2001, Evolving Post-World War II U.S. Inflation Dynamics, NBER Macroeconomics Annual 2001, volume 16, De Wind, J, 2014, Accounting for time-varying and nonlinear relationships in macroeconomic models, dissertation, University of Amsterdam. Available at Lancaster, T., 2004, An introduction to modern Bayesian econometrics, Blackwell Publishing. very nice textbook covering lots of stuff in an understandable way 38 / 38
Dynamic 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 informationGMM, HAC estimators, & Standard Errors for Business Cycle Statistics
GMM, HAC estimators, & Standard Errors for Business Cycle Statistics Wouter J. Den Haan London School of Economics c Wouter J. Den Haan Overview Generic GMM problem Estimation Heteroskedastic and Autocorrelation
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 informationGibbs Sampling in Linear Models #2
Gibbs Sampling in Linear Models #2 Econ 690 Purdue University Outline 1 Linear Regression Model with a Changepoint Example with Temperature Data 2 The Seemingly Unrelated Regressions Model 3 Gibbs sampling
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 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 informationGibbs Sampling in Latent Variable Models #1
Gibbs Sampling in Latent Variable Models #1 Econ 690 Purdue University Outline 1 Data augmentation 2 Probit Model Probit Application A Panel Probit Panel Probit 3 The Tobit Model Example: Female Labor
More informationECO 513 Fall 2008 C.Sims KALMAN FILTER. s t = As t 1 + ε t Measurement equation : y t = Hs t + ν t. u t = r t. u 0 0 t 1 + y t = [ H I ] u t.
ECO 513 Fall 2008 C.Sims KALMAN FILTER Model in the form 1. THE KALMAN FILTER Plant equation : s t = As t 1 + ε t Measurement equation : y t = Hs t + ν t. Var(ε t ) = Ω, Var(ν t ) = Ξ. ε t ν t and (ε t,
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 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 informationA Fully Nonparametric Modeling Approach to. BNP Binary Regression
A Fully Nonparametric Modeling Approach to Binary Regression Maria Department of Applied Mathematics and Statistics University of California, Santa Cruz SBIES, April 27-28, 2012 Outline 1 2 3 Simulation
More informationEcon 423 Lecture Notes: Additional Topics in Time Series 1
Econ 423 Lecture Notes: Additional Topics in Time Series 1 John C. Chao April 25, 2017 1 These notes are based in large part on Chapter 16 of Stock and Watson (2011). They are for instructional purposes
More informationThe Kalman filter, Nonlinear filtering, and Markov Chain Monte Carlo
NBER Summer Institute Minicourse What s New in Econometrics: Time Series Lecture 5 July 5, 2008 The Kalman filter, Nonlinear filtering, and Markov Chain Monte Carlo Lecture 5, July 2, 2008 Outline. Models
More informationMotivation Scale Mixutres of Normals Finite Gaussian Mixtures Skew-Normal Models. Mixture Models. Econ 690. Purdue University
Econ 690 Purdue University In virtually all of the previous lectures, our models have made use of normality assumptions. From a computational point of view, the reason for this assumption is clear: combined
More informationStatistics & Data Sciences: First Year Prelim Exam May 2018
Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book
More informationLearning in Macroeconomic Models
Learning in Macroeconomic Models Wouter J. Den Haan London School of Economics c by Wouter J. Den Haan Overview A bit of history of economic thought How expectations are formed can matter in the long run
More informationParameterized Expectations Algorithm
Parameterized Expectations Algorithm Wouter J. Den Haan London School of Economics c by Wouter J. Den Haan Overview Two PEA algorithms Explaining stochastic simulations PEA Advantages and disadvantages
More informationSteven L. Scott. Presented by Ahmet Engin Ural
Steven L. Scott Presented by Ahmet Engin Ural Overview of HMM Evaluating likelihoods The Likelihood Recursion The Forward-Backward Recursion Sampling HMM DG and FB samplers Autocovariance of samplers Some
More informationLatent Variable Models for Binary Data. Suppose that for a given vector of explanatory variables x, the latent
Latent Variable Models for Binary Data Suppose that for a given vector of explanatory variables x, the latent variable, U, has a continuous cumulative distribution function F (u; x) and that the binary
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationSwitching Regime Estimation
Switching Regime Estimation Series de Tiempo BIrkbeck March 2013 Martin Sola (FE) Markov Switching models 01/13 1 / 52 The economy (the time series) often behaves very different in periods such as booms
More informationOnline Appendix for: A Bounded Model of Time Variation in Trend Inflation, NAIRU and the Phillips Curve
Online Appendix for: A Bounded Model of Time Variation in Trend Inflation, NAIRU and the Phillips Curve Joshua CC Chan Australian National University Gary Koop University of Strathclyde Simon M Potter
More informationTime-Varying Parameters
Kalman Filter and state-space models: time-varying parameter models; models with unobservable variables; basic tool: Kalman filter; implementation is task-specific. y t = x t β t + e t (1) β t = µ + Fβ
More informationMID-TERM EXAM ANSWERS. p t + δ t = Rp t 1 + η t (1.1)
ECO 513 Fall 2005 C.Sims MID-TERM EXAM ANSWERS (1) Suppose a stock price p t and the stock dividend δ t satisfy these equations: p t + δ t = Rp t 1 + η t (1.1) δ t = γδ t 1 + φp t 1 + ε t, (1.2) where
More informationWorking Paper Series
Discussion of Cogley and Sargent s Drifts and Volatilities: Monetary Policy and Outcomes in the Post WWII U.S. Marco Del Negro Working Paper 23-26 October 23 Working Paper Series Federal Reserve Bank of
More informationA Bayesian Mixture Model with Application to Typhoon Rainfall Predictions in Taipei, Taiwan 1
Int. J. Contemp. Math. Sci., Vol. 2, 2007, no. 13, 639-648 A Bayesian Mixture Model with Application to Typhoon Rainfall Predictions in Taipei, Taiwan 1 Tsai-Hung Fan Graduate Institute of Statistics National
More informationECO 513 Fall 2009 C. Sims HIDDEN MARKOV CHAIN MODELS
ECO 513 Fall 2009 C. Sims HIDDEN MARKOV CHAIN MODELS 1. THE CLASS OF MODELS y t {y s, s < t} p(y t θ t, {y s, s < t}) θ t = θ(s t ) P[S t = i S t 1 = j] = h ij. 2. WHAT S HANDY ABOUT IT Evaluating the
More informationState-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Financial Econometrics / 49
State-space Model Eduardo Rossi University of Pavia November 2013 Rossi State-space Model Financial Econometrics - 2013 1 / 49 Outline 1 Introduction 2 The Kalman filter 3 Forecast errors 4 State smoothing
More informationMarginal Specifications and a Gaussian Copula Estimation
Marginal Specifications and a Gaussian Copula Estimation Kazim Azam Abstract Multivariate analysis involving random variables of different type like count, continuous or mixture of both is frequently required
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods Tomas McKelvey and Lennart Svensson Signal Processing Group Department of Signals and Systems Chalmers University of Technology, Sweden November 26, 2012 Today s learning
More information1. Introduction. Hang Qian 1 Iowa State University
Users Guide to the VARDAS Package Hang Qian 1 Iowa State University 1. Introduction The Vector Autoregression (VAR) model is widely used in macroeconomics. However, macroeconomic data are not always observed
More informationAppendix for Monetary Policy Shifts and the Term Structure
Appendix for Monetary Policy Shifts and the Term Structure Andrew Ang Columbia University and NBER Jean Boivin Bank of Canada, CIRPÉE, CIRANO and NBER Sen Dong Ortus Capital Management Rudy Loo-Kung Columbia
More informationin Trend University Gary Koop # 590
A Bounded Model of Time Variation in Trend Inflation, NAIRU and the Phillips Curve Joshua C.C. Chan Research School of Economics, Australian National University Gary Koop University of Strathclyde Simon
More informationForecasting with ARMA
Forecasting with ARMA Eduardo Rossi University of Pavia October 2013 Rossi Forecasting Financial Econometrics - 2013 1 / 32 Mean Squared Error Linear Projection Forecast of Y t+1 based on a set of variables
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 informationBayesian IV: the normal case with multiple endogenous variables
Baesian IV: the normal case with multiple endogenous variables Timoth Cogle and Richard Startz revised Januar 05 Abstract We set out a Gibbs sampler for the linear instrumental-variable model with normal
More informationIntroduction to Numerical Methods
Introduction to Numerical Methods Wouter J. Den Haan London School of Economics c by Wouter J. Den Haan "D", "S", & "GE" Dynamic Stochastic General Equilibrium What is missing in the abbreviation? DSGE
More informationA Mixed Frequency Steady-State Bayesian Vector Autoregression: Forecasting the Macroeconomy
A Mixed Frequency Steady-State Bayesian Vector Autoregression: Forecasting the Macroeconomy By Måns Unosson Department of Statistics Uppsala University Supervisor: Yukai Yang 2016 Abstract This thesis
More information9. Multivariate Linear Time Series (II). MA6622, Ernesto Mordecki, CityU, HK, 2006.
9. Multivariate Linear Time Series (II). MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Introduction to Time Series and Forecasting. P.J. Brockwell and R. A. Davis, Springer Texts
More informationBayesian nonparametric estimation of finite population quantities in absence of design information on nonsampled units
Bayesian nonparametric estimation of finite population quantities in absence of design information on nonsampled units Sahar Z Zangeneh Robert W. Keener Roderick J.A. Little Abstract In Probability proportional
More informationThe Bayesian Approach to Multi-equation Econometric Model Estimation
Journal of Statistical and Econometric Methods, vol.3, no.1, 2014, 85-96 ISSN: 2241-0384 (print), 2241-0376 (online) Scienpress Ltd, 2014 The Bayesian Approach to Multi-equation Econometric Model Estimation
More informationBayesian spatial hierarchical modeling for temperature extremes
Bayesian spatial hierarchical modeling for temperature extremes Indriati Bisono Dr. Andrew Robinson Dr. Aloke Phatak Mathematics and Statistics Department The University of Melbourne Maths, Informatics
More informationResearch 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 informationReducing Dimensions in a Large TVP-VAR
Reducing Dimensions in a Large TVP-VAR Eric Eisenstat University of Queensland Joshua C.C. Chan University of Technology Sydney Rodney W. Strachan University of Queensland March, 8 Corresponding author:
More informationThe Effects of Monetary Policy on Stock Market Bubbles: Some Evidence
The Effects of Monetary Policy on Stock Market Bubbles: Some Evidence Jordi Gali Luca Gambetti ONLINE APPENDIX The appendix describes the estimation of the time-varying coefficients VAR model. The model
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 informationGibbs Sampling in Endogenous Variables Models
Gibbs Sampling in Endogenous Variables Models Econ 690 Purdue University Outline 1 Motivation 2 Identification Issues 3 Posterior Simulation #1 4 Posterior Simulation #2 Motivation In this lecture we take
More informationSTAT 518 Intro Student Presentation
STAT 518 Intro Student Presentation Wen Wei Loh April 11, 2013 Title of paper Radford M. Neal [1999] Bayesian Statistics, 6: 475-501, 1999 What the paper is about Regression and Classification Flexible
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 informationProjection. Wouter J. Den Haan London School of Economics. c by Wouter J. Den Haan
Projection Wouter J. Den Haan London School of Economics c by Wouter J. Den Haan Model [ ] ct ν = E t βct+1 ν αz t+1kt+1 α 1 c t + k t+1 = z t k α t ln(z t+1 ) = ρ ln (z t ) + ε t+1 ε t+1 N(0, σ 2 ) k
More informationProf. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis
Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation
More informationEstimating Unobservable Inflation Expectations in the New Keynesian Phillips Curve
econometrics Article Estimating Unobservable Inflation Expectations in the New Keynesian Phillips Curve Francesca Rondina ID Department of Economics, University of Ottawa, 120 University Private, Ottawa,
More informationBayesian Linear Regression
Bayesian Linear Regression Sudipto Banerjee 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. September 15, 2010 1 Linear regression models: a Bayesian perspective
More informationLong-Run Covariability
Long-Run Covariability Ulrich K. Müller and Mark W. Watson Princeton University October 2016 Motivation Study the long-run covariability/relationship between economic variables great ratios, long-run Phillips
More informationEfficient Bayesian Multivariate Surface Regression
Efficient Bayesian Multivariate Surface Regression Feng Li (joint with Mattias Villani) Department of Statistics, Stockholm University October, 211 Outline of the talk 1 Flexible regression models 2 The
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 informationAssessing Structural Convergence between Romanian Economy and Euro Area: A Bayesian Approach
Vol. 3, No.3, July 2013, pp. 372 383 ISSN: 2225-8329 2013 HRMARS www.hrmars.com Assessing Structural Convergence between Romanian Economy and Euro Area: A Bayesian Approach Alexie ALUPOAIEI 1 Ana-Maria
More informationAn Introduction to Bayesian Linear Regression
An Introduction to Bayesian Linear Regression APPM 5720: Bayesian Computation Fall 2018 A SIMPLE LINEAR MODEL Suppose that we observe explanatory variables x 1, x 2,..., x n and dependent variables y 1,
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationECONOMETRIC METHODS II: TIME SERIES LECTURE NOTES ON THE KALMAN FILTER. The Kalman Filter. We will be concerned with state space systems of the form
ECONOMETRIC METHODS II: TIME SERIES LECTURE NOTES ON THE KALMAN FILTER KRISTOFFER P. NIMARK The Kalman Filter We will be concerned with state space systems of the form X t = A t X t 1 + C t u t 0.1 Z t
More informationLARGE TIME-VARYING PARAMETER VARS
WP 12-11 Gary Koop University of Strathclyde, UK The Rimini Centre for Economic Analysis (RCEA), Italy Dimitris Korobilis University of Glasgow, UK The Rimini Centre for Economic Analysis (RCEA), Italy
More informationLecture 2: Univariate Time Series
Lecture 2: Univariate Time Series Analysis: Conditional and Unconditional Densities, Stationarity, ARMA Processes Prof. Massimo Guidolin 20192 Financial Econometrics Spring/Winter 2017 Overview Motivation:
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 12: Gaussian Belief Propagation, State Space Models and Kalman Filters Guest Kalman Filter Lecture by
More informationBayesian Inference: Probit and Linear Probability Models
Utah State University DigitalCommons@USU All Graduate Plan B and other Reports Graduate Studies 5-1-2014 Bayesian Inference: Probit and Linear Probability Models Nate Rex Reasch Utah State University Follow
More informationMULTILEVEL IMPUTATION 1
MULTILEVEL IMPUTATION 1 Supplement B: MCMC Sampling Steps and Distributions for Two-Level Imputation This document gives technical details of the full conditional distributions used to draw regression
More informationBayesian linear regression
Bayesian linear regression Linear regression is the basis of most statistical modeling. The model is Y i = X T i β + ε i, where Y i is the continuous response X i = (X i1,..., X ip ) T is the corresponding
More informationForecast comparison of principal component regression and principal covariate regression
Forecast comparison of principal component regression and principal covariate regression Christiaan Heij, Patrick J.F. Groenen, Dick J. van Dijk Econometric Institute, Erasmus University Rotterdam Econometric
More informationDefault Priors and Effcient Posterior Computation in Bayesian
Default Priors and Effcient Posterior Computation in Bayesian Factor Analysis January 16, 2010 Presented by Eric Wang, Duke University Background and Motivation A Brief Review of Parameter Expansion Literature
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 informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and
More informationChapter 3 - Temporal processes
STK4150 - Intro 1 Chapter 3 - Temporal processes Odd Kolbjørnsen and Geir Storvik January 23 2017 STK4150 - Intro 2 Temporal processes Data collected over time Past, present, future, change Temporal aspect
More informationStatistical Machine Learning Lecture 8: Markov Chain Monte Carlo Sampling
1 / 27 Statistical Machine Learning Lecture 8: Markov Chain Monte Carlo Sampling Melih Kandemir Özyeğin University, İstanbul, Turkey 2 / 27 Monte Carlo Integration The big question : Evaluate E p(z) [f(z)]
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 informationLecture 4: Dynamic models
linear s Lecture 4: s Hedibert Freitas Lopes The University of Chicago Booth School of Business 5807 South Woodlawn Avenue, Chicago, IL 60637 http://faculty.chicagobooth.edu/hedibert.lopes hlopes@chicagobooth.edu
More informationEcon 623 Econometrics II Topic 2: Stationary Time Series
1 Introduction Econ 623 Econometrics II Topic 2: Stationary Time Series In the regression model we can model the error term as an autoregression AR(1) process. That is, we can use the past value of the
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Basics of System Identification Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont (UBC) EECE574 - Basics of
More informationGraphical Models for Collaborative Filtering
Graphical Models for Collaborative Filtering Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012 Sequence modeling HMM, Kalman Filter, etc.: Similarity: the same graphical model topology,
More informationProblem Set 2: Box-Jenkins methodology
Problem Set : Box-Jenkins methodology 1) For an AR1) process we have: γ0) = σ ε 1 φ σ ε γ0) = 1 φ Hence, For a MA1) process, p lim R = φ γ0) = 1 + θ )σ ε σ ε 1 = γ0) 1 + θ Therefore, p lim R = 1 1 1 +
More informationLecture 5: Spatial probit models. James P. LeSage University of Toledo Department of Economics Toledo, OH
Lecture 5: Spatial probit models James P. LeSage University of Toledo Department of Economics Toledo, OH 43606 jlesage@spatial-econometrics.com March 2004 1 A Bayesian spatial probit model with individual
More informationNoncausal Bayesian Vector Autoregression. Markku Lanne University of Helsinki and CREATES. Jani Luoto* University of Helsinki. August 2014.
Noncausal Bayesian Vector Autoregression Markku Lanne University of Helsinki and CREATES Jani Luoto* University of Helsinki August 2014 Abstract We propose a Bayesian inferential procedure for the noncausal
More informationLarge Vector Autoregressions with asymmetric priors and time varying volatilities
Large Vector Autoregressions with asymmetric priors and time varying volatilities Andrea Carriero Queen Mary, University of London a.carriero@qmul.ac.uk Todd E. Clark Federal Reserve Bank of Cleveland
More informationBayesian Econometrics - Computer section
Bayesian Econometrics - Computer section Leandro Magnusson Department of Economics Brown University Leandro Magnusson@brown.edu http://www.econ.brown.edu/students/leandro Magnusson/ April 26, 2006 Preliminary
More informationPart 6: Multivariate Normal and Linear Models
Part 6: Multivariate Normal and Linear Models 1 Multiple measurements Up until now all of our statistical models have been univariate models models for a single measurement on each member of a sample of
More informationMeasurement Errors and the Kalman Filter: A Unified Exposition
Luiss Lab of European Economics LLEE Working Document no. 45 Measurement Errors and the Kalman Filter: A Unified Exposition Salvatore Nisticò February 2007 Outputs from LLEE research in progress, as well
More informationNew Information Response Functions
New Information Response Functions Caroline JARDET (1) Banque de France Alain MONFORT (2) CNAM, CREST and Banque de France Fulvio PEGORARO (3) Banque de France and CREST First version : February, 2009.
More information(5) Multi-parameter models - Gibbs sampling. ST440/540: Applied Bayesian Analysis
Summarizing a posterior Given the data and prior the posterior is determined Summarizing the posterior gives parameter estimates, intervals, and hypothesis tests Most of these computations are integrals
More informationSequential Monte Carlo Methods (for DSGE Models)
Sequential Monte Carlo Methods (for DSGE Models) Frank Schorfheide University of Pennsylvania, PIER, CEPR, and NBER October 23, 2017 Some References These lectures use material from our joint work: Tempered
More informationFinancial Econometrics
Financial Econometrics Multivariate Time Series Analysis: VAR Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) VAR 01/13 1 / 25 Structural equations Suppose have simultaneous system for supply
More informationStat 451 Lecture Notes Markov Chain Monte Carlo. Ryan Martin UIC
Stat 451 Lecture Notes 07 12 Markov Chain Monte Carlo Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapters 8 9 in Givens & Hoeting, Chapters 25 27 in Lange 2 Updated: April 4, 2016 1 / 42 Outline
More informationAn EM algorithm for Gaussian Markov Random Fields
An EM algorithm for Gaussian Markov Random Fields Will Penny, Wellcome Department of Imaging Neuroscience, University College, London WC1N 3BG. wpenny@fil.ion.ucl.ac.uk October 28, 2002 Abstract Lavine
More informationMCMC Methods: Gibbs and Metropolis
MCMC Methods: Gibbs and Metropolis Patrick Breheny February 28 Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/30 Introduction As we have seen, the ability to sample from the posterior distribution
More informationUniversità di Pavia. Forecasting. Eduardo Rossi
Università di Pavia Forecasting Eduardo Rossi Mean Squared Error Forecast of Y t+1 based on a set of variables observed at date t, X t : Yt+1 t. The loss function MSE(Y t+1 t ) = E[Y t+1 Y t+1 t ]2 The
More informationBayesian modelling of real GDP rate in Romania
Bayesian modelling of real GDP rate in Romania Mihaela SIMIONESCU 1* 1 Romanian Academy, Institute for Economic Forecasting, Bucharest, Calea 13 Septembrie, no. 13, District 5, Bucharest 1 Abstract The
More information1 Class Organization. 2 Introduction
Time Series Analysis, Lecture 1, 2018 1 1 Class Organization Course Description Prerequisite Homework and Grading Readings and Lecture Notes Course Website: http://www.nanlifinance.org/teaching.html wechat
More informationState-Space Methods for Inferring Spike Trains from Calcium Imaging
State-Space Methods for Inferring Spike Trains from Calcium Imaging Joshua Vogelstein Johns Hopkins April 23, 2009 Joshua Vogelstein (Johns Hopkins) State-Space Calcium Imaging April 23, 2009 1 / 78 Outline
More informationLarge time varying parameter VARs for macroeconomic forecasting
Large time varying parameter VARs for macroeconomic forecasting Gianni Amisano (FRB) Domenico Giannone (NY Fed, CEPR and ECARES) Michele Lenza (ECB and ECARES) 1 Invited talk at the Cleveland Fed, May
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationTime-Varying Vector Autoregressive Models with Structural Dynamic Factors
Time-Varying Vector Autoregressive Models with Structural Dynamic Factors Paolo Gorgi, Siem Jan Koopman, Julia Schaumburg http://sjkoopman.net Vrije Universiteit Amsterdam School of Business and Economics
More informationEcon 424 Time Series Concepts
Econ 424 Time Series Concepts Eric Zivot January 20 2015 Time Series Processes Stochastic (Random) Process { 1 2 +1 } = { } = sequence of random variables indexed by time Observed time series of length
More informationSome Customers Would Rather Leave Without Saying Goodbye
Some Customers Would Rather Leave Without Saying Goodbye Eva Ascarza, Oded Netzer, and Bruce G. S. Hardie Web Appendix A: Model Estimation In this appendix we describe the hierarchical Bayesian framework
More information