HPCFinance research project 8
|
|
- Dayna Riley
- 6 years ago
- Views:
Transcription
1 HPCFinance research projec 8 Financial models, volailiy risk, and Bayesian algorihms Hanxue Yang Tampere Universiy of Technology March 14, 2016
2 Research projec 8 12/ /2015, Tampere Universiy of Technology, Finland Supervisor: Professor Juho Kanniainen Projec opics: - Realisic financial models o capure he volailiy risk and jump risk. - Opion valuaion under advanced models. - Bayesian esimaion mehods for models wih differen jump srucures.
3 Resuls and secondmens Deliverable 3.5: Opion Pricing Performance of Differen Opion Pricing Models Docoral hesis: Markov Chain Mone Carlo Esimaion of Sochasic Volailiy Models wih Finie and Infinie Aciviy Lévy Jumps: Evidence for Efficien Models and Algorihms Publicaions: - Kanniainen, J., B. Lin, and H. Yang (2014), Esimaing and Using GARCH Models wih VIX Daa for Opion Valuaion, Journal of Banking and Finance, 43, Marino, L., H. Yang, D. Luengo, J. Kanniainen, J. Corander (2015), A Fas Universal Self-uned Sampler wihin Gibbs Sampling, Digial Signal Processing, 47, Yang, H. and J. Kanniainen (2016), Jump and Volailiy Dynamics for he S&P 500: Evidence for Infinie-Aciviy Jumps wih Non-Affine Volailiy Dynamics from Sock and Opion Marke, forhcoming in Review of Finance Secondmens: - Maasrich Universiy, Bgaor,
4 Robus financial models o capure he volailiy risk and jump risk We sudy sochasic volailiy models wih finie/infinie-aciviy Lévy jumps in index reurns, or in boh reurns and spo variance, and wih affine/non-affine variance dynamics and compare heir empirical performance in erms of goodness of fi and opion pricing errors.
5 Models Under he risk-neural measure Q, he log reurns and spo variance of he S&P 500 index follow dy = (r 12 ) v + φ QJ ( i) d + v dw Y, + dj Y,, dv = κ Q (θ Q v )d + αv β (ρdw Y, + ) 1 ρ 2 dw v, + dj v,. Benchmark model: affine Heson model. Finie-aciviy Lévy jumps: compound Poisson process. Infinie-aciviy Lévy jumps: Variance Gamma and Normal Inverse Gaussian processes. Non-affine variance specificaion: β = 1.
6 VIX index pricing formula According o he model assumpions and he generic payoff expansion resul of Carr e al. (2001), he heoreical value of he squared VIX follows: VIX 2 T 10 4 = ( [ 2 S log T F 0 (T ) EQ log S ] ) T S 0 (1) = 2 T r sds 2 ( T r sds + T (M J + φ Q J T 0 T [ T ] ) 0 2 EQ v sds, 0 = 2(M J + φ Q J ( i)) e κq T v 0 T κ Q + θ P κ P + Mv T κ Q 1 e κq T κ Q. T is he mauriy of he VIX index, i.e. 22/252 using rading day coun convenion. M J is he expeced mean size of he reurn jumps under Q. For he Variance Gamma jump, M J = γ Q, for he Normal Inverse Gaussian jump, M J = γ Q /ν Q, and for he compound Poisson jump, M J = λ(µ Q y + ρ J µ v). M v is he expeced mean size of he variance jumps under Q. We assume ha he heoreical value of he VIX is relaed o is marke value as follows: where VIX Marke,T = VIX Model,T eε, ε +1 = ρ εε + σ εɛ ε,+1, ɛ ε,+1 N(0, 1).
7 Mehods Daa: S&P 500 index reurn and 30-day VIX from January 1996 o December Esimaion mehod: Markov Chain Mone Carlo. According o he Bayes rule, he join poserior of laen variables and parameers is p(θ, v, J vix, Y) p(vix, Y, v, J, Θ) where p(θ) is he produc of priors of parameers and = p(vix Y, v, J, Θ)p(Y, v J, Θ)p(J Θ)p(Θ), [ ( ) ( )] vix Marke +1 vix Marke 2 log T 1 1 vix Model ρ ε log +1 vix Model p(vix Y, v, J, Θ) = exp =0 2πσε 2σε 2. T 1 p(y, v J, Θ) =0 1 αv 1/2+β exp 1 ρ 2 { ɛ2 Y,+1 2ρɛ Y,+1ɛ v,+1 + ɛ 2 v,+1 2(1 ρ 2 ) where ɛ Y, = ( Y Y 1 ( r 1 2 v 1 + φ Q J ( i) ηsv 1 ) J Y, ) / v 1 and ɛ v, = ( v v 1 κ P (θ P v 1 ) J v, ) /αv β 1 },
8 Goodness of fi and opion pricing performance Table : Deviance Informaion Crierion (DIC) and opion VRMSEs for affine (i = 0) and non-affine (i = 1) models. Sample A: Opions on Wednesdays from January 1996 o December Sample B: Opions on Thursdays from January 1996 o December Sample C: Opions on Wednesdays and Thursdays from January 2010 o December The bolded and ialic numbers are he bes and wors resuls, respecively, and he models are ranked according o heir performance in Sample A. Model Jump ype β DIC A B C NIG-SQR Normal Inverse Gaussian Heson-SQR None Baes-SQR compound Poisson VG-SQR Variance Gamma Model Jump ype β DIC A B C NIG-linear Normal Inverse Gaussian VG-linear Variance Gamma Heson-linear None Baes-linear compound Poisson
9 Discussion The inclusion of infinie-aciviy reurn jumps is criical for non-affine (linear) variance specificaion. Non-affine variance specificaions improve model robusness.
10 Convenional MCMC mehods: Gibbs sampling and Meropolis-Hasings algorihms The general idea of MCMC algorihms is o generae a Markov chain by Mone Carlo echniques, whose equilibrium disribuion is a given arge disribuion, ha is, he join poserior p(θ, X Y), where Θ represens he se of model parameers, X represens he sae variable, and Y represens he se of observaions. Suppose we would like o esimae Θ = {κ} and v = {v } in he Heson model from index reurns Y = {Y }: Gibbs sampling: M is he pre-specified number of MCMC runs, Sep 1 Iniialisaion: Assume he iniial sae: (κ (0), v (0) ), Sep 2 Recursion: For i = 1,..., M, 1. Sample v (i) from p(v (i) Y, v (i 1) +1, v(i) 1, κ(i 1) ), for = 1,..., T, 2. Sample κ (i) from p(κ (i) v (i), Y) = N(a, b), Meropolis-Hasings (MH): for = 1,..., T, 1. Sample v (i) 2. Accep v (i) α(v (i) from he proposal densiy q(v (i) wih probabiliy (, v (i 1) ) = min where π( ) is he poserior. q(v (i 1) q(v (i) v (i 1) ) = N(v (i 1), σ 2 ), v (i) )π(v (i) ) ( ) v (i 1) )π(v (i 1), 1 = min ) π(v (i) ) π(v (i 1), 1 ) The success of MCMC algorihms lies in an appropriae choice of he proposal densiy. ),
11 Advanced MCMC mehods Adapive Meropolis (AM), proposed by Haario e al. (1999), updaes he variance of a Gaussian proposal based on pas samples generaed in he MCMC chain: σ i = s cov(x (1), x (2),..., x (i 1) ) + s d ɛ. Paricle MCMC mehods: Paricle Marginal Meropolis-Hasings (PMMH), proposed by Andrieu e al. (2010), and Paricle Gibbs Ancesor Sampling (PGAS) by Lindsen e al. (2014) approximae he performance of he marginal MH or Gibbs sampling algorihm by generaing a paricle se for sae variables o consruc he proposal disribuion used in he MH algorihm.
12 Fas Universal Self-uned Sampler (FUSS) Proposed by Marino, L., H. Yang, D. Luengo, J. Kanniainen, J. Corander (2015), A Fas Universal Self-uned Sampler wihin Gibbs Sampling, Digial Signal Processing, 47, Iniialisaion: Choose a se of suppor poins, S M = {s 1,..., s M }, such ha s 1 < s 2 <... < s M. 2. Pruning: Remove suppor poins according o a pre-specified crierion, aaining a final se S m, wih m < M. 3. Consrucion: Build a proposal funcion p(x S m) given S m, using some appropriae pre-defined mechanism. 4. MCMC algorihm: Perform K seps of he MCMC mehod using p(x S m) as a proposal pdf, hus yielding a se of samples {x 1,..., x K }. V (x) W (x) s5,v(s5) (x) p(x) (x) p(x) x x (a) Example of he proposal consrucion wih 9 suppor poins, in he log-domain. (b) The corresponding unnormalised densiies p(x) = e W (x) and π(x) = e V (x). (c) The corresponding normalised densiies p(x) p(x) and π(x) π(x).
13 Mehods Simulaion sudies: we esimae he Heson, Baes, and Normal Inverse Gaussian (NIG) models using simulaed one-year log asse prices. Only he physical model dynamics is esimaed. Empirical sudies: we esimae he non-affine NIG model using he join informaion of he S&P 500 index reurns and he 30-day VIX index in
14 Seleced resuls of simulaion sudies Table: Esimaion resuls of he Heson, Baes, and NIG models wih differen algorihms. This able shows he mean squared error of spo variance v (MSE v), he log-likelihood of he parameers and laen sae variables, and he compuaion ime. The compuaion ime is he running ime for he whole algorihm, normalised wih respec o he ime required by he MH algorihm wih a chain of 100,000 runs. Heson model Algorihms Chain lengh Number of paricles MSE v Log-likelihood Time MH E AM E FUSS-RC-P E PMMH-PF E PGAS E Baes model Algorihms Chain lengh Number of paricles MSE v Log-likelihood Time MH E AM E FUSS-RC-P E PMMH-PF E PGAS E NIG model Algorihms Chain lengh Number of paricles MSE v Log-likelihood Time MH E AM E FUSS-RC-P E PMMH-PF E PGAS E
15 v 300 v 600 v 300 v 600 Seleced resuls of empirical sudies: Auocorrelaion funcion (ACF) v 400 v v 400 v v 700 v 800 v v 700 v 800 v (a) MH: ACFs of variance a seleced ime poins. (b) AM: ACFs of variance a seleced ime poins. 1 v v v 400 v 500 v 600 v v 400 v 500 v 800 v v 600 v 700 v 800 v (c) FUSS-RC-P4: ACFs of variance a seleced ime poins. (d) PGAS: ACFs of variance a seleced ime poins.
16 Discussion AM algorihm: able o significanly increase he accepance rae of he MH algorihm; less efficien when models become complex. FUSS algorihms: able o generae virually independen samples wih he fases mixing speed, sable regardless of model specificaions; compuaionally expensive. PMCMC mehods: compuaionally cheaper, efficien in mixing; vulnerable o large independen jumps.
Air Traffic Forecast Empirical Research Based on the MCMC Method
Compuer and Informaion Science; Vol. 5, No. 5; 0 ISSN 93-8989 E-ISSN 93-8997 Published by Canadian Cener of Science and Educaion Air Traffic Forecas Empirical Research Based on he MCMC Mehod Jian-bo Wang,
More informationCH Sean Han QF, NTHU, Taiwan BFS2010. (Joint work with T.-Y. Chen and W.-H. Liu)
CH Sean Han QF, NTHU, Taiwan BFS2010 (Join work wih T.-Y. Chen and W.-H. Liu) Risk Managemen in Pracice: Value a Risk (VaR) / Condiional Value a Risk (CVaR) Volailiy Esimaion: Correced Fourier Transform
More informationRecent Developments In Evolutionary Data Assimilation And Model Uncertainty Estimation For Hydrologic Forecasting Hamid Moradkhani
Feb 6-8, 208 Recen Developmens In Evoluionary Daa Assimilaion And Model Uncerainy Esimaion For Hydrologic Forecasing Hamid Moradkhani Cener for Complex Hydrosysems Research Deparmen of Civil, Consrucion
More informationTwo Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017
Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =
More informationMaximum Likelihood Parameter Estimation in State-Space Models
Maximum Likelihood Parameer Esimaion in Sae-Space Models Arnaud Douce Deparmen of Saisics, Oxford Universiy Universiy College London 4 h Ocober 212 A. Douce (UCL Maserclass Oc. 212 4 h Ocober 212 1 / 32
More informationRichard A. Davis Colorado State University Bojan Basrak Eurandom Thomas Mikosch University of Groningen
Mulivariae Regular Variaion wih Applicaion o Financial Time Series Models Richard A. Davis Colorado Sae Universiy Bojan Basrak Eurandom Thomas Mikosch Universiy of Groningen Ouline + Characerisics of some
More informationSample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen
Sample Auocorrelaions for Financial Time Series Models Richard A. Davis Colorado Sae Universiy Thomas Mikosch Universiy of Copenhagen Ouline Characerisics of some financial ime series IBM reurns NZ-USA
More informationInference of Sparse Gene Regulatory Network from RNA-Seq Time Series Data
Inference of Sparse Gene Regulaory Nework from RNA-Seq Time Series Daa Alireza Karbalayghareh and Tao Hu Texas A&M Universiy December 16, 2015 Alireza Karbalayghareh GRN Inference from RNA-Seq Time Series
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationSequential Bayesian Learning for Merton s Jump Model with Stochastic Volatility
Sequenial Bayesian Learning for Meron s Jump Model wih Sochasic Volailiy arxiv:1610.09750v1 [sa.ap] 31 Oc 2016 Eric Jacquier Boson Universiy, Quesrom School of Business Nicholas Polson The Universiy of
More informationObject tracking: Using HMMs to estimate the geographical location of fish
Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
More informationEstimation of Poses with Particle Filters
Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU
More informationAnno accademico 2006/2007. Davide Migliore
Roboica Anno accademico 2006/2007 Davide Migliore migliore@ele.polimi.i Today Eercise session: An Off-side roblem Robo Vision Task Measuring NBA layers erformance robabilisic Roboics Inroducion The Bayesian
More informationWhat Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix
Wha Ties Reurn Volailiies o Price Valuaions and Fundamenals? On-Line Appendix Alexander David Haskayne School of Business, Universiy of Calgary Piero Veronesi Universiy of Chicago Booh School of Business,
More informationCSE-473. A Gentle Introduction to Particle Filters
CSE-473 A Genle Inroducion o Paricle Filers Bayes Filers for Robo Localizaion Dieer Fo 2 Bayes Filers: Framework Given: Sream of observaions z and acion daa u: d Sensor model Pz. = { u, z2, u 1, z 1 Dynamics
More informationOptimal Investment under Dynamic Risk Constraints and Partial Information
Opimal Invesmen under Dynamic Risk Consrains and Parial Informaion Wolfgang Puschögl Johann Radon Insiue for Compuaional and Applied Mahemaics (RICAM) Ausrian Academy of Sciences www.ricam.oeaw.ac.a 2
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationGeneralized Least Squares
Generalized Leas Squares Augus 006 1 Modified Model Original assumpions: 1 Specificaion: y = Xβ + ε (1) Eε =0 3 EX 0 ε =0 4 Eεε 0 = σ I In his secion, we consider relaxing assumpion (4) Insead, assume
More informationSpeaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis
Speaker Adapaion Techniques For Coninuous Speech Using Medium and Small Adapaion Daa Ses Consaninos Boulis Ouline of he Presenaion Inroducion o he speaker adapaion problem Maximum Likelihood Sochasic Transformaions
More informationStochastic Modelling in Finance - Solutions to sheet 8
Sochasic Modelling in Finance - Soluions o shee 8 8.1 The price of a defaulable asse can be modeled as ds S = µ d + σ dw dn where µ, σ are consans, (W ) is a sandard Brownian moion and (N ) is a one jump
More informationHas the Business Cycle Changed? Evidence and Explanations. Appendix
Has he Business Ccle Changed? Evidence and Explanaions Appendix Augus 2003 James H. Sock Deparmen of Economics, Harvard Universi and he Naional Bureau of Economic Research and Mark W. Wason* Woodrow Wilson
More informationExtreme Value GARCH modelling with Bayesian Inference. Xin Zhao, Les Oxley, Carl Scarrott and Marco Reale
Ereme Value GARCH modelling wih Bayesian Inference Xin Zhao, Les Oley, Carl Scarro and Marco Reale No: 5/9 WORKING PAPER No: 5/9 Ereme Value GARCH modelling wih Bayesian inference Zhao, X., L.Oley, C.Scarro,
More informationA new flexible Weibull distribution
Communicaions for Saisical Applicaions and Mehods 2016, Vol. 23, No. 5, 399 409 hp://dx.doi.org/10.5351/csam.2016.23.5.399 Prin ISSN 2287-7843 / Online ISSN 2383-4757 A new flexible Weibull disribuion
More informationChapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull
Chaper 14 Wiener Processes and Iô s Lemma Copyrigh John C. Hull 014 1 Sochasic Processes! Describes he way in which a variable such as a sock price, exchange rae or ineres rae changes hrough ime! Incorporaes
More informationExcel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand
Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338
More informationSpeech and Language Processing
Speech and Language rocessing Lecure 4 Variaional inference and sampling Informaion and Communicaions Engineering Course Takahiro Shinozaki 08//5 Lecure lan (Shinozaki s par) I gives he firs 6 lecures
More informationLinear Gaussian State Space Models
Linear Gaussian Sae Space Models Srucural Time Series Models Level and Trend Models Basic Srucural Model (BSM Dynamic Linear Models Sae Space Model Represenaion Level, Trend, and Seasonal Models Time Varying
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationFORECASTING STOCK MARKET RETURNS VOLATILITY
Universiy of Rhode Island DigialCommons@URI Open Access Maser's Theses 014 FORECASTING STOCK MARKET RETURNS VOLATILITY Yanan Li Universiy of Rhode Island, yanan_li@my.uri.edu Follow his and addiional works
More informationBias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé
Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationZápadočeská Univerzita v Plzni, Czech Republic and Groupe ESIEE Paris, France
ADAPTIVE SIGNAL PROCESSING USING MAXIMUM ENTROPY ON THE MEAN METHOD AND MONTE CARLO ANALYSIS Pavla Holejšovsá, Ing. *), Z. Peroua, Ing. **), J.-F. Bercher, Prof. Assis. ***) Západočesá Univerzia v Plzni,
More informationRobert Kollmann. 6 September 2017
Appendix: Supplemenary maerial for Tracable Likelihood-Based Esimaion of Non- Linear DSGE Models Economics Leers (available online 6 Sepember 207) hp://dx.doi.org/0.06/j.econle.207.08.027 Rober Kollmann
More informationWATER LEVEL TRACKING WITH CONDENSATION ALGORITHM
WATER LEVEL TRACKING WITH CONDENSATION ALGORITHM Shinsuke KOBAYASHI, Shogo MURAMATSU, Hisakazu KIKUCHI, Masahiro IWAHASHI Dep. of Elecrical and Elecronic Eng., Niigaa Universiy, 8050 2-no-cho Igarashi,
More informationRegular Variation and Financial Time Series Models
Regular Variaion and Financial Time Series Models Richard A. Davis Colorado Sae Universiy www.sa.colosae.edu/~rdavis Thomas Mikosch Universiy of Copenhagen Bojan Basrak Eurandom Ouline Characerisics of
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More informationGINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad,
GINI MEAN DIFFEENCE AND EWMA CHATS Muhammad iaz, Deparmen of Saisics, Quaid-e-Azam Universiy Islamabad, Pakisan. E-Mail: riaz76qau@yahoo.com Saddam Akbar Abbasi, Deparmen of Saisics, Quaid-e-Azam Universiy
More informationSMC in Estimation of a State Space Model
SMC in Esimaion of a Sae Space Model Dong-Whan Ko Deparmen of Economics Rugers, he Sae Universiy of New Jersey December 31, 2012 Absrac I briefly summarize procedures for macroeconomic Dynamic Sochasic
More informationSequential Importance Resampling (SIR) Particle Filter
Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle
More informationThe complexity of climate model drifts
The complexiy of climae model drifs Davide Zanchein Angelo Rubino Maeregu Arisido Carlo Gaean Universiy of Venice, Dep. of Environmeal Sc., Informaics and Saisics A conribuion o PREFACE-WP10: (Saisical
More informationProbabilistic Robotics
Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae
More informationGraphical Event Models and Causal Event Models. Chris Meek Microsoft Research
Graphical Even Models and Causal Even Models Chris Meek Microsof Research Graphical Models Defines a join disribuion P X over a se of variables X = X 1,, X n A graphical model M =< G, Θ > G =< X, E > is
More informationThe Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite
American Journal of Operaions Research, 08, 8, 8-9 hp://wwwscirporg/journal/ajor ISSN Online: 60-8849 ISSN Prin: 60-8830 The Opimal Sopping Time for Selling an Asse When I Is Uncerain Wheher he Price Process
More informationSupplementary Document
Saisica Sinica (2013): Preprin 1 Supplemenary Documen for Funcional Linear Model wih Zero-value Coefficien Funcion a Sub-regions Jianhui Zhou, Nae-Yuh Wang, and Naisyin Wang Universiy of Virginia, Johns
More informationFinancial Econometrics Kalman Filter: some applications to Finance University of Evry - Master 2
Financial Economerics Kalman Filer: some applicaions o Finance Universiy of Evry - Maser 2 Eric Bouyé January 27, 2009 Conens 1 Sae-space models 2 2 The Scalar Kalman Filer 2 21 Presenaion 2 22 Summary
More informationDEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND
DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND Asymmery and Leverage in Condiional Volailiy Models Michael McAleer WORKING PAPER
More informationA Bayesian Approach to Spectral Analysis
Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2
More informationEffects of Finite-Sample and Realized Kernels on Standardized Returns on the Tokyo Stock Exchange
Effecs of Finie-Sample and Realized Kernels on Sandardized Reurns on he Toyo Soc Exchange The Third Inernaional Conference High-Frequency Daa Analysis in Financial Mares /7/0 Tesuya Taaishi Hiroshima Universiy
More informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
More informationFiltering Turbulent Signals Using Gaussian and non-gaussian Filters with Model Error
Filering Turbulen Signals Using Gaussian and non-gaussian Filers wih Model Error June 3, 3 Nan Chen Cener for Amosphere Ocean Science (CAOS) Couran Insiue of Sciences New York Universiy / I. Ouline Use
More informationReferences are appeared in the last slide. Last update: (1393/08/19)
SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be
More information6. Stochastic calculus with jump processes
A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio
More informationParameter Estimation for State. Space Models using Sequential. Monte Carlo Algorithms
Parameer Esimaion for Sae Space Models using Sequenial Mone Carlo Algorihms Chrisopher John Nemeh, B.Sc. (Hons.), M.Sc., M.Res. Submied for he degree of Docor of Philosophy a Lancaser Universiy. November
More informationStructural Break Detection in Time Series Models
Srucural Break Deecion in Time Series Models Richard A. Davis Thomas Lee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures) This research suppored in par by an IBM faculy
More informationOption pricing and implied volatilities in a 2-hypergeometric stochastic volatility model
Opion pricing and implied volailiies in a 2-hypergeomeric sochasic volailiy model Nicolas Privaul Qihao She Division of Mahemaical Sciences School of Physical and Mahemaical Sciences Nanyang Technological
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More informationStylized Facts and Discrete Stochastic Volatility Models
The Academy of Economic Sudies The Faculy of Finance, Insurance, Banking and Sock Exchange Docoral School of Finance and Banking Disseraion Paper Sylized Facs and Discree Sochasic Volailiy Models MSc suden
More informationIII. Module 3. Empirical and Theoretical Techniques
III. Module 3. Empirical and Theoreical Techniques Applied Saisical Techniques 3. Auocorrelaion Correcions Persisence affecs sandard errors. The radiional response is o rea he auocorrelaion as a echnical
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationNonparametric Option-Implied Volatility
Nonparameric Opion-Implied Volailiy Vikor odorov May 11, 2017 Absrac We propose a nonparameric esimaor of spo volailiy from shor-daed opions. he esimaor is based on forming porfolios of opions wih differen
More informationFinancial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS
Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If
More informationST2352. Stochastic Processes constructed via Conditional Simulation. 09/02/2014 ST2352 Week 4 1
ST35 Sochasic Processes consruced via Condiional Simulaion 09/0/014 ST35 Week 4 1 Sochasic Processes consruced via Condiional Simulaion Markov Processes Simulaing Random Tex Google Sugges n grams Random
More informationUnderstanding the asymptotic behaviour of empirical Bayes methods
Undersanding he asympoic behaviour of empirical Bayes mehods Boond Szabo, Aad van der Vaar and Harry van Zanen EURANDOM, 11.10.2011. Conens 2/20 Moivaion Nonparameric Bayesian saisics Signal in Whie noise
More informationUsing the Kalman filter Extended Kalman filter
Using he Kalman filer Eended Kalman filer Doz. G. Bleser Prof. Sricker Compuer Vision: Objec and People Tracking SA- Ouline Recap: Kalman filer algorihm Using Kalman filers Eended Kalman filer algorihm
More informationLicenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A
Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In
More informationARCH IN SHORT-TERM INTEREST RATES: CASE STUDY USA
Arch in Shor-Term Ineres Raes: Case Sudy USA ARCH IN SHORT-TERM INTEREST RATES: CASE STUDY USA Adrian Ausin, Universiy of Wes Georgia Swarna (Bashu) Du, Universiy of Wes Georgia ABSTRACT We invesigae ARCH
More informationSimulating models with heterogeneous agents
Simulaing models wih heerogeneous agens Wouer J. Den Haan London School of Economics c by Wouer J. Den Haan Individual agen Subjec o employmen shocks (ε i, {0, 1}) Incomplee markes only way o save is hrough
More informationAnnouncements. Recap: Filtering. Recap: Reasoning Over Time. Example: State Representations for Robot Localization. Particle Filtering
Inroducion o Arificial Inelligence V22.0472-001 Fall 2009 Lecure 18: aricle & Kalman Filering Announcemens Final exam will be a 7pm on Wednesday December 14 h Dae of las class 1.5 hrs long I won ask anyhing
More informationStochastic Modelling of Electricity and Related Markets: Chapter 3
Sochasic Modelling of Elecriciy and Relaed Markes: Chaper 3 Fred Espen Benh, Jūraė Šalyė Benh and Seen Koekebakker Presener: Tony Ware Universiy of Calgary Ocober 14, 2009 The Schwarz model S() = S(0)
More informationDYNAMIC ECONOMETRIC MODELS Vol. 7 Nicolaus Copernicus University Toruń Piotr Fiszeder Nicolaus Copernicus University in Toruń
DYNAMIC ECONOMETRIC MODELS Vol. 7 Nicolaus Copernicus Universiy Toruń 006 Pior Fiszeder Nicolaus Copernicus Universiy in Toruń Consequences of Congruence for GARCH Modelling. Inroducion In 98 Granger formulaed
More informationBayesian estimation of correlation matrices: Application to the Multivariate GARCH model
Bayesian esimaion of correlaion marices: Applicaion o he Mulivariae GARCH model By Pierre Evarise NGUIMKEU NGUEDIA To fulfill he degree of Maser in Economics Advisor: Prof. William McCAUSLAND Deparmen
More informationTom Heskes and Onno Zoeter. Presented by Mark Buller
Tom Heskes and Onno Zoeer Presened by Mark Buller Dynamic Bayesian Neworks Direced graphical models of sochasic processes Represen hidden and observed variables wih differen dependencies Generalize Hidden
More informationIntroduction to Mobile Robotics
Inroducion o Mobile Roboics Bayes Filer Kalman Filer Wolfram Burgard Cyrill Sachniss Giorgio Grisei Maren Bennewiz Chrisian Plagemann Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel
More informationStructural Break Detection in Time Series Models
Srucural Break Deecion in Time Series Models Richard A. Davis Thomas Lee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures) This research suppored in par by: Naional
More informationAuthors. Introduction. Introduction
Auhors Hidden Applied in Agriculural Crops Classificaion Caholic Universiy of Rio de Janeiro (PUC-Rio Paula B. C. Leie Raul Q. Feiosa Gilson A. O. P. Cosa Hidden Applied in Agriculural Crops Classificaion
More informationTransform Techniques. Moment Generating Function
Transform Techniques A convenien way of finding he momens of a random variable is he momen generaing funcion (MGF). Oher ransform echniques are characerisic funcion, z-ransform, and Laplace ransform. Momen
More informationMANY FACET, COMMON LATENT TRAIT POLYTOMOUS IRT MODEL AND EM ALGORITHM. Dimitar Atanasov
Pliska Sud. Mah. Bulgar. 20 (2011), 5 12 STUDIA MATHEMATICA BULGARICA MANY FACET, COMMON LATENT TRAIT POLYTOMOUS IRT MODEL AND EM ALGORITHM Dimiar Aanasov There are many areas of assessmen where he level
More informationStructural Break Detection in Time Series Models
Srucural Break Deecion in Time Series Models Richard A. Davis Thomas Lee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures) This research suppored in par by an IBM faculy
More informationRandom Walk with Anti-Correlated Steps
Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and
More informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
More informationSequential Monte Carlo Methods in the nimble R Package arxiv: v1 [stat.co] 17 Mar 2017
Sequenial Mone Carlo Mehods in he nimble R Package arxiv:1703.06206v1 [sa.co] 17 Mar 2017 Nicholas Michaud 12, Perry de Valpine 2, Daniel Turek 3, Chrisopher J. Paciorek 1 1 Deparmen of Saisics, Universiy
More informationFigure 1. Jaw RMS target-tracker difference for a.9hz sinusoidal target.
Approaches o Bayesian Smooh Unimodal Regression George Woodworh Dec 3, 999 (Draf - please send commens o george-woodworh@uiowa.edu). Background Speech Ariculaion Daa The daa in Figure were obained by asking
More information1. Consider a pure-exchange economy with stochastic endowments. The state of the economy
Answer 4 of he following 5 quesions. 1. Consider a pure-exchange economy wih sochasic endowmens. The sae of he economy in period, 0,1,..., is he hisory of evens s ( s0, s1,..., s ). The iniial sae is given.
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationEnergy Storage Benchmark Problems
Energy Sorage Benchmark Problems Daniel F. Salas 1,3, Warren B. Powell 2,3 1 Deparmen of Chemical & Biological Engineering 2 Deparmen of Operaions Research & Financial Engineering 3 Princeon Laboraory
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More informationStatistics versus mean-field limit for Hawkes process. with Sylvain Delattre (P7)
Saisics versus mean-field limi for Hawkes process wih Sylvain Delare (P7) The model We have individuals. Z i, := number of acions of he i-h individual unil ime. Z i, jumps (is increased by 1) a rae λ i,
More informationRecent Developments in the Unit Root Problem for Moving Averages
Recen Developmens in he Uni Roo Problem for Moving Averages Richard A. Davis Colorado Sae Universiy Mei-Ching Chen Chaoyang Insiue of echnology homas Miosch Universiy of Groningen Non-inverible MA() Model
More informationOutline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests
Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange
More informationStochastic Model for Cancer Cell Growth through Single Forward Mutation
Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com
More informationPure Jump Lévy Processes and Self-decomposability in Financial Modelling
Pure Jump Lévy Processes and Self-decomposabiliy in Financial Modelling Ömer Önalan Faculy of Business Adminisraion and Economics, Marmara Universiy, Beykoz/Isanbul/ Turkey E-mail : omeronalan@marmara.edu.r
More information2017 3rd International Conference on E-commerce and Contemporary Economic Development (ECED 2017) ISBN:
7 3rd Inernaional Conference on E-commerce and Conemporary Economic Developmen (ECED 7) ISBN: 978--6595-446- Fuures Arbirage of Differen Varieies and based on he Coinegraion Which is under he Framework
More informationA Local Vector Autoregressive Framework and its Applications to Multivariate Time Series Monitoring and Forecasting
A Local Vecor Auoregressive Framework and is Applicaions o Mulivariae Time Series Monioring and Forecasing Ying Chen, Bo Li and Linlin Niu Absrac Our proposed local vecor auoregressive (LVAR) model has
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationStructural Break Detection in Time Series Models
Srucural Break Deecion in Time Series Models Richard A. Davis Thomas Lee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures) This research suppored in par by an IBM faculy
More informationA Shooting Method for A Node Generation Algorithm
A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan
More informationRobust critical values for unit root tests for series with conditional heteroscedasticity errors: An application of the simple NoVaS transformation
WORKING PAPER 01: Robus criical values for uni roo ess for series wih condiional heeroscedasiciy errors: An applicaion of he simple NoVaS ransformaion Panagiois Manalos ECONOMETRICS AND STATISTICS ISSN
More information