Stein s Method for Steady-State Approximations: A Toolbox of Techniques
|
|
- Grant Malone
- 5 years ago
- Views:
Transcription
1 Stein s Method for Steady-State Approximations: A Toolbox of Techniques Anton Braverman Based on joint work with Jim Dai (Cornell), Itai Gurvich (Cornell), and Junfei Huang (CUHK). November 17, 2017
2 Outline Toolbox for moment bounds: E X( ) C. Toolbox for gradient bounds: D 2 f h (x), D 3 f h (x). MDPs and control
3 Foster-Lyapunov Condition {X(t)} positive recurrent Markov chain V (x) R + is a Lyapunov function satisfying E x V (X(1)) V (x) cf (x) + b, where b, c, f (x) 0. Implies Ef (X( )) b/c. Problem: finding such a Lyapunov function is a bottleneck for many problems.
4 Join the Shortest Queue (JSQ)
5 JSQ Model Buffer exp(1) Poisson(nλ) arrival process Buffer Buffer exp(1) exp(1). Buffer exp(1) n servers, each has own buffer/queue, service rate µ = 1. Customer joins shortest queue, ties broken uniformly. Halfin-Whitt regime: λ = 1 β/ n.
6 State Space. X 4 (t) X 3 (t) X 2 (t) X 1 (t) X i (t) = number of servers with at least i customers at time t. {(X 1 (t), X 2 (t),...) Z } t 0 is a Markov chain.
7 Dynamics - Arrivals X 3 (t) X 2 (t) X 1 (t)
8 Dynamics - Arrivals X 3 (t) X 2 (t) X 1 (t)
9 Dynamics - Arrivals X 3 (t) X 2 (t) X 1 (t)
10 Dynamics - Arrivals X 3 (t) X 2 (t) X 1 (t)
11 Dynamics - Arrivals X 3 (t) X 2 (t) X 1 (t)
12 Dynamics - Arrivals X 3 (t) X 2 (t) X 1 (t)
13 Dynamics - Arrivals X 3 (t) X 2 (t) X 1 (t) General arrival transition: X i X i + 1 with rate nλ1(x 1 = X 2 =... = X i 1 = n, X i < n).
14 Dynamics - Departures X 3 (t) X 2 (t) X 1 (t)
15 Dynamics - Departures X 3 (t) X 2 (t) X 1 (t)
16 Dynamics - Departures X 3 (t) X 2 (t) X 1 (t) X 1 X 1 1 with rate X 1 X 2 = 2.
17 Dynamics - Departures X 3 (t) X 2 (t) X 1 (t) X 1 X 1 1 with rate X 1 X 2 = 2. X 2 X 2 1 with rate X 2 X 3 = 2.
18 Dynamics - Departures X 3 (t) X 2 (t) X 1 (t) X 1 X 1 1 with rate X 1 X 2 = 2. X 2 X 2 1 with rate X 2 X 3 = 2. X 3 X 3 1 with rate X 3 X 4 = 1.
19 Dynamics - Departures X 3 (t) X 2 (t) X 1 (t) X 1 X 1 1 with rate X 1 X 2 = 2. X 2 X 2 1 with rate X 2 X 3 = 2. X 3 X 3 1 with rate X 3 X 4 = 1. In general: X i X i 1 with rate X i X i+1.
20 Known Results Let X(t) = ( X 1 (t), X 2 (t),...) R, where X 1 (t) = δ(x 1 (t) n), idle servers X i (t) = δx i (t), i 2. Gamarnik & Eschenfeldt (2017) establish: Functional Law of Large Numbers for δ = 1/n. Functional Central Limit Theorem for δ = 1/ n. Steady-state convergence not proved.
21 State Space Collapse (SSC) X 3 (t) X 2 (t) X 1 (t) Above is a picture of an unlikely state: For X 3 to increase, X 1 + X 2 must build up to 2n up rate nλ = n β n vs. down rate n when X 1 = n. In steady-state, expect fluid equilibrium + noise: X 1 ( ) n = β n + O( n), X 2 ( ) = 0 + O( n), X i ( ) = 0 + O(1), i 3. number of idle servers number of non-empty buffers
22 Applying Stein to JSQ
23 Markov Chain Generator Recall that X 1 (t) = δ(x 1 (t) n) and X i (t) = δx i (t) for i 2. Consider a function f (x 1, x 2 ) : R 2 R. For x 1 = δ(q 1 n), and x i = δq i, CTMC generator G X f (x) = nλ1(q 1 < n) ( f (x 1 + δ, x 2 ) f (x) ) + nλ1(q 1 = n, q 2 < n) ( f (x 1, x 2 + δ) f (x) ) + (q 1 q 2 ) ( f (x 1 δ, x 2 ) f (x) ) + (q 2 q 3 ) ( f (x 1, x 2 δ) f (x) ).
24 Markov Chain Generator Recall that X 1 (t) = δ(x 1 (t) n) and X i (t) = δx i (t) for i 2. Consider a function f (x 1, x 2 ) : R 2 R. For x 1 = δ(q 1 n), and x i = δq i, CTMC generator G X f (x) = nλ1(q 1 < n) ( f (x 1 + δ, x 2 ) f (x) ) + nλ1(q 1 = n, q 2 < n) ( f (x 1, x 2 + δ) f (x) ) + (q 1 q 2 ) ( f (x 1 δ, x 2 ) f (x) ) + (q 2 q 3 ) ( f (x 1, x 2 δ) f (x) ). Example third order error term: δ 2 x 3 1 f (x)δ(q 1 n) = δ 2 x 3 1 f (x)x 1. Suitable moment bounds needed on E X 1.
25 Moment Bounds via Fluid Approximation For diffusion approximation (δ = 1/ n), need moment bounds E(n X 1 ( )) n C. Equivalently, can we show that E(n X 1 ( )) n β n C/ n? β/ n is the fluid equilibrium of E(n X 1 ( ))/n.
26 Moment Bounds via Fluid Approximation For diffusion approximation (δ = 1/ n), need moment bounds E(n X 1 ( )) n C. Equivalently, can we show that E(n X 1 ( )) n β n C/ n? β/ n is the fluid equilibrium of E(n X 1 ( ))/n. Stein s method for fluid, or mean-field models (δ = 1/n). Stolyar (2015), Ying (2016), Gast (2017).
27 First Order Taylor Expansion Going forward, δ = 1/n: X 1 (t) = X 1(t) n n [ 1, 0], X 2 (t) = X 2(t) n ( X 1 (t), X 2 (t)) T = {y 1 [ 1, 0], 0 y 2 y 1 + 1}. Fluid equilibrium for ( X 1 ( ), X 2 ( )) is ( β/ n, 0). [0, 1], 1 1
28 The Partial Differential Equation (PDE) G X f (x) ( x 1 + x 2 β/ n ) x1 f (x) x 2 x2 f (x) + λ1(x 1 = 0) ( x1 f (x) + x2 f (x) ). For h(x) with h( β/ n, 0) = 0, consider the first-order PDE G F f h (x) = h(x), x T, x1 f h (x) x2 f h (x) = 0, x 1 = 0, where G F f (x) = ( x 1 + x 2 β/ n ) x1 f (x) x 2 x2 f (x).
29 Taylor Expansion Stein step: Eh( X( )) = EG F f h ( X( )) = EG X f h( X( )) EG F f h ( X( )) The expected generator difference equals 1 [ 2n E x1 x 1 f h (ξ (1) ) ( 1 + λ1( X 1 ( ) < 0) + X 1 ( ) + X 2 ( ) )] + 1 [ 2n E x2 x 2 f h (ξ (2) ) ( λ1( X 1 ( ) = 0, X 2 ( ) < 1) + X 2 ( ) )] λp( X 1 ( ) = 0, X 2 ( ) = 1) x2 f h (0, 1) ] + 1 n [X E 3 ( ) x2 f h (ξ (3) SSC error term. ) With fluid scaling, moment bounds simpler: X i ( ) 1. Still need gradient bounds...
30 Fluid Model for JSQ G F f (x) = ( x 1 + x 2 β/ n ) x1 f (x) x 2 x2 f (x). Reflection condition x1 f h (x) x2 f h (x) = 0 for x 1 = 0. Suggests dynamical system y(t) T, where ẏ(t) = F (y) R 2, and vector field (F 1 (y), F 2 (y)) equals ( y1 + y 2 β/ ) n, y(t) T \ {y y 1 = 0, y 2 [β/ n, 1]} 2 ( ) 0 β/, y(t) {y n 1 = 0, y 2 [β/ n, 1]}.
31 F (y) = ( y 1 + y 2 β/ n, y 2 ) 1 β/ n 1 β/ n
32 F (y) = ( y 1 + y 2 β/ n, y 2 ) 1 β/ n 1 β/ n
33 F (y) = ( y 1 + y 2 β/ n, y 2 ) 1 β/ n 1 β/ n
34 F (y) = (0, β/ n) 1 β/ n 1 β/ n
35 Gradient Bounds Recall the PDE Lemma G F f h (x) = h(x), x T, x1 f (x) x2 f (x) = 0, x 1 = 0. For nice h(x) : T R with h( β/ n, 0) = 0, f h (x) = 0 h(y(t))dt, where y(t) solves ẏ(t) = F (y) with initial condition y(0) = x. Perturbation theory: Ying(2016), Gast (2017) gradient bounds for f h (x) when vector field F (y) is continuous and y(t) exponentially stable.
36 Gradient Bounds X 1 ( ) Lemma Fix κ > β and let h(x) = min(x 1 + κ/ n, 0). Then n xi x i f h (x) κ β, x 2 f h (0, 1) = 0, x2 f h (x) 0, x T. 0 Eh( X 1 ( )) = 1 [ 2n E x1 x 1 f h (ξ (1) ) ( 1 + λ1( X 1 ( ) < 0) + X 1 ( ) + X 2 ( ) )] + 1 [ 2n E x2 x 2 f h (ξ (2) ) ( λ1( X 1 ( ) = 0, X 2 ( ) < 1) + X 2 ( ) )] + 1 ] [X n E 3 ( ) x2 f h (ξ (3) ) 2 n(κ β). As a consequence, E X 1 ( ) n / n C(β).
37 Gradient Bounds X 2 ( ) Conjecture (Work in Progress) Fix κ > β and let h(x) = max(x 2 κ/ n, 0). Then xi f h (x) nc(β), xi x j f h (x) nc(β), x T. Can then show that EX 2 ( )/ ( ) n C(β) 1 + E(X 3 ( )). Conjecture: EX 3 ( ) C(β), n 1.
38 Summary To prove diffusion scale moment bounds: 1. Perform first order Taylor expansion (fluid approximation). Fluid-scale moment bounds are much easier. 2. First order PDE solved by f h (x) = 0 h(y(t))dt, y(0) = x. 3. For gradient bounds, either 3.1 Verify conditions in Ying (2016) or Gast (2017). 3.2 brute-force track the behavior of fluid system like for JSQ; see Stolyar (2015) for another example.
39 Gradient Bounds for Multi-Dimensional Diffusion Approximations
40 The Poisson Equation Consider the Poisson equation G Y f (x) = Eh(Y ( )) h(x), x R d. Goal: bound Df (x), D 2 f (x), and D 3 f (x). Leverage solution form f h (x) = 0 ( Ex h(y (t)) Eh(Y ( )) ) dt. 1. Use elliptic PDE theory (Gilbarg and Trudinger) 2. Coupling arguments 3. x E x Y (t), xx E x Y (t), xxx E x Y (t)?
41 Gradient Bounds Diffusion Setting: {Y (t)} t 0 is a diffusion process on R d. Lipschitz drift (b 1 (x),..., b d (x))r d R d b(x) b(y) c b x y. Diffusion coefficient A = (a ij ) symmetric positive definite. Generator G Y f (x) = d i=1 b i (x) xi f (x) di,j=1 a ij xi x j f (x)
42 Schauder Interior Estimates For x R d, let B x = {y R d : y x 1 1+ x }. Lemma (Gilbarg & Trudinger (2001)) Suppose f h (x) solves G Y f h (x) = Eh(Y ( )) h(x). There exists a constant C depending only on A and c b, such that Df h (x) + D 2 D 2 f h (z) D 2 f h (y) f h (x) + sup y,z B x,y z z y ( C sup f h (y) + Eh(Y ( )) y B x ) h(z) h(y) + sup h(y) + sup (1 + x ) 3. y B x y,z B x,y z z y Gurvich (2014) was first to apply this idea.
43 Using Geometric Ergodicity How to bound f h (x)? Recall f h (x) = Suppose we know that 0 ( Ex h(y (t)) Eh(Y ( )) ) dt. E x h(y (t)) Eh(Y ( )) V (x)e ηt, η > 0. This condition is known as V -geometric ergodicity. Then E x h(y (t)) Eh(Y ( )) dt 0 0 V (x)e ηt dt CV (x). Can combine with Schauder estimates on previous slide.
44 Schauder Estimates Drawbacks Difficulties of using Schauder estimates: General result, not specialized to any system. Geometric ergodicity can be difficult to verify in practice. Non-smooth PDE domains Do we have an alternative?
45 One Idea Coupling For any ε > 0, xi f h (x) f h(x + εe (i) ) f h (x) ε = 1 ( Ex+εe (i)h(y (t)) E x h(y (t)) ) dt. ε 0 where e (i) R d equals 1 in the ith component, and zero elsewhere. Main idea: construct a coupling {Ȳ (t)} t 0 such that 1. Ȳ (0) = x + εe (i) 2. Ȳ (t) Y (t) = O(ε). 3. Mixing time of {Ȳ (t)} and {Y (t)} can be bounded.
46 Examples in the Literature Barbour 1988 Poisson process approximation Barbour 1990 Multivariate Normal approximation Mackey & Gorham 2016 Log-concave distributions Gan, Röllin, Ross 2017 Dirichlet approximation
47 Networks of Single Server Stations
48 A G/G/1 Queue Consider a single-server queue operating under first-come-first-serve discipline. A, A 1, A 2,... i.i.d. inter-arrival times with mean 1/λ > 1. S, S 1, S 2,... i.i.d. service times with mean 1. Traffic intensity ρ = λ < 1. Lindley recursion for waiting times: Lindley recursion for W n the nth customer s waiting time in queue: W n+1 = (W n + S n A n+1 ) +, x + := max(x, 0). A n, S n inter-arrival and service time of nth customer, respectively.
49 Generator Scaled steady-state customer waiting time Markov process generator where W ( ) = (1 ρ)w ( ). [ G W f (w) := E f ( (w + X) +)] f (w), w 0, X d = (1 ρ)(s A), EX = (1 ρ)2. ρ
50 Taylor Expansion Taylor expansion suggests diffusion generator bf (w) σ2 f (w) f (0)EX, where b = EX, and σ 2 = EX 2. Poisson equation G Y f h (w) = Eh(Y ( )) h(w), w 0, f h(0) = 0, where G Y f h (w) := bf (w) σ2 f (w).
51 Taylor Expansion Taylor expansion suggests diffusion generator bf (w) σ2 f (w) f (0)EX, where b = EX, and σ 2 = EX 2. Poisson equation G Y f h (w) = Eh(Y ( )) h(w), w 0, f h(0) = 0, where G Y f h (w) := bf (w) σ2 f (w). Gradient bounds via 1 ε 0 ( Ew+ε h(y (t)) E w h(y (t)) ) dt.
52 Reflected Brownian Motion Reflected Brownian motion on R + Y (t) = Y (0) bt + σb(t) + R(t). Regulator R(t) = sup 0 s t (Y (s)). BM( b, σ 2 ) 0 0 0
53 Coupling Recall f h (w) 1 ε Let 0 E w+ε h(y (t)) = E w h(ȳ ε (t)). Same source of randomness. ( Ew+ε h(y (t)) E w h(y (t)) ) dt. Y (t) =w bt + σb(t) + R(t), Ȳ ε (t) =w + ε bt + σb(t) + R(t) 0 0 0
54 Result Let τ w = inf{t 0 Y (t) = 0 given Y (0) = w}. Eτ w = w b. Standard result about Brownian motion. Lemma f h(w) = f h (w) = f h (w) = 0 + lim ] E w [h (Y (t))1(t τ w ) dt ] E w [h (Y (t))1(t τ w ) 0 1 ε 0 ε lim ε 0 ε dt ] E ε [h (Y (t))1(t τ ε ) dt ] E w [h (Y (t))1(t τ w ) dt 0 ] E ε [h (Y (t))1(t τ ε ) dt
55 Multidimensional SRBMs Consider the M/M/1 /M/1 tandem system, we are interested in the queue lengths. The approximating diffusion process is a two-dimensional semimartingale reflecting Brownian motion (SRBM) Y = {(Y 1 (t), Y 2 (t)) R 2 +, t 0}. See Williams (1995) for a review of SRBMs.
56 2-d SRBM ( ) ( ) Y1 (t) b1 = Y (0) + t + Σ Y 2 (t) b 2 b 1 ( ) ( B1 (t) + B 2 (t) b 1, b 2 > 0, {B(t)} is standard Brownian motion. Denote generator by G Y. R 1 (t) R 1 (t) + R 2 (t) )
57 f h (x) = 0 ( Ex h(y (t)) Eh(Y ( )) ) dt. xi f h (x) 1 ε 0 ( Ex+εe (i)h(y (t)) E x h(y (t)) ) dt.
58 Coupling Visualized x2 f h (x) = [ 0 E x x2 h(y (t))1(t τ x (2) ] (0)) dt. τ x (i) (s) = inf{t s Y i (t) = 0 given Y (0) = x}.
59 Coupling Visualized x1 f h (x) = [ ] E x x1 h(y (t))1(t τ x (1) (0)) dt E x [ x2 h(y (t))1(τ (1) x (0) t τ (2) x (τ (1) x )) ] dt.
60 First Derivatives Lemma x1 f h (x) = x2 f h (x) = [ ] E x x1 h(y (t))1(t τ x (1) (0)) dt E x [ x2 h(y (t))1(τ (1) x (0) t τ (2) x (τ (1) x )) E x [ x2 h(y (t))1(t τ (2) x (0)) ] dt. ] dt Expected hitting time bounds Dupuis & Williams (1994). Can t get explicit form of second derivative (can bound it though). To reach third derivative, one direction is to use higher order finite difference schemes: e.g. Barbour (1988), Gorham & Mackey (2016), Gan et al (2017).
61 Directional Derivatives of SRBM Lipshutz & Ramanan (2016) On directional derivatives of Skorokhod maps in convex polyhedral domains Characterize xi E x Y (t).
62 Conclusion Higher order derivatives tougher. One idea: apply first derivative bounds to Schauder estimates Df h (x) + D 2 D 2 f h (z) D 2 f h (y) f h (x) + sup y,z B x,y z z y ( C sup f h (y) + Eh(Y ( )) y B x ) h(z) h(y) + sup h(y) + sup (1 + x ) 3. y B x z y y,z B x,y z Need to be careful about existence/regularity of PDE solution.
63 Summary of Gradient Bound Techniques Obtaining multi-dimensional gradient bounds is difficult, and problem specific. Coupling-based bounds Schauder interior estimates needs bound on f h (x) Geometric ergodicity. Bounds on Df h (x) from coupling. Sensitivity analysis: x E x Y (t), xx E x Y (t), xxx E x Y (t).
64 Applications to Control: Taylor Expansions of Value Functions Based on joint work with Itai Gurvich and Junfei Huang
65 MDP Overview Controlled discrete time Markov chain {X(t) S} t Z+. S Z d, S < (think S = rectangle). r(x, u) reward for performing action u at state x. E x,u (X(1)) = E ( X(1) X(0) = x, u(x(0)) = u ). U(x) set of feasible actions at state x S. Policy π : S U(x). Inifinite horizon discounted expected reward: V π (x) = α t ( ) E x r(x(t), π(x(t)), x S, t=0 α < 1 is the discount factor.
66 Bellman Equation Goal is to find Bellman equation: V (x) = max π (x) arg max π(x), π (optimal policy) V (x) = max π(x), π (optimal value). u U(x) ( ) r(x, u) + αe x,u V (X(1)), ( 0 = max r(x, u) + α ( E x,u V (X(1)) V (x) ) ) (1 α)v (x) u U(x) Idea: perform Taylor expansion on E x,u V (X(1)) V (x).
67 Taylor Expansion Focus on 1-d case for notational convenience. E x,u V (X(1)) V (x) = y S P (u) x,y V (y) V (x) 2nd order V (x) y S = V (x) E x,u (X(1) x) }{{} b u(x) P x,y (u) (y x) V (x) P x,y (u) (y x) 2 y S V ( ) 2 (x) E x,u X(1) x. }{{} a u(x)
68 Taylored Control Equation Bellman vs. Taylored: vs. ( ) V (x) = max r(x, u) + αe x,u V (X(1)), u U(x) π (x) arg max u U(x) (1 α) V (x) = max u U(x) ˆπ (x) arg max u U(x) ( ) r(x, u) + αe x,u V (X(1)), (r(x, u) + α ( b u (x) V (x) a u(x) 2 V (x) )) ( r(x, u) + α ( b u (x) V (x) a u(x) 2 V (x) )).
69 Main Result Let d π = max π(x) P y x (maximal jump size). x,y >0 Theorem (B, Gurvich, Huang (2017)) Suppose there exists a solution V C 3 (R d ), that V (x) C(1 + x m ), and that dˆπ, d π <. Then, ( V ) (x) V (x) Vˆπ (x) V (x) d 3ˆπ d 3 [ π (E x,ˆπ α n D 3 V ] Xn±d ˆπ n=0 n=0 [ ]) + E x,π α n D 3 V Xn±d π, x S. First guarantees of this type appeared in Gurvich & Huang (2016) for a controlled M/G/1 + G system.
70 So What? Went from MDP on discrete, finite state space S to continuous space control problem.
71 So What? Went from MDP on discrete, finite state space S to continuous space control problem. Computational implications of Tayloring: Up to us to choose discretization level to solve Taylored problem useful when S is very large. Makes search over action space tractable.
72 Policy Iteration (PI) 1. Start with some initial policy π (0) (x). 2. For k = 0, 1, find V (k) (x) that satisfies V (k) (x) = r(x, π (k) (x)) + αe x,π (k) (x)[ V (k) (X(1)) ], x S. This is known as the policy evaluation step. 2.2 Set { π (k+1) (x) = arg max r(x, u) + αe x,u [V (k) (X(1)) ]}, x S. u U(x) This is known as the policy improvement step.
73 Taylored Policy Iteration (TPI) 1. Start with initial policy π (0) (x). 2. For k = 0, 1, For x R, find V (k) (x) satisfying (1 α) V (k) (x) = (r(x, π (k) (x)) + α ( b π (x)(x) V (k) (x) + 1 (k) 2 a π (x)(x) 2 V (k) (x) )) (k) 2.2 Set ˆπ (k+1) (x) equal to ( arg max r(x, u) + α ( b π (x)(x) V (k) (x) + 1 u U(x) (k) 2 a π (x)(x) 2 V (k) (x) )). (k) Not just a heuristic.
74 Dupuis and Kushner (2013) Coarse Grid Solution (1 α) V ( (x) = max r(x, u) + α ( b u (x) V (x) + 1 u U(x) 2 a u(x) 2 V (x) )). Use finite differences V (x) ( V (x + h) V (x h))/2h, which leads to new Bellman equation with discount β α,h = (1 + (1 α)h 2 ) 1 : ( V h (x) = max h 2 ( β α,h r(x, u) + β α,h P h,u V x,x+h h (x + h) u U(x) Choice of h is up to us! Sparse transition matrix + P h,u x,x h V h (x h) + P h,u x,x V h (x) )).
75 Example (based on Dai & Shi 2017) Author: Multi-pool model with overflow 4 Manufacturing & Service Operations Management 00(0), pp , c 0000 INFORMS 1(t) 2(t) 3(t) 4(t) 5(t) GeMed Surg Ortho Card OtMed GeMed Surg Ortho Card OtMed Figure 1 Illustration of the five-class, five-pool system. Pools 1 to 5 are the primary server pools for General Medicine (GeMed), Surgery (Surg), Orthopedic (Ortho), Cardiology (Card), and Other Medicine (OtMed), respectively; λ solid black lines represent the primary bed assignment, dashed red lines represent i = 0.8N i p i λ i = N i p i the preferred overflow assignment, and dashed blue lines represent the secondary overflow assignment. α h TPI rel. error α h TPI rel. error waiting patient to a non-primary ward, and (b) which ward to use 2 if multiple non-primary wards have available beds. A 1primary bed assignment, if possible, is made 1 immediately upon a patient arrival or bed release since 2 this is the most ideal situation and thus, 2 it is excluded from the decisions at each epoch. We model the overflow decisions as a Markov decision process (MDP), where the Table: J = decision 2, N maker aims to minimize the long-run average cost over (infinitely) many epochs. The total 1 = N 2 = 10, M = 10, (p 1, p 2 ) = (0.56, 0.56), (H 1, H 2 ) = (1, 4), and (B 12, B 21 ) = (5, 1). cost in each epoch includes: the holding cost associated with the number of waiting patients, which captures the system Relative congestion; errors displayed are sup x S V (x) Vˆπ (x) the overflow cost associated with each overflow patient, which captures the non-desirableness of placing a patient in a non-primary ward. These two costs represent the key tradeo in the overflow decision and provide us a tuning knob
76 References Barbour (1988), Stein s Method and Poisson Process Convergence, Journal of Applied Probability, Vol. 25, pp Barbour (1990), Stein s method for diffusion approximations, Probability Theory and Related Fields, September 1990, Volume 84, Issue 3, pp Mackey & Gorham (2016), Multivariate Stein Factors for a Class of Strongly Log-concave Distributions. Electronic Communications in Probability 2016, Vol. 21, paper no. 56, Gan, Röllin & Ross (2017), Dirichlet approximation of equilibrium distributions in Cannings models with mutation. Advances in Applied Probability, Volume 49, Issue 3, September 2017, pp
77 references Lipshutz & Ramanan (2016), On directional derivatives of Skorokhod maps in convex polyhedral domains. To appear The Annals of Applied Probability. Kushner & Dupuis (2013), Numerical Methods for Stochastic Control Problems in Continuous Time. Gurvich & Huang (2016), Beyond Heavy-Traffic Regimes: Universal Bounds and Controls for the Single-Server Queue. Dai & Shi (2017), Inpatient Bed Overflow: An Approximate Dynamic Programming Approach.
Stein s Method for Steady-State Approximations: Error Bounds and Engineering Solutions
Stein s Method for Steady-State Approximations: Error Bounds and Engineering Solutions Jim Dai Joint work with Anton Braverman and Jiekun Feng Cornell University Workshop on Congestion Games Institute
More informationLoad Balancing in Distributed Service System: A Survey
Load Balancing in Distributed Service System: A Survey Xingyu Zhou The Ohio State University zhou.2055@osu.edu November 21, 2016 Xingyu Zhou (OSU) Load Balancing November 21, 2016 1 / 29 Introduction and
More informationStein s Method for Steady-State Approximations: Error Bounds and Engineering Solutions
Stein s Method for Steady-State Approximations: Error Bounds and Engineering Solutions Jim Dai Joint work with Anton Braverman, Jiekun Feng Cornell University and Pengyi Shi Purdue University ACM Sigmetrics
More informationClass 11 Non-Parametric Models of a Service System; GI/GI/1, GI/GI/n: Exact & Approximate Analysis.
Service Engineering Class 11 Non-Parametric Models of a Service System; GI/GI/1, GI/GI/n: Exact & Approximate Analysis. G/G/1 Queue: Virtual Waiting Time (Unfinished Work). GI/GI/1: Lindley s Equations
More informationarxiv: v1 [math.pr] 27 Apr 2017
STEIN S METHOD FOR STEADY-STATE DIFFUSION APPROXIMATIONS arxiv:174.8398v1 [math.pr] 7 Apr 17 A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment
More informationOther properties of M M 1
Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected
More informationDynamic Control of a Tandem Queueing System with Abandonments
Dynamic Control of a Tandem Queueing System with Abandonments Gabriel Zayas-Cabán 1 Jungui Xie 2 Linda V. Green 3 Mark E. Lewis 1 1 Cornell University Ithaca, NY 2 University of Science and Technology
More informationOn the Resource/Performance Tradeoff in Large Scale Queueing Systems
On the Resource/Performance Tradeoff in Large Scale Queueing Systems David Gamarnik MIT Joint work with Patrick Eschenfeldt, John Tsitsiklis and Martin Zubeldia (MIT) High level comments High level comments
More informationOverflow Networks: Approximations and Implications to Call-Center Outsourcing
Overflow Networks: Approximations and Implications to Call-Center Outsourcing Itai Gurvich (Northwestern University) Joint work with Ohad Perry (CWI) Call Centers with Overflow λ 1 λ 2 Source of complexity:
More informationLecture 20: Reversible Processes and Queues
Lecture 20: Reversible Processes and Queues 1 Examples of reversible processes 11 Birth-death processes We define two non-negative sequences birth and death rates denoted by {λ n : n N 0 } and {µ n : n
More informationMaximum pressure policies for stochastic processing networks
Maximum pressure policies for stochastic processing networks Jim Dai Joint work with Wuqin Lin at Northwestern Univ. The 2011 Lunteren Conference Jim Dai (Georgia Tech) MPPs January 18, 2011 1 / 55 Outline
More informationPerformance Evaluation of Queuing Systems
Performance Evaluation of Queuing Systems Introduction to Queuing Systems System Performance Measures & Little s Law Equilibrium Solution of Birth-Death Processes Analysis of Single-Station Queuing Systems
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationUniformly Uniformly-ergodic Markov chains and BSDEs
Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,
More informationOPEN MULTICLASS HL QUEUEING NETWORKS: PROGRESS AND SURPRISES OF THE PAST 15 YEARS. w 1. v 2. v 3. Ruth J. Williams University of California, San Diego
OPEN MULTICLASS HL QUEUEING NETWORKS: PROGRESS AND SURPRISES OF THE PAST 15 YEARS v 2 w3 v1 w 2 w 1 v 3 Ruth J. Williams University of California, San Diego 1 PERSPECTIVE MQN SPN Sufficient conditions
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 6
MATH 56A: STOCHASTIC PROCESSES CHAPTER 6 6. Renewal Mathematically, renewal refers to a continuous time stochastic process with states,, 2,. N t {,, 2, 3, } so that you only have jumps from x to x + and
More informationAsymptotic Coupling of an SPDE, with Applications to Many-Server Queues
Asymptotic Coupling of an SPDE, with Applications to Many-Server Queues Mohammadreza Aghajani joint work with Kavita Ramanan Brown University March 2014 Mohammadreza Aghajanijoint work Asymptotic with
More informationElectronic Companion Fluid Models for Overloaded Multi-Class Many-Server Queueing Systems with FCFS Routing
Submitted to Management Science manuscript MS-251-27 Electronic Companion Fluid Models for Overloaded Multi-Class Many-Server Queueing Systems with FCFS Routing Rishi Talreja, Ward Whitt Department of
More informationTime Reversibility and Burke s Theorem
Queuing Analysis: Time Reversibility and Burke s Theorem Hongwei Zhang http://www.cs.wayne.edu/~hzhang Acknowledgement: this lecture is partially based on the slides of Dr. Yannis A. Korilis. Outline Time-Reversal
More informationIEOR 6711, HMWK 5, Professor Sigman
IEOR 6711, HMWK 5, Professor Sigman 1. Semi-Markov processes: Consider an irreducible positive recurrent discrete-time Markov chain {X n } with transition matrix P (P i,j ), i, j S, and finite state space.
More informationA Diffusion Approximation for Stationary Distribution of Many-Server Queueing System In Halfin-Whitt Regime
A Diffusion Approximation for Stationary Distribution of Many-Server Queueing System In Halfin-Whitt Regime Mohammadreza Aghajani joint work with Kavita Ramanan Brown University APS Conference, Istanbul,
More informationIntroduction to Markov Chains, Queuing Theory, and Network Performance
Introduction to Markov Chains, Queuing Theory, and Network Performance Marceau Coupechoux Telecom ParisTech, departement Informatique et Réseaux marceau.coupechoux@telecom-paristech.fr IT.2403 Modélisation
More informationStability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk
Stability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk ANSAPW University of Queensland 8-11 July, 2013 1 Outline (I) Fluid
More information(b) What is the variance of the time until the second customer arrives, starting empty, assuming that we measure time in minutes?
IEOR 3106: Introduction to Operations Research: Stochastic Models Fall 2006, Professor Whitt SOLUTIONS to Final Exam Chapters 4-7 and 10 in Ross, Tuesday, December 19, 4:10pm-7:00pm Open Book: but only
More informationOperations Research Letters. Instability of FIFO in a simple queueing system with arbitrarily low loads
Operations Research Letters 37 (2009) 312 316 Contents lists available at ScienceDirect Operations Research Letters journal homepage: www.elsevier.com/locate/orl Instability of FIFO in a simple queueing
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 informationCentral-limit approach to risk-aware Markov decision processes
Central-limit approach to risk-aware Markov decision processes Jia Yuan Yu Concordia University November 27, 2015 Joint work with Pengqian Yu and Huan Xu. Inventory Management 1 1 Look at current inventory
More informationMarkov processes and queueing networks
Inria September 22, 2015 Outline Poisson processes Markov jump processes Some queueing networks The Poisson distribution (Siméon-Denis Poisson, 1781-1840) { } e λ λ n n! As prevalent as Gaussian distribution
More informationState-dependent Limited Processor Sharing - Approximations and Optimal Control
State-dependent Limited Processor Sharing - Approximations and Optimal Control Varun Gupta University of Chicago Joint work with: Jiheng Zhang (HKUST) State-dependent Limited Processor Sharing Processor
More informationDesigning load balancing and admission control policies: lessons from NDS regime
Designing load balancing and admission control policies: lessons from NDS regime VARUN GUPTA University of Chicago Based on works with : Neil Walton, Jiheng Zhang ρ K θ is a useful regime to study the
More informationStability of the two queue system
Stability of the two queue system Iain M. MacPhee and Lisa J. Müller University of Durham Department of Mathematical Science Durham, DH1 3LE, UK (e-mail: i.m.macphee@durham.ac.uk, l.j.muller@durham.ac.uk)
More informationDynamic Control of Parallel-Server Systems
Dynamic Control of Parallel-Server Systems Jim Dai Georgia Institute of Technology Tolga Tezcan University of Illinois at Urbana-Champaign May 13, 2009 Jim Dai (Georgia Tech) Many-Server Asymptotic Optimality
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 22 12/09/2013. Skorokhod Mapping Theorem. Reflected Brownian Motion
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.7J Fall 213 Lecture 22 12/9/213 Skorokhod Mapping Theorem. Reflected Brownian Motion Content. 1. G/G/1 queueing system 2. One dimensional reflection mapping
More informationThe Transition Probability Function P ij (t)
The Transition Probability Function P ij (t) Consider a continuous time Markov chain {X(t), t 0}. We are interested in the probability that in t time units the process will be in state j, given that it
More informationChapter 5. Continuous-Time Markov Chains. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan
Chapter 5. Continuous-Time Markov Chains Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan Continuous-Time Markov Chains Consider a continuous-time stochastic process
More informationContinuous-Time Markov Chain
Continuous-Time Markov Chain Consider the process {X(t),t 0} with state space {0, 1, 2,...}. The process {X(t),t 0} is a continuous-time Markov chain if for all s, t 0 and nonnegative integers i, j, x(u),
More informationAsymptotics for Polling Models with Limited Service Policies
Asymptotics for Polling Models with Limited Service Policies Woojin Chang School of Industrial and Systems Engineering Georgia Institute of Technology Atlanta, GA 30332-0205 USA Douglas G. Down Department
More informationLecture 7: Simulation of Markov Processes. Pasi Lassila Department of Communications and Networking
Lecture 7: Simulation of Markov Processes Pasi Lassila Department of Communications and Networking Contents Markov processes theory recap Elementary queuing models for data networks Simulation of Markov
More information2905 Queueing Theory and Simulation PART III: HIGHER DIMENSIONAL AND NON-MARKOVIAN QUEUES
295 Queueing Theory and Simulation PART III: HIGHER DIMENSIONAL AND NON-MARKOVIAN QUEUES 16 Queueing Systems with Two Types of Customers In this section, we discuss queueing systems with two types of customers.
More informationStrategic Dynamic Jockeying Between Two Parallel Queues
Strategic Dynamic Jockeying Between Two Parallel Queues Amin Dehghanian 1 and Jeffrey P. Kharoufeh 2 Department of Industrial Engineering University of Pittsburgh 1048 Benedum Hall 3700 O Hara Street Pittsburgh,
More informationM/G/1 and M/G/1/K systems
M/G/1 and M/G/1/K systems Dmitri A. Moltchanov dmitri.moltchanov@tut.fi http://www.cs.tut.fi/kurssit/elt-53606/ OUTLINE: Description of M/G/1 system; Methods of analysis; Residual life approach; Imbedded
More informationOn the fast convergence of random perturbations of the gradient flow.
On the fast convergence of random perturbations of the gradient flow. Wenqing Hu. 1 (Joint work with Chris Junchi Li 2.) 1. Department of Mathematics and Statistics, Missouri S&T. 2. Department of Operations
More informationA Mean Field Approach for Optimization in Discrete Time
oname manuscript o. will be inserted by the editor A Mean Field Approach for Optimization in Discrete Time icolas Gast Bruno Gaujal the date of receipt and acceptance should be inserted later Abstract
More informationON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME
ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME SAUL D. JACKA AND ALEKSANDAR MIJATOVIĆ Abstract. We develop a general approach to the Policy Improvement Algorithm (PIA) for stochastic control problems
More informationDynamic Matching Models
Dynamic Matching Models Ana Bušić Inria Paris - Rocquencourt CS Department of École normale supérieure joint work with Varun Gupta, Jean Mairesse and Sean Meyn 3rd Workshop on Cognition and Control January
More informationNear-Optimal Control of Queueing Systems via Approximate One-Step Policy Improvement
Near-Optimal Control of Queueing Systems via Approximate One-Step Policy Improvement Jefferson Huang March 21, 2018 Reinforcement Learning for Processing Networks Seminar Cornell University Performance
More informationStochastic Gradient Descent in Continuous Time
Stochastic Gradient Descent in Continuous Time Justin Sirignano University of Illinois at Urbana Champaign with Konstantinos Spiliopoulos (Boston University) 1 / 27 We consider a diffusion X t X = R m
More informationData analysis and stochastic modeling
Data analysis and stochastic modeling Lecture 7 An introduction to queueing theory Guillaume Gravier guillaume.gravier@irisa.fr with a lot of help from Paul Jensen s course http://www.me.utexas.edu/ jensen/ormm/instruction/powerpoint/or_models_09/14_queuing.ppt
More informationTechnical Appendix for: When Promotions Meet Operations: Cross-Selling and Its Effect on Call-Center Performance
Technical Appendix for: When Promotions Meet Operations: Cross-Selling and Its Effect on Call-Center Performance In this technical appendix we provide proofs for the various results stated in the manuscript
More informationQ = (c) Assuming that Ricoh has been working continuously for 7 days, what is the probability that it will remain working at least 8 more days?
IEOR 4106: Introduction to Operations Research: Stochastic Models Spring 2005, Professor Whitt, Second Midterm Exam Chapters 5-6 in Ross, Thursday, March 31, 11:00am-1:00pm Open Book: but only the Ross
More informationStability and Asymptotic Optimality of h-maxweight Policies
Stability and Asymptotic Optimality of h-maxweight Policies - A. Rybko, 2006 Sean Meyn Department of Electrical and Computer Engineering University of Illinois & the Coordinated Science Laboratory NSF
More informationEconomy of Scale in Multiserver Service Systems: A Retrospective. Ward Whitt. IEOR Department. Columbia University
Economy of Scale in Multiserver Service Systems: A Retrospective Ward Whitt IEOR Department Columbia University Ancient Relics A. K. Erlang (1924) On the rational determination of the number of circuits.
More informationThe G/GI/N Queue in the Halfin-Whitt Regime I: Infinite Server Queue System Equations
The G/GI/ Queue in the Halfin-Whitt Regime I: Infinite Server Queue System Equations J. E Reed School of Industrial and Systems Engineering Georgia Institute of Technology October 17, 27 Abstract In this
More informationLecture 3: Hamilton-Jacobi-Bellman Equations. Distributional Macroeconomics. Benjamin Moll. Part II of ECON Harvard University, Spring
Lecture 3: Hamilton-Jacobi-Bellman Equations Distributional Macroeconomics Part II of ECON 2149 Benjamin Moll Harvard University, Spring 2018 1 Outline 1. Hamilton-Jacobi-Bellman equations in deterministic
More information1 IEOR 4701: Continuous-Time Markov Chains
Copyright c 2006 by Karl Sigman 1 IEOR 4701: Continuous-Time Markov Chains A Markov chain in discrete time, {X n : n 0}, remains in any state for exactly one unit of time before making a transition (change
More informationPositive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network
Positive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network Yoni Nazarathy a,1, Gideon Weiss a,1 a Department of Statistics, The University of Haifa, Mount Carmel 31905, Israel.
More informationTutorial: Optimal Control of Queueing Networks
Department of Mathematics Tutorial: Optimal Control of Queueing Networks Mike Veatch Presented at INFORMS Austin November 7, 2010 1 Overview Network models MDP formulations: features, efficient formulations
More informationQueues and Queueing Networks
Queues and Queueing Networks Sanjay K. Bose Dept. of EEE, IITG Copyright 2015, Sanjay K. Bose 1 Introduction to Queueing Models and Queueing Analysis Copyright 2015, Sanjay K. Bose 2 Model of a Queue Arrivals
More information1 Markov decision processes
2.997 Decision-Making in Large-Scale Systems February 4 MI, Spring 2004 Handout #1 Lecture Note 1 1 Markov decision processes In this class we will study discrete-time stochastic systems. We can describe
More informationA Robust Queueing Network Analyzer Based on Indices of Dispersion
A Robust Queueing Network Analyzer Based on Indices of Dispersion Wei You (joint work with Ward Whitt) Columbia University INFORMS 2018, Phoenix November 6, 2018 1/20 Motivation Many complex service systems
More informationSolutions to Homework Discrete Stochastic Processes MIT, Spring 2011
Exercise 6.5: Solutions to Homework 0 6.262 Discrete Stochastic Processes MIT, Spring 20 Consider the Markov process illustrated below. The transitions are labelled by the rate q ij at which those transitions
More information5 Lecture 5: Fluid Models
5 Lecture 5: Fluid Models Stability of fluid and stochastic processing networks Stability analysis of some fluid models Optimization of fluid networks. Separated continuous linear programming 5.1 Stability
More informationTCOM 501: Networking Theory & Fundamentals. Lecture 6 February 19, 2003 Prof. Yannis A. Korilis
TCOM 50: Networking Theory & Fundamentals Lecture 6 February 9, 003 Prof. Yannis A. Korilis 6- Topics Time-Reversal of Markov Chains Reversibility Truncating a Reversible Markov Chain Burke s Theorem Queues
More informationStein s Method for Mean-Field Approximations in Light and Heavy Traffic Regimes
Stein s Method for Mean-Field Approximations in Light and Heavy Traffic Regimes LEI YIG, School of Electrical, Computer and Energy Engineering, Arizona State University Mean-field analysis is an analytical
More informationQueueing Theory I Summary! Little s Law! Queueing System Notation! Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K "
Queueing Theory I Summary Little s Law Queueing System Notation Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K " Little s Law a(t): the process that counts the number of arrivals
More informationFluid Models of Parallel Service Systems under FCFS
Fluid Models of Parallel Service Systems under FCFS Hanqin Zhang Business School, National University of Singapore Joint work with Yuval Nov and Gideon Weiss from The University of Haifa, Israel Queueing
More informationConvergence to equilibrium for rough differential equations
Convergence to equilibrium for rough differential equations Samy Tindel Purdue University Barcelona GSE Summer Forum 2017 Joint work with Aurélien Deya (Nancy) and Fabien Panloup (Angers) Samy T. (Purdue)
More informationIEOR 6711: Stochastic Models I, Fall 2003, Professor Whitt. Solutions to Final Exam: Thursday, December 18.
IEOR 6711: Stochastic Models I, Fall 23, Professor Whitt Solutions to Final Exam: Thursday, December 18. Below are six questions with several parts. Do as much as you can. Show your work. 1. Two-Pump Gas
More 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 informationIntroduction to Queuing Networks Solutions to Problem Sheet 3
Introduction to Queuing Networks Solutions to Problem Sheet 3 1. (a) The state space is the whole numbers {, 1, 2,...}. The transition rates are q i,i+1 λ for all i and q i, for all i 1 since, when a bus
More informationA random perturbation approach to some stochastic approximation algorithms in optimization.
A random perturbation approach to some stochastic approximation algorithms in optimization. Wenqing Hu. 1 (Presentation based on joint works with Chris Junchi Li 2, Weijie Su 3, Haoyi Xiong 4.) 1. Department
More informationChapter 2. Poisson Processes. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan
Chapter 2. Poisson Processes Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan Outline Introduction to Poisson Processes Definition of arrival process Definition
More informationTOWARDS BETTER MULTI-CLASS PARAMETRIC-DECOMPOSITION APPROXIMATIONS FOR OPEN QUEUEING NETWORKS
TOWARDS BETTER MULTI-CLASS PARAMETRIC-DECOMPOSITION APPROXIMATIONS FOR OPEN QUEUEING NETWORKS by Ward Whitt AT&T Bell Laboratories Murray Hill, NJ 07974-0636 March 31, 199 Revision: November 9, 199 ABSTRACT
More informationTD Learning. Sean Meyn. Department of Electrical and Computer Engineering University of Illinois & the Coordinated Science Laboratory
TD Learning Sean Meyn Department of Electrical and Computer Engineering University of Illinois & the Coordinated Science Laboratory Joint work with I.-K. Cho, Illinois S. Mannor, McGill V. Tadic, Sheffield
More informationCPSC 531: System Modeling and Simulation. Carey Williamson Department of Computer Science University of Calgary Fall 2017
CPSC 531: System Modeling and Simulation Carey Williamson Department of Computer Science University of Calgary Fall 2017 Motivating Quote for Queueing Models Good things come to those who wait - poet/writer
More informationBeyond heavy-traffic regimes: Universal bounds and controls for the single-server (M/GI/1+GI) queue
1 / 33 Beyond heavy-traffic regimes: Universal bounds and controls for the single-server (M/GI/1+GI) queue Itai Gurvich Northwestern University Junfei Huang Chinese University of Hong Kong Stochastic Networks
More informationSynchronized Queues with Deterministic Arrivals
Synchronized Queues with Deterministic Arrivals Dimitra Pinotsi and Michael A. Zazanis Department of Statistics Athens University of Economics and Business 76 Patission str., Athens 14 34, Greece Abstract
More informationMaximizing throughput in zero-buffer tandem lines with dedicated and flexible servers
Maximizing throughput in zero-buffer tandem lines with dedicated and flexible servers Mohammad H. Yarmand and Douglas G. Down Department of Computing and Software, McMaster University, Hamilton, ON, L8S
More informationNEW FRONTIERS IN APPLIED PROBABILITY
J. Appl. Prob. Spec. Vol. 48A, 209 213 (2011) Applied Probability Trust 2011 NEW FRONTIERS IN APPLIED PROBABILITY A Festschrift for SØREN ASMUSSEN Edited by P. GLYNN, T. MIKOSCH and T. ROLSKI Part 4. Simulation
More informationTechnical Appendix for: When Promotions Meet Operations: Cross-Selling and Its Effect on Call-Center Performance
Technical Appendix for: When Promotions Meet Operations: Cross-Selling and Its Effect on Call-Center Performance In this technical appendix we provide proofs for the various results stated in the manuscript
More informationLatent voter model on random regular graphs
Latent voter model on random regular graphs Shirshendu Chatterjee Cornell University (visiting Duke U.) Work in progress with Rick Durrett April 25, 2011 Outline Definition of voter model and duality with
More informationStatistics 150: Spring 2007
Statistics 150: Spring 2007 April 23, 2008 0-1 1 Limiting Probabilities If the discrete-time Markov chain with transition probabilities p ij is irreducible and positive recurrent; then the limiting probabilities
More informationSeries Expansions in Queues with Server
Series Expansions in Queues with Server Vacation Fazia Rahmoune and Djamil Aïssani Abstract This paper provides series expansions of the stationary distribution of finite Markov chains. The work presented
More informationSingle-Server Service-Station (G/G/1)
Service Engineering July 997 Last Revised January, 006 Single-Server Service-Station (G/G/) arrivals queue 000000000000 000000000000 departures Arrivals A = {A(t), t 0}, counting process, e.g., completely
More informationLecture 9: Deterministic Fluid Models and Many-Server Heavy-Traffic Limits. IEOR 4615: Service Engineering Professor Whitt February 19, 2015
Lecture 9: Deterministic Fluid Models and Many-Server Heavy-Traffic Limits IEOR 4615: Service Engineering Professor Whitt February 19, 2015 Outline Deterministic Fluid Models Directly From Data: Cumulative
More informationQueues with Many Servers and Impatient Customers
MATHEMATICS OF OPERATIOS RESEARCH Vol. 37, o. 1, February 212, pp. 41 65 ISS 364-765X (print) ISS 1526-5471 (online) http://dx.doi.org/1.1287/moor.111.53 212 IFORMS Queues with Many Servers and Impatient
More informationName of the Student:
SUBJECT NAME : Probability & Queueing Theory SUBJECT CODE : MA 6453 MATERIAL NAME : Part A questions REGULATION : R2013 UPDATED ON : November 2017 (Upto N/D 2017 QP) (Scan the above QR code for the direct
More informationStructured Markov Chains
Structured Markov Chains Ivo Adan and Johan van Leeuwaarden Where innovation starts Book on Analysis of structured Markov processes (arxiv:1709.09060) I Basic methods Basic Markov processes Advanced Markov
More informationBlind Fair Routing in Large-Scale Service Systems with Heterogeneous Customers and Servers
OPERATIONS RESEARCH Vol. 6, No., January February 23, pp. 228 243 ISSN 3-364X (print) ISSN 526-5463 (online) http://dx.doi.org/.287/opre.2.29 23 INFORMS Blind Fair Routing in Large-Scale Service Systems
More informationState Space Collapse in Many-Server Diffusion Limits of Parallel Server Systems. J. G. Dai. Tolga Tezcan
State Space Collapse in Many-Server Diffusion imits of Parallel Server Systems J. G. Dai H. Milton Stewart School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia
More informationBRAVO for QED Queues
1 BRAVO for QED Queues Yoni Nazarathy, The University of Queensland Joint work with Daryl J. Daley, The University of Melbourne, Johan van Leeuwaarden, EURANDOM, Eindhoven University of Technology. Applied
More informationarxiv: v1 [cs.pf] 22 Dec 2017
Scalable Load Balancing in etworked Systems: Universality Properties and Stochastic Coupling Methods Mark van der Boor 1, Sem C. Borst 1,2, Johan S.H. van Leeuwaarden 1, and Debankur Mukherjee 1 1 Eindhoven
More informationExercises Stochastic Performance Modelling. Hamilton Institute, Summer 2010
Exercises Stochastic Performance Modelling Hamilton Institute, Summer Instruction Exercise Let X be a non-negative random variable with E[X ]
More informationQueuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem. Wade Trappe
Queuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem Wade Trappe Lecture Overview Network of Queues Introduction Queues in Tandem roduct Form Solutions Burke s Theorem What
More informationDesign, staffing and control of large service systems: The case of a single customer class and multiple server types. April 2, 2004 DRAFT.
Design, staffing and control of large service systems: The case of a single customer class and multiple server types. Mor Armony 1 Avishai Mandelbaum 2 April 2, 2004 DRAFT Abstract Motivated by modern
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationBlind Fair Routing in Large-Scale Service Systems with Heterogeneous Customers and Servers
Blind Fair Routing in Large-Scale Service Systems with Heterogeneous Customers and Servers Mor Armony Amy R. Ward 2 October 6, 2 Abstract In a call center, arriving customers must be routed to available
More informationOPTIMALITY OF RANDOMIZED TRUNK RESERVATION FOR A PROBLEM WITH MULTIPLE CONSTRAINTS
OPTIMALITY OF RANDOMIZED TRUNK RESERVATION FOR A PROBLEM WITH MULTIPLE CONSTRAINTS Xiaofei Fan-Orzechowski Department of Applied Mathematics and Statistics State University of New York at Stony Brook Stony
More informationPath Decomposition of Markov Processes. Götz Kersting. University of Frankfurt/Main
Path Decomposition of Markov Processes Götz Kersting University of Frankfurt/Main joint work with Kaya Memisoglu, Jim Pitman 1 A Brownian path with positive drift 50 40 30 20 10 0 0 200 400 600 800 1000-10
More informationRouting and Staffing in Large-Scale Service Systems: The Case of Homogeneous Impatient Customers and Heterogeneous Servers 1
Routing and Staffing in Large-Scale Service Systems: The Case of Homogeneous Impatient Customers and Heterogeneous Servers 1 Mor Armony 2 Avishai Mandelbaum 3 June 25, 2008 Abstract Motivated by call centers,
More information