Convergence of the long memory Markov switching model to Brownian motion

Size: px
Start display at page:

Download "Convergence of the long memory Markov switching model to Brownian motion"

Transcription

1 Convergence of the long memory Marov switching model to Brownian motion Changryong Bae Sungyunwan University Natércia Fortuna CEF.UP, Universidade do Porto Vladas Pipiras University of North Carolina February 5, 04 Abstract This technical note shows that the partial sum process of the Marov switching model converges to Brownian motion in the sense of finite-dimensional distributions. Definitions The Marov switching MS, in short model is defined as follows. Let s T {s T },...,T be a stationary Marov chain, taing the value 0 or, with the probability transition matrix p00 p P 00.. p p We shall tae below p 00 p 00 T and p p T as functions of T, which explains the use of the index T in the notation s T. Let µ 0, µ R. Consider the series X T µ s T + ϵ,,..., T,. {ϵ } is a sequence of i.i.d. random variables with mean 0 and variance σϵ, and its partial sum process S T t X T EXT, t [0, ]. The term µ s T in X T is not constant. It varies between two values µ 0 and µ. The switching mechanism is dictated by the Marov chain s T,,..., T. If one is in regime j, j 0,, one stays there with probability p jj and switches to the other regime with probability p jj. Suppose now that µ 0 µ.3 and, for j 0,, p jj p jj T c j T δ j,.4 JEL classification: C, C3, C4, C50. Keywords and phrases: Marov switching model, long memory, changes in mean.

2 0 < c j < and 0 < δ 0, δ <..5 Since p jj T, the Marov chain s T tends to stay in the same state as T. Under the assumptions.3.5, Diebold and Inoue 00, Proposition 3, showed that Var S T VarX T XT T CT d+,.6 d min{δ 0, δ } δ 0 δ 0,..7 The behavior.6 of the variance of the partial sum is consistent with that for long memory time series. Recall that a long memory time series X {X } Z is commonly defined as a secondorder stationary time series with autocovariance function CovX 0, X h ch d, as h,.8 c > 0 and d 0, /. Under.8, one can easily show that VarX X T CT d+, as T. The parameter d is nown as the long memory parameter. Convergence result We are interested here not only in the second-order properties of the partial sum process S T t but also in its convergence in distribution. For a long memory time series, a common limit of its partial sum process is fractional Brownian motion with parameter H d + /. Recall that B H {B H t} t R is fractional Brownian motion FBM if it is Gaussian, has stationary increments and is H-self-similar i.e. for any c > 0, the process B H ct and c H B H t have the same finitedimensional distributions. It is called standard if EB H, and if H /, then it is the usual Brownian motion. Does the partial sum process S T t of the Marov switching model converge to FBM with parameter H d + /? The answer to this question is negative, as will be shown below under additional assumptions. In fact, the limit turns out to be the usual Brownian motion. Consequences of this fact for estimation of the long memory parameter in the MS model are discussed in Bae, Fortuna and Pipiras 04. Before providing the convergence result, we state an auxiliary lemma which will be used in the proof of the main result. The result of the lemma is implicit in the proof of Proposition 3 in Diebold and Inoue 00. See Appendix for a proof. Lemma. The term µ s T in the Marov switching model. satisfies: for j, T, Cov µ s, µ T s µ Γ Tm m µ p 00 p p 00 p µ 0 µ λ m,. and µ µ0 µ, Γ j p 00 p p 00 p λ j λ j λ j λ j. λ p 00 + p..3 In the next result, we establish the convergence of the partial sum process to Brownian motion. See Appendix for a proof.

3 Theorem. Suppose that.3.5 hold and that Then, 3 min{δ 0, δ } + δ 0 δ <..4 S T t T /+d fdd σbt,.5 d is given in.7, fdd stands for the convergence of finite-dimensional distributions and B is a standard Brownian motion with and σ c 0c µ 0 µ c 3.6 c 0, if δ 0 < δ, c c, if δ < δ 0, c 0 + c, if δ 0 δ. The condition.4 in Theorem. was used to obtain a normal limit for the partial sum process S T. Note that, when δ δ 0 δ > 0, this condition becomes 0 < δ < /3.7 and hence see.7 d δ 0,. 6 More generally suppose δ 0 δ. Then.7 becomes d δ 0 δ. The condition d > 0 translates to δ 0 > δ /, thus δ / < δ 0 δ. Moreover,.4 becomes δ 0 + δ <, that is, δ 0 < δ /. Two cases need to be distinguished next. First, suppose that δ < δ / or δ < /3. Then, the two conditions δ / < δ 0 δ and δ 0 < δ / reduce to the first one, namely, δ / < δ 0 δ. In particular, d δ 0 δ δ and d 0,, 6 since δ < /3. Second, suppose that δ δ / or δ /3. Then, the two conditions δ / < δ 0 δ and δ 0 < δ / reduce to δ / < δ 0 δ / supposing δ < δ / or δ < /; otherwise, the intervals in the two conditions are disjoint. In particular, d δ 0 δ δ and d 0, 6 since δ /3. What happens when the condition.4 does not hold that is, for other values of δ 0, δ? When t, the convergence to a normal limit follows immediately from the results of Dobrushin 96. In fact, the asymptotic normality is expected for all finite-dimensional distributions. In the proof of Theorem., we used general convergence results for nonstationary Marov chains. In a nearly stationary case, the conditions of those results are expected to be too strong as discussed in the footnote on p. 76 of Dobrushin 956. But to the best of our nowledge, there are currently no results addressing this issue aside from the result corresponding to t. 3,

4 A Technical proofs Proof of Lemma.: Let ξ {s T 0} {s T }, so that µ s T µ ξ, prime denotes the transpose. Observe further that Cov µ s, µ T s Cov µ ξ Tm, µ ξ m µ Covξ, ξ m µ µ Γ m µ, A. Γ j Covξ 0, ξ j Cov{s T 0 {s T 0} j Cov {s T 0}, 0 {s T } j Cov {s T 0 }, {s T 0} j Cov {s T }, 0 {s T } j Ps T 0 0, s T j 0 PsT 0 0PsT j 0 PsT 0 0, st j PsT 0 0PsT j Ps T 0, st j 0 PsT 0 PsT j 0 PsT 0, st j PsT 0 PsT j Denoting by π π 0 π the invariant stationary probabilities of s T, and by p j nm the transition probability from state n to state m in j steps, we have π Γ j 0 p j 00 π 0 π 0 p j 0 π π p j 0 π 0 π p j π. A. To evaluate Γ j, we need expressions for π 0, π and p j nm. One can show that π 0 p p 00 p, π. p 00 p 00 p A.3 and that the probability transition matrix P in. is diagonalizable as π 0 π0 π P 0 : E π 0 0 λ 0 λ F, A.4 λ is given by.3. Note that and λ in A.4 are the eigenvalues of the matrix P and the matrices E and F are made of the right- and left-eigenvectors of P, satisfying F E. The latter fact implies that p 00 p P j π π λ j π0 π p + p 00 λ j p 00 λ j p λ j p 00 + p λ j π0 + π λ j π λ j π 0 λ j π + π 0 λ j p j 00 p j 0 p j 0 p j. A.5 By using A.3 and A.5, we obtain from A. that Γ j established by using A.. is given by. and. is finally Proof of Theorem.: It is enough to show that, for t,..., t n [0, ] and θ,..., θ n R, n S T t T /+d d N 0, σ n A.6 4

5 and that Since X T σ n lim T E n is a stationary series the relation A.7 follows from n S T t T /+d σ E Bt. A.7 ST t E T /+d σ t, A.8 as T. Write S T t µ s T Eµ s T + ϵ : S,T t + S,T t. Then, by the independence of s T and {ϵ }, A.9 and Hence, it is enough to show that Observe that ES,T t Var From Lemma. and by using ST t S,T t S,T t E T /+d E T /+d + E T /+d S,T t E σ ϵ 0. T +d µ s T T /+d µ Γ 0 + S,T t E σ t. A.0 T /+d Var µ ξ µ, jγ j + Γ j µ. j Covξ, ξ µ j λ j λ λ λ, j jλ j λ λ λ λ +, A. we get V 00 p 00 p p 00 p ES,T t V 00 µ 0 µ, λ λ + + λ+ λ λ. A. But by.3 and.4, λ p 00 + p c 0 T δ 0 c T δ. 5

6 Hence, as T, c is defined in.7, and λ λ c 0 T δ 0 + c T δ T min{δ0,δ}, c c 0 T δ c c { } 0 T δ T min{δ exp ct min{δ 0,δ } o, 0,δ } since δ 0, δ <. This shows that Now, V 00 p 00 p λ p 00 p λ p 00 p p 00 p T min{δ 0,δ }. c p 00 p p 00 p c 0T δ0 c T δ c T min{δ 0,δ } c 0c c T δ 0 δ, A.3 since δ 0 δ + min{δ 0, δ } δ 0 δ. Thus, by.4, Relations A. and A.4 now yield V 00 c 0c c 3 T min{δ 0,δ } δ 0 δ c 0c c 3 T d. A.4 ES,T T +d c 0c µ 0 µ t c 3 σ t, which establishes A.0. We now turn to the convergence A.6. In view of the decomposition A.9 and since, by the central limit theorem, as T, n 0 < σ,n <, it is enough to show that σ n is defined in A.7. Write n n f T t x S,T t T / d N 0, σ,n, S,T t T /+d d N 0, σ n, A.5 S,T t T /+d n T t f T t s T t, T /+d {t [T t ]}µ x Eµ s T t. Recall that s T t may be equal to 0 or. We shall apply next a central limit theorem for functionals of Marov chains. Let T : max Ps T t+ A s T t x Ps T t+ A s T t x x,x,a 6

7 be the so-called contraction coefficient for the Marov chain s T. Note first that T max Ps T t+ A s T t 0 Ps T t+ A s T t, A since x, x {0, } and the difference Ps T t+ A st t x Ps T t+ A st t x is zero whenever x x. Note next that the difference Ps T t+ A st t 0 Ps T t+ A st t is zero for A {0, }, and equals p 00 + p for A {0} and A {}. Hence, By.4, as T. Observe that T p 00 + p. α T : T c 0 T δ 0 + c T δ c T min{δ 0,δ }, sup t T sup f T t x C T : C x T /+d. A.6 By Theorem. in Sethuraman and Varadhan 005 see also Theorem in Dobrushin 956, the convergence A.5 follows if C T α 3 T T Since {t [T t ]} for large enough T, note that, for large enough T, Varf T t In view of A.3, this yields Varf T t t s T t n p 00 p T +d Varf T t s T t 0. A.7 s T t n T +d Varµ s T t n T +d µ Γ 0 µ. p 00 p µ 0 µ n µ 0 µ c 0 c c T d δ0 δ Combining this with the asymptotic form of α T given in A.6, we conclude that the condition A.7 holds if T +d T 3 min{δ 0,δ } T d+ δ 0 δ T 3 min{δ 0,δ }+ δ 0 δ 0. The convergence to 0 follows from the assumption.4. References Bae, C., Fortuna, N. and Pipiras, V. 04, Can Marov switching model generate long memory?, Preprint. Available at Diebold, F. X. and Inoue, A. 00, Long memory and regime switching, Journal of Econometrics 05, Dobrushin, R. L. 956, Central limit theorem for non-stationary Marov chains. I, Theory of Probability and its Applications, Dobrushin, R. L. 96, Limit theorems for a Marov chain of two states, in Select. Transl. Math. Statist. and Probability, Vol., Inst. Math. Statist. and Amer. Math. Soc., Providence, R.I., pp T. 7

8 Sethuraman, S. and Varadhan, S. 005, A martingale proof of Dobrushin s theorem for non-homogeneous Marov chains, Electronic Journal of Probability 0. Changryong Bae Natércia Fortuna Dept. of Statistics Faculdade de Economia Sungyunwan University Universidade do Porto 5-, Sungyunwan-ro, Jongno-gu Rua Dr. Roberto Frias Seoul, 0-745, Korea Porto, Portugal crbae@su.edu nfortuna@fep.up.pt Vladas Pipiras Dept. of Statistics and Operations Research UNC at Chapel Hill CB#360, Hanes Hall Chapel Hill, NC 7599, USA pipiras@ .unc.edu 8

On the usefulness of wavelet-based simulation of fractional Brownian motion

On the usefulness of wavelet-based simulation of fractional Brownian motion On the usefulness of wavelet-based simulation of fractional Brownian motion Vladas Pipiras University of North Carolina at Chapel Hill September 16, 2004 Abstract We clarify some ways in which wavelet-based

More information

Brownian survival and Lifshitz tail in perturbed lattice disorder

Brownian survival and Lifshitz tail in perturbed lattice disorder Brownian survival and Lifshitz tail in perturbed lattice disorder Ryoki Fukushima Kyoto niversity Random Processes and Systems February 16, 2009 6 B T 1. Model ) ({B t t 0, P x : standard Brownian motion

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

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

Bernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010

Bernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010 1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes

More information

Gaussian Processes. 1. Basic Notions

Gaussian Processes. 1. Basic Notions Gaussian Processes 1. Basic Notions Let T be a set, and X : {X } T a stochastic process, defined on a suitable probability space (Ω P), that is indexed by T. Definition 1.1. We say that X is a Gaussian

More information

Modeling and testing long memory in random fields

Modeling and testing long memory in random fields Modeling and testing long memory in random fields Frédéric Lavancier lavancier@math.univ-lille1.fr Université Lille 1 LS-CREST Paris 24 janvier 6 1 Introduction Long memory random fields Motivations Previous

More information

Stationary independent increments. 1. Random changes of the form X t+h X t fixed h > 0 are called increments of the process.

Stationary independent increments. 1. Random changes of the form X t+h X t fixed h > 0 are called increments of the process. Stationary independent increments 1. Random changes of the form X t+h X t fixed h > 0 are called increments of the process. 2. If each set of increments, corresponding to non-overlapping collection of

More information

In terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.

In terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3. 1. GAUSSIAN PROCESSES A Gaussian process on a set T is a collection of random variables X =(X t ) t T on a common probability space such that for any n 1 and any t 1,...,t n T, the vector (X(t 1 ),...,X(t

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

Integral representations in models with long memory

Integral representations in models with long memory Integral representations in models with long memory Georgiy Shevchenko, Yuliya Mishura, Esko Valkeila, Lauri Viitasaari, Taras Shalaiko Taras Shevchenko National University of Kyiv 29 September 215, Ulm

More information

Lecture 9. d N(0, 1). Now we fix n and think of a SRW on [0,1]. We take the k th step at time k n. and our increments are ± 1

Lecture 9. d N(0, 1). Now we fix n and think of a SRW on [0,1]. We take the k th step at time k n. and our increments are ± 1 Random Walks and Brownian Motion Tel Aviv University Spring 011 Lecture date: May 0, 011 Lecture 9 Instructor: Ron Peled Scribe: Jonathan Hermon In today s lecture we present the Brownian motion (BM).

More information

A remark on the maximum eigenvalue for circulant matrices

A remark on the maximum eigenvalue for circulant matrices IMS Collections High Dimensional Probability V: The Luminy Volume Vol 5 (009 79 84 c Institute of Mathematical Statistics, 009 DOI: 04/09-IMSCOLL5 A remark on the imum eigenvalue for circulant matrices

More information

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory

More information

CONVERGENCE TO FRACTIONAL BROWNIAN MOTION AND LOSS PROBABILITY. Jin-Chun Kim and Hee-Choon Lee

CONVERGENCE TO FRACTIONAL BROWNIAN MOTION AND LOSS PROBABILITY. Jin-Chun Kim and Hee-Choon Lee Kangweon-Kyungki Math. Jour. (2003), No., pp. 35 43 CONVERGENCE TO FRACTIONAL BROWNIAN MOTION AND LOSS PROBABILITY Jin-Chun Kim and Hee-Choon Lee Abstract. We study the weak convergence to Fractional Brownian

More information

PACKING-DIMENSION PROFILES AND FRACTIONAL BROWNIAN MOTION

PACKING-DIMENSION PROFILES AND FRACTIONAL BROWNIAN MOTION PACKING-DIMENSION PROFILES AND FRACTIONAL BROWNIAN MOTION DAVAR KHOSHNEVISAN AND YIMIN XIAO Abstract. In order to compute the packing dimension of orthogonal projections Falconer and Howroyd 997) introduced

More information

1.1 Limits and Continuity. Precise definition of a limit and limit laws. Squeeze Theorem. Intermediate Value Theorem. Extreme Value Theorem.

1.1 Limits and Continuity. Precise definition of a limit and limit laws. Squeeze Theorem. Intermediate Value Theorem. Extreme Value Theorem. STATE EXAM MATHEMATICS Variant A ANSWERS AND SOLUTIONS 1 1.1 Limits and Continuity. Precise definition of a limit and limit laws. Squeeze Theorem. Intermediate Value Theorem. Extreme Value Theorem. Definition

More information

In Memory of Wenbo V Li s Contributions

In Memory of Wenbo V Li s Contributions In Memory of Wenbo V Li s Contributions Qi-Man Shao The Chinese University of Hong Kong qmshao@cuhk.edu.hk The research is partially supported by Hong Kong RGC GRF 403513 Outline Lower tail probabilities

More information

GENERALIZED CANTOR SETS AND SETS OF SUMS OF CONVERGENT ALTERNATING SERIES

GENERALIZED CANTOR SETS AND SETS OF SUMS OF CONVERGENT ALTERNATING SERIES Journal of Applied Analysis Vol. 7, No. 1 (2001), pp. 131 150 GENERALIZED CANTOR SETS AND SETS OF SUMS OF CONVERGENT ALTERNATING SERIES M. DINDOŠ Received September 7, 2000 and, in revised form, February

More information

Adaptive wavelet decompositions of stochastic processes and some applications

Adaptive wavelet decompositions of stochastic processes and some applications Adaptive wavelet decompositions of stochastic processes and some applications Vladas Pipiras University of North Carolina at Chapel Hill SCAM meeting, June 1, 2012 (joint work with G. Didier, P. Abry)

More information

Central Limit Theorem for Non-stationary Markov Chains

Central Limit Theorem for Non-stationary Markov Chains Central Limit Theorem for Non-stationary Markov Chains Magda Peligrad University of Cincinnati April 2011 (Institute) April 2011 1 / 30 Plan of talk Markov processes with Nonhomogeneous transition probabilities

More information

Gaussian, Markov and stationary processes

Gaussian, Markov and stationary processes Gaussian, Markov and stationary processes Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ November

More information

ON THE CONVERGENCE OF FARIMA SEQUENCE TO FRACTIONAL GAUSSIAN NOISE. Joo-Mok Kim* 1. Introduction

ON THE CONVERGENCE OF FARIMA SEQUENCE TO FRACTIONAL GAUSSIAN NOISE. Joo-Mok Kim* 1. Introduction JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volume 26, No. 2, May 2013 ON THE CONVERGENCE OF FARIMA SEQUENCE TO FRACTIONAL GAUSSIAN NOISE Joo-Mok Kim* Abstract. We consider fractional Gussian noise

More information

. Find E(V ) and var(v ).

. Find E(V ) and var(v ). Math 6382/6383: Probability Models and Mathematical Statistics Sample Preliminary Exam Questions 1. A person tosses a fair coin until she obtains 2 heads in a row. She then tosses a fair die the same number

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

Lecture 19 : Brownian motion: Path properties I

Lecture 19 : Brownian motion: Path properties I Lecture 19 : Brownian motion: Path properties I MATH275B - Winter 2012 Lecturer: Sebastien Roch References: [Dur10, Section 8.1], [Lig10, Section 1.5, 1.6], [MP10, Section 1.1, 1.2]. 1 Invariance We begin

More information

ELEMENTS OF PROBABILITY THEORY

ELEMENTS OF PROBABILITY THEORY ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable

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

Summary of the Ph.D. Thesis

Summary of the Ph.D. Thesis Summary of the Ph.D. Thesis Katalin Varga Institute of Mathematics and Informatics Debrecen University Pf. 2, H 4 Debrecen, Hungary Introduction In the dissertation stable and nearly unstable multidimensional

More information

Asymptotic distribution of the Yule-Walker estimator for INAR(p) processes

Asymptotic distribution of the Yule-Walker estimator for INAR(p) processes Asymptotic distribution of the Yule-Walker estimator for IAR(p) processes Isabel Silva a,b, M Eduarda Silva b,c a Departamento de Engenharia Civil, Faculdade de Engenharia - Universidade do Porto, Rua

More information

Lectures 2 3 : Wigner s semicircle law

Lectures 2 3 : Wigner s semicircle law Fall 009 MATH 833 Random Matrices B. Való Lectures 3 : Wigner s semicircle law Notes prepared by: M. Koyama As we set up last wee, let M n = [X ij ] n i,j=1 be a symmetric n n matrix with Random entries

More information

Jae Gil Choi and Young Seo Park

Jae Gil Choi and Young Seo Park Kangweon-Kyungki Math. Jour. 11 (23), No. 1, pp. 17 3 TRANSLATION THEOREM ON FUNCTION SPACE Jae Gil Choi and Young Seo Park Abstract. In this paper, we use a generalized Brownian motion process to define

More information

Statistics & Data Sciences: First Year Prelim Exam May 2018

Statistics & Data Sciences: First Year Prelim Exam May 2018 Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book

More information

arxiv: v1 [math.co] 17 Dec 2007

arxiv: v1 [math.co] 17 Dec 2007 arxiv:07.79v [math.co] 7 Dec 007 The copies of any permutation pattern are asymptotically normal Milós Bóna Department of Mathematics University of Florida Gainesville FL 36-805 bona@math.ufl.edu Abstract

More information

A Concise Course on Stochastic Partial Differential Equations

A Concise Course on Stochastic Partial Differential Equations A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original

More information

LONG RANGE DEPENDENCE, UNBALANCED HAAR WAVELET TRANSFORMATION AND CHANGE IN LOCAL MEAN LEVEL.

LONG RANGE DEPENDENCE, UNBALANCED HAAR WAVELET TRANSFORMATION AND CHANGE IN LOCAL MEAN LEVEL. International Journal of Wavelets, Multiresolution and Information Processing c World Scientific Publishing Company LONG RANGE DEPENDENCE, UNBALANCED HAAR WAVELET TRANSFORMATION AND CHANGE IN LOCAL MEAN

More information

ON THE ZERO-ONE LAW AND THE LAW OF LARGE NUMBERS FOR RANDOM WALK IN MIXING RAN- DOM ENVIRONMENT

ON THE ZERO-ONE LAW AND THE LAW OF LARGE NUMBERS FOR RANDOM WALK IN MIXING RAN- DOM ENVIRONMENT Elect. Comm. in Probab. 10 (2005), 36 44 ELECTRONIC COMMUNICATIONS in PROBABILITY ON THE ZERO-ONE LAW AND THE LAW OF LARGE NUMBERS FOR RANDOM WALK IN MIXING RAN- DOM ENVIRONMENT FIRAS RASSOUL AGHA Department

More information

{ ϕ(f )Xt, 1 t n p, X t, n p + 1 t n, { 0, t=1, W t = E(W t W s, 1 s t 1), 2 t n.

{ ϕ(f )Xt, 1 t n p, X t, n p + 1 t n, { 0, t=1, W t = E(W t W s, 1 s t 1), 2 t n. Innovations, Sufficient Statistics, And Maximum Lielihood In ARM A Models Georgi N. Boshnaov Institute of Mathematics, Bulgarian Academy of Sciences 1. Introduction. Let X t, t Z} be a zero mean ARMA(p,

More information

Bernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012

Bernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012 1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.

More information

Tilburg University. Strongly Regular Graphs with Maximal Energy Haemers, W. H. Publication date: Link to publication

Tilburg University. Strongly Regular Graphs with Maximal Energy Haemers, W. H. Publication date: Link to publication Tilburg University Strongly Regular Graphs with Maximal Energy Haemers, W. H. Publication date: 2007 Link to publication Citation for published version (APA): Haemers, W. H. (2007). Strongly Regular Graphs

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

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

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

A Functional Central Limit Theorem for an ARMA(p, q) Process with Markov Switching Communications for Statistical Applications and Methods 2013, Vol 20, No 4, 339 345 DOI: http://dxdoiorg/105351/csam2013204339 A Functional Central Limit Theorem for an ARMAp, q) Process with Markov Switching

More information

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title Composition operators on Hilbert spaces of entire functions Author(s) Doan, Minh Luan; Khoi, Le Hai Citation

More information

Yunhi Cho and Young-One Kim

Yunhi Cho and Young-One Kim Bull. Korean Math. Soc. 41 (2004), No. 1, pp. 27 43 ANALYTIC PROPERTIES OF THE LIMITS OF THE EVEN AND ODD HYPERPOWER SEQUENCES Yunhi Cho Young-One Kim Dedicated to the memory of the late professor Eulyong

More information

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3 Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................

More information

Maximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments

Maximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments Austrian Journal of Statistics April 27, Volume 46, 67 78. AJS http://www.ajs.or.at/ doi:.773/ajs.v46i3-4.672 Maximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments Yuliya

More information

Estimation of arrival and service rates for M/M/c queue system

Estimation of arrival and service rates for M/M/c queue system Estimation of arrival and service rates for M/M/c queue system Katarína Starinská starinskak@gmail.com Charles University Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics

More information

University of Toronto Department of Statistics

University of Toronto Department of Statistics Norm Comparisons for Data Augmentation by James P. Hobert Department of Statistics University of Florida and Jeffrey S. Rosenthal Department of Statistics University of Toronto Technical Report No. 0704

More information

Lower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions

Lower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions Lower Tail Probabilities and Normal Comparison Inequalities In Memory of Wenbo V. Li s Contributions Qi-Man Shao The Chinese University of Hong Kong Lower Tail Probabilities and Normal Comparison Inequalities

More information

Nonminimal sets, their projections and integral representations of stable processes

Nonminimal sets, their projections and integral representations of stable processes Stochastic Processes and their Applications 117 (2007) 1285 1302 www.elsevier.com/locate/spa Nonminimal sets, their projections and integral representations of stable processes Vladas Pipiras Department

More information

Documents de Travail du Centre d Economie de la Sorbonne

Documents de Travail du Centre d Economie de la Sorbonne Documents de Travail du Centre d Economie de la Sorbonne A note on self-similarity for discrete time series Dominique GUEGAN, Zhiping LU 2007.55 Maison des Sciences Économiques, 106-112 boulevard de L'Hôpital,

More information

Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet 1

Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet 1 Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet. For ξ, ξ 2, i.i.d. with P(ξ i = ± = /2 define the discrete-time random walk W =, W n = ξ +... + ξ n. (i Formulate and prove the property

More information

11. Further Issues in Using OLS with TS Data

11. Further Issues in Using OLS with TS Data 11. Further Issues in Using OLS with TS Data With TS, including lags of the dependent variable often allow us to fit much better the variation in y Exact distribution theory is rarely available in TS applications,

More information

MODELS FOR COMPUTER NETWORK TRAFFIC

MODELS FOR COMPUTER NETWORK TRAFFIC MODELS FOR COMPUTER NETWORK TRAFFIC Murad S. Taqqu Boston University Joint work with Walter Willinger, Joshua Levy and Vladas Pipiras,... Web Site http://math.bu.edu/people/murad OUTLINE Background: 1)

More information

ODE Final exam - Solutions

ODE Final exam - Solutions ODE Final exam - Solutions May 3, 018 1 Computational questions (30 For all the following ODE s with given initial condition, find the expression of the solution as a function of the time variable t You

More information

Asymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½

Asymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½ University of Pennsylvania ScholarlyCommons Statistics Papers Wharton Faculty Research 1998 Asymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½ Lawrence D. Brown University

More information

MINIMAL POLYNOMIALS AND CHARACTERISTIC POLYNOMIALS OVER RINGS

MINIMAL POLYNOMIALS AND CHARACTERISTIC POLYNOMIALS OVER RINGS JP Journal of Algebra, Number Theory and Applications Volume 0, Number 1, 011, Pages 49-60 Published Online: March, 011 This paper is available online at http://pphmj.com/journals/jpanta.htm 011 Pushpa

More information

Independence of some multiple Poisson stochastic integrals with variable-sign kernels

Independence of some multiple Poisson stochastic integrals with variable-sign kernels Independence of some multiple Poisson stochastic integrals with variable-sign kernels Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological

More information

MOMENT CONVERGENCE RATES OF LIL FOR NEGATIVELY ASSOCIATED SEQUENCES

MOMENT CONVERGENCE RATES OF LIL FOR NEGATIVELY ASSOCIATED SEQUENCES J. Korean Math. Soc. 47 1, No., pp. 63 75 DOI 1.4134/JKMS.1.47..63 MOMENT CONVERGENCE RATES OF LIL FOR NEGATIVELY ASSOCIATED SEQUENCES Ke-Ang Fu Li-Hua Hu Abstract. Let X n ; n 1 be a strictly stationary

More information

Multivariate Time Series: VAR(p) Processes and Models

Multivariate Time Series: VAR(p) Processes and Models Multivariate Time Series: VAR(p) Processes and Models A VAR(p) model, for p > 0 is X t = φ 0 + Φ 1 X t 1 + + Φ p X t p + A t, where X t, φ 0, and X t i are k-vectors, Φ 1,..., Φ p are k k matrices, with

More information

JUSTIN HARTMANN. F n Σ.

JUSTIN HARTMANN. F n Σ. BROWNIAN MOTION JUSTIN HARTMANN Abstract. This paper begins to explore a rigorous introduction to probability theory using ideas from algebra, measure theory, and other areas. We start with a basic explanation

More information

Introduction to Machine Learning CMU-10701

Introduction to Machine Learning CMU-10701 Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos & Aarti Singh Contents Markov Chain Monte Carlo Methods Goal & Motivation Sampling Rejection Importance Markov

More information

From Fractional Brownian Motion to Multifractional Brownian Motion

From Fractional Brownian Motion to Multifractional Brownian Motion From Fractional Brownian Motion to Multifractional Brownian Motion Antoine Ayache USTL (Lille) Antoine.Ayache@math.univ-lille1.fr Cassino December 2010 A.Ayache (USTL) From FBM to MBM Cassino December

More information

A Matrix Theoretic Derivation of the Kalman Filter

A Matrix Theoretic Derivation of the Kalman Filter A Matrix Theoretic Derivation of the Kalman Filter 4 September 2008 Abstract This paper presents a matrix-theoretic derivation of the Kalman filter that is accessible to students with a strong grounding

More information

ON REGULARITY OF FINITE REFLECTION GROUPS. School of Mathematics and Statistics, University of Sydney, NSW 2006, Australia

ON REGULARITY OF FINITE REFLECTION GROUPS. School of Mathematics and Statistics, University of Sydney, NSW 2006, Australia ON REGULARITY OF FINITE REFLECTION GROUPS R. B. Howlett and Jian-yi Shi School of Mathematics and Statistics, University of Sydney, NSW 2006, Australia Department of Mathematics, East China Normal University,

More information

Filtrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition

Filtrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition Filtrations, Markov Processes and Martingales Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition David pplebaum Probability and Statistics Department,

More information

CEF.UP Working Paper ON RANK ESTIMATION IN SEMIDEFINITE MATRICES. Stephen G. Donald Natércia Fortuna Vladas Pipiras

CEF.UP Working Paper ON RANK ESTIMATION IN SEMIDEFINITE MATRICES. Stephen G. Donald Natércia Fortuna Vladas Pipiras CEF.UP Working Paper 2010-02 ON RANK ESTIMATION IN SEMIDEFINITE MATRICES Stephen G. Donald Natércia Fortuna Vladas Pipiras On rank estimation in semidefinite matrices Stephen G. Donald University of Texas

More information

On a class of stochastic differential equations in a financial network model

On a class of stochastic differential equations in a financial network model 1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University

More information

ON THE DIFFUSION OPERATOR IN POPULATION GENETICS

ON THE DIFFUSION OPERATOR IN POPULATION GENETICS J. Appl. Math. & Informatics Vol. 30(2012), No. 3-4, pp. 677-683 Website: http://www.kcam.biz ON THE DIFFUSION OPERATOR IN POPULATION GENETICS WON CHOI Abstract. W.Choi([1]) obtains a complete description

More information

Stein s Method and Characteristic Functions

Stein s Method and Characteristic Functions Stein s Method and Characteristic Functions Alexander Tikhomirov Komi Science Center of Ural Division of RAS, Syktyvkar, Russia; Singapore, NUS, 18-29 May 2015 Workshop New Directions in Stein s method

More information

A COUNTEREXAMPLE TO AN ENDPOINT BILINEAR STRICHARTZ INEQUALITY TERENCE TAO. t L x (R R2 ) f L 2 x (R2 )

A COUNTEREXAMPLE TO AN ENDPOINT BILINEAR STRICHARTZ INEQUALITY TERENCE TAO. t L x (R R2 ) f L 2 x (R2 ) Electronic Journal of Differential Equations, Vol. 2006(2006), No. 5, pp. 6. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A COUNTEREXAMPLE

More information

On the estimation of the heavy tail exponent in time series using the max spectrum. Stilian A. Stoev

On the estimation of the heavy tail exponent in time series using the max spectrum. Stilian A. Stoev On the estimation of the heavy tail exponent in time series using the max spectrum Stilian A. Stoev (sstoev@umich.edu) University of Michigan, Ann Arbor, U.S.A. JSM, Salt Lake City, 007 joint work with:

More information

1 Simulating normal (Gaussian) rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions

1 Simulating normal (Gaussian) rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions Copyright c 2007 by Karl Sigman 1 Simulating normal Gaussian rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions Fundamental to many applications

More information

IEOR 4701: Stochastic Models in Financial Engineering. Summer 2007, Professor Whitt. SOLUTIONS to Homework Assignment 9: Brownian motion

IEOR 4701: Stochastic Models in Financial Engineering. Summer 2007, Professor Whitt. SOLUTIONS to Homework Assignment 9: Brownian motion IEOR 471: Stochastic Models in Financial Engineering Summer 27, Professor Whitt SOLUTIONS to Homework Assignment 9: Brownian motion In Ross, read Sections 1.1-1.3 and 1.6. (The total required reading there

More information

Asymptotically sufficient statistics in nonparametric regression experiments with correlated noise

Asymptotically sufficient statistics in nonparametric regression experiments with correlated noise Vol. 0 0000 1 0 Asymptotically sufficient statistics in nonparametric regression experiments with correlated noise Andrew V Carter University of California, Santa Barbara Santa Barbara, CA 93106-3110 e-mail:

More information

Joint Parameter Estimation of the Ornstein-Uhlenbeck SDE driven by Fractional Brownian Motion

Joint Parameter Estimation of the Ornstein-Uhlenbeck SDE driven by Fractional Brownian Motion Joint Parameter Estimation of the Ornstein-Uhlenbeck SDE driven by Fractional Brownian Motion Luis Barboza October 23, 2012 Department of Statistics, Purdue University () Probability Seminar 1 / 59 Introduction

More information

Supermodular ordering of Poisson arrays

Supermodular ordering of Poisson arrays Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore

More information

Spectral representations and ergodic theorems for stationary stochastic processes

Spectral representations and ergodic theorems for stationary stochastic processes AMS 263 Stochastic Processes (Fall 2005) Instructor: Athanasios Kottas Spectral representations and ergodic theorems for stationary stochastic processes Stationary stochastic processes Theory and methods

More information

Synthesis of Gaussian and non-gaussian stationary time series using circulant matrix embedding

Synthesis of Gaussian and non-gaussian stationary time series using circulant matrix embedding Synthesis of Gaussian and non-gaussian stationary time series using circulant matrix embedding Vladas Pipiras University of North Carolina at Chapel Hill UNC Graduate Seminar, November 10, 2010 (joint

More information

Hessenberg Pairs of Linear Transformations

Hessenberg Pairs of Linear Transformations Hessenberg Pairs of Linear Transformations Ali Godjali November 21, 2008 arxiv:0812.0019v1 [math.ra] 28 Nov 2008 Abstract Let K denote a field and V denote a nonzero finite-dimensional vector space over

More information

Introduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued

Introduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued Introduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and

More information

Ergodic Theorems. Samy Tindel. Purdue University. Probability Theory 2 - MA 539. Taken from Probability: Theory and examples by R.

Ergodic Theorems. Samy Tindel. Purdue University. Probability Theory 2 - MA 539. Taken from Probability: Theory and examples by R. Ergodic Theorems Samy Tindel Purdue University Probability Theory 2 - MA 539 Taken from Probability: Theory and examples by R. Durrett Samy T. Ergodic theorems Probability Theory 1 / 92 Outline 1 Definitions

More information

AN ASYMPTOTICALLY UNBIASED MOMENT ESTIMATOR OF A NEGATIVE EXTREME VALUE INDEX. Departamento de Matemática. Abstract

AN ASYMPTOTICALLY UNBIASED MOMENT ESTIMATOR OF A NEGATIVE EXTREME VALUE INDEX. Departamento de Matemática. Abstract AN ASYMPTOTICALLY UNBIASED ENT ESTIMATOR OF A NEGATIVE EXTREME VALUE INDEX Frederico Caeiro Departamento de Matemática Faculdade de Ciências e Tecnologia Universidade Nova de Lisboa 2829 516 Caparica,

More information

Multi-operator scaling random fields

Multi-operator scaling random fields Available online at www.sciencedirect.com Stochastic Processes and their Applications 121 (2011) 2642 2677 www.elsevier.com/locate/spa Multi-operator scaling random fields Hermine Biermé a, Céline Lacaux

More information

Estimation of matrix rank: historical overview and more recent developments

Estimation of matrix rank: historical overview and more recent developments Estimation of matrix rank: historical overview and more recent developments Vladas Pipiras CEMAT, IST and University of North Carolina, Chapel Hill UTL Probability and Statistics Seminar in Lisbon November

More information

Superconcentration inequalities for centered Gaussian stationnary processes

Superconcentration inequalities for centered Gaussian stationnary processes Superconcentration inequalities for centered Gaussian stationnary processes Kevin Tanguy Toulouse University June 21, 2016 1 / 22 Outline What is superconcentration? Convergence of extremes (Gaussian case).

More information

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T 1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.

More information

SpringerBriefs in Probability and Mathematical Statistics

SpringerBriefs in Probability and Mathematical Statistics SpringerBriefs in Probability and Mathematical Statistics Editor-in-chief Mark Podolskij, Aarhus C, Denmark Series editors Nina Gantert, Münster, Germany Richard Nickl, Cambridge, UK Sandrine Péché, Paris,

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

arxiv: v2 [math.st] 20 Feb 2013

arxiv: v2 [math.st] 20 Feb 2013 Exact bounds on the closeness between the Student and standard normal distributions arxiv:1101.3328v2 [math.st] 20 Feb 2013 Contents Iosif Pinelis Department of Mathematical Sciences Michigan Technological

More information

On Stochastic Adaptive Control & its Applications. Bozenna Pasik-Duncan University of Kansas, USA

On Stochastic Adaptive Control & its Applications. Bozenna Pasik-Duncan University of Kansas, USA On Stochastic Adaptive Control & its Applications Bozenna Pasik-Duncan University of Kansas, USA ASEAS Workshop, AFOSR, 23-24 March, 2009 1. Motivation: Work in the 1970's 2. Adaptive Control of Continuous

More information

Monte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan

Monte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan Monte-Carlo MMD-MA, Université Paris-Dauphine Xiaolu Tan tan@ceremade.dauphine.fr Septembre 2015 Contents 1 Introduction 1 1.1 The principle.................................. 1 1.2 The error analysis

More information

Packing-Dimension Profiles and Fractional Brownian Motion

Packing-Dimension Profiles and Fractional Brownian Motion Under consideration for publication in Math. Proc. Camb. Phil. Soc. 1 Packing-Dimension Profiles and Fractional Brownian Motion By DAVAR KHOSHNEVISAN Department of Mathematics, 155 S. 1400 E., JWB 233,

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

GAUSSIAN MEASURE OF SECTIONS OF DILATES AND TRANSLATIONS OF CONVEX BODIES. 2π) n

GAUSSIAN MEASURE OF SECTIONS OF DILATES AND TRANSLATIONS OF CONVEX BODIES. 2π) n GAUSSIAN MEASURE OF SECTIONS OF DILATES AND TRANSLATIONS OF CONVEX BODIES. A. ZVAVITCH Abstract. In this paper we give a solution for the Gaussian version of the Busemann-Petty problem with additional

More information

Extremes and ruin of Gaussian processes

Extremes and ruin of Gaussian processes International Conference on Mathematical and Statistical Modeling in Honor of Enrique Castillo. June 28-30, 2006 Extremes and ruin of Gaussian processes Jürg Hüsler Department of Math. Statistics, University

More information

EnKF and Catastrophic filter divergence

EnKF and Catastrophic filter divergence EnKF and Catastrophic filter divergence David Kelly Kody Law Andrew Stuart Mathematics Department University of North Carolina Chapel Hill NC bkelly.com June 4, 2014 SIAM UQ 2014 Savannah. David Kelly

More information

K-ANTITHETIC VARIATES IN MONTE CARLO SIMULATION ISSN k-antithetic Variates in Monte Carlo Simulation Abdelaziz Nasroallah, pp.

K-ANTITHETIC VARIATES IN MONTE CARLO SIMULATION ISSN k-antithetic Variates in Monte Carlo Simulation Abdelaziz Nasroallah, pp. K-ANTITHETIC VARIATES IN MONTE CARLO SIMULATION ABDELAZIZ NASROALLAH Abstract. Standard Monte Carlo simulation needs prohibitive time to achieve reasonable estimations. for untractable integrals (i.e.

More information

Department of Mathematics, Nanjing University Nanjing , People s Republic of China

Department of Mathematics, Nanjing University Nanjing , People s Republic of China Proc Amer Math Soc 1382010, no 1, 37 46 SOME CONGRUENCES FOR THE SECOND-ORDER CATALAN NUMBERS Li-Lu Zhao, Hao Pan Zhi-Wei Sun Department of Mathematics, Naning University Naning 210093, People s Republic

More information