2008 Hotelling Lectures
|
|
- Conrad Hancock
- 5 years ago
- Views:
Transcription
1 First Prev Next Go To Go Back Full Screen Close Quit 1 28 Hotelling Lectures 1. Stochastic models for chemical reactions 2. Identifying separated time scales in stochastic models of reaction networks 3. Averaging fast subsystems
2 First Prev Next Go To Go Back Full Screen Close Quit 2 Indeed the use of differential equations supplies the statistician with a powerful tool, replacing the purely empirical fitting of arbitrary curves by a reasonable resultant of general considerations with particular data. But this growing statistical use of differential equations must inevitably face the fact that our a priori knowledge can never supply us with a definite relation between a variable and its rate of change, but only with a correlation. In astronomy, physics and chemistry nearly all correlations obtained in good work are very close either to zero or to perfection. Consequently by discarding the first set and identifying the other with complete causation it has been possible and natural to avoid consideration of correlated but not rigidly connected variables. In this way it has come about that, for those fields in which large random fluctuations and merely correlated variables must be dealt with, the existing theory of differential equations is inadequate and needs to be supplemented by a new theory involving probability. Harold Hotelling JASA 1927 [11]
3 First Prev Next Go To Go Back Full Screen Close Quit 3 Stochastic models for chemical reactions Application of the LLN and CLT to Poisson processes Formulating Markov models Reaction networks Classical scaling and the law of mass action Central limit theorem Diffusion approximations Multiple scales Example: Michaelis-Menten equation Example: Model of a viral infection References Abstract Collaboration with David Anderson, Karen Ball, George Craciun, Hye-Won Kang, Lea Popovic, Greg Rempala
4 First Prev Next Go To Go Back Full Screen Close Quit 4 Poisson processes A Poisson process is a model for a series for random observations occurring in time. x x x x x x x x t Let Y (t) denote the number of observations by time t. In the figure above, Y (t) = 6. Note that for t < s, Y (s) Y (t) is the number of observations in the time interval (t, s]. We make the following assumptions about the model. 1) Observations occur one at a time. 2) Numbers of observations in disjoint time intervals are independent random variables, i.e., if t < t 1 < < t m, then Y (t k ) Y (t k 1 ), k = 1,..., m are independent random variables. 3) The distribution of Y (t + a) Y (t) does not depend on t.
5 First Prev Next Go To Go Back Full Screen Close Quit 5 Characterization of a Poisson process Theorem 1 Under assumptions 1), 2), and 3), there is a constant λ > such that, for t < s, Y (s) Y (t) is Poisson distributed with parameter λ(s t), that is, P {Y (s) Y (t) = k} = (λ(s t))k e λ(s t). k! If λ = 1, then Y is a unit (or rate one) Poisson process. If Y is a unit Poisson process and Y λ (t) Y (λt), then Y λ is a Poisson process with parameter λ.
6 First Prev Next Go To Go Back Full Screen Close Quit 6 Application of the LLN and CLT to Poisson processes Theorem 2 If Y is a unit Poisson process, then for each u >, lim sup Y (Ku) u = a.s. K u u K Proof. For fixed u, by the independent increments assumption, the result is just the ordinary law of large numbers. The uniformity follows by monotonicity. The central limit theorem suggests that for large K Y (Ku) Ku K W (u), where W is standard Brownian motion. Y (Ku) K u + 1 K W (u)
7 First Prev Next Go To Go Back Full Screen Close Quit 7 Formulating Markov models Suppose Y λ (t) = Y (λt) and F t represents the information obtained by observing Y λ (s), for s t. P {Y λ (t+ t) Y λ (t) = 1 F t } = P {Y λ (t+ t) Y λ (t) = 1} = 1 e λ t λ t A continuous time Markov chain X taking values in Z d is specified by giving its transition intensities that determine P {X(t + t) X(t) = l F X t } β l (X(t)) t, l Z d.
8 First Prev Next Go To Go Back Full Screen Close Quit 8 Counting process representation If we write X(t) = X() + ln l (t) l where N l (t) is the number of jumps of l at or before time t, then P {N l (t + t) N l (t) = 1 Ft X } β l (X(t)) t, l Z d. N l is a counting process with intensity (propensity in the chemical literature) β l (X(t)) and we can write N l (t) = Y l ( β l (X(s))ds), where the Y l are independent, unit Poisson processes. Consequently, X(t) = X() + ly l ( l β l (X(s))ds).
9 First Prev Next Go To Go Back Full Screen Close Quit 9 Random jump equation Alternatively, setting β(k) = l β l(k), and N(t) = Y ( X(t) = X() + β(x(s))ds) F (X(s ), ξ N(s ) )dn(s) where Y is a unit Poisson process, {ξ i } are iid uniform [, 1], and P {F (k, ξ) = l} = β l(k) β(k).
10 First Prev Next Go To Go Back Full Screen Close Quit 1 Connections to simulation schemes Simulating the random-jump equation gives Gillespie s [7, 8] direct method (the stochastic simulation algorithm SSA). Simulating the time-change equation gives the next jump) method as defined by Gibson and Bruck [6]. reaction (next For = τ (x) < τ 1 (x) <, satisfying τ k (x) = τ k (x τk ), where { x τk x(s) s < τ (s) = k (x) x(τ k (x) ) s τ k (x) (typically, τ k+1 (x) = τ k (x) + g k+1 (x(τ k ))), simulation of ˆX(t) = X() + ( ) ly l β l ( ˆX(τ k ))(τ k+1 t τ k t) l k gives Gillespie s [9] τ-leap method
11 First Prev Next Go To Go Back Full Screen Close Quit 11 Bilingual dictionary Chemistry propensity master equation nonlinear diffusion approximation Langevin approximation Van Kampen approximation quasi steady state/partial equilibrium Probability intensity forward equation diffusion approximation central limit theorem averaging
12 First Prev Next Go To Go Back Full Screen Close Quit 12 Reaction networks Standard notation for chemical reactions A + B k C is interpreted as a molecule of A combines with a molecule of B to give a molecule of C. A + B C means that the reaction can go in either direction, that is, a molecule of C can dissociate into a molecule of A and a molecule of B. We consider a network of reactions involving m chemical species, A 1,..., A m. m m ν ik A i i=1 i=1 ν ika i where the ν ik and ν ik are nonnegative integers.
13 First Prev Next Go To Go Back Full Screen Close Quit 13 Markov chain models X(t) number of molecules of each species in the system at time t. ν k number of molecules of each chemical species consumed in the kth reaction. ν k number of molecules of each species created by the kth reaction. λ k (x) rate at which the kth reaction occurs. (The propensity/intensity.) If the kth reaction occurs at time t, the new state becomes X(t) = X(t ) + ν k ν k. The number of times that the kth reaction occurs by time t is given by the counting process satisfying R k (t) = Y k ( λ k (X(s))ds), where the Y k are independent unit Poisson processes.
14 First Prev Next Go To Go Back Full Screen Close Quit 14 Equations for the system state The state of the system satisfies X(t) = X() + k = X() + k R k (t)(ν k ν k ) Y k ( λ k (X(s))ds)(ν k ν k ) = (ν ν)r(t) ν is the matrix with columns given by the ν k. ν is the matrix with columns given by the ν k. R(t) is the vector with components R k (t).
15 First Prev Next Go To Go Back Full Screen Close Quit 15 Rates for the law of mass action For a binary reaction A 1 + A 2 A 3 or A 1 + A 2 A 3 + A 4 λ k (x) = κ k x 1 x 2 For A 1 A 2 or A 1 A 2 + A 3, λ k (x) = κ k x 1. For 2A 1 A 2, λ k (x) = κ k x 1 (x 1 1). For a binary reaction A 1 +A 2 A 3, the rate should vary inversely with volume, so it would be better to write λ N k (x) = κ k N 1 x 1 x 2 = Nκ k z 1 z 2, where classically, N is taken to be the volume of the system times Avogadro s number and z i = N 1 x i is the concentration in moles per unit volume. Note that unary reaction rates also satisfy λ k (x) = κ k x i = Nκ k z i.
16 First Prev Next Go To Go Back Full Screen Close Quit 16 General form for classical scaling All the rates naturally satisfy λ N k (x) Nκ k i z ν ik i N λ k (z).
17 First Prev Next Go To Go Back Full Screen Close Quit 17 First scaling limit Setting C N (t) = N 1 X(t) C N (t) = C N () + k C N () + k N 1 Y k ( N 1 Y k (N λ N k (X(s))ds)(ν k ν k ) λ k (C N (s))ds)(ν k ν k ) The law of large numbers for the Poisson process implies N 1 Y (Nu) u, C N (t) C N () + κ k Ci N (s) ν ik (ν k ν k )ds, k i which in the large volume limit gives the classical deterministic law of mass action Ċ(t) = κ k C i (t) ν ik (ν k ν k ) F (C(t)). k i
18 First Prev Next Go To Go Back Full Screen Close Quit 18 Central limit theorem/van Kampen approximation V N (t) N(C N (t) C(t)) V N () + N( 1 t N Y k(n λ k (C N (s))ds)(ν k ν k ) k = V N () + 1 t Ỹ k (N λ k (C N (s))ds)(ν k ν k ) k N + N(F (CN (s)) F (C(s)))ds V N () + W k ( λ k (C(s))ds)(ν k ν k ) k + F (C(s)))V N (s)ds F (C(s))ds)
19 First Prev Next Go To Go Back Full Screen Close Quit 19 Gaussian limit V N converges to the solution of V (t) = V () + k W k ( λ k (C(s))ds)(ν k ν k ) + F (C(s)))V (s)ds C N (t) C(t) + 1 N V (t)
20 First Prev Next Go To Go Back Full Screen Close Quit 2 Diffusion approximation where C N (t) = C N () + k N 1 Y k ( C N () + N 1/2 W k ( k + F (C N (s))ds, F (c) = k λ k (X N (s))ds)(ν k ν k ) λ k (c)(ν k ν k ). λ k (C N (s))ds)(ν k ν k ) The diffusion approximation is given by the equation C N (t) = C N ()+ k N 1/2 W k ( λ k ( C N (s))ds)(ν k ν k )+ F ( C N (s))ds.
21 First Prev Next Go To Go Back Full Screen Close Quit 21 Itô formulation The time-change formulation is equivalent to the Itô equation C N (t) = N 1/2 λ k ( C N (s))d W k (s)(ν k ν k ) C N () + k = C N () + k + F ( C N (s))ds N 1/2 σ( C N (s))d W (s) + where σ(c) is the matrix with columns λ k (c)(ν k ν k). F ( C N (s))ds, See Kurtz [12], Ethier and Kurtz [4], Chapter 1, Gardiner [5], Chapter 7, and Van Kampen [14].
22 First Prev Next Go To Go Back Full Screen Close Quit 22 Multiple scales Take N to be of the order of magnitude of the abundance of the most abundant species in the system. For each species i, define the normalized abundances (or simply, the abundances) by Z i (t) = N α i X i (t), where α i 1 should be selected so that Z i = O(1). Note that the abundance may be the species number (α i = ) or the species concentration or something else. The rate constants may also vary over several orders of magnitude κ k = κ kn β k, so for a binary reaction κ kx i x j = N β k+α i +α j κ k z i z j
23 First Prev Next Go To Go Back Full Screen Close Quit 23 A parameterized family of models Let Z N i (t) = Z i () + k N α i Y k ( N β k+ν k α λ k (Z N (s))ds)(ν ik ν ik ). Then the true model is Z = Z N.
24 First Prev Next Go To Go Back Full Screen Close Quit 24 Example: Michaelis-Menten kinetics Consider the reaction system A + E AE B + E modeled as a continuous time Markov chain satisfying X A (t) = X A () Y 1 ( X E (t) = X E () Y 1 ( +Y 3 ( κ 1X A (s)x E (s)ds) + Y 2 ( κ 1X A (s)x E (s)ds) + Y 2 ( κ 3X AE (s)ds) X B (t) = Y 3 ( κ 3X AE (s)ds) κ 2, κ 3 >> κ 1 κ 2X AE (s)ds) κ 2X AE (s)ds)
25 First Prev Next Go To Go Back Full Screen Close Quit 25 Scaling Note that M = X AE (t) + X E (t) is constant. Let N = O(X A ) >> M. Setting β 2 = β 3 = 1, α A = 1, α E = α AE =, κ 1 = κ 1, κ 2 = κ 2N 1, κ 3 = κ 3N 1 V E (t) = M 1 X E (s)ds, Z A (t) = N 1 X A (t) Z A (t) = Z A () N 1 Y 1 (N X E (t) = X E () Y 1 (N +Y 3 (Nκ 3 M(t V E (t))) X B (t) = Y 3 (Nκ 3 M(t V E (t))) κ 1 MZ A (s)dv E (s)) + N 1 Y 2 (Nκ 2 M(t V E (t))) κ 1 MZ A (s)dv E (s)) + Y 2 (Nκ 2 M(t V E (t)))
26 First Prev Next Go To Go Back Full Screen Close Quit 26 Analysis Divide X E (t) = X E () Y 1 (N by N so that and lim N lim N +Y 3 (Nκ 3 M(t V E (t))) ( (κ 2 + κ 3 )M(t V E (t))) κ 1 MZ A (s)dv E (s)) + Y 2 (Nκ 2 M(t V E (t))) ) κ 1 MZ A (s)dv E (s) = ( ) Z A (t) Z A () + κ 1 MZ A (s)dv E (s) κ 2 M(t V E (t)) =
27 First Prev Next Go To Go Back Full Screen Close Quit 27 Derivation of Michaelis-Menten equation Theorem 3 (Darden [2, 3]) Assume that N and ZA N() = X A()/N x A (). Then (ZA N, V E N) converges to (x A(t), v E (t)) satisfying x A (t) = x A () = and hence v E (s) = κ 1 Mx A (s) v E (s)ds + κ 1 x A (s) v E (s)ds + κ 2 +κ 3 κ 2 +κ 3 +κ 1 x A (s) and ẋ A (t) = Mκ 1κ 3 x A (t) κ 2 + κ 3 + κ 1 x A (s). κ 2 M(1 v E (s))ds(1) (κ 2 + κ 3 )(1 v E (s))ds,
28 First Prev Next Go To Go Back Full Screen Close Quit 28 Example: Model of a viral infection Srivastava, You, Summers, and Yin [13], Haseltine and Rawlings [1], Ball, Kurtz, Popovic, and Rampala [1] Three time-varying species, the viral template, the viral genome, and the viral structural protein (indexed, 1, 2, 3 respectively). The model involves six reactions, T + stuff κ 1 T + G G κ 2 T T + stuff κ 3 T + S T κ 4 S κ 5 G + S κ 6 V
29 First Prev Next Go To Go Back Full Screen Close Quit 29 Stochastic system X 1 (t) = X 1 () + Y b ( X 2 (t) = X 2 () + Y a ( X 3 (t) = X 3 () + Y c ( κ 2X 2 (s)ds) Y d ( κ 1X 1 (s)ds) Y b ( Y f ( κ 6X 2 (s)x 3 (s)ds) κ 3X 1 (s)ds) Y e ( Y f ( κ 6X 2 (s)x 3 (s)ds) κ 4X 1 (s)ds) κ 2X 2 (s)ds) κ 5X 3 (s)ds)
30 Figure 1: Simulation (Haseltine and Rawlings 22) First Prev Next Go To Go Back Full Screen Close Quit 3
31 First Prev Next Go To Go Back Full Screen Close Quit 31 Scaling parameters Each X i is scaled according to its abundance in the system. For N = 1, X 1 = O(N ), X 2 = O(N 2/3 ), and X 3 = O(N ) and we take Z 1 = X 1, Z 2 = X 2 N 2/3, and Z 3 = X 3 N 1. Expressing the rate constants in terms of N = 1 κ κ N 2/3 κ 3 1 N κ κ κ N 5/3
32 First Prev Next Go To Go Back Full Screen Close Quit 32 Normalized system With the scaled rate constants, we have Z N 1 (t) = Z N 1 () + Y b ( Z N 2 (t) = Z N 2 () + N 2/3 Y a ( 2.5Z N 2 (s)ds) Y d ( N 2/3 Y f ( Z3 N (t) = Z3 N () + N 1 Y c ( N 1 Y f ( Z N 1 (s)ds) N 2/3 Y b (.25Z N 1 (s)ds).75z N 2 (s)z N 3 (s)ds) NZ N 1 (s)ds) N 1 Y e (.75Z N 2 (s)z N 3 (s)ds), 2.5Z N 2 (s)ds) 2NZ N 3 (s)ds)
33 First Prev Next Go To Go Back Full Screen Close Quit 33 Limiting system With the scaled rate constants, we have Z 1 (t) = Z 1 () + Y b ( Z 2 (t) = Z 2 () Z 3 (t) = Z 3 () + 2.5Z 2 (s)ds) Y d ( Z 1 (s)ds 2Z 3 (s)ds.25z 1 (s)ds)
34 First Prev Next Go To Go Back Full Screen Close Quit 34 Fast time scale Define V N i (t) = Z i (N 2/3 t). V N 1 (t) = V N 1 () + Y b ( V N 2 (t) = V N 2 () + N 2/3 Y a ( V3 N (t) = V3 N () + N 1 Y c ( 2.5N 2/3 V2 N (s)ds) Y d ( N 2/3 V1 N (s)ds) N 2/3 Y b ( N 2/3 Y f (N 2/3 2.5N 2/3 V N 2 (s)ds) N 5/3 V1 N (s)ds) N 1 Y e ( N 1 Y f (.25N 2/3 V N 1 (s)ds).75v2 N (s)v3 N (s)ds).75n 2/3 V N 2 (s)v N 3 (s)ds) 2N 5/3 V N 3 (s)ds)
35 First Prev Next Go To Go Back Full Screen Close Quit 35 Averaging As N, dividing the equations for V1 N and V3 N by N 2/3 shows that V N 1 (s)ds 1 V N 3 (s)ds 5 V N 2 (s)ds V N 2 (s)ds. The assertion for V3 N and the fact that V2 N is asymptotically regular imply V2 N (s)v3 N (s)ds 5 V2 N (s) 2 ds. It follows that V2 N converges to the solution of (2).
36 First Prev Next Go To Go Back Full Screen Close Quit 36 Law of large numbers Theorem 4 For each δ > and t >, where V 2 is the solution of lim P { sup V2 N (s) V 2 (s) δ} =, N s t V 2 (t) = V 2 () + 7.5V 2 (s)ds) 3.75V 2 (s) 2 ds. (2)
37 First Prev Next Go To Go Back Full Screen Close Quit The Whole System The Reduced System when γ= 3 X1:Viral Template X1:Viral Template 2.5 The Number of X The Number of X Time Time The Whole System The Reduced System when γ= 12 X3:Viral Structural Protein 12 X3:Viral Structural Protein The Number of X The Number of X Time Time
38 First Prev Next Go To Go Back Full Screen Close Quit The Whole System The Reduced System when γ= 3 X2:Viral Genome X2:Viral Genome 25 The Number of X The Number of X Time Time
39 First Prev Next Go To Go Back Full Screen Close Quit The Whole System The Reduced System when γ=2/3 25 X2:Viral Genome X2:Viral Genome The Number of X The Number of X Time Time
40 First Prev Next Go To Go Back Full Screen Close Quit The Whole System The Reduced System when γ=2/3 35 X1:Viral Template X1:Viral Template 3 The Number of X The Number of X Time Time 16 The Whole System X3:Viral Structural Protein 16 The Reduced System when γ=2/3 X3:Viral Structural Protein The Number of X The Number of X Time Time
41 First Prev Next Go To Go Back Full Screen Close Quit 41 References [1] Karen Ball, Thomas G. Kurtz, Lea Popovic, and Greg Rempala. Asymptotic analysis of multiscale approximations to reaction networks. Ann. Appl. Probab., 16(4): , 26. [2] Thomas Darden. A pseudo-steady state approximation for stochastic chemical kinetics. Rocky Mountain J. Math., 9(1):51 71, Conference on Deterministic Differential Equations and Stochastic Processes Models for Biological Systems (San Cristobal, N.M., 1977). [3] Thomas A. Darden. Enzyme kinetics: stochastic vs. deterministic models. In Instabilities, bifurcations, and fluctuations in chemical systems (Austin, Tex., 198), pages Univ. Texas Press, Austin, TX, [4] Stewart N. Ethier and Thomas G. Kurtz. Markov processes. Wiley Series in Probability and Mathematical Statistics: Probability and Mathematical Statistics. John Wiley & Sons Inc., New York, Characterization and convergence. [5] C. W. Gardiner. Handbook of stochastic methods for physics, chemistry and the natural sciences, volume 13 of Springer Series in Synergetics. Springer-Verlag, Berlin, third edition, 24. [6] M. A. Gibson and Bruck J. Efficient exact simulation of chemical systems with many species and many channels. J. Phys. Chem. A, 14(9): , 2. [7] Daniel T. Gillespie. A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J. Computational Phys., 22(4):43 434, [8] Daniel T. Gillespie. Exact stochastic simulation of coupled chemical reactions. J. Phys. Chem., 81:234 61, [9] Daniel T. Gillespie. Approximate accelerated stochastic simulation of chemically reacting systems. The Journal of Chemical Physics, 115(4): , 21. [1] Eric L. Haseltine and James B. Rawlings. Approximate simulation of coupled fast and slow reactions for stochastic chemical kinetics. J. Chem. Phys., 117(15): , 22.
42 First Prev Next Go To Go Back Full Screen Close Quit 42 [11] Harold Hotelling. Differential equations subject to error, and population estimates. Journal of the American Statistical Association, 22(159): , [12] Thomas G. Kurtz. Strong approximation theorems for density dependent Markov chains. Stochastic Processes Appl., 6(3):223 24, 1977/78. [13] R. Srivastava, L. You, J. Summers, and J. Yin. Stochastic vs. deterministic modeling of intracellular viral kinetics. J. Theoret. Biol., 218(3):39 321, 22. [14] N. G. van Kampen. Stochastic processes in physics and chemistry. North-Holland Publishing Co., Amsterdam, Lecture Notes in Mathematics, 888.
43 First Prev Next Go To Go Back Full Screen Close Quit 43 Abstract Stochastic models for chemical reactions Attempts to model chemical reactions within biological cells have led to renewed interest in stochastic models for these systems. The classical stochastic models for chemical reaction networks will be reviewed, and multiscale methods for model reduction will be described. The methods will be illustrated with derivation of the Michaelis-Menten model for enzyme reactions and a simple model of viral infection of a cell.
Stochastic models for chemical reactions
First Prev Next Go To Go Back Full Screen Close Quit 1 Stochastic models for chemical reactions Formulating Markov models Two stochastic equations Simulation schemes Reaction Networks Scaling limit Central
More informationStochastic models of biochemical systems
Stochastic models of biochemical systems David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison University of Amsterdam November 14th, 212 Stochastic models
More informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationComputational methods for continuous time Markov chains with applications to biological processes
Computational methods for continuous time Markov chains with applications to biological processes David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Penn.
More informationLongtime behavior of stochastically modeled biochemical reaction networks
Longtime behavior of stochastically modeled biochemical reaction networks David F. Anderson Department of Mathematics University of Wisconsin - Madison ASU Math Biology Seminar February 17th, 217 Overview
More informationStochastic and deterministic models of biochemical reaction networks
Stochastic and deterministic models of biochemical reaction networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Case Western Reserve November 9th, 215 Goals of the talk
More informationEfficient Leaping Methods for Stochastic Chemical Systems
Efficient Leaping Methods for Stochastic Chemical Systems Ioana Cipcigan Muruhan Rathinam November 18, 28 Abstract. Well stirred chemical reaction systems which involve small numbers of molecules for some
More informationExtending the Tools of Chemical Reaction Engineering to the Molecular Scale
Extending the Tools of Chemical Reaction Engineering to the Molecular Scale Multiple-time-scale order reduction for stochastic kinetics James B. Rawlings Department of Chemical and Biological Engineering
More informationfor all f satisfying E[ f(x) ] <.
. Let (Ω, F, P ) be a probability space and D be a sub-σ-algebra of F. An (H, H)-valued random variable X is independent of D if and only if P ({X Γ} D) = P {X Γ}P (D) for all Γ H and D D. Prove that if
More informationChemical reaction network theory for stochastic and deterministic models of biochemical reaction systems
Chemical reaction network theory for stochastic and deterministic models of biochemical reaction systems University of Wisconsin at Madison anderson@math.wisc.edu MBI Workshop for Young Researchers in
More informationStochastic analysis of biochemical reaction networks with absolute concentration robustness
Stochastic analysis of biochemical reaction networks with absolute concentration robustness David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Boston University
More informationPersistence and Stationary Distributions of Biochemical Reaction Networks
Persistence and Stationary Distributions of Biochemical Reaction Networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Discrete Models in Systems Biology SAMSI December
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationFast Probability Generating Function Method for Stochastic Chemical Reaction Networks
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. 71 (2014) 57-69 ISSN 0340-6253 Fast Probability Generating Function Method for Stochastic Chemical Reaction
More informationSTOCHASTIC CHEMICAL KINETICS
STOCHASTIC CHEICAL KINETICS Dan Gillespie GillespieDT@mailaps.org Current Support: Caltech (NIGS) Caltech (NIH) University of California at Santa Barbara (NIH) Past Support: Caltech (DARPA/AFOSR, Beckman/BNC))
More informationarxiv: v2 [q-bio.mn] 31 Aug 2007
A modified Next Reaction Method for simulating chemical systems with time dependent propensities and delays David F. Anderson Department of Mathematics, University of Wisconsin-Madison, Madison, Wi 53706
More informationMean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations
Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations Ram Sharan Adhikari Assistant Professor Of Mathematics Rogers State University Mathematical
More informationLAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC
LAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC R. G. DOLGOARSHINNYKH Abstract. We establish law of large numbers for SIRS stochastic epidemic processes: as the population size increases the paths of SIRS epidemic
More informationSlow Scale Tau-leaping Method
Slow Scale Tau-leaping Method Yang Cao Linda Petzold Abstract For chemical systems involving both fast and slow scales, stiffness presents challenges for efficient stochastic simulation. Two different
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 informationSupporting Information for The Stochastic Quasi-Steady-State Assumption: Reducing the Model but Not the Noise
TWCCC Texas Wisconsin California Control Consortium Technical report number 2011-03 Supporting Information for The Stochastic Quasi-Steady-State Assumption: Reducing the Model but Not the Noise Rishi Srivastava
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 informationThe Adaptive Explicit-Implicit Tau-Leaping Method with Automatic Tau Selection
The Adaptive Explicit-Implicit Tau-Leaping Method with Automatic Tau Selection Yang Cao Department of Computer Science, 660 McBryde Hall, Virginia Tech, Blacksburg, VA 24061 Daniel T. Gillespie Dan T.
More informationStochastic Simulation of Biochemical Reactions
1 / 75 Stochastic Simulation of Biochemical Reactions Jorge Júlvez University of Zaragoza 2 / 75 Outline 1 Biochemical Kinetics 2 Reaction Rate Equation 3 Chemical Master Equation 4 Stochastic Simulation
More informationModelling in Systems Biology
Modelling in Systems Biology Maria Grazia Vigliotti thanks to my students Anton Stefanek, Ahmed Guecioueur Imperial College Formal representation of chemical reactions precise qualitative and quantitative
More informationManual for SOA Exam MLC.
Chapter 10. Poisson processes. Section 10.5. Nonhomogenous Poisson processes Extract from: Arcones Fall 2009 Edition, available at http://www.actexmadriver.com/ 1/14 Nonhomogenous Poisson processes Definition
More informationOn the numerical solution of the chemical master equation with sums of rank one tensors
ANZIAM J. 52 (CTAC21) pp.c628 C643, 211 C628 On the numerical solution of the chemical master equation with sums of rank one tensors Markus Hegland 1 Jochen Garcke 2 (Received 2 January 211; revised 4
More informationNested Stochastic Simulation Algorithm for Chemical Kinetic Systems with Disparate Rates. Abstract
Nested Stochastic Simulation Algorithm for Chemical Kinetic Systems with Disparate Rates Weinan E Department of Mathematics and PACM, Princeton University, Princeton, NJ 08544, USA Di Liu and Eric Vanden-Eijnden
More informationCellular Automata Approaches to Enzymatic Reaction Networks
Cellular Automata Approaches to Enzymatic Reaction Networks Jörg R. Weimar Institute of Scientific Computing, Technical University Braunschweig, D-38092 Braunschweig, Germany J.Weimar@tu-bs.de, http://www.jweimar.de
More informationDiscrepancies between extinction events and boundary equilibria in reaction networks
Discrepancies between extinction events and boundary equilibria in reaction networks David F. nderson Daniele Cappelletti September 2, 208 bstract Reaction networks are mathematical models of interacting
More informationRegular Variation and Extreme Events for Stochastic Processes
1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for
More informationMultivariate Risk Processes with Interacting Intensities
Multivariate Risk Processes with Interacting Intensities Nicole Bäuerle (joint work with Rudolf Grübel) Luminy, April 2010 Outline Multivariate pure birth processes Multivariate Risk Processes Fluid Limits
More informationApplied Mathematics Letters. Stationary distribution, ergodicity and extinction of a stochastic generalized logistic system
Applied Mathematics Letters 5 (1) 198 1985 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Stationary distribution, ergodicity
More informatione - c o m p a n i o n
OPERATIONS RESEARCH http://dx.doi.org/1.1287/opre.111.13ec e - c o m p a n i o n ONLY AVAILABLE IN ELECTRONIC FORM 212 INFORMS Electronic Companion A Diffusion Regime with Nondegenerate Slowdown by Rami
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationDerivations for order reduction of the chemical master equation
2 TWMCC Texas-Wisconsin Modeling and Control Consortium 1 Technical report number 2006-02 Derivations for order reduction of the chemical master equation Ethan A. Mastny, Eric L. Haseltine, and James B.
More informationApproximation of BSDEs using least-squares regression and Malliavin weights
Approximation of BSDEs using least-squares regression and Malliavin weights Plamen Turkedjiev (turkedji@math.hu-berlin.de) 3rd July, 2012 Joint work with Prof. Emmanuel Gobet (E cole Polytechnique) Plamen
More informationSTAT 331. Martingale Central Limit Theorem and Related Results
STAT 331 Martingale Central Limit Theorem and Related Results In this unit we discuss a version of the martingale central limit theorem, which states that under certain conditions, a sum of orthogonal
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationAn Introduction to Stochastic Simulation
Stephen Gilmore Laboratory for Foundations of Computer Science School of Informatics University of Edinburgh PASTA workshop, London, 29th June 2006 Background The modelling of chemical reactions using
More information1. Introduction to Chemical Kinetics
1. Introduction to Chemical Kinetics objectives of chemical kinetics 1) Determine empirical rate laws H 2 + I 2 2HI How does the concentration of H 2, I 2, and HI change with time? 2) Determine the mechanism
More informationSelf-intersection local time for Gaussian processes
Self-intersection local Olga Izyumtseva, olaizyumtseva@yahoo.com (in collaboration with Andrey Dorogovtsev, adoro@imath.kiev.ua) Department of theory of random processes Institute of mathematics Ukrainian
More information1 Basic concepts from probability theory
Basic concepts from probability theory This chapter is devoted to some basic concepts from probability theory.. Random variable Random variables are denoted by capitals, X, Y, etc. The expected value or
More informationPoint Process Control
Point Process Control The following note is based on Chapters I, II and VII in Brémaud s book Point Processes and Queues (1981). 1 Basic Definitions Consider some probability space (Ω, F, P). A real-valued
More informationIntroduction to asymptotic techniques for stochastic systems with multiple time-scales
Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x
More informationBrownian Motion and Conditional Probability
Math 561: Theory of Probability (Spring 2018) Week 10 Brownian Motion and Conditional Probability 10.1 Standard Brownian Motion (SBM) Brownian motion is a stochastic process with both practical and theoretical
More informationMA 777: Topics in Mathematical Biology
MA 777: Topics in Mathematical Biology David Murrugarra Department of Mathematics, University of Kentucky http://www.math.uky.edu/~dmu228/ma777/ Spring 2018 David Murrugarra (University of Kentucky) Lecture
More information6. Brownian Motion. Q(A) = P [ ω : x(, ω) A )
6. Brownian Motion. stochastic process can be thought of in one of many equivalent ways. We can begin with an underlying probability space (Ω, Σ, P) and a real valued stochastic process can be defined
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 informationMathematical Modeling I
Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA February 14, 2017 Introduction When we choose to model
More informationKinetic Monte Carlo. Heiko Rieger. Theoretical Physics Saarland University Saarbrücken, Germany
Kinetic Monte Carlo Heiko Rieger Theoretical Physics Saarland University Saarbrücken, Germany DPG school on Efficient Algorithms in Computational Physics, 10.-14.9.2012, Bad Honnef Intro Kinetic Monte
More informationc 2018 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 78, No. 5, pp. 2692 2713 c 2018 Society for Industrial and Applied Mathematics SOME NETWORK CONDITIONS FOR POSITIVE RECURRENCE OF STOCHASTICALLY MODELED REACTION NETWORKS DAVID
More informationJump-type Levy Processes
Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More informationA Central Limit Theorem for Fleming-Viot Particle Systems Application to the Adaptive Multilevel Splitting Algorithm
A Central Limit Theorem for Fleming-Viot Particle Systems Application to the Algorithm F. Cérou 1,2 B. Delyon 2 A. Guyader 3 M. Rousset 1,2 1 Inria Rennes Bretagne Atlantique 2 IRMAR, Université de Rennes
More informationLet (Ω, F) be a measureable space. A filtration in discrete time is a sequence of. F s F t
2.2 Filtrations Let (Ω, F) be a measureable space. A filtration in discrete time is a sequence of σ algebras {F t } such that F t F and F t F t+1 for all t = 0, 1,.... In continuous time, the second condition
More informationStochastic approaches for modelling in vivo reactions
Stochastic approaches for modelling in vivo reactions T. E. Turner a, S. Schnell b,c, and K. Burrage d. a Oxford Centre for Industrial and Applied Mathematics, Mathematical Institute, 24-29 St. Giles,
More informationGeneralized Chemical Langevin Equation for Chemical Reaction Systems
ATCH Communications in athematical and in Computer Chemistry ATCH Commun. ath. Comput. Chem. 78 (207) 28-305 ISSN 0340-6253 Generalized Chemical Langevin Equation for Chemical Reaction Systems Tao Li a,
More informationMaximum Process Problems in Optimal Control Theory
J. Appl. Math. Stochastic Anal. Vol. 25, No., 25, (77-88) Research Report No. 423, 2, Dept. Theoret. Statist. Aarhus (2 pp) Maximum Process Problems in Optimal Control Theory GORAN PESKIR 3 Given a standard
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. XI Stochastic Stability - H.J. Kushner
STOCHASTIC STABILITY H.J. Kushner Applied Mathematics, Brown University, Providence, RI, USA. Keywords: stability, stochastic stability, random perturbations, Markov systems, robustness, perturbed systems,
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
More informationLIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974
LIMITS FOR QUEUES AS THE WAITING ROOM GROWS by Daniel P. Heyman Ward Whitt Bell Communications Research AT&T Bell Laboratories Red Bank, NJ 07701 Murray Hill, NJ 07974 May 11, 1988 ABSTRACT We study the
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationNested Stochastic Simulation Algorithm for KMC with Multiple Time-Scales
Nested Stochastic Simulation Algorithm for KMC with Multiple Time-Scales Eric Vanden-Eijnden Courant Institute Join work with Weinan E and Di Liu Chemical Systems: Evolution of an isothermal, spatially
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationSome Properties of NSFDEs
Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline
More informationOn the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem
On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family
More informationEULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS
Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show
More informationTyler Hofmeister. University of Calgary Mathematical and Computational Finance Laboratory
JUMP PROCESSES GENERALIZING STOCHASTIC INTEGRALS WITH JUMPS Tyler Hofmeister University of Calgary Mathematical and Computational Finance Laboratory Overview 1. General Method 2. Poisson Processes 3. Diffusion
More informationThe multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications
The multidimensional Ito Integral and the multidimensional Ito Formula Eric Mu ller June 1, 215 Seminar on Stochastic Geometry and its applications page 2 Seminar on Stochastic Geometry and its applications
More information1 Brownian Local Time
1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =
More informationIntroduction to stochastic multiscale modelling in tumour growth
Introduction to stochastic multiscale modelling in tumour growth Tomás Alarcón ICREA & Centre de Recerca Matemàtica T. Alarcón (ICREA & CRM, Barcelona, Spain) Lecture 3 CIMPA, Santiago de Cuba, June 2016
More informationAPPLICATION OF THE KALMAN-BUCY FILTER IN THE STOCHASTIC DIFFERENTIAL EQUATIONS FOR THE MODELING OF RL CIRCUIT
Int. J. Nonlinear Anal. Appl. 2 (211) No.1, 35 41 ISSN: 28-6822 (electronic) http://www.ijnaa.com APPICATION OF THE KAMAN-BUCY FITER IN THE STOCHASTIC DIFFERENTIA EQUATIONS FOR THE MODEING OF R CIRCUIT
More informationInvariant Manifold Reductions for Markovian Ion Channel Dynamics
Invariant Manifold Reductions for Markovian Ion Channel Dynamics James P. Keener Department of Mathematics, University of Utah, Salt Lake City, UT 84112 June 4, 2008 Abstract We show that many Markov models
More informationHandbook of Stochastic Methods
Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei
More informationMath 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA Molecular biology and Biochemical kinetics Molecular biology is one of
More informationStationary distribution and pathwise estimation of n-species mutualism system with stochastic perturbation
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 6), 936 93 Research Article Stationary distribution and pathwise estimation of n-species mutualism system with stochastic perturbation Weiwei
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
More informationA Note on the Central Limit Theorem for a Class of Linear Systems 1
A Note on the Central Limit Theorem for a Class of Linear Systems 1 Contents Yukio Nagahata Department of Mathematics, Graduate School of Engineering Science Osaka University, Toyonaka 560-8531, Japan.
More informationA connection between the stochastic heat equation and fractional Brownian motion, and a simple proof of a result of Talagrand
A connection between the stochastic heat equation and fractional Brownian motion, and a simple proof of a result of Talagrand Carl Mueller 1 and Zhixin Wu Abstract We give a new representation of fractional
More informationBrownian motion and the Central Limit Theorem
Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall
More informationStochastic Processes at Single-molecule and Single-cell levels
Stochastic Processes at Single-molecule and Single-cell levels Hao Ge haoge@pu.edu.cn Beijing International Center for Mathematical Research 2 Biodynamic Optical Imaging Center Peing University, China
More informationLie Symmetry of Ito Stochastic Differential Equation Driven by Poisson Process
American Review of Mathematics Statistics June 016, Vol. 4, No. 1, pp. 17-30 ISSN: 374-348 (Print), 374-356 (Online) Copyright The Author(s).All Rights Reserved. Published by American Research Institute
More informationGENERATION OF COLORED NOISE
International Journal of Modern Physics C, Vol. 12, No. 6 (2001) 851 855 c World Scientific Publishing Company GENERATION OF COLORED NOISE LORENZ BARTOSCH Institut für Theoretische Physik, Johann Wolfgang
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationSTOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY
STOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY Third edition N.G. VAN KAMPEN Institute for Theoretical Physics of the University at Utrecht ELSEVIER Amsterdam Boston Heidelberg London New York Oxford Paris
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 informationMath 735: Stochastic Analysis
First Prev Next Go To Go Back Full Screen Close Quit 1 Math 735: Stochastic Analysis 1. Introduction and review 2. Notions of convergence 3. Continuous time stochastic processes 4. Information and conditional
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationCentral limit theorems for ergodic continuous-time Markov chains with applications to single birth processes
Front. Math. China 215, 1(4): 933 947 DOI 1.17/s11464-15-488-5 Central limit theorems for ergodic continuous-time Markov chains with applications to single birth processes Yuanyuan LIU 1, Yuhui ZHANG 2
More informationAccurate hybrid stochastic simulation of a system of coupled chemical or biochemical reactions
THE JOURNAL OF CHEMICAL PHYSICS 122, 054103 2005 Accurate hybrid schastic simulation of a system of coupled chemical or biochemical reactions Howard Salis and Yiannis Kaznessis a) Department of Chemical
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationOne-Parameter Processes, Usually Functions of Time
Chapter 4 One-Parameter Processes, Usually Functions of Time Section 4.1 defines one-parameter processes, and their variations (discrete or continuous parameter, one- or two- sided parameter), including
More informationTransport of single molecules along the periodic parallel lattices with coupling
THE JOURNAL OF CHEMICAL PHYSICS 124 204901 2006 Transport of single molecules along the periodic parallel lattices with coupling Evgeny B. Stukalin The James Franck Institute The University of Chicago
More informationPositive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network. Haifa Statistics Seminar May 5, 2008
Positive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network Yoni Nazarathy Gideon Weiss Haifa Statistics Seminar May 5, 2008 1 Outline 1 Preview of Results 2 Introduction Queueing
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 informationAvailable online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing
More informationUpper and lower bounds for ruin probability
Upper and lower bounds for ruin probability E. Pancheva,Z.Volkovich and L.Morozensky 3 Institute of Mathematics and Informatics, the Bulgarian Academy of Sciences, 3 Sofia, Bulgaria pancheva@math.bas.bg
More informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationQuasi-Stationary Simulation: the Subcritical Contact Process
Brazilian Journal of Physics, vol. 36, no. 3A, September, 6 685 Quasi-Stationary Simulation: the Subcritical Contact Process Marcelo Martins de Oliveira and Ronald Dickman Departamento de Física, ICEx,
More information