Limit theorems for discrete-time metapopulation models
|
|
- Clara Quinn
- 5 years ago
- Views:
Transcription
1 MASCOS Insiu de Mahémaiques de Toulouse, June Page 1 Limi heorems for discree-ime meapopulaion models Phil Polle Deparmen of Mahemaics The Universiy of Queensland hp:// pkp AUSTRALIAN RESEARCH COUNCIL Cenre of Excellence for Mahemaics and Saisics of Complex Sysems
2 Meapopulaions MASCOS Insiu de Mahémaiques de Toulouse, June Page 2
3 MASCOS Insiu de Mahémaiques de Toulouse, June Page 3 Meapopulaions Colonizaion
4 Meapopulaions MASCOS Insiu de Mahémaiques de Toulouse, June Page 4
5 MASCOS Insiu de Mahémaiques de Toulouse, June Page 5 Meapopulaions Local Exincion
6 Meapopulaions MASCOS Insiu de Mahémaiques de Toulouse, June Page 6
7 Meapopulaions MASCOS Insiu de Mahémaiques de Toulouse, June Page 7
8 Meapopulaions MASCOS Insiu de Mahémaiques de Toulouse, June Page 8
9 MASCOS Insiu de Mahémaiques de Toulouse, June Page 9 Meapopulaions Toal Exincion
10 Meapopulaions MASCOS Insiu de Mahémaiques de Toulouse, June Page 10
11 Mainland-island configuraion MASCOS Insiu de Mahémaiques de Toulouse, June Page 11
12 MASCOS Insiu de Mahémaiques de Toulouse, June Page 12 Mainland-island configuraion Colonizaion from he mainland
13 Meapopulaions MASCOS Insiu de Mahémaiques de Toulouse, June Page 13
14 MASCOS Insiu de Mahémaiques de Toulouse, June Page 14 Pach-occupancy models We record he number n of occupied paches a each ime. A ypical approach is o suppose ha (n, 0) is Markovian.
15 MASCOS Insiu de Mahémaiques de Toulouse, June Page 14 Pach-occupancy models We record he number n of occupied paches a each ime. A ypical approach is o suppose ha (n, 0) is Markovian. Suppose ha here are N paches. Each occupied pach becomes empy a rae e (he local exincion rae), colonizaion of empy paches occurs a rae c/n for each suiable pair (c is he colonizaion rae) and immigraion from he mainland occurs ha rae v (he immigraion rae).
16 MASCOS Insiu de Mahémaiques de Toulouse, June Page 15 A coninuous-ime sochasic model The sae space of he Markov chain (n, 0) is S = {0, 1,...,N} and he ransiions are: n n + 1 a rae n n 1 a rae en (ν + c ) N n (N n)
17 MASCOS Insiu de Mahémaiques de Toulouse, June Page 15 A coninuous-ime sochasic model The sae space of he Markov chain (n, 0) is S = {0, 1,...,N} and he ransiions are: n n + 1 a rae n n 1 a rae en (ν + c ) N n (N n) This an example of Feller s sochasic logisic (SL) model, sudied in deail by J.V. Ross. Ross, J.V. (2006) Sochasic models for mainland-island meapopulaions in saic and dynamic landscapes. Bullein of Mahemaical Biology 68, Feller, W. (1939) Die grundlagen der volerraschen heorie des kampfes ums dasein in wahrscheinlichkeiseoreischer behandlung. Aca Bioheoreica 5,
18 MASCOS Insiu de Mahémaiques de Toulouse, June Page 16 Accouning for life cycle Many species have life cycles (ofen annual) ha consis of disinc phases, and he propensiy for colonizaion and local exincion is differen in each phase.
19 MASCOS Insiu de Mahémaiques de Toulouse, June Page 16 Accouning for life cycle Many species have life cycles (ofen annual) ha consis of disinc phases, and he propensiy for colonizaion and local exincion is differen in each phase. Examples: The Vernal pool fairy shrimp (Branchineca lynchi) and he California linderiella (Linderiella occidenalis), boh lised under he Endangered Species Ac (USA) The Jasper Ridge populaion of Bay checkerspo buerfly (Euphydryas ediha bayensis), now exinc
20 MASCOS Insiu de Mahémaiques de Toulouse, June Page 17 Colonizaion and exincion phases For he buerfly, colonizaion is resriced o he adul phase and here is a greaer propensiy for local exincion in he non-adul phases.
21 MASCOS Insiu de Mahémaiques de Toulouse, June Page 17 Colonizaion and exincion phases For he buerfly, colonizaion is resriced o he adul phase and here is a greaer propensiy for local exincion in he non-adul phases. We will assume ha ha colonizaion (C) and exincion (E) occur in separae disinc phases.
22 MASCOS Insiu de Mahémaiques de Toulouse, June Page 17 Colonizaion and exincion phases For he buerfly, colonizaion is resriced o he adul phase and here is a greaer propensiy for local exincion in he non-adul phases. We will assume ha ha colonizaion (C) and exincion (E) occur in separae disinc phases. There are several ways o model his: A quasi-birh-deah process wih wo phases A non-homogeneous coninuous-ime Markov chain (cycle beween wo ses of ransiion raes) A discree-ime Markov chain
23 MASCOS Insiu de Mahémaiques de Toulouse, June Page 18 Colonizaion and exincion phases For he buerfly, colonizaion is resriced o he adul phase and here is a greaer propensiy for local exincion in he non-adul phases. We will assume ha ha colonizaion (C) and exincion (E) occur in separae disinc phases. There are several ways o model his: A quasi-birh-deah process wih wo phases A non-homogeneous coninuous-ime Markov chain (cycle beween wo ses of ransiion raes) A discree-ime Markov chain
24 MASCOS Insiu de Mahémaiques de Toulouse, June Page 19 A discree-ime Markovian model Recall ha here are N paches and ha n is he number of occupied paches a ime. We suppose ha (n, = 0, 1,...) is a discree-ime Markov chain aking values in S = {0, 1,...,N} wih a 1-sep ransiion marix P = (p ij ) consruced as follows.
25 MASCOS Insiu de Mahémaiques de Toulouse, June Page 19 A discree-ime Markovian model Recall ha here are N paches and ha n is he number of occupied paches a ime. We suppose ha (n, = 0, 1,...) is a discree-ime Markov chain aking values in S = {0, 1,...,N} wih a 1-sep ransiion marix P = (p ij ) consruced as follows. The exincion and colonizaion phases are governed by heir own ransiion marices, E = (e ij ) and C = (c ij ). We le P = EC if he census is aken afer he colonizaion phase or P = CE if he census is aken afer he exincion phase.
26 MASCOS Insiu de Mahémaiques de Toulouse, June Page 20 EC versus CE P = EC { P = CE {
27 MASCOS Insiu de Mahémaiques de Toulouse, June Page 21 Assumpions The number of exincions when here are i paches occupied follows a Bin(i,e) law (0 < e < 1): e i,i k = ( ) i e k (1 e) i k k (k = 0, 1,...,i). (e ij = 0 if j > i.) The number of colonizaions when here are i paches occupied follows a Bin(N i,c i ) law: c i,i+k = (c ij = 0 if j < i.) ( ) N i c k i (1 c i ) N i k (k = 0, 1,...,N i). k
28 MASCOS Insiu de Mahémaiques de Toulouse, June Page 22 Chain-binomial srucure Thus, we have he following chain-binomial srucure: n +1 = ñ + Bin(N ñ,cñ ) ñ = n Bin(n,e) (EC) n +1 = ñ Bin(ñ,e) ñ = n + Bin(N n,c n ). (CE)
29 MASCOS Insiu de Mahémaiques de Toulouse, June Page 22 Chain-binomial srucure Thus, we have he following chain-binomial srucure: n +1 = ñ + Bin(N ñ,cñ ) ñ = n Bin(n,e) (EC) n +1 = ñ Bin(ñ,e) ñ = n + Bin(N n,c n ). (CE) For he CE model (only) i is easy o show ha n +1 has he same disribuion as he sum of wo independen binomial random variables: n +1 D = Bin(n, 1 e) + Bin(N n, (1 e)c n ).
30 MASCOS Insiu de Mahémaiques de Toulouse, June Page 22 Chain-binomial srucure Thus, we have he following chain-binomial srucure: n +1 = ñ + Bin(N ñ,cñ ) ñ = n Bin(n,e) (EC) n +1 = ñ Bin(ñ,e) ñ = n + Bin(N n,c n ). (CE) For he CE model (only) i is easy o show ha n +1 has he same disribuion as he sum of wo independen binomial random variables: n +1 D = Bin(n, 1 e) + Bin(N n, (1 e)c n ). So, (1 e)c i is he effecive colonisaion probabiliy when here are i occupied paches.
31 MASCOS Insiu de Mahémaiques de Toulouse, June Page 23 Examples of c i c i = (i/n)c, where c (0, 1] is he maximum colonizaion poenial. (This enails c 0j = δ 0j, so ha 0 is an absorbing sae and {1,...,N} is a communicaing class.)
32 MASCOS Insiu de Mahémaiques de Toulouse, June Page 23 Examples of c i c i = (i/n)c, where c (0, 1] is he maximum colonizaion poenial. (This enails c 0j = δ 0j, so ha 0 is an absorbing sae and {1,...,N} is a communicaing class.) c i = c, where c (0, 1] is a fixed colonizaion poenial mainland colonizaion dominan. (Now {0, 1,...,N} is irreducible.)
33 MASCOS Insiu de Mahémaiques de Toulouse, June Page 23 Examples of c i c i = (i/n)c, where c (0, 1] is he maximum colonizaion poenial. (This enails c 0j = δ 0j, so ha 0 is an absorbing sae and {1,...,N} is a communicaing class.) c i = c, where c (0, 1] is a fixed colonizaion poenial mainland colonizaion dominan. (Now {0, 1,...,N} is irreducible.) Oher possibiliies include c i = c 0 (1 (1 c 1 /c 0 ) i ), c i = 1 exp( iβ/n) and c i = c 0 + (i/n)c, where c 0 + c (0, 1] (mainland and island colonizaion).
34 MASCOS Insiu de Mahémaiques de Toulouse, June Page 24 The proporion of occupied paches Henceforh we shall be concerned wih X (N) he proporion of occupied paches a ime. = n /N,
35 MASCOS Insiu de Mahémaiques de Toulouse, June Page 25 Simulaion: EC Model wih c i = c 1 Mainland-Island simulaion P = EC (N=100, x 0 =0.05, e =0.01, c =0.05) X (N)
36 MASCOS Insiu de Mahémaiques de Toulouse, June Page 26 The proporion of occupied paches Henceforh we shall be concerned wih X (N) he proporion of occupied paches a ime. = n /N,
37 MASCOS Insiu de Mahémaiques de Toulouse, June Page 26 The proporion of occupied paches Henceforh we shall be concerned wih X (N) he proporion of occupied paches a ime. = n /N, In he mainland-island case c i = c, he disribuion of n can be evaluaed explicily, and we have esablished large-n deerminisic and Gaussian approximaions for (X (N) ). Buckley, F.M. and Polle, P.K. (2009) Analyical mehods for a sochasic mainlandisland meapopulaion model. Ecological Modelling. In press (acceped 24/02/10).
38 MASCOS Insiu de Mahémaiques de Toulouse, June Page 27 Mainland-Island c i = c (Summary) Le p = 1 e(1 c) q = c (EC model) p = 1 e q = (1 e)c. (CE model) and define sequences (p ) and (q ) by q = q (1 a ) and p = q + a ( 0), where a = p q = (1 e)(1 c) (he same for boh EC and CE) and q = q/(1 a).
39 MASCOS Insiu de Mahémaiques de Toulouse, June Page 27 Mainland-Island c i = c (Summary) Le p = 1 e(1 c) q = c (EC model) p = 1 e q = (1 e)c. (CE model) and define sequences (p ) and (q ) by q = q (1 a ) and p = q + a ( 0), where a = p q = (1 e)(1 c) (he same for boh EC and CE) and q = q/(1 a). Then, n D = Bin(n0,p ) + Bin(N n 0,q ) (independen binomial random variables).
40 MASCOS Insiu de Mahémaiques de Toulouse, June Page 28 Mainland-Island c i = c (Summary) Le p = 1 e(1 c) q = c (EC model) p = 1 e q = (1 e)c. (CE model) and define sequences (p ) and (q ) by q = q (1 a ) and p = q + a ( 0), where a = p q = (1 e)(1 c) (he same for boh EC and CE) and q = q/(1 a). Then, n D = Bin(n0,p ) + Bin(N n 0,q ) ( D Bin(N,q ) ) (independen binomial random variables).
41 MASCOS Insiu de Mahémaiques de Toulouse, June Page 29 Mainland-Island c i = c (Summary) Le X (N) = n /N be he proporion occupied a ime. If X (N) 0 P x 0, as N, hen X (N) x = x 0 p + (1 x 0 )q. P x, where
42 MASCOS Insiu de Mahémaiques de Toulouse, June Page 30 Simulaion: EC Model wih c i = c 1 Mainland-Island simulaion P = EC (N=100, x 0 =0.05, e =0.01, c =0.05) X (N)
43 MASCOS Insiu de Mahémaiques de Toulouse, June Page 31 Simulaion: EC Model (Deerminisic pah) 1 Mainland-Island simulaion P = EC (N=100, x 0 =0.05, e =0.01, c =0.05) X (N) Deerminisic pah
44 MASCOS Insiu de Mahémaiques de Toulouse, June Page 32 Mainland-Island c i = c (Summary) Le X (N) = n /N be he proporion occupied a ime. If X (N) 0 P x 0, as N, hen X (N) x = x 0 p + (1 x 0 )q. P x, where
45 MASCOS Insiu de Mahémaiques de Toulouse, June Page 32 Mainland-Island c i = c (Summary) Le X (N) = n /N be he proporion occupied a ime. If X (N) 0 P x 0, as N, hen X (N) x = x 0 p + (1 x 0 )q. P x, where Now pu Z (N) := N(X (N) x ).
46 MASCOS Insiu de Mahémaiques de Toulouse, June Page 32 Mainland-Island c i = c (Summary) Le X (N) = n /N be he proporion occupied a ime. If X (N) 0 P x 0, as N, hen X (N) x = x 0 p + (1 x 0 )q. P x, where Now pu Z (N) Z (N) D := N(X (N) N(a z 0,V ), where x ). Then, if Z (N) 0 D z 0, V = x 0 p (1 p ) + (1 x 0 )q (1 q ).
47 MASCOS Insiu de Mahémaiques de Toulouse, June Page 33 Simulaion: EC Model (Gaussian approx.) 1 Mainland-Island simulaion P = EC (N=100, x 0 =0.05, e =0.01, c =0.05) X (N) Deerminisic pah ± wo sandard deviaions
48 MASCOS Insiu de Mahémaiques de Toulouse, June Page 34 Gaussian approximaions Can we esablish deerminisic and Gaussian approximaions for he basic N-pach models (where he disribuion of n is no known explicily)?
49 MASCOS Insiu de Mahémaiques de Toulouse, June Page 35 Simulaion: EC Model wih c i = (i/n)c 100 Meapopulaion simulaion P = EC (N=100, n 0 =95, e =0.3, c =0.8) n
50 MASCOS Insiu de Mahémaiques de Toulouse, June Page 36 Sim. & qsd: EC Model wih c i = (i/n)c 100 Meapopulaion simulaion P = EC (N=100, n 0 =95, e =0.3, c =0.8) n
51 MASCOS Insiu de Mahémaiques de Toulouse, June Page 37 Gaussian approximaions Can we esablish deerminisic and Gaussian approximaions for he basic N-pach models (where he disribuion of n is no known explicily)?
52 MASCOS Insiu de Mahémaiques de Toulouse, June Page 37 Gaussian approximaions Can we esablish deerminisic and Gaussian approximaions for he basic N-pach models (where he disribuion of n is no known explicily)? Is here a general heory of convergence for discree-ime Markov chains ha share he salien feaures of he pach-occupancy models presened here?
53 MASCOS Insiu de Mahémaiques de Toulouse, June Page 38 General srucure: densiy dependence We have a sequence of Markov chains (n (N) ) indexed by N, ogeher wih funcions (f ) such ha E(n (N) +1 n(n) ) = Nf (n (N) /N).
54 MASCOS Insiu de Mahémaiques de Toulouse, June Page 39 General srucure: densiy dependence We have a sequence of Markov chains (n (N) ) indexed by N, ogeher wih funcions (f ) such ha E(n (N) +1 n(n) ) = Nf (n (N) /N). We hen define (X (N) ) by X (N) = n (N) /N.
55 MASCOS Insiu de Mahémaiques de Toulouse, June Page 40 General srucure: densiy dependence We have a sequence of Markov chains (n (N) ) indexed by N, ogeher wih funcions (f ) such ha E(X (N) +1 X(N) ) = f (X (N) ).
56 MASCOS Insiu de Mahémaiques de Toulouse, June Page 41 General srucure: densiy dependence We have a sequence of Markov chains (n (N) ) indexed by N, ogeher wih funcions (f ) such ha E(n (N) +1 n(n) ) = Nf (n (N) /N). We hen define (X (N) D ) by X (N) = n (N) /N. We hope ha ) FDD (x ), where (x ) if X (N) 0 x 0 as N, hen (X (N) saisfies x +1 = f (x ) (he limiing deerminisic model).
57 MASCOS Insiu de Mahémaiques de Toulouse, June Page 42 General srucure: densiy dependence Nex we suppose ha here are funcions (s ) such ha Var(n (N) +1 n(n) ) = Ns(n (N) /N).
58 MASCOS Insiu de Mahémaiques de Toulouse, June Page 43 General srucure: densiy dependence Nex we suppose ha here are funcions (s ) such ha N Var(X (N) +1 X(N) ) = s(x (N) ).
59 MASCOS Insiu de Mahémaiques de Toulouse, June Page 44 General srucure: densiy dependence Nex we suppose ha here are funcions (s ) such ha Var(n (N) +1 n(n) ) = Ns (n (N) /N). We hen define (Z (N) ) by Z (N) = N(X (N) x ).
60 MASCOS Insiu de Mahémaiques de Toulouse, June Page 45 General srucure: densiy dependence Nex we suppose ha here are funcions (s ) such ha Var(Z (N) +1 X(N) ) = s (X (N) ). We hen define (Z (N) ) by Z (N) = N(X (N) x ).
61 MASCOS Insiu de Mahémaiques de Toulouse, June Page 46 General srucure: densiy dependence Nex we suppose ha here are funcions (s ) such ha Var(n (N) +1 n(n) ) = Ns (n (N) /N). We hen define (Z (N) ) by Z (N) hope ha if N(X (N) = N(X (N) x ). We ) FDD (Z ), 0 x 0 ) D z 0, hen (Z (N) where (Z ) is a Gaussian Markov chain wih Z 0 = z 0.
62 MASCOS Insiu de Mahémaiques de Toulouse, June Page 47 General srucure: densiy dependence Wha will be he form of his chain?
63 MASCOS Insiu de Mahémaiques de Toulouse, June Page 47 General srucure: densiy dependence Wha will be he form of his chain? Consider he ime-homogeneous case, f = f and s = s.
64 MASCOS Insiu de Mahémaiques de Toulouse, June Page 47 General srucure: densiy dependence Wha will be he form of his chain? Consider he ime-homogeneous case, f = f and s = s. Formally, by Taylor s heorem, f(x (N) ) f(x ) = (X (N) x )f (x ) + and so, since E(X (N) +1 X(N) ) = f(x (N) ) and x +1 = f(x ), E(Z (N) +1 ) = N (E(X (N) +1 ) f(x )) = f (x ) E(Z (N) ) +, suggesing ha E(Z +1 ) = a E(Z ), where a = f (x ).
65 MASCOS Insiu de Mahémaiques de Toulouse, June Page 48 General srucure: densiy dependence We have Var(X (N) +1 ) = Var(E(X(N) So, since N Var(X (N) +1 X(N) Var(Z (N) +1 ) = N Var(X(N) +1 X(N) )) + E(Var(X (N) +1 X(N) )). ) = s(x (N) ), +1 ) = N Var(f(X(N) )) + E(s(X (N) )) a 2 N Var(X (N) ) + E(s(X (N) )) (where a = f (x )) = a 2 Var(Z (N) ) + E(s(X (N) )), suggesing ha Var(Z +1 ) = a 2 Var(Z ) + s(x ).
66 MASCOS Insiu de Mahémaiques de Toulouse, June Page 48 General srucure: densiy dependence We have Var(X (N) +1 ) = Var(E(X(N) So, since N Var(X (N) +1 X(N) Var(Z (N) +1 ) = N Var(X(N) +1 X(N) )) + E(Var(X (N) +1 X(N) )). ) = s(x (N) ), +1 ) = N Var(f(X(N) )) + E(s(X (N) )) a 2 N Var(X (N) ) + E(s(X (N) )) (where a = f (x )) = a 2 Var(Z (N) ) + E(s(X (N) )), suggesing ha Var(Z +1 ) = a 2 Var(Z ) + s(x ). And, since (Z ) will be Markovian,...
67 MASCOS Insiu de Mahémaiques de Toulouse, June Page 49 General srucure: densiy dependence And, since (Z ) will be Markovian, we migh hope ha Z +1 = a Z + E (Z 0 = z 0 ), where a = f (x ) and E ( = 0, 1,...) are independen Gaussian random variables wih E N(0,s(x )).
68 MASCOS Insiu de Mahémaiques de Toulouse, June Page 49 General srucure: densiy dependence And, since (Z ) will be Markovian, we migh hope ha Z +1 = a Z + E (Z 0 = z 0 ), where a = f (x ) and E ( = 0, 1,...) are independen Gaussian random variables wih E N(0,s(x )). If x eq is a fixed poin of f, and N(X (N) 0 x eq ) z 0, hen we migh hope ha (Z (N) (Z ), where (Z ) is he AR-1 process defined by Z +1 = az + E, Z 0 = z 0, where a = f (x eq ) and E ( = 0, 1,...) are iid Gaussian N(0,s(x eq )) random variables. ) FDD
69 MASCOS Insiu de Mahémaiques de Toulouse, June Page 50 Convergence of Markov chains We can adap resuls of Alan Karr for our purpose. Karr, A.F. (1975) Weak convergence of a sequence of Markov chains. Probabiliy Theory and Relaed Fields 33, He considered a sequence of ime-homogeneous Markov chains (X (n) ) on a general sae space (Ω, F) = (E, E) N wih ransiion kernels (K n (x,a), x E,A E) and iniial disribuions (π n (A),A E). He proved ha if (i) π n π and (ii) x n x in E implies K n (x n, ) K(x, ), hen he corresponding probabiliy measures (P π n n ) on (Ω, F) also converge: P π n n P π.
70 MASCOS Insiu de Mahémaiques de Toulouse, June Page 51 N-pach models: convergence Theorem For he N-pach models wih c i = (i/n)c, if D x 0 as N, hen X (N) 0 (X (N) 1,X (N) 2,...,X (N) n ) D (x 1,x 2,...,x n ), for any finie sequence of imes 1, 2,..., n, where (x ) is defined by he recursion x +1 = f(x ) wih f(x) = (1 e)(1 + c c(1 e)x)x f(x) = (1 e)(1 + c cx)x (EC model) (CE model)
71 MASCOS Insiu de Mahémaiques de Toulouse, June Page 52 N-pach models: convergence Theorem If, addiionally, N(X (N) 0 x 0 ) D z 0, hen (Z (N) ) FDD (Z ), where (Z ) is he Gaussian Markov chain defined by Z +1 = f (x )Z + E (Z 0 = z 0 ), where E ( = 0, 1,...) are independen Gaussian random variables wih E N(0,s(x )) and s(x) = (1 e)[c(1 (1 e)x)(1 c(1 e)x) + e(1 + c 2c(1 e)x) 2 ]x (EC model) s(x) = (1 e)[e + c(1 x)(1 c(1 e)x)]x (CE model)
72 MASCOS Insiu de Mahémaiques de Toulouse, June Page 53 Simulaion: EC Model 1 Meapopulaion simulaion P = EC (N=100, x 0 =0.95, e =0.4, c =0.8) X (N)
73 MASCOS Insiu de Mahémaiques de Toulouse, June Page 54 Simulaion: EC Model (Deerminisic pah) 1 Meapopulaion simulaion P = EC (N=100, x 0 =0.95, e =0.4, c =0.8) Deerminisic pah X (N)
74 MASCOS Insiu de Mahémaiques de Toulouse, June Page 55 Simulaion: EC Model (Gaussian approx.) 1 Meapopulaion simulaion P = EC (N=100, x 0 =0.95, e =0.4, c =0.8) Deerminisic pah ± wo sandard deviaions X (N)
75 MASCOS Insiu de Mahémaiques de Toulouse, June Page 56 N-pach models: convergence In boh cases (EC and CE) he deerminisic model has wo equilibria, x = 0 and x = x, given by x = 1 1 e x = 1 ( 1 e c(1 e) e c(1 e) ) (EC model) (CE model)
76 MASCOS Insiu de Mahémaiques de Toulouse, June Page 56 N-pach models: convergence In boh cases (EC and CE) he deerminisic model has wo equilibria, x = 0 and x = x, given by x = 1 1 e x = 1 ( 1 e c(1 e) e c(1 e) ) (EC model) (CE model) Indeed, we may wrie f(x) = x (1 + r (1 x/x )), r = c(1 e) e for boh models (he form of he discree-ime logisic model), and we obain he condiion c > e/(1 e) for x o be posiive and hen sable.
77 MASCOS Insiu de Mahémaiques de Toulouse, June Page 57 N-pach models: convergence Corollary If c > e/(1 e), so ha x given above is sable, and N(X (N) 0 x ) D z 0, hen (Z (N) where (Z ) is he AR-1 process defined by ) FDD (Z ), Z +1 = (1 + e c(1 e))z + E (Z 0 = z 0 ), where E ( = 0, 1,...) are independen Gaussian N(0,σ 2 ) random variables wih σ 2 = (1 e)[c(1 (1 e)x )(1 c(1 e)x ) + e(1 + c 2c(1 e)x ) 2 ]x (EC model) σ 2 = (1 e)[e + c(1 x )(1 c(1 e)x )]x (CE model)
78 MASCOS Insiu de Mahémaiques de Toulouse, June Page 58 Simulaion: EC Model 1 Meapopulaion simulaion P = EC (N=100, x 0 =0.95, e =0.3, c =0.8) X (N) x =
79 MASCOS Insiu de Mahémaiques de Toulouse, June Page 59 Simulaion: EC Model (AR-1 approx.) 1 Meapopulaion simulaion P = EC (N=100, x 0 =0.95, e =0.3, c =0.8) X (N) x =
80 MASCOS Insiu de Mahémaiques de Toulouse, June Page 60 AR-1 Simulaion: EC Model 1 AR-1 simulaion P = EC (N=100, x 0 = , e =0.3, c =0.8) X (N) x =
81 MASCOS Insiu de Mahémaiques de Toulouse, June Page 61 Recen developmens Buckley, F.M. and Polle, P.K. (2010) Limi heorems for discree-ime meapopulaion models. Probabiliy Surveys 7,
82 MASCOS Insiu de Mahémaiques de Toulouse, June Page 61 Recen developmens Buckley, F.M. and Polle, P.K. (2010) Limi heorems for discree-ime meapopulaion models. Probabiliy Surveys 7, A general heory of convergence for sequences of ime-inhomogeneous densiy-dependen Markov chains.
83 MASCOS Insiu de Mahémaiques de Toulouse, June Page 61 Recen developmens Buckley, F.M. and Polle, P.K. (2010) Limi heorems for discree-ime meapopulaion models. Probabiliy Surveys 7, A general heory of convergence for sequences of ime-inhomogeneous densiy-dependen Markov chains. Analysis of he scheme n +1 = ñ + Bin(N ñ,c(ñ /N)) ñ = n Bin(n,e) (EC) n +1 = ñ Bin(ñ,e) ñ = n + Bin(N n,c(n /N)), (CE) where c is coninuous, increasing and concave, wih c(0) 0 and c(x) 1.
84 MASCOS Insiu de Mahémaiques de Toulouse, June Page 62 Recen developmens Sabiliy analysis of he limiing deerminisic model: (i) Saionariy: c(0) > 0. (ii) Evanescence: c(0) = 0 and c (0) e/(1 e). (iii) Quasi saionariy: c(0) = 0 and c (0) > e/(1 e).
85 MASCOS Insiu de Mahémaiques de Toulouse, June Page 62 Recen developmens Sabiliy analysis of he limiing deerminisic model: (i) Saionariy: c(0) > 0. (ii) Evanescence: c(0) = 0 and c (0) e/(1 e). (iii) Quasi saionariy: c(0) = 0 and c (0) > e/(1 e). Infinie-pach models. If c(0) = 0 and c(x) has a coninuous second derivaive near 0, hen Bin(N n,c(n/n)) D Poi(mn) as N, where m = c (0).
86 MASCOS Insiu de Mahémaiques de Toulouse, June Page 62 Recen developmens Sabiliy analysis of he limiing deerminisic model: (i) Saionariy: c(0) > 0. (ii) Evanescence: c(0) = 0 and c (0) e/(1 e). (iii) Quasi saionariy: c(0) = 0 and c (0) > e/(1 e). Infinie-pach models. If c(0) = 0 and c(x) has a coninuous second derivaive near 0, hen Bin(N n,c(n/n)) D Poi(mn) as N, where m = c (0). This leads o he scheme n +1 = ñ + Poi(mñ ) ñ = n Bin(n,e) (EC) n +1 = ñ Bin(ñ,e) ñ = n + Poi(mn ), (CE)
87 MASCOS Insiu de Mahémaiques de Toulouse, June Page 63 Recen developmens...which urns ou o be a (Galon-Wason) branching process.
88 MASCOS Insiu de Mahémaiques de Toulouse, June Page 63 Recen developmens...which urns ou o be a (Galon-Wason) branching process. Analysis of he more general scheme n +1 = ñ + Poi(m(ñ )) ñ = n Bin(n,e) (EC) n +1 = ñ Bin(ñ,e) ñ = n + Poi(m(n )), (CE) assuming m(n) = n 0 µ(n/n 0 ).
89 MASCOS Insiu de Mahémaiques de Toulouse, June Page 63 Recen developmens...which urns ou o be a (Galon-Wason) branching process. Analysis of he more general scheme n +1 = ñ + Poi(m(ñ )) ñ = n Bin(n,e) (EC) n +1 = ñ Bin(ñ,e) ñ = n + Poi(m(n )), (CE) assuming m(n) = n 0 µ(n/n 0 ). In he limi as n 0 X (N) := n /n 0 has a deerminisic approximaion ha can exhibi he full range of dynamic behaviour (including chaos).
90 MASCOS Insiu de Mahémaiques de Toulouse, June Page 64 Ricker dynamics: µ(x) = x exp(r(1-x)) 1 (a) 1 (b) x x (c) 3 (d) x 1 x
Limit theorems for discrete-time metapopulation models
MASCOS AustMS Meeting, October 2008 - Page 1 Limit theorems for discrete-time metapopulation models Phil Pollett Department of Mathematics The University of Queensland http://www.maths.uq.edu.au/ pkp AUSTRALIAN
More informationMetapopulations with infinitely many patches
Metapopulations with infinitely many patches Phil. Pollett The University of Queensland UQ ACEMS Research Group Meeting 10th September 2018 Phil. Pollett (The University of Queensland) Infinite-patch metapopulations
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationStochastic Model for Cancer Cell Growth through Single Forward Mutation
Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More information5. Stochastic processes (1)
Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly
More informationProblem Set on Differential Equations
Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()
More informationDiscrete Markov Processes. 1. Introduction
Discree Markov Processes 1. Inroducion 1. Probabiliy Spaces and Random Variables Sample space. A model for evens: is a family of subses of such ha c (1) if A, hen A, (2) if A 1, A 2,..., hen A1 A 2...,
More informationReliability of Technical Systems
eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationBifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays
Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi
More informationMath 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.
Mah 36. Rumbos Spring 1 1 Soluions o Assignmen #6 1. Suppose he growh of a populaion is governed by he differenial equaion where k is a posiive consan. d d = k (a Explain why his model predics ha he populaion
More informationnon -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.
LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis
More informationLECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS
LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN
Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.
More information. Now define y j = log x j, and solve the iteration.
Problem 1: (Disribued Resource Allocaion (ALOHA!)) (Adaped from M& U, Problem 5.11) In his problem, we sudy a simple disribued proocol for allocaing agens o shared resources, wherein agens conend for resources
More informationOn Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction
On Mulicomponen Sysem Reliabiliy wih Microshocks - Microdamages Type of Componens Ineracion Jerzy K. Filus, and Lidia Z. Filus Absrac Consider a wo componen parallel sysem. The defined new sochasic dependences
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More information20. Applications of the Genetic-Drift Model
0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0
More informationStochastic models and their distributions
Sochasic models and heir disribuions Couning cusomers Suppose ha n cusomers arrive a a grocery a imes, say T 1,, T n, each of which akes any real number in he inerval (, ) equally likely The values T 1,,
More informationMean-square Stability Control for Networked Systems with Stochastic Time Delay
JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455
More informationAn Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.
1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard
More informationRight tail. Survival function
Densiy fi (con.) Lecure 4 The aim of his lecure is o improve our abiliy of densiy fi and knowledge of relaed opics. Main issues relaed o his lecure are: logarihmic plos, survival funcion, HS-fi mixures,
More informationModeling Economic Time Series with Stochastic Linear Difference Equations
A. Thiemer, SLDG.mcd, 6..6 FH-Kiel Universiy of Applied Sciences Prof. Dr. Andreas Thiemer e-mail: andreas.hiemer@fh-kiel.de Modeling Economic Time Series wih Sochasic Linear Difference Equaions Summary:
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationTransform Techniques. Moment Generating Function
Transform Techniques A convenien way of finding he momens of a random variable is he momen generaing funcion (MGF). Oher ransform echniques are characerisic funcion, z-ransform, and Laplace ransform. Momen
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationStable approximations of optimal filters
Sable approximaions of opimal filers Joaquin Miguez Deparmen of Signal Theory & Communicaions, Universidad Carlos III de Madrid. E-mail: joaquin.miguez@uc3m.es Join work wih Dan Crisan (Imperial College
More informationCash Flow Valuation Mode Lin Discrete Time
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics
More informationStochastic Structural Dynamics. Lecture-6
Sochasic Srucural Dynamics Lecure-6 Random processes- Dr C S Manohar Deparmen of Civil Engineering Professor of Srucural Engineering Indian Insiue of Science Bangalore 560 0 India manohar@civil.iisc.erne.in
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationA Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients
mahemaics Aricle A Noe on he Equivalence of Fracional Relaxaion Equaions o Differenial Equaions wih Varying Coefficiens Francesco Mainardi Deparmen of Physics and Asronomy, Universiy of Bologna, and he
More informationFamilies with no matchings of size s
Families wih no machings of size s Peer Franl Andrey Kupavsii Absrac Le 2, s 2 be posiive inegers. Le be an n-elemen se, n s. Subses of 2 are called families. If F ( ), hen i is called - uniform. Wha is
More informationExponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits
DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,
More informationInstitute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler
MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion
More information= ( ) ) or a system of differential equations with continuous parametrization (T = R
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationExcel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand
Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338
More informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
More informationSchool and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011
2229-12 School and Workshop on Marke Microsrucure: Design, Efficiency and Saisical Regulariies 21-25 March 2011 Some mahemaical properies of order book models Frederic ABERGEL Ecole Cenrale Paris Grande
More informationTHE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University
THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationOn Gronwall s Type Integral Inequalities with Singular Kernels
Filoma 31:4 (217), 141 149 DOI 1.2298/FIL17441A Published by Faculy of Sciences and Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma On Gronwall s Type Inegral Inequaliies
More information6. Stochastic calculus with jump processes
A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio
More information( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationThe Arcsine Distribution
The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we
More informationKeywords: competition models; density-dependence; ecology; population dynamics; predation models; stochastic models UNESCO EOLSS
POPULATIO MODELS Michael B. Bonsall Deparmen of Zoology, Universiy of Oxford, Oxford, UK Keywords: compeiion models; densiy-dependence; ecology; populaion dynamics; predaion models; sochasic models Conens.
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More informationOn a Fractional Stochastic Landau-Ginzburg Equation
Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn
More informationSample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen
Sample Auocorrelaions for Financial Time Series Models Richard A. Davis Colorado Sae Universiy Thomas Mikosch Universiy of Copenhagen Ouline Characerisics of some financial ime series IBM reurns NZ-USA
More informationLecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT
Lecure 9: Advanced DFT conceps: The Exchange-correlaion funcional and ime-dependen DFT Marie Curie Tuorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dep. of Chemisry and Couran Insiue
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationMANY FACET, COMMON LATENT TRAIT POLYTOMOUS IRT MODEL AND EM ALGORITHM. Dimitar Atanasov
Pliska Sud. Mah. Bulgar. 20 (2011), 5 12 STUDIA MATHEMATICA BULGARICA MANY FACET, COMMON LATENT TRAIT POLYTOMOUS IRT MODEL AND EM ALGORITHM Dimiar Aanasov There are many areas of assessmen where he level
More informationarxiv: v1 [math.ca] 15 Nov 2016
arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy
More informationMARKOV STOCHASTIC PROCESSES IN BIOLOGY AND MATHEMATICS THE SAME, AND YET DIFFERENT
Journal of Advances in Mahemaics MARKOV STOCHASTIC PROCESSES IN BIOLOGY AND MATHEMATICS THE SAME, AND YET DIFFERENT MI los lawa SOKÓ l Absrac. Virually every biological model uilising a random number generaor
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More informationMost Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation
Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,
More informationCosmic String Loop Distribution with a Gravitational Wave Cutoff
Ouline Inroducion Towards he Cosmological Aracor Cosmic Sring Loop Disribuion wih a Graviaional Wave Cuoff Larissa Lorenz Insiue of Mahemaics and Physics Cenre for Cosmology, Paricle Physics and Phenomenology
More informationST2352. Stochastic Processes constructed via Conditional Simulation. 09/02/2014 ST2352 Week 4 1
ST35 Sochasic Processes consruced via Condiional Simulaion 09/0/014 ST35 Week 4 1 Sochasic Processes consruced via Condiional Simulaion Markov Processes Simulaing Random Tex Google Sugges n grams Random
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
More informationA New Perturbative Approach in Nonlinear Singularity Analysis
Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The
More informationIntroduction to Probability and Statistics Slides 4 Chapter 4
Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random
More informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationModule 4: Time Response of discrete time systems Lecture Note 2
Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model
More informationAppendix to Creating Work Breaks From Available Idleness
Appendix o Creaing Work Breaks From Available Idleness Xu Sun and Ward Whi Deparmen of Indusrial Engineering and Operaions Research, Columbia Universiy, New York, NY, 127; {xs2235,ww24}@columbia.edu Sepember
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
More informationNature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.
Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike
More informationOptima and Equilibria for Traffic Flow on a Network
Opima and Equilibria for Traffic Flow on a Nework Albero Bressan Deparmen of Mahemaics, Penn Sae Universiy bressan@mah.psu.edu Albero Bressan (Penn Sae) Opima and equilibria for raffic flow 1 / 1 A Traffic
More informationResearch Article Dual Synchronization of Fractional-Order Chaotic Systems via a Linear Controller
The Scienific World Journal Volume 213, Aricle ID 159194, 6 pages hp://dx.doi.org/1155/213/159194 Research Aricle Dual Synchronizaion of Fracional-Order Chaoic Sysems via a Linear Conroller Jian Xiao,
More informationUnit Root Time Series. Univariate random walk
Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he
More informationTesting the Random Walk Model. i.i.d. ( ) r
he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp
More informationI. Return Calculations (20 pts, 4 points each)
Universiy of Washingon Spring 015 Deparmen of Economics Eric Zivo Econ 44 Miderm Exam Soluions This is a closed book and closed noe exam. However, you are allowed one page of noes (8.5 by 11 or A4 double-sided)
More informationComparison between the Discrete and Continuous Time Models
Comparison beween e Discree and Coninuous Time Models D. Sulsky June 21, 2012 1 Discree o Coninuous Recall e discree ime model Î = AIS Ŝ = S Î. Tese equaions ell us ow e populaion canges from one day o
More informationCHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence
CHEMICL KINETICS: Rae Order Rae law Rae consan Half-life Temperaure Dependence Chemical Reacions Kineics Chemical ineics is he sudy of ime dependence of he change in he concenraion of reacans and producs.
More informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationPersistence and non-persistence of a stochastic food chain model with finite delay
Available online a www.isr-publicaions.com/jnsa J. Nonlinear Sci. Appl., (7), 7 87 Research Aricle Journal Homepage: www.jnsa.com - www.isr-publicaions.com/jnsa Persisence non-persisence of a sochasic
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More information1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC
This documen was generaed a :37 PM, 1/11/018 Copyrigh 018 Richard T. Woodward 1. An inroducion o dynamic opimiaion -- Opimal Conrol and Dynamic Programming AGEC 64-018 I. Overview of opimiaion Opimiaion
More informationComparing Means: t-tests for One Sample & Two Related Samples
Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion
More informationProperties Of Solutions To A Generalized Liénard Equation With Forcing Term
Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick
More informationON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES WITH IMMIGRATION WITH FAMILY SIZES WITHIN RANDOM INTERVAL
ON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES ITH IMMIGRATION ITH FAMILY SIZES ITHIN RANDOM INTERVAL Husna Hasan School of Mahemaical Sciences Universii Sains Malaysia, 8 Minden, Pulau Pinang, Malaysia
More informationMulti-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More informationComparison of Approximation Schemes in Stochastic Simulation Methods for Stiff Chemical Systems
Comparison of Approximaion Schemes in Sochasic Simulaion Mehods for Siff Chemical Sysems by Chad Richard Wells A hesis presened o he Universiy of Waerloo in fulfillmen of he hesis requiremen for he degree
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationGeneralized Snell envelope and BSDE With Two general Reflecting Barriers
1/22 Generalized Snell envelope and BSDE Wih Two general Reflecing Barriers EL HASSAN ESSAKY Cadi ayyad Universiy Poly-disciplinary Faculy Safi Work in progress wih : M. Hassani and Y. Ouknine Iasi, July
More informationComparing Theoretical and Practical Solution of the First Order First Degree Ordinary Differential Equation of Population Model
Open Access Journal of Mahemaical and Theoreical Physics Comparing Theoreical and Pracical Soluion of he Firs Order Firs Degree Ordinary Differenial Equaion of Populaion Model Absrac Populaion dynamics
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationBBP-type formulas, in general bases, for arctangents of real numbers
Noes on Number Theory and Discree Mahemaics Vol. 19, 13, No. 3, 33 54 BBP-ype formulas, in general bases, for arcangens of real numbers Kunle Adegoke 1 and Olawanle Layeni 2 1 Deparmen of Physics, Obafemi
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. XI Control of Stochastic Systems - P.R. Kumar
CONROL OF SOCHASIC SYSEMS P.R. Kumar Deparmen of Elecrical and Compuer Engineering, and Coordinaed Science Laboraory, Universiy of Illinois, Urbana-Champaign, USA. Keywords: Markov chains, ransiion probabiliies,
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More informationPhysics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution
Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his
More informationOn Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems
MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 2, May 2013, pp. 209 227 ISSN 0364-765X (prin) ISSN 1526-5471 (online) hp://dx.doi.org/10.1287/moor.1120.0562 2013 INFORMS On Boundedness of Q-Learning Ieraes
More informationarxiv: v1 [math.fa] 9 Dec 2018
AN INVERSE FUNCTION THEOREM CONVERSE arxiv:1812.03561v1 [mah.fa] 9 Dec 2018 JIMMIE LAWSON Absrac. We esablish he following converse of he well-known inverse funcion heorem. Le g : U V and f : V U be inverse
More information1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population
Scienific Comuing I Module : Poulaion Modelling Coninuous Models Michael Bader Par I ODE Models Lehrsuhl Informaik V Winer 7/ Discree vs. Coniuous Models d d = F,,...) ) =? discree model: coninuous model:
More informationBiol. 356 Lab 8. Mortality, Recruitment, and Migration Rates
Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese
More information