A Functional Central Limit Theorem for an ARMA(p, q) Process with Markov Switching

Size: px
Start display at page:

Download "A Functional Central Limit Theorem for an ARMA(p, q) Process with Markov Switching"

Transcription

1 Communications for Statistical Applications and Methods 2013, Vol 20, No 4, DOI: A Functional Central Limit Theorem for an ARMAp, q) Process with Markov Switching Oesook Lee 1,a a Department of Statistics, Ewha Womans University Abstract In this paper, we give a tractable sufficient condition for functional central limit theorem to hold in Markov switching ARMA p, q) model Keywords: Functional central limit theorem, Markov switching ARMA p, q) process, φ-mixing 1 Introduction Since the seminal work of Hamilton 1989), nonlinear time series models subject to Markov switching are widely used for modelling dynamics in different fields, especially in econometric studies Krolzig, 1997; Hamilton and Raj, 2002) Markov switching autoregressive moving averagemsarma)p, q) model, in which a hidden Markov process governs the behavior of an observable time series, exhibits structural breaks and local linearity When we consider a time series model as a data generating process, one of the important properties to show is the functional) central limit theorem Functional central limit theoremfclt) is applied for statistical inference in time series to establish the asymptotics of various statistics concerning, for example a test for stability such as CUCUM or MOSUM and unit root testing Probabilistic properties of MSARMAp, q) models have been studied in, eg, Francq and Zakoïan 2001), Yao and Attali 2000), Yang 2000), Lee 2005), and Stelzer 2009) The purpose of this paper is to find a sufficient condition under which FCLT holds for the partial sums processes of the given MSARMAp, q) models The typical approach to obtaining the FCLT for nonlinear time series models is to show that a specific dependence property such as various mixing condition, L p -NEDnear-epoch dependent), association or θ, L or ψ-weak dependence holds However, restrictive conditions such as distributional assumptions on errors and higher order moment are required in order to prove such dependence properties Ango Nze and Doukhan, 2004; De Jong and Davidson, 2000; Davidson, 2002; Dedecker et al, 2007; Doukhan and Wintenberger, 2007; Harrndorf, 1984) Our proof is based on Theorem 211 of Billingsley 1968) that is an extension of the results of Ibragimov 1962) The proof is short and relies on a second moment assumption This research was supported by Basic Science Research Program through the NRF funded by the Ministry of Education, Science and Technology2012R1A1A and 2013R1A1A ) 1 Professor, Department of Statistics, Ewha Womans University, Ewhayeodaegil 52, Seoul, Korea oslee@ewhaackr

2 340 Oesook Lee 2 Results We consider the following Markov switching autoregressive moving average MSARMA)p, q) model, where ARMA coefficients are allowed to change over time according to a Markov chain: x t = p ϕ i u t 1 )x t i + q θ j u t 1 )e t j + e t, t Z, 21) j=1 where p 1, q 0, {u t } is an irreducible aperiodic Markov chain on a finite state space E with n-step transition probability matrix P n) = p n) u,v) u,v E and the stationary distribution π, ϕ i u) and θ j u)i = 1, 2,, p, j = 1, 2,, q, u E) are constants and {e t } is a sequence of independent and identically distributed random variables with mean 0 and finite variance σ 2 We assume that {u t } and {e t } are independent Assuming without loss of generality that p = q and letting X t = x t, x t 1,, x t p+1 ) and ε t = e t, e t 1,, e t p+1 ), define p p matrix Φ t = Φu t ) = ϕ 1 u t ) ϕ 2 u t ) ϕ p 1 u t ) ϕ p u t ) and Θ t = Θu t ) = θ 1 u t ) θ 2 u t ) θ p 1 u t ) θ p u t ) Then X t = Φ t 1 X t 1 + Θ t 1 ε t 1 + ε t 22) and x t = c X t, c = 1, 0,, 0) Applying recursively the Equation 22), after m-step we have that m 1 X t = m j=1 Φ t jx t m + j=1 Φ t jφ t k + Θ t k )ε t k + m 1 j=1 Φ t jθ t m ε t m + ε t 23) assume 0 j=1 Φ t j = I) We make the assumptions A1) ρ := p α i < 1 with α i = sup u E Eϕ 2 i u 1 ) u 0 = u) ) 1 2 Note that the assumption A1) depends only on the autoregressive coefficients ϕ i u t ) of the Equation 21) Next theorem is for the existence of a strictly stationary solution to the Equation 21)

3 A Functional Central Limit Theorem for an ARMAp, q) Process with Markov Switching 341 Theorem 1 Suppose the assumption A1) holds Then there is a unique non-anticipative strictly stationary solution x t to 21) with Ex 2 t ) < Proof: Define X t = j=1 Φ t j) Φt k + Θ t k ) ε t k + ε t 24) Then according to Minkowski inequality, we have that Ec Xt ) 2) 1 2 = c Xt 2 c j=1 Φ t jφ t k + Θ t k )ε 2 t k + σ 2 Let β i := E π ϕ i u t ) + θ i u t )) 2 ) 1/2 < Apply the independence of {u t } and {e t } and Minkowski inequality and use the assumption A1) to obtain, after some tedious calculations, that c j=1 Φ t j Φ t k + Θ t k ) ε 2 t k Ck 1), k = 1, 2, 3,, where p C := C0) = σ β i, C1) = α 1 C ρc, C2) = α α 2) ρc, Cp) = α 1 Cp 1) + α 2 Cp 2) + + α p C ρc Moreover, In general, for 1 r p, m 0, Therefore, Cp + 1) = α 1 Cp) + α 2 Cp 1) + + α p C1) ρ 2 C, Cp + 2) = α 1 Cp + 1) + α 2 Cp) + + α p C2) ρ 2 C, Cmp + r) = α 1 Cmp + r 1) + α 2 Cmp + r 2) + + α p Cm 1)p + r) ρ m+1 C c j=1 Φ t jφ t k + Θ t k )ε 2 t k C ρ [ ) p ] +1 < 25) Hence we have that Ec Xt ) 2 < and xt = c Xt strictly stationary sequence satisfying 22) < ae One can easily verify that X t in 24) is a

4 342 Oesook Lee To prove the uniqueness, assume that ˆX t is any non-anticipative, strictly stationary solution to 22) with Ec ˆX t ) 2 < Obviously ˆX t satisfies 23) for every m 1 and thus we obtain that X t ˆX t = m j=1 Φ t j)x t m ˆX t m ) E x t ˆx t = E c ) Xt ˆX t = E c m j=1 Φ ) t j X t m ˆX t m p ) E x t m i+1 ˆx t m i+1 c m j=1 Φ t j) νi p xt m i+1 ˆx t m i+1 c 2 m j=1 Φ 2 t jν i, where ν i = 0,, 0, 1, 0,, 0), 0 except the i th entry which is 1 Adopt the similar method used to derive the inequality 25), then we have that c m j=1 Φ t jν i 2 Mρ [m 1)/p] 0 as m 0, where M = p E π ϕ i u t ) < Thus we have that E xt ˆx t = 0, and xt = ˆx t ae Take x t = xt to get the conclusion A stationary process x t ) t Z is called φ- mixing if φ n = sup { PE 2 E 1 ) PE 2 ) : E 1 F k, E 2 F k+n} 0 as n, where for s < t, F t s := σx s, x s+1,, x t ) Let [x] denote the largest integer not exceeding x and let D denote convergence in distribution Let B denote standard Brownian motion on [0, 1] Following theorem is our main result Theorem 2 Suppose the assumption A1) holds Then for a strictly stationary solution x t of 21) with Ex 2 t ) <, τ 2 = Varx 0 ) k< Covx 0, x t ) is convergent and as n, 1 [nξ] D x t τbξ), 0 ξ 1 n t=1 Proof: Note that v t := u t 1, e t, e t 1,, e t p ) is an aperiodic E R p+1 -valued Markov chain Since {u t } is a positive recurrent Markov chain with finite state space, {u t } is exponentially φ-mixing By independence of {u t } and {e t }, {v t } is also exponentially φ mixing with φ 1/2 n < By the above Theorem 1, and take x t = c x t,l = c l j=1 Φ t j) Φt k + Θ t k ) ε t k + c ε t, j=1 Φ t j) Φt k + Θ t k ) ε t k + c ε t

5 A Functional Central Limit Theorem for an ARMAp, q) Process with Markov Switching 343 {x t,l } is stationary for each l For some measurable functions f and f l, we can write x t = f v t, v t 1, v t 2, ), and x t,l = f l v t, v t 1, v t 2,, v t l+1 ) Then we have that Thus, x t x t,l 2 = l=1 c C k=l+1 j=1 Φ ) t j Φt k + Θ t k ) ε t k 2 c j=1 Φ 2 t j) Φt k + Θ t k ) ε t k k=l+1 ρ k=l+1 [ ) ] p +1 x t x 2 t,l C Cp 2 ρ l=1 k=l+1 l=1 i=0 = Cp 2 ρ 1 ρ) 2 < [ ) ] p +1 ρ i+l Applying Theorem 211 in Billingsley 1968) yields the conclusion Remark 1 For x t in 21), define a process Y t = Au t )Y t 1 + η t, with Y t = x t 1,, x t p, e t,, e t q+1 ) for properly defined p + q) p + q) matrix Au t ) and p + q) 1 vector η t If we take W t = u t, Y t ), then W t is a Markov chain One way to prove the FCLT is to use the Markovian structure of the model to obtain mixing properties If W t is ϕ-irreducible weak Feller and top Lyapunov exponent of Au t ) is negative, then the process Y t is geometrically ergodic and β-mixing and the FCLT for x t is obtained Lee, 2005) The stationarity and geometric ergodicity for vector valued MSARMAp, q) process with a general state space parameter chain are examined by Stelzer 2009) Remark 2 It is known that if the largest modulus of the eigenvalues λ of the set Φu), u E is less than 1, then the FCLT holds for x t Davidson, 2002, p256) Following is a simple example, which shows that λ > 1, but the FCLT holds due to Theorem 2 Example 1 Consider a MSARMA3, 2) process given by x t = 3 ϕ i u t 1 )x t i + 2 θ j u t 1 )e t j + e t, 26) j=1

6 344 Oesook Lee where {u t } is a Markov chain with state space E = {1, 2, 3} and one step transition probability matrix P is given by P = Take ϕ 1 1), ϕ 2 1), ϕ 3 1)) = 01, 03, 001), ϕ 1 2), ϕ 2 2), ϕ 3 2)) = 11, 02, 02), ϕ 1 3), ϕ 2 3), ϕ 3 3)) = 02, 003, 02), and Φu) = ϕ 1 u) ϕ 2 u) ϕ 3 u) , u E For the above example, the largest modulus of the eigenvalues λ of the set Φ1), Φ2) and Φ3) is larger than 1 but 3 sup u E Eϕ2 i u 1) u 0 = u)) 1/2 < 1 Hence Theorem 1 and 2 guarantee the existence of a strictly stationary solution x t of 26) and the FCLT for {x t } References Ango Nze, P and Doukhan, P 2004) Weak dependence: Models and applications to econometrics, Econometric Theory, 20, Billingsley, P 1968) Convergence of Probability Measures, Wiley, New York Davidson, J 2002) Establishing conditions for functional central limit theorem in nonlinear and semiparametric time series processes, Journal of Econometrics, 106, Dedecker, J, Doukhan, P, Lang, G, Leon, J R, Louhichi, S and Prieur, C 2007) Weak dependence, Examples and Applications In Springer Lecture Notes in Statistics, 190 De Jong, R M and Davidson, J 2000) The functional central limit theorem and weak convergence to stochastic integrals I: weakly dependent processes, Econometric Theory, 16, Doukhan, P and Wintenberger, O 2007) An invariance principle for weakly dependent stationary general models, Probability and Mathematical Statistics, 27, Francq, C and Zakoïan, J M 2001) Stationarity of multivariate Markov-switching ARMA models, Journal of Econometrics, 102, Hamilton, J D 1989) A new approach to the economic analysis of nonstationary time series and business cycle, Econometrica, 57, Hamilton, J D and Raj, Beds) 2002) Advances in Markov- switching Models - Applications in Business Cycle Research and Finance, Physica-Verlag, Heidelberg Harrndorf, N 1984) A functional central limit theorem for weakly dependent sequences of random variables, The Annals of Probability, 12, Ibragimov, I A 1962) Some limit theorems for stationary processes, Theory of Probability and Its Applications, 7, Krolzig, H M 1997) Markov-Switching Vector Autoregressions: Lecture Notes in Economics and Mathematical Systems, 454, Springer, Berlin

7 A Functional Central Limit Theorem for an ARMAp, q) Process with Markov Switching 345 Lee, O 2005) Probabilistic properties of a nonlinear ARMA process with Markov switching, Communications in Statistics: Theory and Methods, 34, Stelzer, R 2009) On Markov-switching ARMA processes-stationarity, existence of moments and geometric ergodicity, Econometric Theory, 29, Yang, M 2000) Some properties of vector autoregressive processes with Markov switching coefficients, Econometric Theory, 16, Yao, J and Attali, J G 2000) On stability of nonlinear AR processes with Markov switching, Advances in Applied Probability, 32, Received June 12, 2013; Revised July 8, 2013; Accepted July 8, 2013

Multivariate Markov-switching ARMA processes with regularly varying noise

Multivariate Markov-switching ARMA processes with regularly varying noise Multivariate Markov-switching ARMA processes with regularly varying noise Robert Stelzer 26th February 2007 Abstract The tail behaviour of stationary R d -valued Markov-Switching ARMA processes driven

More information

Multivariate Markov-switching ARMA processes with regularly varying noise

Multivariate Markov-switching ARMA processes with regularly varying noise Multivariate Markov-switching ARMA processes with regularly varying noise Robert Stelzer 23rd June 2006 Abstract The tail behaviour of stationary R d -valued Markov-Switching ARMA processes driven by a

More information

Continuous Time Approximations to GARCH(1, 1)-Family Models and Their Limiting Properties

Continuous Time Approximations to GARCH(1, 1)-Family Models and Their Limiting Properties Communications for Statistical Applications and Methods 214, Vol. 21, No. 4, 327 334 DOI: http://dx.doi.org/1.5351/csam.214.21.4.327 Print ISSN 2287-7843 / Online ISSN 2383-4757 Continuous Time Approximations

More information

Maddalena Cavicchioli

Maddalena Cavicchioli SPECTRAL REPRESENTATION AND AUTOCOVARIANCE STRUCTURE OF MARKOV SWITCHING DSGE MODELS Maddalena Cavicchioli Abstract. In this paper we investigate the L 2 -structure of Markov switching Dynamic Stochastic

More information

ARTICLE IN PRESS Statistics and Probability Letters ( )

ARTICLE IN PRESS Statistics and Probability Letters ( ) Statistics and Probability Letters ( ) Contents lists available at ScienceDirect Statistics and Probability Letters journal homepage: www.elsevier.com/locate/stapro The functional central limit theorem

More information

Practical conditions on Markov chains for weak convergence of tail empirical processes

Practical conditions on Markov chains for weak convergence of tail empirical processes Practical conditions on Markov chains for weak convergence of tail empirical processes Olivier Wintenberger University of Copenhagen and Paris VI Joint work with Rafa l Kulik and Philippe Soulier Toronto,

More information

KALMAN-TYPE RECURSIONS FOR TIME-VARYING ARMA MODELS AND THEIR IMPLICATION FOR LEAST SQUARES PROCEDURE ANTONY G AU T I E R (LILLE)

KALMAN-TYPE RECURSIONS FOR TIME-VARYING ARMA MODELS AND THEIR IMPLICATION FOR LEAST SQUARES PROCEDURE ANTONY G AU T I E R (LILLE) PROBABILITY AND MATHEMATICAL STATISTICS Vol 29, Fasc 1 (29), pp 169 18 KALMAN-TYPE RECURSIONS FOR TIME-VARYING ARMA MODELS AND THEIR IMPLICATION FOR LEAST SQUARES PROCEDURE BY ANTONY G AU T I E R (LILLE)

More information

Consistency of Quasi-Maximum Likelihood Estimators for the Regime-Switching GARCH Models

Consistency of Quasi-Maximum Likelihood Estimators for the Regime-Switching GARCH Models Consistency of Quasi-Maximum Likelihood Estimators for the Regime-Switching GARCH Models Yingfu Xie Research Report Centre of Biostochastics Swedish University of Report 2005:3 Agricultural Sciences ISSN

More information

GARCH( 1, 1) processes are near epoch dependent

GARCH( 1, 1) processes are near epoch dependent Economics Letters 36 (1991) 181-186 North-Holland 181 GARCH( 1, 1) processes are near epoch dependent Bruce E. Hansen 1Jnrrwsity of Rochester, Rochester, NY 14627, USA Keceived 22 October 1990 Accepted

More information

QED. Queen s Economics Department Working Paper No. 1244

QED. Queen s Economics Department Working Paper No. 1244 QED Queen s Economics Department Working Paper No. 1244 A necessary moment condition for the fractional functional central limit theorem Søren Johansen University of Copenhagen and CREATES Morten Ørregaard

More information

Asymptotic inference for a nonstationary double ar(1) model

Asymptotic inference for a nonstationary double ar(1) model Asymptotic inference for a nonstationary double ar() model By SHIQING LING and DONG LI Department of Mathematics, Hong Kong University of Science and Technology, Hong Kong maling@ust.hk malidong@ust.hk

More information

Goodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach

Goodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach Goodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach By Shiqing Ling Department of Mathematics Hong Kong University of Science and Technology Let {y t : t = 0, ±1, ±2,

More information

Establishing Stationarity of Time Series Models via Drift Criteria

Establishing Stationarity of Time Series Models via Drift Criteria Establishing Stationarity of Time Series Models via Drift Criteria Dawn B. Woodard David S. Matteson Shane G. Henderson School of Operations Research and Information Engineering Cornell University January

More information

Weak invariance principles for sums of dependent random functions

Weak invariance principles for sums of dependent random functions Available online at www.sciencedirect.com Stochastic Processes and their Applications 13 (013) 385 403 www.elsevier.com/locate/spa Weak invariance principles for sums of dependent random functions István

More information

Discrete time processes

Discrete time processes Discrete time processes Predictions are difficult. Especially about the future Mark Twain. Florian Herzog 2013 Modeling observed data When we model observed (realized) data, we encounter usually the following

More information

Research Article Least Squares Estimators for Unit Root Processes with Locally Stationary Disturbance

Research Article Least Squares Estimators for Unit Root Processes with Locally Stationary Disturbance Advances in Decision Sciences Volume, Article ID 893497, 6 pages doi:.55//893497 Research Article Least Squares Estimators for Unit Root Processes with Locally Stationary Disturbance Junichi Hirukawa and

More information

Time Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY

Time Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY PREFACE xiii 1 Difference Equations 1.1. First-Order Difference Equations 1 1.2. pth-order Difference Equations 7

More information

Module 9: Stationary Processes

Module 9: Stationary Processes Module 9: Stationary Processes Lecture 1 Stationary Processes 1 Introduction A stationary process is a stochastic process whose joint probability distribution does not change when shifted in time or space.

More information

CREATES Research Paper A necessary moment condition for the fractional functional central limit theorem

CREATES Research Paper A necessary moment condition for the fractional functional central limit theorem CREATES Research Paper 2010-70 A necessary moment condition for the fractional functional central limit theorem Søren Johansen and Morten Ørregaard Nielsen School of Economics and Management Aarhus University

More information

On the Central Limit Theorem for an ergodic Markov chain

On the Central Limit Theorem for an ergodic Markov chain Stochastic Processes and their Applications 47 ( 1993) 113-117 North-Holland 113 On the Central Limit Theorem for an ergodic Markov chain K.S. Chan Department of Statistics and Actuarial Science, The University

More information

Switching Regime Estimation

Switching 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 information

Strong Mixing Conditions. Richard C. Bradley Department of Mathematics, Indiana University, Bloomington, Indiana, USA

Strong Mixing Conditions. Richard C. Bradley Department of Mathematics, Indiana University, Bloomington, Indiana, USA Strong Mixing Conditions Richard C. Bradley Department of Mathematics, Indiana University, Bloomington, Indiana, USA There has been much research on stochastic models that have a well defined, specific

More information

Weak dependence for infinite ARCH-type bilinear models. hal , version 1-3 May 2007

Weak dependence for infinite ARCH-type bilinear models. hal , version 1-3 May 2007 Author manuscript, published in "Statistics A Journal of Theoretical and Applied Statistics 41, 1 (2007) 31-45" DOI : 10.1080/02331880601107064 Statistics: A Journal of Theoretical and Applied Statistics

More information

3. ARMA Modeling. Now: Important class of stationary processes

3. ARMA Modeling. Now: Important class of stationary processes 3. ARMA Modeling Now: Important class of stationary processes Definition 3.1: (ARMA(p, q) process) Let {ɛ t } t Z WN(0, σ 2 ) be a white noise process. The process {X t } t Z is called AutoRegressive-Moving-Average

More information

Lectures on Stochastic Stability. Sergey FOSS. Heriot-Watt University. Lecture 4. Coupling and Harris Processes

Lectures on Stochastic Stability. Sergey FOSS. Heriot-Watt University. Lecture 4. Coupling and Harris Processes Lectures on Stochastic Stability Sergey FOSS Heriot-Watt University Lecture 4 Coupling and Harris Processes 1 A simple example Consider a Markov chain X n in a countable state space S with transition probabilities

More information

1 Linear Difference Equations

1 Linear Difference Equations ARMA Handout Jialin Yu 1 Linear Difference Equations First order systems Let {ε t } t=1 denote an input sequence and {y t} t=1 sequence generated by denote an output y t = φy t 1 + ε t t = 1, 2,... with

More information

Time Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY

Time Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY & Contents PREFACE xiii 1 1.1. 1.2. Difference Equations First-Order Difference Equations 1 /?th-order Difference

More information

Convolution Based Unit Root Processes: a Simulation Approach

Convolution Based Unit Root Processes: a Simulation Approach International Journal of Statistics and Probability; Vol., No. 6; November 26 ISSN 927-732 E-ISSN 927-74 Published by Canadian Center of Science and Education Convolution Based Unit Root Processes: a Simulation

More information

Econometric Forecasting

Econometric Forecasting Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 1, 2014 Outline Introduction Model-free extrapolation Univariate time-series models Trend

More information

M-estimators for augmented GARCH(1,1) processes

M-estimators for augmented GARCH(1,1) processes M-estimators for augmented GARCH(1,1) processes Freiburg, DAGStat 2013 Fabian Tinkl 19.03.2013 Chair of Statistics and Econometrics FAU Erlangen-Nuremberg Outline Introduction The augmented GARCH(1,1)

More information

Mi-Hwa Ko. t=1 Z t is true. j=0

Mi-Hwa Ko. t=1 Z t is true. j=0 Commun. Korean Math. Soc. 21 (2006), No. 4, pp. 779 786 FUNCTIONAL CENTRAL LIMIT THEOREMS FOR MULTIVARIATE LINEAR PROCESSES GENERATED BY DEPENDENT RANDOM VECTORS Mi-Hwa Ko Abstract. Let X t be an m-dimensional

More information

Almost Sure Convergence of Two Time-Scale Stochastic Approximation Algorithms

Almost Sure Convergence of Two Time-Scale Stochastic Approximation Algorithms Almost Sure Convergence of Two Time-Scale Stochastic Approximation Algorithms Vladislav B. Tadić Abstract The almost sure convergence of two time-scale stochastic approximation algorithms is analyzed under

More information

Nonparametric regression with martingale increment errors

Nonparametric regression with martingale increment errors S. Gaïffas (LSTA - Paris 6) joint work with S. Delattre (LPMA - Paris 7) work in progress Motivations Some facts: Theoretical study of statistical algorithms requires stationary and ergodicity. Concentration

More information

Department of Econometrics and Business Statistics

Department of Econometrics and Business Statistics ISSN 440-77X Australia Department of Econometrics and Business Statistics http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/ Nonlinear Regression with Harris Recurrent Markov Chains Degui Li, Dag

More information

Stochastic process for macro

Stochastic process for macro Stochastic process for macro Tianxiao Zheng SAIF 1. Stochastic process The state of a system {X t } evolves probabilistically in time. The joint probability distribution is given by Pr(X t1, t 1 ; X t2,

More information

The largest eigenvalues of the sample covariance matrix. in the heavy-tail case

The largest eigenvalues of the sample covariance matrix. in the heavy-tail case The largest eigenvalues of the sample covariance matrix 1 in the heavy-tail case Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia NY), Johannes Heiny (Aarhus University)

More information

A Study on "Spurious Long Memory in Nonlinear Time Series Models"

A Study on Spurious Long Memory in Nonlinear Time Series Models A Study on "Spurious Long Memory in Nonlinear Time Series Models" Heri Kuswanto and Philipp Sibbertsen Leibniz University Hanover, Department of Economics, Germany Discussion Paper No. 410 November 2008

More information

Geometric ρ-mixing property of the interarrival times of a stationary Markovian Arrival Process

Geometric ρ-mixing property of the interarrival times of a stationary Markovian Arrival Process Author manuscript, published in "Journal of Applied Probability 50, 2 (2013) 598-601" Geometric ρ-mixing property of the interarrival times of a stationary Markovian Arrival Process L. Hervé and J. Ledoux

More information

SUPPLEMENT TO TESTING FOR REGIME SWITCHING: A COMMENT (Econometrica, Vol. 80, No. 4, July 2012, )

SUPPLEMENT TO TESTING FOR REGIME SWITCHING: A COMMENT (Econometrica, Vol. 80, No. 4, July 2012, ) Econometrica Supplementary Material SUPPLEMENT TO TESTING FOR REGIME SWITCHING: A COMMENT Econometrica, Vol. 80, No. 4, July 01, 1809 181 BY ANDREW V. CARTER ANDDOUGLAS G. STEIGERWALD We present proof

More information

AR, MA and ARMA models

AR, MA and ARMA models AR, MA and AR by Hedibert Lopes P Based on Tsay s Analysis of Financial Time Series (3rd edition) P 1 Stationarity 2 3 4 5 6 7 P 8 9 10 11 Outline P Linear Time Series Analysis and Its Applications For

More information

Tail bound inequalities and empirical likelihood for the mean

Tail bound inequalities and empirical likelihood for the mean Tail bound inequalities and empirical likelihood for the mean Sandra Vucane 1 1 University of Latvia, Riga 29 th of September, 2011 Sandra Vucane (LU) Tail bound inequalities and EL for the mean 29.09.2011

More information

Central limit theorems for ergodic continuous-time Markov chains with applications to single birth processes

Central limit theorems for ergodic continuous-time Markov chains with applications to single birth processes Front. Math. China 215, 1(4): 933 947 DOI 1.17/s11464-15-488-5 Central limit theorems for ergodic continuous-time Markov chains with applications to single birth processes Yuanyuan LIU 1, Yuhui ZHANG 2

More information

Nonlinear time series

Nonlinear time series Based on the book by Fan/Yao: Nonlinear Time Series Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 27, 2009 Outline Characteristics of

More information

Selected Exercises on Expectations and Some Probability Inequalities

Selected Exercises on Expectations and Some Probability Inequalities Selected Exercises on Expectations and Some Probability Inequalities # If E(X 2 ) = and E X a > 0, then P( X λa) ( λ) 2 a 2 for 0 < λ

More information

A regeneration proof of the central limit theorem for uniformly ergodic Markov chains

A regeneration proof of the central limit theorem for uniformly ergodic Markov chains A regeneration proof of the central limit theorem for uniformly ergodic Markov chains By AJAY JASRA Department of Mathematics, Imperial College London, SW7 2AZ, London, UK and CHAO YANG Department of Mathematics,

More information

Lesson 9: Autoregressive-Moving Average (ARMA) models

Lesson 9: Autoregressive-Moving Average (ARMA) models Lesson 9: Autoregressive-Moving Average (ARMA) models Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@ec.univaq.it Introduction We have seen

More information

Ph.D. Seminar Series in Advanced Mathematical Methods in Economics and Finance UNIT ROOT DISTRIBUTION THEORY I. Roderick McCrorie

Ph.D. Seminar Series in Advanced Mathematical Methods in Economics and Finance UNIT ROOT DISTRIBUTION THEORY I. Roderick McCrorie Ph.D. Seminar Series in Advanced Mathematical Methods in Economics and Finance UNIT ROOT DISTRIBUTION THEORY I Roderick McCrorie School of Economics and Finance University of St Andrews 23 April 2009 One

More information

Change detection problems in branching processes

Change detection problems in branching processes Change detection problems in branching processes Outline of Ph.D. thesis by Tamás T. Szabó Thesis advisor: Professor Gyula Pap Doctoral School of Mathematics and Computer Science Bolyai Institute, University

More information

Ch. 14 Stationary ARMA Process

Ch. 14 Stationary ARMA Process Ch. 14 Stationary ARMA Process A general linear stochastic model is described that suppose a time series to be generated by a linear aggregation of random shock. For practical representation it is desirable

More information

UNIT ROOT TESTING FOR FUNCTIONALS OF LINEAR PROCESSES

UNIT ROOT TESTING FOR FUNCTIONALS OF LINEAR PROCESSES Econometric Theory, 22, 2006, 1 14+ Printed in the United States of America+ DOI: 10+10170S0266466606060014 UNIT ROOT TESTING FOR FUNCTIONALS OF LINEAR PROCESSES WEI BIAO WU University of Chicago We consider

More information

Consistency of the maximum likelihood estimator for general hidden Markov models

Consistency of the maximum likelihood estimator for general hidden Markov models Consistency of the maximum likelihood estimator for general hidden Markov models Jimmy Olsson Centre for Mathematical Sciences Lund University Nordstat 2012 Umeå, Sweden Collaborators Hidden Markov models

More information

Strong approximation for additive functionals of geometrically ergodic Markov chains

Strong approximation for additive functionals of geometrically ergodic Markov chains Strong approximation for additive functionals of geometrically ergodic Markov chains Florence Merlevède Joint work with E. Rio Université Paris-Est-Marne-La-Vallée (UPEM) Cincinnati Symposium on Probability

More information

Some functional (Hölderian) limit theorems and their applications (II)

Some functional (Hölderian) limit theorems and their applications (II) Some functional (Hölderian) limit theorems and their applications (II) Alfredas Račkauskas Vilnius University Outils Statistiques et Probabilistes pour la Finance Université de Rouen June 1 5, Rouen (Rouen

More information

FE570 Financial Markets and Trading. Stevens Institute of Technology

FE570 Financial Markets and Trading. Stevens Institute of Technology FE570 Financial Markets and Trading Lecture 5. Linear Time Series Analysis and Its Applications (Ref. Joel Hasbrouck - Empirical Market Microstructure ) Steve Yang Stevens Institute of Technology 9/25/2012

More information

The Cusum Test for Parameter Change in GARCH(1,1) Models

The Cusum Test for Parameter Change in GARCH(1,1) Models The Cusum Test for Parameter Change in GARCH(,) Models Sangyeol Lee, Yasuyoshi Tokutsu and Koichi Maekawa 3 Department of Statistics, Seoul National University, Seoul, 5-74, Korea Graduate School of Social

More information

3 Theory of stationary random processes

3 Theory of stationary random processes 3 Theory of stationary random processes 3.1 Linear filters and the General linear process A filter is a transformation of one random sequence {U t } into another, {Y t }. A linear filter is a transformation

More information

ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications

ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications Yongmiao Hong Department of Economics & Department of Statistical Sciences Cornell University Spring 2019 Time and uncertainty

More information

Basic concepts and terminology: AR, MA and ARMA processes

Basic concepts and terminology: AR, MA and ARMA processes ECON 5101 ADVANCED ECONOMETRICS TIME SERIES Lecture note no. 1 (EB) Erik Biørn, Department of Economics Version of February 1, 2011 Basic concepts and terminology: AR, MA and ARMA processes This lecture

More information

Analytical derivates of the APARCH model

Analytical derivates of the APARCH model Analytical derivates of the APARCH model Sébastien Laurent Forthcoming in Computational Economics October 24, 2003 Abstract his paper derives analytical expressions for the score of the APARCH model of

More information

Online Appendix. j=1. φ T (ω j ) vec (EI T (ω j ) f θ0 (ω j )). vec (EI T (ω) f θ0 (ω)) = O T β+1/2) = o(1), M 1. M T (s) exp ( isω)

Online Appendix. j=1. φ T (ω j ) vec (EI T (ω j ) f θ0 (ω j )). vec (EI T (ω) f θ0 (ω)) = O T β+1/2) = o(1), M 1. M T (s) exp ( isω) Online Appendix Proof of Lemma A.. he proof uses similar arguments as in Dunsmuir 979), but allowing for weak identification and selecting a subset of frequencies using W ω). It consists of two steps.

More information

Pure Random process Pure Random Process or White Noise Process: is a random process {X t, t 0} which has: { σ 2 if k = 0 0 if k 0

Pure Random process Pure Random Process or White Noise Process: is a random process {X t, t 0} which has: { σ 2 if k = 0 0 if k 0 MODULE 9: STATIONARY PROCESSES 7 Lecture 2 Autoregressive Processes 1 Moving Average Process Pure Random process Pure Random Process or White Noise Process: is a random process X t, t 0} which has: E[X

More information

Time Series Models and Inference. James L. Powell Department of Economics University of California, Berkeley

Time Series Models and Inference. James L. Powell Department of Economics University of California, Berkeley Time Series Models and Inference James L. Powell Department of Economics University of California, Berkeley Overview In contrast to the classical linear regression model, in which the components of the

More information

GARCH processes probabilistic properties (Part 1)

GARCH processes probabilistic properties (Part 1) GARCH processes probabilistic properties (Part 1) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/

More information

Lecture 10. Theorem 1.1 [Ergodicity and extremality] A probability measure µ on (Ω, F) is ergodic for T if and only if it is an extremal point in M.

Lecture 10. Theorem 1.1 [Ergodicity and extremality] A probability measure µ on (Ω, F) is ergodic for T if and only if it is an extremal point in M. Lecture 10 1 Ergodic decomposition of invariant measures Let T : (Ω, F) (Ω, F) be measurable, and let M denote the space of T -invariant probability measures on (Ω, F). Then M is a convex set, although

More information

Test for Parameter Change in ARIMA Models

Test for Parameter Change in ARIMA Models Test for Parameter Change in ARIMA Models Sangyeol Lee 1 Siyun Park 2 Koichi Maekawa 3 and Ken-ichi Kawai 4 Abstract In this paper we consider the problem of testing for parameter changes in ARIMA models

More information

On Reparametrization and the Gibbs Sampler

On Reparametrization and the Gibbs Sampler On Reparametrization and the Gibbs Sampler Jorge Carlos Román Department of Mathematics Vanderbilt University James P. Hobert Department of Statistics University of Florida March 2014 Brett Presnell Department

More information

Generalized Method of Moments: I. Chapter 9, R. Davidson and J.G. MacKinnon, Econometric Theory and Methods, 2004, Oxford.

Generalized Method of Moments: I. Chapter 9, R. Davidson and J.G. MacKinnon, Econometric Theory and Methods, 2004, Oxford. Generalized Method of Moments: I References Chapter 9, R. Davidson and J.G. MacKinnon, Econometric heory and Methods, 2004, Oxford. Chapter 5, B. E. Hansen, Econometrics, 2006. http://www.ssc.wisc.edu/~bhansen/notes/notes.htm

More information

LIST OF PUBLICATIONS. 1. J.-P. Kreiss and E. Paparoditis, Bootstrap for Time Series: Theory and Applications, Springer-Verlag, New York, To appear.

LIST OF PUBLICATIONS. 1. J.-P. Kreiss and E. Paparoditis, Bootstrap for Time Series: Theory and Applications, Springer-Verlag, New York, To appear. LIST OF PUBLICATIONS BOOKS 1. J.-P. Kreiss and E. Paparoditis, Bootstrap for Time Series: Theory and Applications, Springer-Verlag, New York, To appear. JOURNAL PAPERS 61. D. Pilavakis, E. Paparoditis

More information

Time Series Analysis -- An Introduction -- AMS 586

Time Series Analysis -- An Introduction -- AMS 586 Time Series Analysis -- An Introduction -- AMS 586 1 Objectives of time series analysis Data description Data interpretation Modeling Control Prediction & Forecasting 2 Time-Series Data Numerical data

More information

The Functional Central Limit Theorem and Testing for Time Varying Parameters

The Functional Central Limit Theorem and Testing for Time Varying Parameters NBER Summer Institute Minicourse What s New in Econometrics: ime Series Lecture : July 4, 008 he Functional Central Limit heorem and esting for ime Varying Parameters Lecture -, July, 008 Outline. FCL.

More information

Introduction to Stochastic processes

Introduction to Stochastic processes Università di Pavia Introduction to Stochastic processes Eduardo Rossi Stochastic Process Stochastic Process: A stochastic process is an ordered sequence of random variables defined on a probability space

More information

THE GEOMETRICAL ERGODICITY OF NONLINEAR AUTOREGRESSIVE MODELS

THE GEOMETRICAL ERGODICITY OF NONLINEAR AUTOREGRESSIVE MODELS Statistica Sinica 6(1996), 943-956 THE GEOMETRICAL ERGODICITY OF NONLINEAR AUTOREGRESSIVE MODELS H. Z. An and F. C. Huang Academia Sinica Abstract: Consider the nonlinear autoregressive model x t = φ(x

More information

Stochastic Equicontinuity in Nonlinear Time Series Models. Andreas Hagemann University of Notre Dame. June 10, 2013

Stochastic Equicontinuity in Nonlinear Time Series Models. Andreas Hagemann University of Notre Dame. June 10, 2013 Stochastic Equicontinuity in Nonlinear Time Series Models Andreas Hagemann University of Notre Dame June 10, 2013 In this paper I provide simple and easily verifiable conditions under which a strong form

More information

Asymptotics for the Arc Length of a Multivariate Time Series and Its Applications as a Measure of Risk

Asymptotics for the Arc Length of a Multivariate Time Series and Its Applications as a Measure of Risk Clemson University TigerPrints All Dissertations Dissertations 12-2012 Asymptotics for the Arc Length of a Multivariate Time Series and Its Applications as a Measure of Risk Tharanga Wickramarachchi Clemson

More information

SIMPLE RANDOM WALKS ALONG ORBITS OF ANOSOV DIFFEOMORPHISMS

SIMPLE RANDOM WALKS ALONG ORBITS OF ANOSOV DIFFEOMORPHISMS SIMPLE RANDOM WALKS ALONG ORBITS OF ANOSOV DIFFEOMORPHISMS V. YU. KALOSHIN, YA. G. SINAI 1. Introduction Consider an automorphism S of probability space (M, M, µ). Any measurable function p : M (0, 1)

More information

STOCHASTIC PROCESSES Basic notions

STOCHASTIC PROCESSES Basic notions J. Virtamo 38.3143 Queueing Theory / Stochastic processes 1 STOCHASTIC PROCESSES Basic notions Often the systems we consider evolve in time and we are interested in their dynamic behaviour, usually involving

More information

Inference in VARs with Conditional Heteroskedasticity of Unknown Form

Inference in VARs with Conditional Heteroskedasticity of Unknown Form Inference in VARs with Conditional Heteroskedasticity of Unknown Form Ralf Brüggemann a Carsten Jentsch b Carsten Trenkler c University of Konstanz University of Mannheim University of Mannheim IAB Nuremberg

More information

Time Series Analysis Fall 2008

Time Series Analysis Fall 2008 MIT OpenCourseWare http://ocw.mit.edu 14.384 Time Series Analysis Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Introduction 1 14.384 Time

More information

Chapter 2. Some basic tools. 2.1 Time series: Theory Stochastic processes

Chapter 2. Some basic tools. 2.1 Time series: Theory Stochastic processes Chapter 2 Some basic tools 2.1 Time series: Theory 2.1.1 Stochastic processes A stochastic process is a sequence of random variables..., x 0, x 1, x 2,.... In this class, the subscript always means time.

More information

SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES

SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES RUTH J. WILLIAMS October 2, 2017 Department of Mathematics, University of California, San Diego, 9500 Gilman Drive,

More information

Stochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno

Stochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.

More information

Chapter 9: Forecasting

Chapter 9: Forecasting Chapter 9: Forecasting One of the critical goals of time series analysis is to forecast (predict) the values of the time series at times in the future. When forecasting, we ideally should evaluate the

More information

High-dimensional Markov Chain Models for Categorical Data Sequences with Applications Wai-Ki CHING AMACL, Department of Mathematics HKU 19 March 2013

High-dimensional Markov Chain Models for Categorical Data Sequences with Applications Wai-Ki CHING AMACL, Department of Mathematics HKU 19 March 2013 High-dimensional Markov Chain Models for Categorical Data Sequences with Applications Wai-Ki CHING AMACL, Department of Mathematics HKU 19 March 2013 Abstract: Markov chains are popular models for a modelling

More information

A time series is called strictly stationary if the joint distribution of every collection (Y t

A time series is called strictly stationary if the joint distribution of every collection (Y t 5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a

More information

Krzysztof Burdzy University of Washington. = X(Y (t)), t 0}

Krzysztof Burdzy University of Washington. = X(Y (t)), t 0} VARIATION OF ITERATED BROWNIAN MOTION Krzysztof Burdzy University of Washington 1. Introduction and main results. Suppose that X 1, X 2 and Y are independent standard Brownian motions starting from 0 and

More information

CUSUM TEST FOR PARAMETER CHANGE IN TIME SERIES MODELS. Sangyeol Lee

CUSUM TEST FOR PARAMETER CHANGE IN TIME SERIES MODELS. Sangyeol Lee CUSUM TEST FOR PARAMETER CHANGE IN TIME SERIES MODELS Sangyeol Lee 1 Contents 1. Introduction of the CUSUM test 2. Test for variance change in AR(p) model 3. Test for Parameter Change in Regression Models

More information

CURRICULUM VITAE. December Robert M. de Jong Tel: (614) Ohio State University Fax: (614)

CURRICULUM VITAE. December Robert M. de Jong Tel: (614) Ohio State University Fax: (614) CURRICULUM VITAE December 2011 Robert M. de Jong Tel: (614) 292-2051 Ohio State University Fax: (614) 292-3906 Department of Economics Email: de-jong.8@osu.edu 444 Arps Hall Columbus, Ohio 43210, USA PERSONAL

More information

STAT 443 Final Exam Review. 1 Basic Definitions. 2 Statistical Tests. L A TEXer: W. Kong

STAT 443 Final Exam Review. 1 Basic Definitions. 2 Statistical Tests. L A TEXer: W. Kong STAT 443 Final Exam Review L A TEXer: W Kong 1 Basic Definitions Definition 11 The time series {X t } with E[X 2 t ] < is said to be weakly stationary if: 1 µ X (t) = E[X t ] is independent of t 2 γ X

More information

Applied time-series analysis

Applied time-series analysis Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 18, 2011 Outline Introduction and overview Econometric Time-Series Analysis In principle,

More information

On Spatial Processes and Asymptotic Inference under. Near-Epoch Dependence

On Spatial Processes and Asymptotic Inference under. Near-Epoch Dependence On Spatial Processes and Asymptotic Inference under Near-Epoch Dependence Nazgul Jenish 1 and Ingmar R. Prucha 2 May 8, 2012 Long version of paper forthcoming in the Journal of Econometrics. This version

More information

SMOOTHED BLOCK EMPIRICAL LIKELIHOOD FOR QUANTILES OF WEAKLY DEPENDENT PROCESSES

SMOOTHED BLOCK EMPIRICAL LIKELIHOOD FOR QUANTILES OF WEAKLY DEPENDENT PROCESSES Statistica Sinica 19 (2009), 71-81 SMOOTHED BLOCK EMPIRICAL LIKELIHOOD FOR QUANTILES OF WEAKLY DEPENDENT PROCESSES Song Xi Chen 1,2 and Chiu Min Wong 3 1 Iowa State University, 2 Peking University and

More information

Testing for non-stationarity

Testing for non-stationarity 20 November, 2009 Overview The tests for investigating the non-stationary of a time series falls into four types: 1 Check the null that there is a unit root against stationarity. Within these, there are

More information

Evgeny Spodarev WIAS, Berlin. Limit theorems for excursion sets of stationary random fields

Evgeny Spodarev WIAS, Berlin. Limit theorems for excursion sets of stationary random fields Evgeny Spodarev 23.01.2013 WIAS, Berlin Limit theorems for excursion sets of stationary random fields page 2 LT for excursion sets of stationary random fields Overview 23.01.2013 Overview Motivation Excursion

More information

Prof. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis

Prof. 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 information

Stochastic relations of random variables and processes

Stochastic relations of random variables and processes Stochastic relations of random variables and processes Lasse Leskelä Helsinki University of Technology 7th World Congress in Probability and Statistics Singapore, 18 July 2008 Fundamental problem of applied

More information

Supplement to Quantile-Based Nonparametric Inference for First-Price Auctions

Supplement to Quantile-Based Nonparametric Inference for First-Price Auctions Supplement to Quantile-Based Nonparametric Inference for First-Price Auctions Vadim Marmer University of British Columbia Artyom Shneyerov CIRANO, CIREQ, and Concordia University August 30, 2010 Abstract

More information

Econ 623 Econometrics II Topic 2: Stationary Time Series

Econ 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 information

If we want to analyze experimental or simulated data we might encounter the following tasks:

If we want to analyze experimental or simulated data we might encounter the following tasks: Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction

More information

Markov Chains, Random Walks on Graphs, and the Laplacian

Markov Chains, Random Walks on Graphs, and the Laplacian Markov Chains, Random Walks on Graphs, and the Laplacian CMPSCI 791BB: Advanced ML Sridhar Mahadevan Random Walks! There is significant interest in the problem of random walks! Markov chain analysis! Computer

More information