Stochastic Volatility (SV) Models
|
|
- Blaze Leonard
- 6 years ago
- Views:
Transcription
1 1 Motivations Stochastic Volatility (SV) Models Jun Yu Some stylised facts about financial asset return distributions: 1. Distribution is leptokurtic 2. Volatility clustering 3. Volatility responds to return news asymmetrically 4. Returnsarecrossdependent 5. Volatilities are cross dependent 6. Often a lower dimensional factor structure explains most of the correlation 7. Time varying correlations 1
2 Series : sp9303daily ACF Lag ACF for typical financial asset returns If an ARMA model has to be used, this model suggests an ARMA(0,0). The ACF for squared return suggests that ARMA-type models are NOT suitable for financial returns. 2
3 Series : sp9303daily^2 ACF Lag 3
4 2 Univariate SV Models Model 1: The basic SV model (Taylor, 1982) Suppose y t is the return for a financial asset. The basic SV model is yt = σ t ε t =exp(0.5h t )ε t,ε t N(0, 1) h t = λ + αh t 1 + v t,v t N(0,σ 2 v), where {ε t } and {v t } are uncorrelated. Rewriting: yt = σ exp(0.5h t )ε t,ε t N(0, 1) h t = αh t 1 + v t,v t N(0,σ 2 v), Volatility evolves according to an AR(1) process. Condition on I t 1,h t is not deterministic any more since v t is not realized. The assumption could be more realistic than ARCH-type models. In the SV model α captures volatility clustering. Kurtosis of y t is 3exp(σ 2 v/(1 α 2 )) which is larger than 3 as long as σ v 6= 0. So the unconditional distribution of y t has fatter tails than the normal distribution. The tail thickness is due to the presence of the second noise. E(y t y s )=0, t 6= s. E(y 2 t y 2 t i) =exp(2μ h + σ 2 h (1 + αi )), whereμ h = λ 1 α and σ2 h = We can linearize the model, xt =ln(yt 2 )=h t +ln(ε 2 t ). h t = λ + αh t 1 + v t σ2 v. 1 α 2 4
5 Model 2: The t-sv model: yt = σ t ε t =exp(0.5h t )ε t,ε t t (k) h t = λ + αh t 1 + v t,v t N(0,σ 2 v). Model 3: Asymmetric SV Model (Yu, 2005): yt = σ t ε t =exp(0.5h t )ε t,ε t N(0, 1) h t = λ + αh t 1 + v t,v t N(0,σ 2 v), (2.1) where corr( t,v t+1 )=ρ. Yu (2005) showed that E(h t+1 X t )=λ + λφ µ y 1 φ + ρσ 2 v exp σv 4 4(1 φ 2 ) + σv 2λ 2 (1 φ 2 )(1 φ) y t. This is a linear function in y t and implies that if ρ<0, and everything else is held constant, a fall in the stock price/return leads to an increase of E(ln σ 2 t+1 y t ). So this model allows volatility responds asymmetrically to good/bad news 5
6 3 Multivariate SV Models For the purpose of illustration, we focus on the bivariate multivariate SV models. Let the observed log-returns at time t be denoted by y t =(y 1t,y 2t ) 0 for t = 1,...,T. Let ² t =( 1t, 2t ) 0, η t =(η 1t,η 2t ) 0, μ = (μ 1,μ 2 ) 0, h t =(h 1,t,h 2,t ) 0, Ω t = diag(exp(h t /2)), and Φ = µ µ φ11 φ 12 1 ρ, Σ φ 21 φ = 22 ρ 1 µ, Σ η = ση 2 1 ρ η σ η1 σ η2 ρ η σ η1 σ η2 ση 2 2. Model 1 (Basic-MSV or MSV): y t = Ω t ² t, ² t N(0,I), h t+1 = μ + diag(φ 11,φ 22 )(h t μ)+η t, η t N(0, diag(σ 2 η 1,σ 2 η 2 )), with h 0 = μ. This model is equivalent to stacking two basic univariate SV models together. Clearly, this specification does not allow for correlation across the returns or across the volatilities. However, it does allow for leptokurtic return distributions and volatility clustering. 6
7 Model 2 (Constant Correlation-MSV or CC-MSV): y t = Ω t ² t, ² t N(0, Σ ), h t+1 = μ + diag(φ 11,φ 22 )(h t μ)+η t, η t N(0, diag(σ 2 η 1,σ 2 η 2 )), with h 0 = μ. In this model, the return shocks are allowed to be correlated. As a result, the returns are cross dependent. Model 3 (MSV with Granger Causality or GC-MSV): y t = Ω t ² t, ² t N(0, Σ ), h t+1 = μ + Φ(h t μ)+η t, η t N(0, diag(σ 2 η 1,σ 2 η 2 )), with h 0 = μ and φ 12 =0.Sinceφ 21 can be different from zero, the volatility of the second asset is allowed to be Granger caused by the volatility of the first asset. Consequently, both the returns and volatilities are cross dependent. However, the cross-dependence of volatilities are realized via Granger causality and volatility clustering jointly. Furthermore, when both φ 12 and φ 21 are nonzero, a bilateral Granger causality in volatility between the two assets is allowed. 7
8 Model 4 (Generalized CC-MSV or GCC-MSV): y t = Ω t ² t, ² t N(0, Σ ), h t+1 = μ + diag(φ 11,φ 22 )(h t μ)+η t, η t N(0, Σ η ), with h 0 = μ. In this model both returns and volatilities are cross dependent. Obviously, both GC-MSV and GCC-MSV can generate cross dependence in volatilities. Model 5 (Dynamic Correlation-MSV or DC-MSV): y t = Ω t ² t, ² t Ω t N(0, Σ,t ), µ 1 ρt Σ,t =, ρ t 1 h t+1 = μ + diag(φ 11,φ 22 )(h t μ)+η t, η t N(0, diag(ση 2 1,ση 2 2 )), q t+1 = ψ 0 + ψ(q t ψ 0 )+σ ρ v t,v t N(0, 1), ρ t = exp(q t) 1 exp(q t )+1. with h 0 = μ,q 0 = ψ 0. In this model, not only volatilities but also correlation coefficientsaretimevarying. Ofcourse,ρ t has to be bounded by 1 and 1 for Σ to be a well-defined correlation matrix. This constraint is achieved by using the Fisher transformation. 8
9 Model 6 (Heavy-tailed MSV or t-msv): y t = Ω t ² t, ² t t(0, Σ,ν), h t+1 = μ + diag(φ 11,φ 22 )(h t μ)+η t, η t N(0, diag(σ 2 η 1,σ 2 η 2 )), with h 0 = μ. In this model, a heavy-tailed multivariate Student t distribution for the return shock is used and hence extra excess kurtosis is allowed. 9
10 Model 7 (Additive Factor-MSV or AFactor-MSV): y t = Df t + ² t, ² t N(0, diag(σ 1,σ 2 2)) 2 f t = exp(h t /2)u t,u t N(0, 1), h t+1 = μ + φ(h t μ)+σ η η t,η t N(0, 1), with h 0 =0.Thefirst component in the return equation has a smaller number of factors which capture the information relevant to the pricing of all assets while the second one is idiosyncratic noise which captures the asset specific information. Like the univariate SV model, the AFactor-MSV model allows for excess kurtosis and volatility clustering. Clearly, it also allows for cross dependence in both returns and volatilities. Note that in this model that will be introduced below, h t represents the log-volatility of the common factor, f t. The conditional correlation coefficient between y 1t and y 2t is given by: d exp(h t ) p (exp(ht )+σ 2 1)(d 2 exp(h t )+σ 2 2) = d p (1 + σ 2 1 exp( h t ))(d 2 + σ 2 2 exp( h t )). Unless σ 1 2 = σ 2 2 =0, the correlation coefficientsaretimevaryingbutthedynamicsofthecorrelationsdependonthedynamicsofh t. Moreover, correlation is an increasing function of h t, implying that the higher the volatility of the common factor, the higher the correlation in returns. 10
11 4 Estimation of SV Models 4.1 Likelihood function The ML method is more difficult to use for the SV models than for the ARCHtype models since the likelihood function is numerically more difficult to evaluate for the SV models. Let θ =(α, λ, σ v ) be the parameters of interest in the basic SV model. We wish to estimate θ from y = {y 1,y 2,,y T }. Denote the vector of log-volatilities by h = {h 1,h 2,,h T }. The likelihood function of the parameter vector θ canbewrittenas Z L(θ y) =pdf(y; θ) = pdf(y, h; θ)dh (4.2) The integral cannot be solved analytically. It is of T dimensional. Its dimension cannot be reduced either. 11
12 4.2 QMLEwithKalmanFilter The basic SV model can be represented by a linearized version without losing any information, yt = ln(rt 2 )=h t +ln(ε 2 t )= h t + μ t h t = λ + αh t 1 + v t, (4.3) with E(μ t )=0,Var(μ t )=π 2 /2. If we approximate the distribution of μ t by a normal distribution with mean 0 and variance π 2 /2, the linearized SV model is approximated by a linear Gaussian State-Space model. We here follow the standard notations. yt = A 0 x t + H 0 ξ t + w t ξ t+1 = Fξ t + v t+1, (4.4) with A 0 = λ 1 α, x t =1, H 0 =1, ξ t = h t λ 1 α, F = α, Q = σ2 v, R = π 2 /2. 12
13 Based on the linear Gassuain State-Space approximation, Kalman filter can be applied as, Initialisation: ˆξ1 0 = 0 Σ 1 0 = σ 2 v/(1 α 2 ), (4.5) Sequential updating: ˆξt t = ˆξ t t 1 + Σ t t 1 (Σ t t 1 + π 2 /2) 1 (y t λ 1 α ˆξ t t 1 ) Σ t t = Σ t t 1 Σ t t 1 (Σ t t 1 + π 2 /2) 1, (4.6) Σ t t 1 In-sample sequential prediction: ( ˆξt+1 t = αˆξ t t 1 + α(1 + π2 Σ t+1 t = α 2 Σ t t + σ 2 v 2Σ t t 1 ) 1 (y t λ 1 α ˆξ t t 1 ), (4.7) ŷt+1 t = λ 1 α + ˆξ t+1 t E[(y t+1 ŷ t+1 t )(y t+1 ŷ t+1 t ) 0 ]=Σ t+1 t + π 2 /2, (4.8) Out-of-sample forecasting: ( ˆξT +h T = α h Ê(ξ T I T )=α h ˆξT T ŷ T +h T = λ 1 α αh ˆξT T, (4.9) Smoothing: ˆξ t T = ˆξ t t + J t [ˆξ t+1 T ˆξ t+1 t ] Σ t+1 T = Σ t t + J t ( Σ t+1 T + Σ t+1 t )Jt 0 J t = Σ t t ασ 1 t+1 t, (4.10) with t = T 1,T 2,, 1. 13
14 The quasi-likelihood is computed by ln L(α, λ, σv)= 2 1 X log(σt t 1 + π 2 /2) 1 X (y t λ ˆξ 1 α t t 1 ) Σ t t 1 + π 2 /2 The h-day ahead forecast is computed by (4.9) with the QML estimates plugged in. 14
15 4.3 Simulated Maximum Likelihood QML estimators are inefficient To do the full ML estimation, we need to calculate the likelihood function, ie, to evaluate the multi-dimensional integral numerically. One way to do this is to usingamontecarlomethod. Here we introduce an class of Monte Carlo methods via importance sampling technique. We rewrite the multi-dimensional integral by Z Z Z pdf(y, h; θ) pdf(y, h; θ) L(θ y) = pdf(y, h; θ)dh = q(h)dh = dq(h) q(h) q(h) (4.11) where q(h) is an importance density and Q(h) is an importance distribution function. The idea of the simulated ML method is to draw sample h (1),...,h (S) from q so that we can approximate L(θ y) by 1 S P S i=1 pdf(y,h (s) ;θ) q(h (s) ). 15
16 The key of the simulated ML method is to match pdf(y, h; θ) and q(h) as closely as possible while ensuring that it is easy to simulate from q. Todothat,wepropose to base the importance sampler on the Laplace approximation to pdf(y, h; θ). In particular, we choose the mean of q to be the mode of ln pdf(y, h; θ) with respect of h, and the variance of q to be the negative of inverse of the second derivative of ln pdf(y,h; θ) with respect of h evaluated at the mode. h (s) N(h, Ω 1 ) where That is h =argmaxln pdf(y,h; θ) (4.12) h and Ω = 2 ln pdf(y, h ; θ) h h 0 (4.13) 16
17 4.4 MCMC Bayesian inference is then based on the posterior distribution of the unobservables given the data. In the sequel, we will denote the probability density function of a random variable θ by p(θ). By successive conditioning, the joint prior density is p(σ 2,α,σ 2 v,h 0,h 1,...,h T )=p(σ 2,α,σ 2 v)p(h 0 σ 2 v) Q n t=1 p(h t h t 1,α,σ 2 v). We assume h 0 N(0,σ 2 v/(1 α 2 )), prior independence of the parameters σ 2,α, and σv, 2 and use the standard noninformative priors, i.e. p(σ 2 ) 1, α 1, σ 2 and p(σv 2) 1. p(h σv 2 t h t 1,α,σv 2 ) is defined through the variance equation. The likelihood p(y 1,...,y T σ 2,α,σv,h 2 0,...,h T ) is defined through the return equation and the conditional independence assumption: p(y 1,...,y T σ 2,α,σ 2 v,h 0,...,h T )= TY p(y t h t,σ 2 ). (4.14) Then, by Bayes theorem, the joint posterior distribution of the unobservables given the data is proportional to the prior times likelihood, i.e. t=1 p(σ 2,α,σv,h 2 0,...,h T y 1,...,y T ) p(σ 2 )p(α)p(τ 2 )p(h 0 σv) Q 2 T t=1 p(h t h t 1,α,σ Q v) 2 T t=1 p(y t h t,σ 2 ). (4.15) 17
Statistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 14th February 2006 Part VII Session 7: Volatility Modelling Session 7: Volatility Modelling
More informationConditional Heteroscedasticity
1 Conditional Heteroscedasticity May 30, 2010 Junhui Qian 1 Introduction ARMA(p,q) models dictate that the conditional mean of a time series depends on past observations of the time series and the past
More informationIdentifying Long-Run Risks: A Bayesian Mixed-Frequency Approach
Identifying : A Bayesian Mixed-Frequency Approach Frank Schorfheide University of Pennsylvania CEPR and NBER Dongho Song University of Pennsylvania Amir Yaron University of Pennsylvania NBER February 12,
More informationImpact of Calendar Effects. in the Volatility of Vale Shares
Journal of Finance and Investment Analysis, vol.2, no.3, 2013, 1-16 ISSN: 2241-0988 (print version), 2241-0996 (online) Scienpress Ltd, 2013 Impact of Calendar Effects in the Volatility of Vale Shares
More informationChapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29
Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 29 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting
More informationEstimation of Stochastic Volatility Models : An Approximation to the Nonlinear State Space Representation
Estimation of Stochastic Volatility Models : An Approximation to the Nonlinear State Space Representation Junji Shimada and Yoshihiko Tsukuda March, 2004 Keywords : Stochastic volatility, Nonlinear state
More informationThailand Statistician January 2016; 14(1): Contributed paper
Thailand Statistician January 016; 141: 1-14 http://statassoc.or.th Contributed paper Stochastic Volatility Model with Burr Distribution Error: Evidence from Australian Stock Returns Gopalan Nair [a] and
More informationModeling skewness and kurtosis in Stochastic Volatility Models
Modeling skewness and kurtosis in Stochastic Volatility Models Georgios Tsiotas University of Crete, Department of Economics, GR December 19, 2006 Abstract Stochastic volatility models have been seen as
More informationARCH and GARCH models
ARCH and GARCH models Fulvio Corsi SNS Pisa 5 Dic 2011 Fulvio Corsi ARCH and () GARCH models SNS Pisa 5 Dic 2011 1 / 21 Asset prices S&P 500 index from 1982 to 2009 1600 1400 1200 1000 800 600 400 200
More informationST440/550: Applied Bayesian Analysis. (5) Multi-parameter models - Summarizing the posterior
(5) Multi-parameter models - Summarizing the posterior Models with more than one parameter Thus far we have studied single-parameter models, but most analyses have several parameters For example, consider
More information12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.
12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. Autoregressive Conditional Heteroscedasticity with Estimates of Variance
More informationCourse information FN3142 Quantitative finance
Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken
More informationBivariate Birnbaum-Saunders Distribution
Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators
More informationGRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS
GRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS Patrick GAGLIARDINI and Christian GOURIÉROUX INTRODUCTION Risk measures such as Value-at-Risk (VaR) Expected
More informationBayesian Analysis of a Stochastic Volatility Model
U.U.D.M. Project Report 2009:1 Bayesian Analysis of a Stochastic Volatility Model Yu Meng Examensarbete i matematik, 30 hp Handledare och examinator: Johan Tysk Februari 2009 Department of Mathematics
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (42 pts) Answer briefly the following questions. 1. Questions
More informationFinancial Econometrics
Financial Econometrics Volatility Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) Volatility 01/13 1 / 37 Squared log returns for CRSP daily GPD (TCD) Volatility 01/13 2 / 37 Absolute value
More informationLecture Note 9 of Bus 41914, Spring Multivariate Volatility Models ChicagoBooth
Lecture Note 9 of Bus 41914, Spring 2017. Multivariate Volatility Models ChicagoBooth Reference: Chapter 7 of the textbook Estimation: use the MTS package with commands: EWMAvol, marchtest, BEKK11, dccpre,
More informationChapter 6 Forecasting Volatility using Stochastic Volatility Model
Chapter 6 Forecasting Volatility using Stochastic Volatility Model Chapter 6 Forecasting Volatility using SV Model In this chapter, the empirical performance of GARCH(1,1), GARCH-KF and SV models from
More informationDYNAMIC ECONOMETRIC MODELS Vol. 8 Nicolaus Copernicus University Toruń Mateusz Pipień Cracow University of Economics
DYNAMIC ECONOMETRIC MODELS Vol. 8 Nicolaus Copernicus University Toruń 2008 Mateusz Pipień Cracow University of Economics On the Use of the Family of Beta Distributions in Testing Tradeoff Between Risk
More informationForecasting volatility in the New Zealand stock market
Forecasting volatility in the New Zealand stock market JUN YU Department of Economics, University of Auckland, Private Bag 92019, Auckland, New Zealand. This paper evaluates the performance of nine alternative
More informationA potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples
1.3 Regime switching models A potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples (or regimes). If the dates, the
More informationChapter 7: Estimation Sections
1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:
More informationHigh-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5]
1 High-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5] High-frequency data have some unique characteristics that do not appear in lower frequencies. At this class we have: Nonsynchronous
More informationLecture 5a: ARCH Models
Lecture 5a: ARCH Models 1 2 Big Picture 1. We use ARMA model for the conditional mean 2. We use ARCH model for the conditional variance 3. ARMA and ARCH model can be used together to describe both conditional
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationBayesian Hierarchical/ Multilevel and Latent-Variable (Random-Effects) Modeling
Bayesian Hierarchical/ Multilevel and Latent-Variable (Random-Effects) Modeling 1: Formulation of Bayesian models and fitting them with MCMC in WinBUGS David Draper Department of Applied Mathematics and
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Consider
More informationRisk Management and Time Series
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Risk Management and Time Series Time series models are often employed in risk management applications. They can be used to estimate
More informationModelling Joint Distribution of Returns. Dr. Sawsan Hilal space
Modelling Joint Distribution of Returns Dr. Sawsan Hilal space Maths Department - University of Bahrain space October 2011 REWARD Asset Allocation Problem PORTFOLIO w 1 w 2 w 3 ASSET 1 ASSET 2 R 1 R 2
More informationFinancial Econometrics Jeffrey R. Russell. Midterm 2014 Suggested Solutions. TA: B. B. Deng
Financial Econometrics Jeffrey R. Russell Midterm 2014 Suggested Solutions TA: B. B. Deng Unless otherwise stated, e t is iid N(0,s 2 ) 1. (12 points) Consider the three series y1, y2, y3, and y4. Match
More informationPosterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties
Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where
More informationAmath 546/Econ 589 Univariate GARCH Models
Amath 546/Econ 589 Univariate GARCH Models Eric Zivot April 24, 2013 Lecture Outline Conditional vs. Unconditional Risk Measures Empirical regularities of asset returns Engle s ARCH model Testing for ARCH
More informationBayesian Estimation of the Markov-Switching GARCH(1,1) Model with Student-t Innovations
Bayesian Estimation of the Markov-Switching GARCH(1,1) Model with Student-t Innovations Department of Quantitative Economics, Switzerland david.ardia@unifr.ch R/Rmetrics User and Developer Workshop, Meielisalp,
More informationROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices
ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices Bachelier Finance Society Meeting Toronto 2010 Henley Business School at Reading Contact Author : d.ledermann@icmacentre.ac.uk Alexander
More informationEstimation of dynamic term structure models
Estimation of dynamic term structure models Greg Duffee Haas School of Business, UC-Berkeley Joint with Richard Stanton, Haas School Presentation at IMA Workshop, May 2004 (full paper at http://faculty.haas.berkeley.edu/duffee)
More informationBooth School of Business, University of Chicago Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay. Solutions to Midterm
Booth School of Business, University of Chicago Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay Solutions to Midterm Problem A: (30 pts) Answer briefly the following questions. Each question has
More informationModel Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16
Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint
More informationDecomposing the Volatility Structure of Inflation
Decomposing the Volatility Structure of Inflation Mengheng Li 1 Siem Jan Koopman 1,2 1 Department of Econometrics, Vrije Universiteit Amsterdam, The Netherlands 2 CREATES, Aarhus University, Denmark November
More informationModelling financial data with stochastic processes
Modelling financial data with stochastic processes Vlad Ardelean, Fabian Tinkl 01.08.2012 Chair of statistics and econometrics FAU Erlangen-Nuremberg Outline Introduction Stochastic processes Volatility
More informationBayesian analysis of GARCH and stochastic volatility: modeling leverage, jumps and heavy-tails for financial time series
Bayesian analysis of GARCH and stochastic volatility: modeling leverage, jumps and heavy-tails for financial time series Jouchi Nakajima Department of Statistical Science, Duke University, Durham 2775,
More informationIs the Potential for International Diversification Disappearing? A Dynamic Copula Approach
Is the Potential for International Diversification Disappearing? A Dynamic Copula Approach Peter Christoffersen University of Toronto Vihang Errunza McGill University Kris Jacobs University of Houston
More informationKey Moments in the Rouwenhorst Method
Key Moments in the Rouwenhorst Method Damba Lkhagvasuren Concordia University CIREQ September 14, 2012 Abstract This note characterizes the underlying structure of the autoregressive process generated
More informationChapter 7: Estimation Sections
1 / 31 : Estimation Sections 7.1 Statistical Inference Bayesian Methods: 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods: 7.5 Maximum Likelihood
More informationBayesian Inference for Stochastic Volatility Models
Bayesian Inference for Stochastic Volatility Models by Zhongxian Men A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Doctor of Philosophy in
More informationMEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL
MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL Isariya Suttakulpiboon MSc in Risk Management and Insurance Georgia State University, 30303 Atlanta, Georgia Email: suttakul.i@gmail.com,
More informationChapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 59
Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 59 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationEstimating a Dynamic Oligopolistic Game with Serially Correlated Unobserved Production Costs. SS223B-Empirical IO
Estimating a Dynamic Oligopolistic Game with Serially Correlated Unobserved Production Costs SS223B-Empirical IO Motivation There have been substantial recent developments in the empirical literature on
More informationFinancial Econometrics Notes. Kevin Sheppard University of Oxford
Financial Econometrics Notes Kevin Sheppard University of Oxford Monday 15 th January, 2018 2 This version: 22:52, Monday 15 th January, 2018 2018 Kevin Sheppard ii Contents 1 Probability, Random Variables
More informationAnalytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model
Analytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model Advances in Computational Economics and Finance Univerity of Zürich, Switzerland Matthias Thul 1 Ally Quan
More informationFinancial Risk Forecasting Chapter 9 Extreme Value Theory
Financial Risk Forecasting Chapter 9 Extreme Value Theory Jon Danielsson 2017 London School of Economics To accompany Financial Risk Forecasting www.financialriskforecasting.com Published by Wiley 2011
More informationHighly Persistent Finite-State Markov Chains with Non-Zero Skewness and Excess Kurtosis
Highly Persistent Finite-State Markov Chains with Non-Zero Skewness Excess Kurtosis Damba Lkhagvasuren Concordia University CIREQ February 1, 2018 Abstract Finite-state Markov chain approximation methods
More information1 Explaining Labor Market Volatility
Christiano Economics 416 Advanced Macroeconomics Take home midterm exam. 1 Explaining Labor Market Volatility The purpose of this question is to explore a labor market puzzle that has bedeviled business
More informationFinancial Econometrics Lecture 5: Modelling Volatility and Correlation
Financial Econometrics Lecture 5: Modelling Volatility and Correlation Dayong Zhang Research Institute of Economics and Management Autumn, 2011 Learning Outcomes Discuss the special features of financial
More informationBayesian Linear Model: Gory Details
Bayesian Linear Model: Gory Details Pubh7440 Notes By Sudipto Banerjee Let y y i ] n i be an n vector of independent observations on a dependent variable (or response) from n experimental units. Associated
More informationAnalysis of the Bitcoin Exchange Using Particle MCMC Methods
Analysis of the Bitcoin Exchange Using Particle MCMC Methods by Michael Johnson M.Sc., University of British Columbia, 2013 B.Sc., University of Winnipeg, 2011 Project Submitted in Partial Fulfillment
More informationExtracting Information from the Markets: A Bayesian Approach
Extracting Information from the Markets: A Bayesian Approach Daniel Waggoner The Federal Reserve Bank of Atlanta Florida State University, February 29, 2008 Disclaimer: The views expressed are the author
More informationStochastic Volatility Models. Hedibert Freitas Lopes
Stochastic Volatility Models Hedibert Freitas Lopes SV-AR(1) model Nonlinear dynamic model Normal approximation R package stochvol Other SV models STAR-SVAR(1) model MSSV-SVAR(1) model Volume-volatility
More informationDiscussion Paper No. DP 07/05
SCHOOL OF ACCOUNTING, FINANCE AND MANAGEMENT Essex Finance Centre A Stochastic Variance Factor Model for Large Datasets and an Application to S&P data A. Cipollini University of Essex G. Kapetanios Queen
More informationModelling Returns: the CER and the CAPM
Modelling Returns: the CER and the CAPM Carlo Favero Favero () Modelling Returns: the CER and the CAPM 1 / 20 Econometric Modelling of Financial Returns Financial data are mostly observational data: they
More informationLecture Note of Bus 41202, Spring 2008: More Volatility Models. Mr. Ruey Tsay
Lecture Note of Bus 41202, Spring 2008: More Volatility Models. Mr. Ruey Tsay The EGARCH model Asymmetry in responses to + & returns: g(ɛ t ) = θɛ t + γ[ ɛ t E( ɛ t )], with E[g(ɛ t )] = 0. To see asymmetry
More informationMultivariate time series models for asset prices
Multivariate time series models for asset prices Christian M. Hafner 1 and Hans Manner 2 1 Institut de statistique and CORE, Université catholique de Louvain, Voie du Roman Pays 20, B-1348 Louvain-la-Neuve,
More informationModeling Co-movements and Tail Dependency in the International Stock Market via Copulae
Modeling Co-movements and Tail Dependency in the International Stock Market via Copulae Katja Ignatieva, Eckhard Platen Bachelier Finance Society World Congress 22-26 June 2010, Toronto K. Ignatieva, E.
More informationSTOR Lecture 15. Jointly distributed Random Variables - III
STOR 435.001 Lecture 15 Jointly distributed Random Variables - III Jan Hannig UNC Chapel Hill 1 / 17 Before we dive in Contents of this lecture 1. Conditional pmf/pdf: definition and simple properties.
More informationMultilevel quasi-monte Carlo path simulation
Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More informationIdiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective
Idiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective Alisdair McKay Boston University June 2013 Microeconomic evidence on insurance - Consumption responds to idiosyncratic
More informationMarket Risk Analysis Volume II. Practical Financial Econometrics
Market Risk Analysis Volume II Practical Financial Econometrics Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume II xiii xvii xx xxii xxvi
More informationdiscussion Papers Some Flexible Parametric Models for Partially Adaptive Estimators of Econometric Models
discussion Papers Discussion Paper 2007-13 March 26, 2007 Some Flexible Parametric Models for Partially Adaptive Estimators of Econometric Models Christian B. Hansen Graduate School of Business at the
More informationStatistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography
Statistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography Aku Seppänen Inverse Problems Group Department of Applied Physics University of Eastern Finland
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 5 Sequential Monte Carlo methods I January
More information(5) Multi-parameter models - Summarizing the posterior
(5) Multi-parameter models - Summarizing the posterior Spring, 2017 Models with more than one parameter Thus far we have studied single-parameter models, but most analyses have several parameters For example,
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationTime-varying Combinations of Bayesian Dynamic Models and Equity Momentum Strategies
TI 2016-099/III Tinbergen Institute Discussion Paper Time-varying Combinations of Bayesian Dynamic Models and Equity Momentum Strategies Nalan Basturk 1 Stefano Grassi 2 Lennart Hoogerheide 3,5 Herman
More informationNon-informative Priors Multiparameter Models
Non-informative Priors Multiparameter Models Statistics 220 Spring 2005 Copyright c 2005 by Mark E. Irwin Prior Types Informative vs Non-informative There has been a desire for a prior distributions that
More information2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises
96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with
More informationCopulas? What copulas? R. Chicheportiche & J.P. Bouchaud, CFM
Copulas? What copulas? R. Chicheportiche & J.P. Bouchaud, CFM Multivariate linear correlations Standard tool in risk management/portfolio optimisation: the covariance matrix R ij = r i r j Find the portfolio
More informationSome Simple Stochastic Models for Analyzing Investment Guarantees p. 1/36
Some Simple Stochastic Models for Analyzing Investment Guarantees Wai-Sum Chan Department of Statistics & Actuarial Science The University of Hong Kong Some Simple Stochastic Models for Analyzing Investment
More informationStochastic Volatility (SV) Models Lecture 9. Morettin & Toloi, 2006, Section 14.6 Tsay, 2010, Section 3.12 Tsay, 2013, Section 4.
Stochastic Volatility (SV) Models Lecture 9 Morettin & Toloi, 2006, Section 14.6 Tsay, 2010, Section 3.12 Tsay, 2013, Section 4.13 Stochastic volatility model The canonical stochastic volatility model
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationY t )+υ t. +φ ( Y t. Y t ) Y t. α ( r t. + ρ +θ π ( π t. + ρ
Macroeconomics ECON 2204 Prof. Murphy Problem Set 6 Answers Chapter 15 #1, 3, 4, 6, 7, 8, and 9 (on pages 462-63) 1. The five equations that make up the dynamic aggregate demand aggregate supply model
More informationProperties of the estimated five-factor model
Informationin(andnotin)thetermstructure Appendix. Additional results Greg Duffee Johns Hopkins This draft: October 8, Properties of the estimated five-factor model No stationary term structure model is
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam.
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (32 pts) Answer briefly the following questions. 1. Suppose
More informationFinancial Econometrics and Volatility Models Stochastic Volatility
Financial Econometrics and Volatility Models Stochastic Volatility Eric Zivot April 26, 2010 Outline Stochastic Volatility and Stylized Facts for Returns Log-Normal Stochastic Volatility (SV) Model SV
More informationAsymmetric Information: Walrasian Equilibria, and Rational Expectations Equilibria
Asymmetric Information: Walrasian Equilibria and Rational Expectations Equilibria 1 Basic Setup Two periods: 0 and 1 One riskless asset with interest rate r One risky asset which pays a normally distributed
More informationWeight Smoothing with Laplace Prior and Its Application in GLM Model
Weight Smoothing with Laplace Prior and Its Application in GLM Model Xi Xia 1 Michael Elliott 1,2 1 Department of Biostatistics, 2 Survey Methodology Program, University of Michigan National Cancer Institute
More informationAnalyzing Oil Futures with a Dynamic Nelson-Siegel Model
Analyzing Oil Futures with a Dynamic Nelson-Siegel Model NIELS STRANGE HANSEN & ASGER LUNDE DEPARTMENT OF ECONOMICS AND BUSINESS, BUSINESS AND SOCIAL SCIENCES, AARHUS UNIVERSITY AND CENTER FOR RESEARCH
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationDiscussion The Changing Relationship Between Commodity Prices and Prices of Other Assets with Global Market Integration by Barbara Rossi
Discussion The Changing Relationship Between Commodity Prices and Prices of Other Assets with Global Market Integration by Barbara Rossi Domenico Giannone Université libre de Bruxelles, ECARES and CEPR
More informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More informationFinancial Time Series and Their Characterictics
Financial Time Series and Their Characterictics Mei-Yuan Chen Department of Finance National Chung Hsing University Feb. 22, 2013 Contents 1 Introduction 1 1.1 Asset Returns..............................
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Describe
More informationObjective Bayesian Analysis for Heteroscedastic Regression
Analysis for Heteroscedastic Regression & Esther Salazar Universidade Federal do Rio de Janeiro Colóquio Inter-institucional: Modelos Estocásticos e Aplicações 2009 Collaborators: Marco Ferreira and Thais
More informationMultistep prediction error decomposition in DSGE models: estimation and forecast performance
Multistep prediction error decomposition in DSGE models: estimation and forecast performance George Kapetanios Simon Price Kings College, University of London Essex Business School Konstantinos Theodoridis
More informationStatistical Methods in Financial Risk Management
Statistical Methods in Financial Risk Management Lecture 1: Mapping Risks to Risk Factors Alexander J. McNeil Maxwell Institute of Mathematical Sciences Heriot-Watt University Edinburgh 2nd Workshop on
More informationLecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.
Lecture 7 Overture to continuous models Before rigorously deriving the acclaimed Black-Scholes pricing formula for the value of a European option, we developed a substantial body of material, in continuous
More informationExercises on the New-Keynesian Model
Advanced Macroeconomics II Professor Lorenza Rossi/Jordi Gali T.A. Daniël van Schoot, daniel.vanschoot@upf.edu Exercises on the New-Keynesian Model Schedule: 28th of May (seminar 4): Exercises 1, 2 and
More informationTechnical Appendix: Policy Uncertainty and Aggregate Fluctuations.
Technical Appendix: Policy Uncertainty and Aggregate Fluctuations. Haroon Mumtaz Paolo Surico July 18, 2017 1 The Gibbs sampling algorithm Prior Distributions and starting values Consider the model to
More informationExtend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty
Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for
More information