Department of Applied Mathematics Faculty of EEMCS. University of Twente. Memorandum No Birth-death processes with killing
|
|
- Elfrieda Gordon
- 5 years ago
- Views:
Transcription
1 Department of Applied Mathematics Faculty of EEMCS t University of Twente The Netherlands P.O. Box AE Enschede The Netherlands Phone: Fax: memo@math.utwente.nl Memorandum No. 75 Birth-death processes with killing E.A. van Doorn and A.I. Zeifman March, 24 ISSN Vologda State Pedagogical University, and Vologda Scientific Coordinate Centre of CEMI RAS, S. Orlova 6, Vologda, Russia
2 BIRTH-DEATH PROCESSES WITH KILLING Erik A. van Doorn and Alexander I. Zeifman Department of Applied Mathematics University of Twente P.O. Box 27, 75 AE Enschede, The Netherlands Vologda State Pedagogical University and Vologda Scientific Coordinate Centre of CEMI RAS S. Orlova 6, Vologda, Russia March 7, 24 Abstract. The purpose of this note is to point out that Karlin and McGregor s integral representation for the transition probabilities of a birth-death process on a semi-infinite lattice with an absorbing bottom state remains valid if one allows the possibility of absorption into the bottom state from any other state. Conditions for uniqueness of the minimal transition function are also given. Keywords and phrases: Karlin-McGregor representation, orthogonal polynomials, transition function, transition probabilities, state-dependent killing rate, total catastrophe 2 Mathematics Subject Classification: Primary 6J8
3 Introduction We consider a Markov chain X {X(t), t } taking values in S {,,...} with q-matrix Q (q ij, i, j S) given by q i,i+ = λ i, q i+,i = µ i+, q ii = (λ i + µ i + γ i ), i, q ij =, i j >, () where λ i > and γ i for i, µ i > for i >, and µ =. The parameters λ i and µ i are the birth and death rates in state i, while γ i may be regarded as the rate of absorption, or killing, into a fictitious state, say. A transition to the absorbing state is sometimes referred to in the literature as a total catastrophe (see for example Chao and Zheng (23)). We will assume that the (standard) transition function P (.) {p ij (.), i, j S}, where p ij (t) Pr{X(t) = j X() = i}, t, i, j S, is the minimal Q-function, that is, the minimal transition function with q- matrix Q. As a consequence P (.) satisfies the system P (t) = QP (t) = P (t)q, t, (2) of backward and forward equations (see for example Anderson (99)). We note parenthetically that imposing the forward equations is equivalent to imposing continuity at infinity, in the sense that almost surely the only discontinuities of the process are simple discontinuities with saltus ±, even if the process goes to infinity in finite time and returns from there (see Karlin and McGregor (959)). A prominent role in what follows will be played by the polynomials {R n } which are uniquely determined by the transition rates of X through the recurrence relation λ n R n+ (x) = (λ n + µ n + γ n x)r n (x) µ n R n (x), n, λ R (x) = λ + γ x, R (x) =. (3) Karlin and McGregor (957) have shown that if the killing rates γ i are all zero except γ (in their notation γ i = for all i, but µ, and the
4 absorbing state is labeled ), then the transition probabilities p ij (t) may be represented in the form p ij (t) = π j e xt R i (x)r j (x)ψ(dx), t, i, j S. (4) Here π n are constants given by π and π n λ λ... λ n µ µ 2... µ n, n >, (5) and ψ is a Borel measure of total mass on [, ) with respect to which the polynomials {R n } are orthogonal. The orthogonalizing measure for {R n } is in virtual all practical cases unique, and in many cases known explicitly. The usefulness of the integral representation (4) derives from the monotonic properties of e xt, and from the fact that the dependence on t, i and j is factored in the integrand. If we allow for killing rates γ i in each state i S, the polynomials {R n } of (3) still constitute an orthogonal polynomial sequence with respect to a Borel measure ψ on [, ). The main purpose of this note is to point out that also the representation (4) remains valid in this case. The remainder of this paper is organized as follows. In Section 2 we will gather some preliminary results, which are needed in Section 3 to prove our claim. In Section 4 we will investigate under which condition on Q the minimal Q-function P (.) is the only Q-function satisfying both backward and forward equations. We conclude in Section 5 with an example. 2 An associated birth-death process We shall have use for the quantities λ i and µ i, i, recurrently defined by µ =, λ = λ + γ µ i = (λ i / λ i )µ i, λi = λ i + γ i + µ i µ i, i >. (6) It is easily seen by induction that λ i λ i + γ i λ i > and µ i+ >, i, 2
5 so that λ i and µ i may be interpreted as the birth and death rates, respectively, of a (conservative) birth-death process X { X(t), t } on S with q-matrix Q ( q ij ) given by q i,i+ = λ i, q i+,i = µ i+, q ii = ( λ i + µ i ), i, q ij =, i j >. (7) We let π and π n λ λ... λ n µ µ 2... µ n, n >, (8) and observe from (), (6) and (7) that Π /2 Q Π /2 = Π /2 QΠ /2, (9) where Π /2 and Π /2 denote the diagonal matrices with entries π n and π n, respectively, on the diagonals. Evidently, X need not be uniquely determined by Q, but in what follows we will assume that X is the minimal Q-process, represented by the minimal Q-function P (.) { p ij (.), i, j S}, which therefore satisfies the system P (t) = Q P (t) = P (t) Q, t, () of backward and forward equations. Also, Karlin and McGregor s integral representation applies to p ij (t), that is, p ij (t) = π j e xt Ri (x) R j (x) ψ(dx), t, i, j S, () where { R n } are the polynomials satisfying the recurrence λ n Rn+ (x) = ( λ n + µ n x) R n (x) µ n Rn (x), n, λ R (x) = λ x, R (x) =, (2) and ψ is a Borel measure (of total mass ) on [, ) with respect to which the { R n } are orthogonal. Actually, if Q is such that ( π n + ) =, (3) λ n π n n= then the Stieltjes moment problem associated with { R n } is determined, which means that ψ is the unique orthogonalizing measure on [, ) for { R n }. In this 3
6 case P (.) is the unique Q-function satisfying the system () of backward and forward equations. If the series in (3) converges, then the Stieltjes moment problem is indeterminate, that is, there are infinitely many orthogonalizing measures for { R n }. In this case there are also infinitely many Q-functions satisfying (). The measure ψ corresponding to the minimal Q-function may now be characterized as the one which is supported by the zeros of the (entire) function R (x) lim n Rn (x); we will refer to ψ as the minimal measure for { R n } (see Karlin and McGregor (957) for more details). It is sometimes desirable for X to be non-explosive, which requires Q to be such that n= λ n π n n π i =, (4) i= which is stronger than (3). If (and only if) condition (4) is satisfied, then the minimal Q-function P (.) is in fact the unique Q-function satisfying just the backward equations P (t) = Q P (t). We conclude this section with the observation that, apart from a multiplicative constant, the polynomials R n and R n are identical. Indeed, it follows readily by induction from (6) and the recurrence relations (3) and (2) that R n (x) = λ λ... λ n πn λ λ... λ R n (x) = R n (x), n. (5) n π n So, the polynomials {R n } of (3) are orthogonal with respect to any measure 3 Representation Our main result is expressed in the corollary to the following theorem. Theorem The minimal Q-function P (.) {p ij (.)} and the minimal which is an orthogonalizing measure for { R n }, and with respect to ψ in particular. Qfunction P (.) { p ij (.)} satisfy P (t) = Π /2 Π/2 P (t) Π /2 Π /2, t. (6) 4
7 Proof Denoting the right-hand side of (6) by F (t), it follows from (9) and () that F (.) satisfies the system of backward and forward equations (2). F (.) is also the minimal non-negative solution of (2) since the opposite would imply that P (.) is not the minimal non-negative solution of (), and therefore (see for example Anderson, 99, Theorem 2.2.2) not the minimal Q-function. Hence, by the same theorem in Anderson (99), F (.) is the minimal Q-function, that is, F (t) = P (t), t. Corollary 2 The minimal Q-function P (.) {p ij (.)} of the process X may be represented in the form (4), where {R n } are the polynomials defined in (3) and ψ is the orthogonalizing measure (of total mass ) for {R n } on [, ) which, if it is not unique, is the minimal measure. Proof By (6) and () we have π i π j πi p ij (t) = p ij (t) = π j π j π i π j π i e xt Ri (x) R j (x) ψ(dx). Substituting (5) and noting, in view of (5), that {R n } is orthogonal with respect to ψ = ψ yields the result. So, as announced, we may conclude that the representation (4) remains valid whether the killing rates γ i, i >, are zero or not. As in Karlin and McGregor (957), certain non-minimal Q-functions (if there are any) may also be represented in the form (4), with ψ replaced by an appropriate (non-minimal) orthogonalizing measure for {R n }. A necessary and sufficient condition for P (.) to be the only Q-function satisfying (2) will be derived in the next section, but otherwise we will not pursue this issue. As an aside we note that the concept of similarity for birth-death processes introduced in Lenin et al. (2) may be extended to birth-death processes with killing. Namely, in view of (6) the transition probability functions p ij (.) and p ij (.) differ only by a multiplicative constant. Hence the process X is similar (in the sense of Lenin et al. (2)) to the pure birth-death process X. 5
8 4 Uniqueness If condition (3) is satisfied then the minimal Q-function P (.) is the unique Q-function satisfying (), and, as we shall see in this section, P (.) is the unique Q-function satisfying (2). On the other hand, (3) is not necessary for uniqueness of P (.). Before giving a necessary and sufficient condition we need some preliminary results. First, we observe from (5) and (6) that for all x λ n π n Rn (x) R n+ (x) = λ n π n R n (x)r n+ (x), n. (7) It is easily seen from (2) that R n () = for all n, so, as a consequence of (5) and (7), we have π n = π n Rn(), 2 n, (8) and λ n π n = λ n π n R n ()R n+ (), n. (9) It will also be useful to note from the recurrence relation (3) and the fact that λ n π n = µ n π n that, for n, λ n π n (R n+ () R n ()) = λ n π n (R n () R n ()) + γ n π n R n (), while λ π (R () R ()) = γ π R (), so that and n λ n π n (R n+ () R n ()) = γ k π k R k (), n, (2) n R n+ () = + λ j π j j= k= j γ k π k R k (), n. (2) k= Our final preliminary result is the following lemma, which is the generalization to the setting at hand of (part of) Lemma 6 (on p. 526) of Karlin and McGregor (957). Recall that ψ (= ψ) is either the unique or else the minimal orthogonalizing measure for { R n }, and hence for {R n }. 6
9 Lemma 3 We have γ j π j j= x R j (x)ψ(dx) = lim n R n (). (22) Proof From Karlin and McGregor, 957, Lemma 6 we know that x ψ(dx) = n= λ n π n, where both members may be infinite. It subsequently follows by induction from the recurrence relation (2) that x Rj (x) ψ(dx) = n=j λ n π n, j. (23) Hence, with the help of (5), (9), (8), and (2) we can write γ j π j x R j (x)ψ(dx) = j= = = γ j πj π j x Rj (x) ψ(dx) j= γ j πj π j j= j= n=j n=j λ n π n γ j π j R j () λ n π n R n ()R n+ () which yields the required result. = = n= n= λ n π n R n ()R n+ () n γ j π j R j () j= ( ) R n (), R n+ () We can now give a necessary and sufficient condition for P (.), the minimal Q-function, to be the unique Q-function satisfying (2). Theorem 4 In order that there be only one Q-function P (.) satisfying the combined system of backward and forward equations it is necessary and sufficient that at least one of the two conditions lim R n() = (24) n 7
10 or n= be satisfied. ( π n + ) = (25) λ n π n Proof As in the proof of Karlin and McGregor, 957, Theorem 5, we consider the Laplace transform D(s) (D ij (s)) of the difference of any two Q-functions satisfying (2) and find that D ij (s) = π j R i ( s)r j ( s)d (s), s. Since the zeros of Rn (x), and hence R n (x), are all positive (see Karlin and McGregor (957)), we have R n ( s) R n (), s, for all n. Hence, if (24) is satisfied then R i ( s) as i. But since D ij (s) is bounded we must have D (s) =, and hence D ij (s) = for all i, j. On the other hand, if (24) is not satisfied then j (λ jπ j ) < as a consequence of (2). So if, at the same time, (25) is satisfied, we must have j π j =, and hence j π jr j ( s) = for s. But, since j D ij(s) must be bounded, it follows again that D (s) =, and hence D ij (s) = for all i, j. Now suppose neither (24) nor (25) are satisfied. Then, in view of (8) and (9), also (3) is not satisfied. As a result there is a number ξ > and a one-parameter family {ψ ξ, ξ ξ } of so-called extremal solutions to the Stieltjes moment problem associated with { R n }. The index ξ in ψ ξ denotes the smallest point in the support of the measure, and the minimal measure ψ = ψ should be identified with ψ ξ (987) for more details). Next letting φ(ξ) (see Karlin and McGregor (957) or van Doorn γ j π j x R j (x)ψ ξ (dx), ξ ξ, j= we note that, by Lemma 3, the boundedness of R n () amounts to φ(ξ ) <. The proof of Karlin and McGregor, 957, Theorem 5 can now be copied to show that φ(ξ) is continuous, so that there must be a number ξ such that 8
11 φ(ξ) < for ξ (ξ, ξ ]. As in Karlin and McGregor (957) the extremal measures ψ ξ, ξ < ξ ξ, can subsequently be used to construct infinitely many distinct Q-functions satisfying (2). If only finitely many γ i are positive then, by (2), statement (24) is equivalent to n= (λ nπ n ) =, so that (25) alone is necessary and sufficient for uniqueness. Thus we have regained the necessary and sufficient condition for uniqueness found by Karlin and McGregor, 957, Theorem 5, in the case that killing is possible in state only. In general, the first condition in the above theorem does not imply the second one. Indeed, from (9) we obtain R n+ () = R n () + γ n R n () + n γ k π k R k (). λ n λ n π n Hence, we may choose λ n and µ n such that the series in (25) converges, and subsequently γ n successively so large that R n+ () R n () +. As a result (24) does hold. We also note that if R n () is bounded then (3) and (25) are equivalent, in view of (8) and (9). So, as we have claimed at the beginning of this section, if P (.) is the unique Q-function satisfying (), then P (.) is the unique Q-function k= satisfying (2). (The reverse does not necessarily hold true.) The integral in (22) seems enigmatic, but has in fact a clear probabilistic interpretation. Namely, let T denote the killing time, that is, the (possibly defective) random variable representing the time at which absorption in the absorbing state occurs. Since, by the forward equations, t Pr{T t X() = } = γ j p j (u)du, we have Pr{T < X() = } = j= γ j p j (u)du, j= so that, by substituting the integral representation (4) and interchanging the integrals, we get Pr{T < X() = } = γ j π j x R j (x)ψ(dx). (26) j= 9
12 So Lemma 3 tells us that absorption at from state (and hence from any state) is certain if and only if R n () is unbounded. More information about the killing time T is given in van Doorn and Zeifman (24). 5 Example If the killing rates are constant, γ i = γ, say, then, by conditioning, the transition probabilities p ij (.) of the process with killing can simply be expressed in terms of the transition probabilities p ij (.) of the process with the same birth and death rates, but zero killing rates. Namely, p ij (t) = e γt p ij (t), i, j S, t, so interesting cases arise when the killing rates are state dependent. As an example we will consider the process X with linear birth, death and killing rates, namely, λ i = λi + θ, µ i = iµ, γ i = iγ, i S, with λ, θ, µ, γ >. It is easy to see that (25) is satisfied, so X is uniquely determined by its rates. Karlin and Tavaré (982) have analysed the process by adroitly relating it to an honest birth-death process with known transition probabilities. We shall see that a direct approach based on the integral representation (4) yields the same result. Indeed, the recurrence relation (3) becomes (λn + θ)r n+ (x) = (θ + (λ + µ + γ)n x)r n (x) µnr n (x), n >, θr (x) = θ x, R (x) =. Now writing and β θ λ, ρ (λ + µ + γ) 2 4λµ, κ θ ( ) 2λ n P n (x) (β) n R n (ρx + κ), n, λ + µ + γ ρ ( ) 2µ λ + µ + γ + ρ where (β) n Γ(β + n)/γ(β), we see after a little algebra that {P n } satisfies the recurrence cp n+ (x) = ((c )x + (c + )n + cβ)p n (x) n(n + β )P n (x),
13 where c = λ + µ + γ ρ λ + µ + γ + ρ. The polynomials P n can now be identified (see Chihara, 978, Section VI.3) with the Meixner polynomials of the first kind, which are orthogonal with respect to a discrete measure with masses at,,.... Specifically, ( c) β x= P m (x)p n (x) cx (β) x x! = δ m,n (β) n n! c n, m, n. As a consequence the orthogonality relation for the polynomials {R n } may be given as ( c) β x= R m (ρx + κ)r n (ρx + κ) cx (β) x x! µ n n! = δ m,n, m, n. (β) n λn By elementary substitution of these findings in the integral representation (4) we regain the result obtained earlier in Karlin and Tavaré (982). References Anderson, W.J. (99), Continuous-Time Markov Chains (Springer, New York). Chao, X. and Zheng, Y. (23), Transient analysis of immigration birth-death processes with total catastrophes, Probab. Engrg. Inform. Sci. 7, Chihara, T.S. (978), An Introduction to Orthogonal Polynomials (Gordon and Breach, New York). van Doorn, E.A. (987), The indeterminate rate problem for birth-death processes, Pacific J. Math. 3, van Doorn, E.A. and Zeifman, A.I. (24), Extinction probability in a birthdeath process with killing, in preparation. Karlin, S. and McGregor, J.L. (957), The differential equations of birth-anddeath processes, and the Stieltjes moment problem, Trans. Amer. Math. Soc. 85, Karlin, S. and McGregor, J.L. (959), A characterization of birth and death processes, Proc. Natl. Acad. Sci. USA 45,
14 Karlin, S. and Tavaré, S. (982), Linear birth and death processes with killing, J. Appl. Probab. 9, Lenin, R.B., Parthasarathy, P.R., Scheinhardt, W.R.W. and van Doorn, E.A. (2), Families of birth-death processes with similar time-dependent behaviour. J. Appl. Probab. 37,
Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 2, pp. 401 412 (2013) http://campus.mst.edu/adsa Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes
More informationEXTINCTION PROBABILITY IN A BIRTH-DEATH PROCESS WITH KILLING
EXTINCTION PROBABILITY IN A BIRTH-DEATH PROCESS WITH KILLING Erik A. van Doorn and Alexander I. Zeifman Department of Applied Mathematics University of Twente P.O. Box 217, 7500 AE Enschede, The Netherlands
More informationErik A. van Doorn Department of Applied Mathematics University of Twente Enschede, The Netherlands
REPRESENTATIONS FOR THE DECAY PARAMETER OF A BIRTH-DEATH PROCESS Erik A. van Doorn Department of Applied Mathematics University of Twente Enschede, The Netherlands University of Leeds 30 October 2014 Acknowledgment:
More informationQUASI-STATIONARY DISTRIBUTIONS FOR BIRTH-DEATH PROCESSES WITH KILLING
QUASI-STATIONARY DISTRIBUTIONS FOR BIRTH-DEATH PROCESSES WITH KILLING PAULINE COOLEN-SCHRIJNER AND ERIK A. VAN DOORN Received 9 January 2006; Revised 9 June 2006; Accepted 28 July 2006 The Karlin-McGregor
More informationQuasi-Stationary Distributions in Linear Birth and Death Processes
Quasi-Stationary Distributions in Linear Birth and Death Processes Wenbo Liu & Hanjun Zhang School of Mathematics and Computational Science, Xiangtan University Hunan 411105, China E-mail: liuwenbo10111011@yahoo.com.cn
More informationAnalysis of birth-death fluid queues
Analysis of birth-death fluid queues E.A. van Doorn and W.R.W. Scheinhardt Faculty of Applied Mathematics University of Twente P.O. Box 217, 7500 AE Enschede, The Netherlands Abstract We present a survey
More informationc 2010 Society for Industrial and Applied Mathematics
THEORY PROBAB APPL Vol 54, No 1, pp 97 113 c 2010 Society for Industrial Applied Mathematics Translated from Russian Journal BOUNDS AND ASYMPTOTICS FOR THE RATE OF CONVERGENCE OF BIRTH-DEATH PROCESSES
More informationUniversity of Twente. Faculty of Mathematical Sciences. The deviation matrix of a continuous-time Markov chain
Faculty of Mathematical Sciences University of Twente University for Technical and Social Sciences P.O. Box 217 75 AE Enschede The Netherlands Phone: +31-53-48934 Fax: +31-53-4893114 Email: memo@math.utwente.nl
More informationIrregular Birth-Death process: stationarity and quasi-stationarity
Irregular Birth-Death process: stationarity and quasi-stationarity MAO Yong-Hua May 8-12, 2017 @ BNU orks with W-J Gao and C Zhang) CONTENTS 1 Stationarity and quasi-stationarity 2 birth-death process
More informationA Birth and Death Process Related to the Rogers-Ramanujan Continued Fraction
A Birth and Death Process Related to the Rogers-Ramanujan Continued Fraction P. R. Parthasarathy, R. B. Lenin Department of Mathematics Indian Institute of Technology, Madras Chennai - 600 036, INDIA W.
More informationStatistics 150: Spring 2007
Statistics 150: Spring 2007 April 23, 2008 0-1 1 Limiting Probabilities If the discrete-time Markov chain with transition probabilities p ij is irreducible and positive recurrent; then the limiting probabilities
More informationBirth and Death Processes. Birth and Death Processes. Linear Growth with Immigration. Limiting Behaviour for Birth and Death Processes
DTU Informatics 247 Stochastic Processes 6, October 27 Today: Limiting behaviour of birth and death processes Birth and death processes with absorbing states Finite state continuous time Markov chains
More informationIntegrals for Continuous-time Markov chains
Integrals for Continuous-time Markov chains P.K. Pollett Abstract This paper presents a method of evaluating the expected value of a path integral for a general Markov chain on a countable state space.
More informationEXTINCTION TIMES FOR A GENERAL BIRTH, DEATH AND CATASTROPHE PROCESS
(February 25, 2004) EXTINCTION TIMES FOR A GENERAL BIRTH, DEATH AND CATASTROPHE PROCESS BEN CAIRNS, University of Queensland PHIL POLLETT, University of Queensland Abstract The birth, death and catastrophe
More informationarxiv: v1 [math.pr] 4 Nov 2016
Perturbations of continuous-time Markov chains arxiv:1611.1254v1 [math.pr 4 Nov 216 Pei-Sen Li School of Mathematical Sciences, Beijing Normal University, Beijing 1875, China E-ma: peisenli@ma.bnu.edu.cn
More informationSIMILAR MARKOV CHAINS
SIMILAR MARKOV CHAINS by Phil Pollett The University of Queensland MAIN REFERENCES Convergence of Markov transition probabilities and their spectral properties 1. Vere-Jones, D. Geometric ergodicity in
More informationTMA4265 Stochastic processes ST2101 Stochastic simulation and modelling
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 7 English Contact during examination: Øyvind Bakke Telephone: 73 9 8 26, 99 4 673 TMA426 Stochastic processes
More informationEXTINCTION TIMES FOR A GENERAL BIRTH, DEATH AND CATASTROPHE PROCESS
J. Appl. Prob. 41, 1211 1218 (2004) Printed in Israel Applied Probability Trust 2004 EXTINCTION TIMES FOR A GENERAL BIRTH, DEATH AND CATASTROPHE PROCESS BEN CAIRNS and P. K. POLLETT, University of Queensland
More informationSMSTC (2007/08) Probability.
SMSTC (27/8) Probability www.smstc.ac.uk Contents 12 Markov chains in continuous time 12 1 12.1 Markov property and the Kolmogorov equations.................... 12 2 12.1.1 Finite state space.................................
More informationMATH 56A SPRING 2008 STOCHASTIC PROCESSES 65
MATH 56A SPRING 2008 STOCHASTIC PROCESSES 65 2.2.5. proof of extinction lemma. The proof of Lemma 2.3 is just like the proof of the lemma I did on Wednesday. It goes like this. Suppose that â is the smallest
More informationA note on the GI/GI/ system with identical service and interarrival-time distributions
A note on the GI/GI/ system with identical service and interarrival-time distributions E.A. van Doorn and A.A. Jagers Department of Applied Mathematics University of Twente P.O. Box 7, 75 AE Enschede,
More informationSTAT STOCHASTIC PROCESSES. Contents
STAT 3911 - STOCHASTIC PROCESSES ANDREW TULLOCH Contents 1. Stochastic Processes 2 2. Classification of states 2 3. Limit theorems for Markov chains 4 4. First step analysis 5 5. Branching processes 5
More informationDifference Equations for Multiple Charlier and Meixner Polynomials 1
Difference Equations for Multiple Charlier and Meixner Polynomials 1 WALTER VAN ASSCHE Department of Mathematics Katholieke Universiteit Leuven B-3001 Leuven, Belgium E-mail: walter@wis.kuleuven.ac.be
More informationBIRTH-DEATH PROCESSES AND q-continued FRACTIONS
BIRTH-DEATH PROCESSES AND q-continued FRACTIONS TONY FENG, RACHEL KIRSCH, ELISE MCCALL, AND MATT WAGE Abstract. In the 997 paper of Parthasarathy, Lenin, Schoutens, and Van Assche [6], the authors study
More informationON COMPOUND POISSON POPULATION MODELS
ON COMPOUND POISSON POPULATION MODELS Martin Möhle, University of Tübingen (joint work with Thierry Huillet, Université de Cergy-Pontoise) Workshop on Probability, Population Genetics and Evolution Centre
More informationDepartment of Applied Mathematics Faculty of EEMCS. University of Twente. Memorandum No. 1796
Department of Applied Mathematics Faculty of EEMCS t University of Twente The Netherlands P.O. Box 217 7500 AE Enschede The Netherlands Phone: +31-53-4893400 Fax: +31-53-4893114 Email: memo@math.utwente.nl
More informationON LINEAR COMBINATIONS OF
Física Teórica, Julio Abad, 1 7 (2008) ON LINEAR COMBINATIONS OF ORTHOGONAL POLYNOMIALS Manuel Alfaro, Ana Peña and María Luisa Rezola Departamento de Matemáticas and IUMA, Facultad de Ciencias, Universidad
More informationThe SIS and SIR stochastic epidemic models revisited
The SIS and SIR stochastic epidemic models revisited Jesús Artalejo Faculty of Mathematics, University Complutense of Madrid Madrid, Spain jesus_artalejomat.ucm.es BCAM Talk, June 16, 2011 Talk Schedule
More informationClassification of Countable State Markov Chains
Classification of Countable State Markov Chains Friday, March 21, 2014 2:01 PM How can we determine whether a communication class in a countable state Markov chain is: transient null recurrent positive
More informationOn hitting times and fastest strong stationary times for skip-free and more general chains
On hitting times and fastest strong stationary times for skip-free and more general chains James Allen Fill 1 Department of Applied Mathematics and Statistics The Johns Hopkins University jimfill@jhu.edu
More informationDepartment of Applied Mathematics Faculty of EEMCS. University of Twente. Memorandum No Owen coalitional value without additivity axiom
Department of Applied Mathematics Faculty of EEMCS t University of Twente The Netherlands P.O. Box 217 7500 AE Enschede The Netherlands Phone: +31-53-4893400 Fax: +31-53-4893114 Email: memo@math.utwente.nl
More information6 Continuous-Time Birth and Death Chains
6 Continuous-Time Birth and Death Chains Angela Peace Biomathematics II MATH 5355 Spring 2017 Lecture notes follow: Allen, Linda JS. An introduction to stochastic processes with applications to biology.
More informationUniversity of Twente. Faculty of Mathematical Sciences. Scheduling split-jobs on parallel machines. University for Technical and Social Sciences
Faculty of Mathematical Sciences University of Twente University for Technical and Social Sciences P.O. Box 217 7500 AE Enschede The Netherlands Phone: +31-53-4893400 Fax: +31-53-4893114 Email: memo@math.utwente.nl
More informationDeterminants Chapter 3 of Lay
Determinants Chapter of Lay Dr. Doreen De Leon Math 152, Fall 201 1 Introduction to Determinants Section.1 of Lay Given a square matrix A = [a ij, the determinant of A is denoted by det A or a 11 a 1j
More informationPrime numbers and Gaussian random walks
Prime numbers and Gaussian random walks K. Bruce Erickson Department of Mathematics University of Washington Seattle, WA 9895-4350 March 24, 205 Introduction Consider a symmetric aperiodic random walk
More informationMarkov Chains, Stochastic Processes, and Matrix Decompositions
Markov Chains, Stochastic Processes, and Matrix Decompositions 5 May 2014 Outline 1 Markov Chains Outline 1 Markov Chains 2 Introduction Perron-Frobenius Matrix Decompositions and Markov Chains Spectral
More informationPOLYNOMIAL SOLUTIONS OF PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS
J. London Math. Soc. 67 (2003) 16 28 C 2003 London Mathematical Society DOI: 10.1112/S002461070200371X POLYNOMIAL SOLUTIONS OF PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS J. MCLAUGHLIN
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Centre for Research on Inner City Health St Michael s Hospital Toronto On leave from Department of Mathematics University of Manitoba Julien Arino@umanitoba.ca
More informationA trigonometric orthogonality with respect to a nonnegative Borel measure
Filomat 6:4 01), 689 696 DOI 10.98/FIL104689M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A trigonometric orthogonality with
More informationMath Linear Algebra Final Exam Review Sheet
Math 15-1 Linear Algebra Final Exam Review Sheet Vector Operations Vector addition is a component-wise operation. Two vectors v and w may be added together as long as they contain the same number n of
More informationOperators with Compatible Ranges
Filomat : (7), 579 585 https://doiorg/98/fil7579d Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://wwwpmfniacrs/filomat Operators with Compatible Ranges
More informationCOMPLETELY INVARIANT JULIA SETS OF POLYNOMIAL SEMIGROUPS
Series Logo Volume 00, Number 00, Xxxx 19xx COMPLETELY INVARIANT JULIA SETS OF POLYNOMIAL SEMIGROUPS RICH STANKEWITZ Abstract. Let G be a semigroup of rational functions of degree at least two, under composition
More informationA NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS. Masaru Nagisa. Received May 19, 2014 ; revised April 10, (Ax, x) 0 for all x C n.
Scientiae Mathematicae Japonicae Online, e-014, 145 15 145 A NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS Masaru Nagisa Received May 19, 014 ; revised April 10, 014 Abstract. Let f be oeprator monotone
More informationarxiv: v2 [math.pr] 4 Sep 2017
arxiv:1708.08576v2 [math.pr] 4 Sep 2017 On the Speed of an Excited Asymmetric Random Walk Mike Cinkoske, Joe Jackson, Claire Plunkett September 5, 2017 Abstract An excited random walk is a non-markovian
More informationSOLUTIONS OF SEMILINEAR WAVE EQUATION VIA STOCHASTIC CASCADES
Communications on Stochastic Analysis Vol. 4, No. 3 010) 45-431 Serials Publications www.serialspublications.com SOLUTIONS OF SEMILINEAR WAVE EQUATION VIA STOCHASTIC CASCADES YURI BAKHTIN* AND CARL MUELLER
More informationRandom Bernstein-Markov factors
Random Bernstein-Markov factors Igor Pritsker and Koushik Ramachandran October 20, 208 Abstract For a polynomial P n of degree n, Bernstein s inequality states that P n n P n for all L p norms on the unit
More informationDiscrete Orthogonal Polynomials on Equidistant Nodes
International Mathematical Forum, 2, 2007, no. 21, 1007-1020 Discrete Orthogonal Polynomials on Equidistant Nodes Alfredo Eisinberg and Giuseppe Fedele Dip. di Elettronica, Informatica e Sistemistica Università
More informationMulti Stage Queuing Model in Level Dependent Quasi Birth Death Process
International Journal of Statistics and Systems ISSN 973-2675 Volume 12, Number 2 (217, pp. 293-31 Research India Publications http://www.ripublication.com Multi Stage Queuing Model in Level Dependent
More informationIEOR 6711: Stochastic Models I Professor Whitt, Thursday, November 29, Weirdness in CTMC s
IEOR 6711: Stochastic Models I Professor Whitt, Thursday, November 29, 2012 Weirdness in CTMC s Where s your will to be weird? Jim Morrison, The Doors We are all a little weird. And life is a little weird.
More informationQuasi-stationary distributions
Quasi-stationary distributions Erik A. van Doorn a and Philip K. Pollett b a Department of Applied Mathematics, University of Twente P.O. Box 217, 7500 AE Enschede, The Netherlands E-mail: e.a.vandoorn@utwente.nl
More informationMoore Penrose inverses and commuting elements of C -algebras
Moore Penrose inverses and commuting elements of C -algebras Julio Benítez Abstract Let a be an element of a C -algebra A satisfying aa = a a, where a is the Moore Penrose inverse of a and let b A. We
More informationMarkov chains. 1 Discrete time Markov chains. c A. J. Ganesh, University of Bristol, 2015
Markov chains c A. J. Ganesh, University of Bristol, 2015 1 Discrete time Markov chains Example: A drunkard is walking home from the pub. There are n lampposts between the pub and his home, at each of
More informationA Method for Evaluating the Distribution of the Total Cost of a Random Process over its Lifetime
A Method for Evaluating the Distribution of the Total Cost of a Random Process over its Lifetime P.K. Pollett a and V.T. Stefanov b a Department of Mathematics, The University of Queensland Queensland
More informationGAUSSIAN 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 informationConstrained Leja points and the numerical solution of the constrained energy problem
Journal of Computational and Applied Mathematics 131 (2001) 427 444 www.elsevier.nl/locate/cam Constrained Leja points and the numerical solution of the constrained energy problem Dan I. Coroian, Peter
More informationMulti-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
International Journal of Difference Equations ISSN 0973-6069, Volume 0, Number, pp. 9 06 205 http://campus.mst.edu/ijde Multi-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 2
MATH 56A: STOCHASTIC PROCESSES CHAPTER 2 2. Countable Markov Chains I started Chapter 2 which talks about Markov chains with a countably infinite number of states. I did my favorite example which is on
More information88 CONTINUOUS MARKOV CHAINS
88 CONTINUOUS MARKOV CHAINS 3.4. birth-death. Continuous birth-death Markov chains are very similar to countable Markov chains. One new concept is explosion which means that an infinite number of state
More informationON THE SIMILARITY OF CENTERED OPERATORS TO CONTRACTIONS. Srdjan Petrović
ON THE SIMILARITY OF CENTERED OPERATORS TO CONTRACTIONS Srdjan Petrović Abstract. In this paper we show that every power bounded operator weighted shift with commuting normal weights is similar to a contraction.
More informationarxiv: v1 [math.ac] 19 Sep 2014
EXISTENCE OF SPECIAL PRIMARY DECOMPOSITIONS IN MULTIGRADED MODULES arxiv:149.5518v1 [math.ac] 19 Sep 214 DIPANKAR GHOSH Abstract. Let R = n N t Rn be a commutative Noetherian Nt -graded ring, and L = n
More informationRELATION BETWEEN SMALL FUNCTIONS WITH DIFFERENTIAL POLYNOMIALS GENERATED BY MEROMORPHIC SOLUTIONS OF HIGHER ORDER LINEAR DIFFERENTIAL EQUATIONS
Kragujevac Journal of Mathematics Volume 38(1) (2014), Pages 147 161. RELATION BETWEEN SMALL FUNCTIONS WITH DIFFERENTIAL POLYNOMIALS GENERATED BY MEROMORPHIC SOLUTIONS OF HIGHER ORDER LINEAR DIFFERENTIAL
More informationConvergence of generalized entropy minimizers in sequences of convex problems
Proceedings IEEE ISIT 206, Barcelona, Spain, 2609 263 Convergence of generalized entropy minimizers in sequences of convex problems Imre Csiszár A Rényi Institute of Mathematics Hungarian Academy of Sciences
More information2. Transience and Recurrence
Virtual Laboratories > 15. Markov Chains > 1 2 3 4 5 6 7 8 9 10 11 12 2. Transience and Recurrence The study of Markov chains, particularly the limiting behavior, depends critically on the random times
More informationCommuting nilpotent matrices and pairs of partitions
Commuting nilpotent matrices and pairs of partitions Roberta Basili Algebraic Combinatorics Meets Inverse Systems Montréal, January 19-21, 2007 We will explain some results on commuting n n matrices and
More informationThe passage time distribution for a birth-and-death chain: Strong stationary duality gives a first stochastic proof
The passage time distribution for a birth-and-death chain: Strong stationary duality gives a first stochastic proof James Allen Fill 1 Department of Applied Mathematics and Statistics The Johns Hopkins
More informationBirth-death chain models (countable state)
Countable State Birth-Death Chains and Branching Processes Tuesday, March 25, 2014 1:59 PM Homework 3 posted, due Friday, April 18. Birth-death chain models (countable state) S = We'll characterize the
More informationOPSF, Random Matrices and Riemann-Hilbert problems
OPSF, Random Matrices and Riemann-Hilbert problems School on Orthogonal Polynomials in Approximation Theory and Mathematical Physics, ICMAT 23 27 October, 2017 Plan of the course lecture 1: Orthogonal
More informationPositivity of Turán determinants for orthogonal polynomials
Positivity of Turán determinants for orthogonal polynomials Ryszard Szwarc Abstract The orthogonal polynomials p n satisfy Turán s inequality if p 2 n (x) p n 1 (x)p n+1 (x) 0 for n 1 and for all x in
More informationProbability Distributions
Lecture 1: Background in Probability Theory Probability Distributions The probability mass function (pmf) or probability density functions (pdf), mean, µ, variance, σ 2, and moment generating function
More informationSTAT 380 Continuous Time Markov Chains
STAT 380 Continuous Time Markov Chains Richard Lockhart Simon Fraser University Spring 2018 Richard Lockhart (Simon Fraser University)STAT 380 Continuous Time Markov Chains Spring 2018 1 / 35 Continuous
More informationCHAPTER 3 Further properties of splines and B-splines
CHAPTER 3 Further properties of splines and B-splines In Chapter 2 we established some of the most elementary properties of B-splines. In this chapter our focus is on the question What kind of functions
More informationLinear-fractional branching processes with countably many types
Branching processes and and their applications Badajoz, April 11-13, 2012 Serik Sagitov Chalmers University and University of Gothenburg Linear-fractional branching processes with countably many types
More informationA SUFFICIENT CONDITION FOR STRICT TOTAL POSITIVITY OF A MATRIX. Thomas Craven and George Csordas
A SUFFICIENT CONDITION FOR STRICT TOTAL POSITIVITY OF A MATRIX Thomas Craven and George Csordas Abstract. We establish a sufficient condition for strict total positivity of a matrix. In particular, we
More informationA NOTE ON THE COMPLETE MOMENT CONVERGENCE FOR ARRAYS OF B-VALUED RANDOM VARIABLES
Bull. Korean Math. Soc. 52 (205), No. 3, pp. 825 836 http://dx.doi.org/0.434/bkms.205.52.3.825 A NOTE ON THE COMPLETE MOMENT CONVERGENCE FOR ARRAYS OF B-VALUED RANDOM VARIABLES Yongfeng Wu and Mingzhu
More informationINTERTWINING OF BIRTH-AND-DEATH PROCESSES
K Y BERNETIKA VOLUM E 47 2011, NUMBER 1, P AGES 1 14 INTERTWINING OF BIRTH-AND-DEATH PROCESSES Jan M. Swart It has been known for a long time that for birth-and-death processes started in zero the first
More informationLimiting behaviour of moving average processes under ρ-mixing assumption
Note di Matematica ISSN 1123-2536, e-issn 1590-0932 Note Mat. 30 (2010) no. 1, 17 23. doi:10.1285/i15900932v30n1p17 Limiting behaviour of moving average processes under ρ-mixing assumption Pingyan Chen
More information16 1 Basic Facts from Functional Analysis and Banach Lattices
16 1 Basic Facts from Functional Analysis and Banach Lattices 1.2.3 Banach Steinhaus Theorem Another fundamental theorem of functional analysis is the Banach Steinhaus theorem, or the Uniform Boundedness
More informationRegularity of the density for the stochastic heat equation
Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department
More informationIEOR 6711, HMWK 5, Professor Sigman
IEOR 6711, HMWK 5, Professor Sigman 1. Semi-Markov processes: Consider an irreducible positive recurrent discrete-time Markov chain {X n } with transition matrix P (P i,j ), i, j S, and finite state space.
More informationDecompositions Of Modules And Matrices
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Publications, Department of Mathematics Mathematics, Department of 1973 Decompositions Of Modules And Matrices Thomas
More informationTRUNCATED TOEPLITZ OPERATORS ON FINITE DIMENSIONAL SPACES
TRUNCATED TOEPLITZ OPERATORS ON FINITE DIMENSIONAL SPACES JOSEPH A. CIMA, WILLIAM T. ROSS, AND WARREN R. WOGEN Abstract. In this paper, we study the matrix representations of compressions of Toeplitz operators
More informationA Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit
A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit David Vener Department of Mathematics, MIT May 5, 3 Introduction In 8.366, we discussed the relationship
More informationUniversity of Twente. Faculty of Mathematical Sciences. A short proof of a conjecture on the T r -choice number of even cycles
Faculty of Mathematical Sciences University of Twente University for Technical and Social Sciences P.O. Box 217 7500 AE Enschede The Netherlands Phone: +31-53-4893400 Fax: +31-53-4893114 Email: memo@math.utwente.nl
More informationGENERALIZED 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 informationMarkov Chains CK eqns Classes Hitting times Rec./trans. Strong Markov Stat. distr. Reversibility * Markov Chains
Markov Chains A random process X is a family {X t : t T } of random variables indexed by some set T. When T = {0, 1, 2,... } one speaks about a discrete-time process, for T = R or T = [0, ) one has a continuous-time
More informationSTATS 3U03. Sang Woo Park. March 29, Textbook: Inroduction to stochastic processes. Requirement: 5 assignments, 2 tests, and 1 final
STATS 3U03 Sang Woo Park March 29, 2017 Course Outline Textbook: Inroduction to stochastic processes Requirement: 5 assignments, 2 tests, and 1 final Test 1: Friday, February 10th Test 2: Friday, March
More informationAARMS Homework Exercises
1 For the gamma distribution, AARMS Homework Exercises (a) Show that the mgf is M(t) = (1 βt) α for t < 1/β (b) Use the mgf to find the mean and variance of the gamma distribution 2 A well-known inequality
More informationLARGE DEVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILED DEPENDENT RANDOM VECTORS*
LARGE EVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILE EPENENT RANOM VECTORS* Adam Jakubowski Alexander V. Nagaev Alexander Zaigraev Nicholas Copernicus University Faculty of Mathematics and Computer Science
More informationGenerating All Circular Shifts by Context-Free Grammars in Chomsky Normal Form
Generating All Circular Shifts by Context-Free Grammars in Chomsky Normal Form Peter R.J. Asveld Department of Computer Science, Twente University of Technology P.O. Box 217, 7500 AE Enschede, the Netherlands
More informationIndependence 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 informationAn Arnoldi Gram-Schmidt process and Hessenberg matrices for Orthonormal Polynomials
[ 1 ] University of Cyprus An Arnoldi Gram-Schmidt process and Hessenberg matrices for Orthonormal Polynomials Nikos Stylianopoulos, University of Cyprus New Perspectives in Univariate and Multivariate
More informationTHE BAILEY TRANSFORM AND FALSE THETA FUNCTIONS
THE BAILEY TRANSFORM AND FALSE THETA FUNCTIONS GEORGE E ANDREWS 1 AND S OLE WARNAAR 2 Abstract An empirical exploration of five of Ramanujan s intriguing false theta function identities leads to unexpected
More informationACI-matrices all of whose completions have the same rank
ACI-matrices all of whose completions have the same rank Zejun Huang, Xingzhi Zhan Department of Mathematics East China Normal University Shanghai 200241, China Abstract We characterize the ACI-matrices
More informationPerov s fixed point theorem for multivalued mappings in generalized Kasahara spaces
Stud. Univ. Babeş-Bolyai Math. 56(2011), No. 3, 19 28 Perov s fixed point theorem for multivalued mappings in generalized Kasahara spaces Alexandru-Darius Filip Abstract. In this paper we give some corresponding
More informationarxiv:math/ v1 [math.pr] 24 Mar 2005
The Annals of Applied Probability 2004, Vol. 14, No. 4, 2057 2089 DOI: 10.1214/105051604000000477 c Institute of Mathematical Statistics, 2004 arxiv:math/0503555v1 [math.pr] 24 Mar 2005 SPECTRAL PROPERTIES
More informationRECENT RESULTS FOR SUPERCRITICAL CONTROLLED BRANCHING PROCESSES WITH CONTROL RANDOM FUNCTIONS
Pliska Stud. Math. Bulgar. 16 (2004), 43-54 STUDIA MATHEMATICA BULGARICA RECENT RESULTS FOR SUPERCRITICAL CONTROLLED BRANCHING PROCESSES WITH CONTROL RANDOM FUNCTIONS Miguel González, Manuel Molina, Inés
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 1
MATH 56A: STOCHASTIC PROCESSES CHAPTER. Finite Markov chains For the sake of completeness of these notes I decided to write a summary of the basic concepts of finite Markov chains. The topics in this chapter
More informationSelf-similar fractals as boundaries of networks
Self-similar fractals as boundaries of networks Erin P. J. Pearse ep@ou.edu University of Oklahoma 4th Cornell Conference on Analysis, Probability, and Mathematical Physics on Fractals AMS Eastern Sectional
More informationSpectral analysis of two doubly infinite Jacobi operators
Spectral analysis of two doubly infinite Jacobi operators František Štampach jointly with Mourad E. H. Ismail Stockholm University Spectral Theory and Applications conference in memory of Boris Pavlov
More informationCONTINUITY OF CP-SEMIGROUPS IN THE POINT-STRONG OPERATOR TOPOLOGY
J. OPERATOR THEORY 64:1(21), 149 154 Copyright by THETA, 21 CONTINUITY OF CP-SEMIGROUPS IN THE POINT-STRONG OPERATOR TOPOLOGY DANIEL MARKIEWICZ and ORR MOSHE SHALIT Communicated by William Arveson ABSTRACT.
More information