Full characterization of the fractional Poisson process
|
|
- Francine Bryan
- 5 years ago
- Views:
Transcription
1 epl draft Full characterization of the fractional Poisson process Mauro Politi,3, Taisei Kaizoji and Enrico Scalas 2,3 SSRI & Department of Economics and Business, International Christian University, 3--2 Osawa, Mitaka, Tokyo, Japan. 2 Dipartimento di Scienze e Tecnologie Avanzate,Università del Piemonte Orientale Amedeo Avogadro, Viale T. Michel, 52 Alessandria, Italy. 3 Basque Center for Applied Mathematics, Bizkaia Technology Park, Building 5, E486, Derio, Spain PACS PACS 2.5.Ey Stochastic processes 5..Ln Monte Carlo methods Abstract The fractional Poisson process (FPP) is a counting process with independent and identically distributed inter-event times following the Mittag-Leffler distribution. This process is very useful in several fields of applied and theoretical physics including models for anomalous diffusion. Contrary to the well-known Poisson process, the fractional Poisson process does not have stationary and independent increments. It is not a Lévy process and it is not a Markov process. In this letter, we present formulae for its finite-dimensional distribution functions, fully characterizing the process. These exact analytical results are compared to Monte Carlo simulations. From a loose mathematical point of view, counting processes N(A) are stochastic processes that count the random number of points in a set A. They are used in many fields of physics and other applied sciences. In this letter, we will consider one-dimensional real sets with the physical meaning of time intervals. The points will be incoming events whose duration is much smaller than the inter-event or inter-arrival waiting time. For instance, counts from a Geiger-Müller counter can be described in this way. The number of counts, N( t), in a given time interval t is known to follow the Poisson distribution P(N( t) = n) = exp( λ t) (λ t)n, () n! where λ is the constant rate of arrival of ionizing particles. Together with the assumption of independent and stationary increments, Eq. () is sufficient to ine the homogeneous Poisson process. Curiously, one of the first occurrences of this process in the scientific literature was connected to the number of casualties by horse kicks in the Prussian army cavalry []. The Poisson process is strictly related to the exponential distribution. The inter-arrival times τ i identically follow the exponential distribution and are independent random variables. This means that the Poisson process is a prototypical renewal process. A justification for the ubiquity of the Poisson process has to do with its relationship with the binomial distribution. Suppose that the time interval of interest (t, t + t) is divided into n equally spaced sub-intervals. Further assume that a counting event appears in such a sub-interval with probability p and does not appear with probability p. Then, P(N( t) = k) = Bin(k; p, n) is a binomial distribution of parameters p and n and the expected number of events in the time interval is given by E[N( t)] = np. If this expected number is kept constant for n, the binomial distribution converges to the Poisson distribution of parameter λ = E[N( t)]/ t, while, in the meantime, p. However, it can be shown that many counting processes with non-stationary increments converge to the Poisson process after a transient period. It is sufficient to require that they are renewal process (i.e. they have independent and identically distributed (iid) inter-arrival times) and that E(τ i ) <. In other words, many counting processes with non-independent and non-stationary increments behave as the Poisson process if observed long after the transient period. In recent times, it has been shown that heavy-tailed distributed inter-arrival times (for which E(τ i ) = ) do play a role in many phenomena such as blinking nano-dots [2, 3], human dynamics [4, 5] and the related inter-trade times in financial markets [6, 7]. Among the counting processes with non-stationary increments, the so-called fractional Poisson process [8], N β (t), is particularly important because it is the thinning limit of counting processes related to renewal prop-
2 Mauro Politi,3 Taisei Kaizoji Enrico Scalas2,3 cesses with power-law distributed inter-arrival times [9, ]. Moreover, it can be used to approximate anomalous diffusion ruled by space-time fractional diffusion equations [9, 6]. It is a straightforward generalization of the Poisson process ined as follows. Let {τ i } i be a sequence of independent and identically distributed positive random variables with the meaning of inter-arrival times and let their common cumulative distribution function (cdf) be F τ (t) = P(τ t) = E β ( t β ), (2) where E β ( t β ) is the one-parameter Mittag-Leffler function, E β (z), ined in the complex plane as E β (z) = n= z n Γ(nβ + ) evaluated in the point z = t β and with the prescription < β. In equation (3), Γ( ) is Euler s Gamma function. The sequence of the epochs, {T n } n, is given by the sums of the inter-arrival times n T n = τ i. (4) i (3) The epochs represent the times in which events arrive or occur. Let f τ (t) = df τ (t)/dt denote the probability density function (pdf) of the inter-arrival times, then the probability density function of the n-th epoch is simply given by the n-fold convolution of f τ (t), written as fτ n (t). In Ref. [], it is shown that f Tn (t) = fτ n (t) = β tnβ (n )! E(n) β ( tβ ), (5) where E (n) β ( tβ ) is the n-th derivative of E β (z) evaluated in z = t β. The counting process N β (t) counts the number of epochs (events) up to time t, assuming that T = is an epoch as well, or, in other words, that the process begins from a renewal point. This assumption will be used all over this paper. N β (t) is given by N β (t) = max{n : T n t}. (6) In Ref. [9], the fractional Poisson distribution is derived and it is given by P(N β (t) = n) = tβn n! E(n) β ( tβ ). (7) Eq. (7) coincides with the Poisson distribution of parameter λ = for β =. In principle, equations (3) and (7) can be directly used to derive the fractional Poisson distribution, but convergence of the series is slow. Fortunately, in a recent paper, Beghin and Orsingher proved that E (n) β ( tβ ) = n! [ exp( u)u n t βn F Sβ (t; u) (n )! ] exp( u)un du, n! (8) P(N β (t )=n ) β=.5, t β=.8, t n Figure : P (N β (T ) = n ) as function of n for three different values of β. The crosses are estimations obtained from 5 Monte Carlo samples and the lines are given to guide the eye. where F Sβ (t; u) is the cdf of a stable random variable S β (ν, γ, δ) with index β, skewness parameter ν =, scale parameter γ = (u cos πβ/2) /β and location δ = [7]. The integral in equation (8) can be evaluated numerically and Fig. shows P(N β (t) = n) for three different values of β. The Monte Carlo simulation of the fractional Poisson process is based on the algorithm presented in equation (2) of Ref. [4]. As a consequence of Kolmogorov s extension theorem, in order to fully characterize the stochastic process N β (t), one has to derive its finite dimensional distributions. A further requirement on the process paths uniquely determines the process, namely that they are right-continuous step functions with left limits [8]. The finite-dimensional distributions are the multivariate probability distribution functions P(N β (t ) = n, N β (t 2 ) = n 2,..., N β (t k ) = n k ) with t < t 2 <... < t k and n n 2... n k. We have already given the formula for the one-point functions in Eq. (7). The general finite dimensional distribution can be computed observing that the event {N β (t ) = n, N β (t 2 ) = n 2,..., N β (t k ) = n k } is equivalent to { < T n < t, T n+ > t, t < T n2 < t 2, T n2+ > p-2
3 Full characterization of the fractional Poisson process P(N β (t ) = n, N β (t 2 ) = n 2 ) = P(N β (t ) = n, N β (t 2 ) N β (t ) = n 2 n ) and, as a consequence of the inition of conditional probability Y n Y n,n2 Y n,n2,n3 T T T 2 T 3 T 4 T 5 t t 2 t 3 τ τ 2 τ 3 τ 4 τ 5 t P(N β (t ) = n, N β (t 2 ) N β (t ) = n 2 n ) = P(N β (t 2 ) N β (t ) = n 2 n N β (t ) = n ) P(N β (t ) = n ). () For β =, when the fractional Poisson process coincides with the standard Poisson process, the increments are iid random variables and one has Figure 2: (Color online) Pictorial illustration of the random variables used in the text. The light blue dots represent the observation points t, t 2 and t 3. The red squares are the epochs T =, T,..., T 5. The conditional residual life-time is the time elapsed between t i and the next epoch T ni +. It depends on previous values of n i, this is the number of events between and t i, with the event at t = T = not considered. Here, we have n =, n 2 = 2 and n 3 = 4. All the equations in this paper can be derived by analyzing this figure. t 2,..., t k < T nk < t k, T nk + > t k }. Therefore, we find P(N β (t ) = n, N β (t 2 ) = n 2,..., N β (t k ) = n k ) = P( < T n < t, T n+ > t, t < T n2 < t 2, T n2+ > t 2, t2 u u 2..., t k < T nk < t k, T nk + > t k ) = t du f n τ (u ) t u u 2 tk 2k 2 i ui... du 3 f (n2 n ) τ (u 3 ) t k 2k 2 t u du 2 f τ (u 2 ) t 2 u u 2 u 3 du 4 f τ (u 4 ) du 2k fτ (nk nk ) (u 2k ) i ui ( )] 2k [ F τ t k u i. (9) i For instance, the two point function is given by P(N β (t ) = n, N β (t 2 ) = n 2 ) = P( < T n < t, T n+ > t, t < T n2 < t 2, T n2+ > t 2 ) = t du f n τ (u ) t2 u u 2 t u u 2 t u du 2 f τ (u 2 ) du 3 f (n2 n ) τ (u 3 ) [ F τ (t 2 u u 2 u 3 )]. () Let us focus on the two-point case for the sake of illustration. As N β (t) is a counting process, one has P(N (t 2 ) N (t ) = n 2 n N (t ) = n ) = P(N (t 2 ) N (t ) = n 2 n ) = exp( (t 2 t )) (t 2 t ) (n2 n). (2) (n 2 n )! On the contrary, for < β <, the increment N β (t 2 ) N β (t ) and N β (t ) are not independent. Note that N β (t ) can be seen as an increment as N β () = by inition. However from Eq. (), the conditional probability of having n 2 n epochs in the interval (t, t 2 ) conditional on the observation of n epochs in the interval (, t ) can be written as a ratio of two finite dimensional distribution: P(N β (t 2 ) N β (t ) = n 2 n N β (t ) = n ) = P(N β (t ) = n, N β (t 2 ) = n 2 ). (3) P(N β (t ) = n ) This probability can be evaluated by means of an alternative method, more appealing for a direct and practical understanding of the dependence structure. Let Y n = [T n+ t N β (t ) = n ] (4) denote the residual lifetime at time t (that is the time to the next epoch or renewal) conditional on N β (t ) = n. With reference to Fig. 2, one can see that the conditional probability P(N β (t 2 ) N β (t ) = n 2 n N β (t ) = n ) is given by the following convolution integral for n 2 n P(N β (t 2 ) N β (t ) = n 2 n N β (t ) = n ) = t2 t P(N β (t 2 t y) = n 2 n )f Yn (y) dy, (5) where f Yn (t) is the pdf of Y n. In the case n 2 n =, one has P(N β (t 2 ) N β (t ) = N β (t ) = n ) = F Yn (t 2 t ) (6) where F Yn (y) is the cdf of Y n. The distribution of the conditional residual lifetime Y n can be evaluated in several ways. For instance, one can notice that it can be decomposed as follows Y n = τ n+ + U n (7) p-3
4 Mauro Politi,3 Taisei Kaizoji Enrico Scalas2,3 where U n is ined as and the predictive probabilities can be evaluated as U n = [T n N β (t ) = n ], (8) and is the position of the last epoch before t conditional on N β (t ) = n, and τ n+ = [τ n+ t T n+ > t ] (9) is the difference between τ n+ and t conditional on T n+ > t. The pdf of U n is given by the following chain of equalities f Un (t)dt =P(T n dt N β (t ) = n ) =P(T n dt T n < t, T n + τ n+ > t ) =P(T n dt T n < t, τ n+ > t T n ) = P(T n dt) t P(τ t n + dw) P(T n < t, τ n+ > t T n ) = f τ n (t)[ F τ (t t)]dt n dufτ (u)[ F τ (t u)], t (2) where we used the independence between T n and τ n+ ( ) and f Tn (x) = fτ n (x) ( ). The pdf of τ n+ is f τn + (t U n )dt =P(τ n+ t dt T n+ > t ) = P(τ n + dt + t ) P(τ n+ > t U n ) = f τ (t + t )dt F τ (t U n ). From Eq. (7), one can write that f Yn (t) = t and this equation leads to f Yn (t) = (2) f τn + (t u u)f U n (u)du (22) t duf n t duf n τ τ (u)f τ (t + t u) (u)[ F τ (t u)] (23) that, together with Eq. (7), gives us the probability of the conditional increments in Eq. (5). Notice that, for n =, one has fτ (u) = δ(u) and Eq. (23) reduces to the familiar equation for the residual life-time pdf in the absence of previous renewals f Y (t) = f τ (t + t ) F τ (t ). (24) This method can be applied to the general multidimensional case. As in Eq. () we can write P(N β (t ) = n,..., N β (t k ) = n k, N β (t k+ ) = n k+ ) = P(N β (t k+ ) N β (t k ) = n k+ n k N β (t ) = n,..., N β (t k ) = n k ) P(N β (t ) = n,..., N β (t k ) = n k ) (25) P(N β (t k+ ) N β (t k ) = n k+ n k... tk+ t k where we ined N β (t ) = n,..., N β (t k ) = n k ) = P(N β (t k+ t k y) = n k+ n k ) f Yn,...,n k (y)dy, (26) Y n,...,n k = [T nk + t k N β (t ) = n,..., N β (t k ) = n k ]. (27) Again, we can use a decomposition of Y n,...,n k where and U nk Y n,...,n k = τ nk + + U nk, (28) = [T nk N β (t ) = n,..., N β (t k ) = n k ], (29) τ nk + = [τ nk + t k T nk + > t k ]. (3) The difference with the two-point case is that U n = [T n N β (t ) = n ] = [ n i τ i N β (t ) = n ] must be replaced by U nk = t k + Y n,...,n k + n k i=n k + τ i N β (t k ) = n k. (3) The time between t k and the next renewal epoch is Y n,...,n k and it is independent from n k i=n k + τ i. Therefore, the convolution q(n,..., n k ; t) = f Yn,...,n k f (n k n k ) τ (t) (32) replaces fτ n (t) in Eq. (2). This leads to f Unk (z) = q(n,..., n k ; t + t k )[ F τ (t k t)] tk t k q(n,..., n k ; u + t k )[ F τ (t k u)]du. (33) On the other hand, f τnk +(t) has the same functional form as f τn (t) given in Eq. (2) with U + n k replacing U n. Therefore, Y n,...,n k has the following pdf f Yn,...,n k (t) = tk t k du q(n,..., n k ; u + t k )f τ (t + t k u) tk t k du q(n,..., n k ; u + t k )[ F τ (t k u)]. (34) In practice, the random variable Y n,...,n k carries the memory of the observations made at times t,..., t k ; the knowledge of f Yn,...,n k allows the computation of f Yn,...,n k, and, via Eqs. (25) and (26), the k + - dimensional distribution can be derived as well. p-4
5 Full characterization of the fractional Poisson process β=.5, t β=.5, t β=.8, t.9 β=.5, t β=.5, t, n.9 β=.5, t β=.7, t =7, n.9 β=.8, t β=.8, t, n.8 β=.5, t, n β=.5, t.8 β=.5, t =7, n β=.5, t.8 β=.8, t, n β=.8, t Probability density function of Un Probability density function of Un Probability density function of Un β=.8, t.9 β=.8, t β=.8, t β=.8, t β=.8, t.8.8 Probability density function of Un Probability density function of Un Probability density function of Un Figure 3: (Color online) Pdf of the random variable U n as given in Eq. (2) (solid black lines) compared to Monte Carlo simulations (colored step lines) for three values of β and two different values of t. 7 different paths were simulated for each value of β and the bin width is.5. is in arbitrary units. p-5
6 Mauro Politi,3 Taisei Kaizoji Enrico Scalas2,3.4 β=.5, t =.4 β=.5, t = β=.8, t =.35 β=.5, t β=.5, t, n.35 β=.5, t β=.5, t =7, n.9 β=.8, t β=.8, t, n Probability density function of Yn β=.5, t β=.5, t, n β=.5, t Probability density function of Yn β=.5, t β=.5, t =7, n β=.5, t Probability density function of Yn β=.8, t β=.8, t, n β=.8, t β=.8, t = = =.9 β=.8, t β=.8, t β=.8, t β=.8, t.8.8 Probability density function of Yn β=.8, t Probability density function of Yn Probability density function of Yn Figure 4: (Color online) Pdf of the random variable Y n as given in Eqs. (23) and (24) (solid black lines) compared to Monte Carlo simulations (colored step lines) for three values of β and two different values of t. 7 different paths were simulated for each value of β and the bin width is.. is in arbitrary units. p-6
7 Full characterization of the fractional Poisson process Figs. 3 and 4 compare the theoretical results of Eqs. (2), (23) and (24) with those of a Monte Carlo simulation based on the algorithm presented in equation (2) of Ref. [4]. The Japanese Society for the Promotion of Science (grant N. PE943) supported MP during his stay at the International Christian University in Tokyo, Japan. References [] Bortkiewicz L. J., Das Gesetz der kleinen Zahlen (B. G. Teubner) 898. [2] Margolin G. and Barkai E., Physical Review Letters, 94 (25) 86. [3] Margolin G., Protasenko V., Kuno M. and Barkai E., Journal of Physical Chemistry B, (26) 953. [4] Barabási A.-L., Nature, 435 (25) 27. [5] Barabási A.-L., Bursts: The Hidden Pattern Behind Everything We Do (Dutton) 2. [6] Scalas E., Gorenflo R., Luckock H., Mainardi F., Mantelli M. and Raberto M., Quantitative Finance, 4 (24) 695. [7] Scalas E., Kaizoji T., Kirchler M., Huber J. and Tedeschi A., Physica A, 366 (26) 463. [8] Laskin N., Communications in Nonlinear Science and Numerical Simulation, 8 (23) 2. [9] Scalas E., Gorenflo R. and Mainardi F., Physical Review E, 69 (24) 7. [] Mainardi F., Gorenflo R. and Scalas E., Vietnam Journal of Mathematics, 32 (24) 53. [] Metzler R. and Klafter J., Journal of Physics A: Mathematical and General, 37 (24) R6. [2] Heinsalu E., Patriarca M., Goychuk I., Schmid G. and Hänggi P., Physical Review E, 73 (26) [3] Magdziarz M., Weron A. and Weron K., Physical Review E, 75 (27) 678. [4] Fulger D., Scalas E. and Germano G., Physical Review E, 77 (28) 222. [5] Germano G., Politi M., Scalas E. and Schilling R. L., Physical Review E, 79 (29) 662. [6] Meerschaert M. M., Nane E. and Vellaisamy P., (2). [7] Beghin L. and Orsingher E., Electronic Journal of Probability, 4 (29) 79. [8] Billingsley P., Probability and Measure 2nd Edition (John Wiley & Sons) 986. p-7
Covariance structure of continuous time random walk limit processes
Covariance structure of continuous time random walk limit processes Alla Sikorskii Department of Statistics and Probability Michigan State University Fractional Calculus, Probability and Non-local Operators:
More informationVolatility and Returns in Korean Futures Exchange Markets
Volatility and Returns in Korean Futures Exchange Markets Kyungsik Kim*,, Seong-Min Yoon and Jum Soo Choi Department of Physics, Pukyong National University, Pusan 608-737, Korea Division of Economics,
More informationSummary of Results on Markov Chains. Abstract
Summary of Results on Markov Chains Enrico Scalas 1, 1 Laboratory on Complex Systems. Dipartimento di Scienze e Tecnologie Avanzate, Università del Piemonte Orientale Amedeo Avogadro, Via Bellini 25 G,
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationON THE CONVERGENCE OF QUADRATIC VARIATION FOR COMPOUND FRACTIONAL POISSON PROCESSES IN FCAA JOURNAL
ON THE CONVERGENCE OF QUADRATIC VARIATION FOR COMPOUND FRACTIONAL POISSON PROCESSES IN FCAA JOURNAL Enrico Scalas 1, Noèlia Viles 2 Abstract The relationship between quadratic variation for compound renewal
More informationarxiv: v1 [cond-mat.stat-mech] 23 Apr 2010
Fractional Fokker-Planck Equations for Subdiffusion with Space-and-Time-Dependent Forces. B.. Henry Department of Applied Mathematics, School of Mathematics and Statistics, University of New South Wales,
More informationCRRAO Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS) Research Report. B. L. S. Prakasa Rao
CRRAO Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS) Research Report Author (s): B. L. S. Prakasa Rao Title of the Report: FILTERED FRACTIONAL POISSON ROCESSES Research Report
More informationPoisson Processes. Stochastic Processes. Feb UC3M
Poisson Processes Stochastic Processes UC3M Feb. 2012 Exponential random variables A random variable T has exponential distribution with rate λ > 0 if its probability density function can been written
More informationPart I Stochastic variables and Markov chains
Part I Stochastic variables and Markov chains Random variables describe the behaviour of a phenomenon independent of any specific sample space Distribution function (cdf, cumulative distribution function)
More informationCoupled continuous time random walks in finance
Physica A 37 (26) 114 118 www.elsevier.com/locate/physa Coupled continuous time random walks in finance Mark M. Meerschaert a,,1, Enrico Scalas b,2 a Department of Mathematics and Statistics, University
More informationLecture 4a: Continuous-Time Markov Chain Models
Lecture 4a: Continuous-Time Markov Chain Models Continuous-time Markov chains are stochastic processes whose time is continuous, t [0, ), but the random variables are discrete. Prominent examples of continuous-time
More informationApplied Probability and Stochastic Processes
Applied Probability and Stochastic Processes In Engineering and Physical Sciences MICHEL K. OCHI University of Florida A Wiley-Interscience Publication JOHN WILEY & SONS New York - Chichester Brisbane
More informationSome Generalizations of Poisson Processes
Journal of Statistical Theory and Applications Volume 11, Number 3, 212, pp. 225-235 ISSN 1538-7887 Some Generalizations of Poisson Processes Seetha Lekshmi, V. Dept.of Statistics, Nirmala College Muvattupuzha,
More informationRENEWAL PROCESSES AND POISSON PROCESSES
1 RENEWAL PROCESSES AND POISSON PROCESSES Andrea Bobbio Anno Accademico 1997-1998 Renewal and Poisson Processes 2 Renewal Processes A renewal process is a point process characterized by the fact that the
More informationPoissonprocessandderivationofBellmanequations
APPENDIX B PoissonprocessandderivationofBellmanequations 1 Poisson process Let us first define the exponential distribution Definition B1 A continuous random variable X is said to have an exponential distribution
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics. Probability Distributions Prof. Dr. Klaus Reygers (lectures) Dr. Sebastian Neubert (tutorials) Heidelberg University WS 07/8 Gaussian g(x; µ, )= p exp (x µ) https://en.wikipedia.org/wiki/normal_distribution
More informationWeak Ergodicity Breaking WCHAOS 2011
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Bel, Burov, Korabel, Margolin, Rebenshtok WCHAOS 211 Outline Single molecule experiments exhibit weak ergodicity breaking. Blinking quantum dots,
More informationRecap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks
Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution
More informationCDA6530: Performance Models of Computers and Networks. Chapter 3: Review of Practical Stochastic Processes
CDA6530: Performance Models of Computers and Networks Chapter 3: Review of Practical Stochastic Processes Definition Stochastic process X = {X(t), t2 T} is a collection of random variables (rvs); one rv
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationPhysica A. Extremal behavior of a coupled continuous time random walk
Physica A 39 (211) 55 511 Contents lists available at ScienceDirect Physica A journal homepage: www.elsevier.com/locate/physa Extremal behavior of a coupled continuous time random walk Rina Schumer a,,
More informationExponential Distribution and Poisson Process
Exponential Distribution and Poisson Process Stochastic Processes - Lecture Notes Fatih Cavdur to accompany Introduction to Probability Models by Sheldon M. Ross Fall 215 Outline Introduction Exponential
More informationCauchy Problems in Bounded Domains and Iterated Processes
Cauchy Problems in Bounded Domains and Iterated Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J
More informationLecture 1: Brief Review on Stochastic Processes
Lecture 1: Brief Review on Stochastic Processes A stochastic process is a collection of random variables {X t (s) : t T, s S}, where T is some index set and S is the common sample space of the random variables.
More informationScientific Measurement
Scientific Measurement SPA-4103 Dr Alston J Misquitta Lecture 5 - The Binomial Probability Distribution Binomial Distribution Probability of k successes in n trials. Gaussian/Normal Distribution Poisson
More informationMarkov Chains. X(t) is a Markov Process if, for arbitrary times t 1 < t 2 <... < t k < t k+1. If X(t) is discrete-valued. If X(t) is continuous-valued
Markov Chains X(t) is a Markov Process if, for arbitrary times t 1 < t 2
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationFrom Renewal Theory To High-Frequency Finance
From Renewal Theory To High-Frequency Finance Mauro Politi * Enrico Scalas ** * Universitá degli studi di Milano Dipartimento di Fisica, via Celoria 16-2133 Milano ITALY mauro.politi@unimi.it ** Universitá
More informationarxiv:cond-mat/ v2 [cond-mat.dis-nn] 11 Oct 2002
arxiv:cond-mat/21166v2 [cond-mat.dis-nn] 11 Oct 22 REVISITING THE DERIVATION OF THE FRACTIONAL DIFFUSION EQUATION ENRICO SCALAS Dipartimento di Scienze e Tecnologie Avanzate, Università del Piemonte Orientale
More informationWeek 1 Quantitative Analysis of Financial Markets Distributions A
Week 1 Quantitative Analysis of Financial Markets Distributions A Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
More informationChapter 2. Poisson Processes. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan
Chapter 2. Poisson Processes Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan Outline Introduction to Poisson Processes Definition of arrival process Definition
More informationSUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION
Vol. 4 (29) ACTA PHYSICA POLONICA B No 5 SUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION Joanna Janczura, Agnieszka Wyłomańska Institute of Mathematics and Computer Science Hugo Steinhaus
More informationContinuous-time Markov Chains
Continuous-time Markov Chains Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ October 23, 2017
More informationCONTINUOUS TIME RANDOM WALK WITH CORRELATED WAITING TIMES
CONTINUOUS TIME RANDOM WALK WITH CORRELATED WAITING TIMES Aleksei V. Chechkin, 1,2 Michael Hofmann 3, and Igor M. Sokolov 3 1 School of Chemistry, Tel Aviv University, Ramat Aviv, 69978 Tel Aviv, Israel
More informationModeling point processes by the stochastic differential equation
Modeling point processes by the stochastic differential equation Vygintas Gontis, Bronislovas Kaulakys, Miglius Alaburda and Julius Ruseckas Institute of Theoretical Physics and Astronomy of Vilnius University,
More informationPerformance Modelling of Computer Systems
Performance Modelling of Computer Systems Mirco Tribastone Institut für Informatik Ludwig-Maximilians-Universität München Fundamentals of Queueing Theory Tribastone (IFI LMU) Performance Modelling of Computer
More informationAn Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems
An Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems M.P. Kaminskiy and V.V. Krivtsov Abstract This paper introduces a simple index that helps to assess the degree of aging
More informationChapter 2 SOME ANALYTICAL TOOLS USED IN THE THESIS
Chapter 2 SOME ANALYTICAL TOOLS USED IN THE THESIS 63 2.1 Introduction In this chapter we describe the analytical tools used in this thesis. They are Markov Decision Processes(MDP), Markov Renewal process
More informationContinuous Time Random Walk Limits in Bounded Domains
Continuous Time Random Walk Limits in Bounded Domains Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY May 16, 2012 Acknowledgement Boris Baeumer, University of Otago, New Zealand.
More informationIntroduction to Rare Event Simulation
Introduction to Rare Event Simulation Brown University: Summer School on Rare Event Simulation Jose Blanchet Columbia University. Department of Statistics, Department of IEOR. Blanchet (Columbia) 1 / 31
More informationLecture 25: Large Steps and Long Waiting Times
Lecture 25: Large Steps and Long Waiting Times Scribe: Geraint Jones (and Martin Z. Bazant) Department of Economics, MIT Proofreader: Sahand Jamal Rahi Department of Physics, MIT scribed: May 10, 2005,
More informationBasics of Stochastic Modeling: Part II
Basics of Stochastic Modeling: Part II Continuous Random Variables 1 Sandip Chakraborty Department of Computer Science and Engineering, INDIAN INSTITUTE OF TECHNOLOGY KHARAGPUR August 10, 2016 1 Reference
More informationChapter 3 sections. SKIP: 3.10 Markov Chains. SKIP: pages Chapter 3 - continued
Chapter 3 sections Chapter 3 - continued 3.1 Random Variables and Discrete Distributions 3.2 Continuous Distributions 3.3 The Cumulative Distribution Function 3.4 Bivariate Distributions 3.5 Marginal Distributions
More informationChapter 3 sections. SKIP: 3.10 Markov Chains. SKIP: pages Chapter 3 - continued
Chapter 3 sections 3.1 Random Variables and Discrete Distributions 3.2 Continuous Distributions 3.3 The Cumulative Distribution Function 3.4 Bivariate Distributions 3.5 Marginal Distributions 3.6 Conditional
More informationThe exponential distribution and the Poisson process
The exponential distribution and the Poisson process 1-1 Exponential Distribution: Basic Facts PDF f(t) = { λe λt, t 0 0, t < 0 CDF Pr{T t) = 0 t λe λu du = 1 e λt (t 0) Mean E[T] = 1 λ Variance Var[T]
More informationCDA5530: Performance Models of Computers and Networks. Chapter 3: Review of Practical
CDA5530: Performance Models of Computers and Networks Chapter 3: Review of Practical Stochastic Processes Definition Stochastic ti process X = {X(t), t T} is a collection of random variables (rvs); one
More informationChapter 4 Continuous Random Variables and Probability Distributions
Chapter 4 Continuous Random Variables and Probability Distributions Part 3: The Exponential Distribution and the Poisson process Section 4.8 The Exponential Distribution 1 / 21 Exponential Distribution
More informationHigh-frequency data modelling using Hawkes processes
High-frequency data modelling using Hawkes processes Valérie Chavez-Demoulin 1 joint work J.A McGill 1 Faculty of Business and Economics, University of Lausanne, Switzerland Boulder, April 2016 Boulder,
More informationPROBABILITY DISTRIBUTIONS
Review of PROBABILITY DISTRIBUTIONS Hideaki Shimazaki, Ph.D. http://goo.gl/visng Poisson process 1 Probability distribution Probability that a (continuous) random variable X is in (x,x+dx). ( ) P x < X
More informationDepartment of Mathematics
Department of Mathematics Ma 3/103 KC Border Introduction to Probability and Statistics Winter 2017 Supplement 2: Review Your Distributions Relevant textbook passages: Pitman [10]: pages 476 487. Larsen
More informationHigh-frequency data modelling using Hawkes processes
Valérie Chavez-Demoulin joint work with High-frequency A.C. Davison data modelling and using A.J. Hawkes McNeil processes(2005), J.A EVT2013 McGill 1 /(201 High-frequency data modelling using Hawkes processes
More informationAsymptotic distribution of the sample average value-at-risk
Asymptotic distribution of the sample average value-at-risk Stoyan V. Stoyanov Svetlozar T. Rachev September 3, 7 Abstract In this paper, we prove a result for the asymptotic distribution of the sample
More informationComposite Fractional Time Evolutions of Cumulants in Fractional Generalized Birth Processes
Adv. Studies Theor. Phys., Vol. 8, 214, no. 5, 195-213 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/astp.214.312143 Composite Fractional Time Evolutions of Cumulants in Fractional Generalized
More informationModelling the risk process
Modelling the risk process Krzysztof Burnecki Hugo Steinhaus Center Wroc law University of Technology www.im.pwr.wroc.pl/ hugo Modelling the risk process 1 Risk process If (Ω, F, P) is a probability space
More informationEli Barkai. Wrocklaw (2015)
1/f Noise and the Low-Frequency Cutoff Paradox Eli Barkai Bar-Ilan University Niemann, Kantz, EB PRL (213) Theory Sadegh, EB, Krapf NJP (214) Experiment Stefani, Hoogenboom, EB Physics Today (29) Wrocklaw
More informationBMIR Lecture Series on Probability and Statistics Fall, 2015 Uniform Distribution
Lecture #5 BMIR Lecture Series on Probability and Statistics Fall, 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University s 5.1 Definition ( ) A continuous random
More informationUniversal examples. Chapter The Bernoulli process
Chapter 1 Universal examples 1.1 The Bernoulli process First description: Bernoulli random variables Y i for i = 1, 2, 3,... independent with P [Y i = 1] = p and P [Y i = ] = 1 p. Second description: Binomial
More informationDepartment of Mathematics
Department of Mathematics Ma 3/103 KC Border Introduction to Probability and Statistics Winter 2018 Supplement 2: Review Your Distributions Relevant textbook passages: Pitman [10]: pages 476 487. Larsen
More informationOn the Nonstationarity of the Exchange Rate Process
JSPS Grants-in-Aid for Creative Scientific Research Understanding Inflation Dynamics of the Japanese Economy Working Paper Series No.65 On the Nonstationarity of the Exchange Rate Process Takaaki Ohnishi
More informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More informationCauchy Problems Solved by Running Subordinate Processes
Cauchy Problems Solved by Running Subordinate Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J 4
More informationarxiv: v1 [q-fin.st] 18 Sep 2017
A new approach to the modeling of financial volumes arxiv:179.583v1 [q-fin.st] 18 Sep 17 Guglielmo D Amico Department of Pharmacy, University G. d Annunzio of Chieti-Pescara e-mail: g.damico@unich.it Filippo
More informationLecturer: Olga Galinina
Renewal models Lecturer: Olga Galinina E-mail: olga.galinina@tut.fi Outline Reminder. Exponential models definition of renewal processes exponential interval distribution Erlang distribution hyperexponential
More informationStochastic Renewal Processes in Structural Reliability Analysis:
Stochastic Renewal Processes in Structural Reliability Analysis: An Overview of Models and Applications Professor and Industrial Research Chair Department of Civil and Environmental Engineering University
More informationLecture 2. Distributions and Random Variables
Lecture 2. Distributions and Random Variables Igor Rychlik Chalmers Department of Mathematical Sciences Probability, Statistics and Risk, MVE300 Chalmers March 2013. Click on red text for extra material.
More informationarxiv: v1 [math.pr] 18 Nov 2018
Integral Equation Approach to Stationary Stochastic Counting Process with Independent Increments arxiv:1811.7262v1 [math.pr] 18 Nov 218 Enzhi Li Suning R & D Center, Palo Alto, USA November 2, 218 Abstract:
More informationWeak Ergodicity Breaking
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University 215 Ergodicity Ergodicity: time averages = ensemble averages. x = lim t t x(τ)dτ t. x = xp eq (x)dx. Ergodicity out of equilibrium δ 2 (, t) = t
More informationarxiv: v1 [math.st] 31 Jul 2013
A FRACTIONAL GENERALIZATION OF THE POISSON PROCESSES AND SOME OF ITS PROPERTIES Nicy Sebastian Indian Statistical Institute, Chennai Centre, Taramani, Chennai - 600113, India nicy@isichennai.res.in Rudolf
More informationChing-Han Hsu, BMES, National Tsing Hua University c 2015 by Ching-Han Hsu, Ph.D., BMIR Lab. = a + b 2. b a. x a b a = 12
Lecture 5 Continuous Random Variables BMIR Lecture Series in Probability and Statistics Ching-Han Hsu, BMES, National Tsing Hua University c 215 by Ching-Han Hsu, Ph.D., BMIR Lab 5.1 1 Uniform Distribution
More information1 Delayed Renewal Processes: Exploiting Laplace Transforms
IEOR 6711: Stochastic Models I Professor Whitt, Tuesday, October 22, 213 Renewal Theory: Proof of Blackwell s theorem 1 Delayed Renewal Processes: Exploiting Laplace Transforms The proof of Blackwell s
More informationComplex Systems Methods 11. Power laws an indicator of complexity?
Complex Systems Methods 11. Power laws an indicator of complexity? Eckehard Olbrich e.olbrich@gmx.de http://personal-homepages.mis.mpg.de/olbrich/complex systems.html Potsdam WS 2007/08 Olbrich (Leipzig)
More informationResearch Reports on Mathematical and Computing Sciences
ISSN 1342-2804 Research Reports on Mathematical and Computing Sciences Long-tailed degree distribution of a random geometric graph constructed by the Boolean model with spherical grains Naoto Miyoshi,
More informationConvergence Rates for Renewal Sequences
Convergence Rates for Renewal Sequences M. C. Spruill School of Mathematics Georgia Institute of Technology Atlanta, Ga. USA January 2002 ABSTRACT The precise rate of geometric convergence of nonhomogeneous
More informationMonotonicity and Aging Properties of Random Sums
Monotonicity and Aging Properties of Random Sums Jun Cai and Gordon E. Willmot Department of Statistics and Actuarial Science University of Waterloo Waterloo, Ontario Canada N2L 3G1 E-mail: jcai@uwaterloo.ca,
More information1.1 Review of Probability Theory
1.1 Review of Probability Theory Angela Peace Biomathemtics II MATH 5355 Spring 2017 Lecture notes follow: Allen, Linda JS. An introduction to stochastic processes with applications to biology. CRC Press,
More informationBirth-Death Processes
Birth-Death Processes Birth-Death Processes: Transient Solution Poisson Process: State Distribution Poisson Process: Inter-arrival Times Dr Conor McArdle EE414 - Birth-Death Processes 1/17 Birth-Death
More informationAN INTEGRAL MEASURE OF AGING/REJUVENATION FOR REPAIRABLE AND NON REPAIRABLE SYSTEMS
R&RAA # 1 (Vol.1) 8, March AN INEGRAL MEASURE OF AGING/REJUVENAION FOR REPAIRABLE AND NON REPAIRABLE SYSEMS M.P. Kaminskiy and V.V. Krivtsov Abstract his paper introduces a simple index that helps to assess
More informationarxiv: v1 [math.pr] 19 Nov 2017
Article A note on Hadamard fractional differential equations with varying coefficients and their applications in probability Roberto Garra 1,, Enzo Orsingher 1 and Federico Polito 2 arxiv:1711.07016v1
More informationOrnstein-Uhlenbeck processes for geophysical data analysis
Ornstein-Uhlenbeck processes for geophysical data analysis Semere Habtemicael Department of Mathematics North Dakota State University May 19, 2014 Outline 1 Introduction 2 Model 3 Characteristic function
More informationLecture 5. If we interpret the index n 0 as time, then a Markov chain simply requires that the future depends only on the present and not on the past.
1 Markov chain: definition Lecture 5 Definition 1.1 Markov chain] A sequence of random variables (X n ) n 0 taking values in a measurable state space (S, S) is called a (discrete time) Markov chain, if
More informationPoisson Processes for Neuroscientists
Poisson Processes for Neuroscientists Thibaud Taillefumier This note is an introduction to the key properties of Poisson processes, which are extensively used to simulate spike trains. For being mathematical
More informationFractional Cauchy Problems for time-changed Processes
Fractional Cauchy Problems for time-changed Processes Erkan Nane Department of Mathematics and Statistics Auburn University February 17, 211 Erkan Nane (Auburn University) Fractional Cauchy Problems February
More informationT. Liggett Mathematics 171 Final Exam June 8, 2011
T. Liggett Mathematics 171 Final Exam June 8, 2011 1. The continuous time renewal chain X t has state space S = {0, 1, 2,...} and transition rates (i.e., Q matrix) given by q(n, n 1) = δ n and q(0, n)
More informationRandomness at the root of things 2: Poisson sequences
SPECIAL FEATURE: FUZZY PHYSICS www.iop.org/journals/physed Randomness at the root of things 2: Poisson sequences Jon Ogborn 1, Simon Collins 2 and Mick Brown 3 1 Institute of Education, University of London,
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More information2905 Queueing Theory and Simulation PART III: HIGHER DIMENSIONAL AND NON-MARKOVIAN QUEUES
295 Queueing Theory and Simulation PART III: HIGHER DIMENSIONAL AND NON-MARKOVIAN QUEUES 16 Queueing Systems with Two Types of Customers In this section, we discuss queueing systems with two types of customers.
More informationIEOR 6711: Stochastic Models I, Fall 2003, Professor Whitt. Solutions to Final Exam: Thursday, December 18.
IEOR 6711: Stochastic Models I, Fall 23, Professor Whitt Solutions to Final Exam: Thursday, December 18. Below are six questions with several parts. Do as much as you can. Show your work. 1. Two-Pump Gas
More informationIrreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1
Irreducibility Irreducible every state can be reached from every other state For any i,j, exist an m 0, such that i,j are communicate, if the above condition is valid Irreducible: all states are communicate
More informationContinuous-time random walks, fractional calculus and stochastic integrals
Continuous-time random walks, fractional calculus and stochastic integrals A model for high-frequency financial time series G. Germano 1 M. Politi 1,2 E. Scalas 3 R. L. Schilling 4 1 Philipps-Universität
More informationEE126: Probability and Random Processes
EE126: Probability and Random Processes Lecture 18: Poisson Process Abhay Parekh UC Berkeley March 17, 2011 1 1 Review 2 Poisson Process 2 Bernoulli Process An arrival process comprised of a sequence of
More informationTCOM 501: Networking Theory & Fundamentals. Lecture 6 February 19, 2003 Prof. Yannis A. Korilis
TCOM 50: Networking Theory & Fundamentals Lecture 6 February 9, 003 Prof. Yannis A. Korilis 6- Topics Time-Reversal of Markov Chains Reversibility Truncating a Reversible Markov Chain Burke s Theorem Queues
More informationarxiv: v1 [math.pr] 4 Feb 2016
Mittag-Leffler Lévy Processes Arun Kumar an N. S. Upahye *Inian Statistical Institute, Chennai Center, Taramani, Chennai-636, Inia an **Department of Mathematics, Inian Institute of Technology Maras, Chennai
More informationLECTURE #6 BIRTH-DEATH PROCESS
LECTURE #6 BIRTH-DEATH PROCESS 204528 Queueing Theory and Applications in Networks Assoc. Prof., Ph.D. (รศ.ดร. อน นต ผลเพ ม) Computer Engineering Department, Kasetsart University Outline 2 Birth-Death
More information(b) What is the variance of the time until the second customer arrives, starting empty, assuming that we measure time in minutes?
IEOR 3106: Introduction to Operations Research: Stochastic Models Fall 2006, Professor Whitt SOLUTIONS to Final Exam Chapters 4-7 and 10 in Ross, Tuesday, December 19, 4:10pm-7:00pm Open Book: but only
More informationHANDBOOK OF APPLICABLE MATHEMATICS
HANDBOOK OF APPLICABLE MATHEMATICS Chief Editor: Walter Ledermann Volume II: Probability Emlyn Lloyd University oflancaster A Wiley-Interscience Publication JOHN WILEY & SONS Chichester - New York - Brisbane
More informationStat 516, Homework 1
Stat 516, Homework 1 Due date: October 7 1. Consider an urn with n distinct balls numbered 1,..., n. We sample balls from the urn with replacement. Let N be the number of draws until we encounter a ball
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 6
MATH 56A: STOCHASTIC PROCESSES CHAPTER 6 6. Renewal Mathematically, renewal refers to a continuous time stochastic process with states,, 2,. N t {,, 2, 3, } so that you only have jumps from x to x + and
More informationBrief Review of Probability
Maura Department of Economics and Finance Università Tor Vergata Outline 1 Distribution Functions Quantiles and Modes of a Distribution 2 Example 3 Example 4 Distributions Outline Distribution Functions
More informationLecture 10: Semi-Markov Type Processes
Lecture 1: Semi-Markov Type Processes 1. Semi-Markov processes (SMP) 1.1 Definition of SMP 1.2 Transition probabilities for SMP 1.3 Hitting times and semi-markov renewal equations 2. Processes with semi-markov
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More information