Estimation of arrival and service rates for M/M/c queue system
|
|
- Clinton Carpenter
- 6 years ago
- Views:
Transcription
1 Estimation of arrival and service rates for M/M/c queue system Katarína Starinská Charles University Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics Prague, Czech Republic Katarína Starinská Arrival and Service Rates / 22
2 Content 1 Notation 2 Estimation for queues from queue length data OU approximation for M/M/c queue Simulation results 3 MLE for single server queues from waiting time data The likelihood function based on the waiting times The consistency and the asymptotic normality of the MLE Application to M/M/1 Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
3 Notation M/M/c The M/M/c queue system is a Markov process {X t, t 0} on the state space S = {0, 1,...} with non-zero transition rates q(m, m + 1) = λ, q(m, m 1) = µ min(m, c), where m S is the state of the process at time t. Poisson arrivals at rate λ independent exponentially distributed service times c servers each operating at rate µ c = 1 - estimating from waiting time data large c(> 40) - estimating from queue length data unknown parameters θ = (λ, µ) Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
4 Estimation for queues from queue length data Estimation for queues from queue length data {m c ( )}...family of Markov processes S c Z...state space of m c ( ) Q c = {q c (m, n), m, n S c }...transition rates of m c ( ) c...the number of servers in system (c > 40) Definition Suppose that there exists an open set E R k and a family {f c, c > 0} of continuous functions, with f c : Z k R, such that q c (m, m + l) = cf c ( m, l), l 0. c Then, the family of Markov chains is asymptotically density dependent if, additionally, there exists a function F : E R such that {F c } converges to F pointwise on E, where F c (x) = l lf c(x, l), x E. The transition rates of the corresponding density process X c (t) = mc(t), t > 0, depend on the present state m only through c Katarína Starinská m/c. (starinskak@gmail.com) Arrival and Service Rates / 22
5 Estimation for queues from queue length data Law of Large Numbers f c (, l) is bounded for each l and c F is Lipschitz continuous on E F c converges uniformly to F on E Then, if lim c X c (0) = x 0, the density process X c ( ) X (, x) uniformly in probability on [0, t], where X (, x) is the unique trajectory satisfying X (0, x) = x, X (s, x) E, s [0, t] and X (s, x) = F (X (s, x)). s Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
6 Estimation for queues from queue length data Central Limit Law f c (, l) is bounded F is Lipschitz continuous and has uniformly continuous first derivative on E lim c sup x E c Fc (x) F (x) = 0 G c (x) = l l 2 f c (x, l), x E converges uniformly to G, where G is uniformly continuous on E Let x 0 E. Then, if lim c c Xc (0) x 0 = z, the family of processes {Z c ( )}, Z c (s) = c(x c (s) X (s x 0 )), 0 s t, converges weakly in D[0, 1] to a Gaussian diffusion Z( ) with initial value Z(0) = z and with mean and variance given by µ s = EZ(s) = M s z = z exp( σ 2 s = VarZ(s) = M 2 s s 0 s 0 B u du), B s = F (X (s, x 0 )) Mu 2 G(X (u, x 0 ))du. Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
7 Estimation for queues from queue length data X c (s) has an approximate normal distribution with VarX c (s) σ 2 s /c. If we take x 0 = X c (0), then EX c (s) X (s, x 0 ). Corollary If F (x 0 ) = 0 then, under the conditions of Central Limit Law, the family {Z c (s)}, defined by Z c (s) = c(x c (s) x 0 ), 0 s t, converges weakly in D[0, t] to an Ornstein-Uhlenbeck process Z( ) with initial value Z(0) = z, local drift B = F (x 0 ) and local variance V = G(x 0 ). In particular, ( Z(s) N µ s = ze Bs, σs 2 = V ) 2B (e2bs 1) Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
8 Estimation for queues from queue length data For large c, X c (s) has an approximate normal distribution with EX c (s) x 0 + e Bs (X c (0) x 0 ) VarX c (s) σ 2 s /c In case of queueing models, we have B < 0 (x 0 asymptotically stable). OU process is strongly stationary if we start in equilibrium where σ 2 = V / 2B. Therefore Z(0) N(0, σ 2 ) X c (0) N(x 0, σ 2 /c) We may approximate Cov(X c (s), X c (s + t)) by c(t) = 1 σ2 Cov(Z(s), Z(s + t)) = c c exp(b t ) Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
9 Estimation for queues from queue length data The likelihood function of a Gaussian vector (X c (t 1 ), X c (t 2 ),..., X c (t n )) { 1 f (x) = (2π) n C exp 1 } 2 (x m)c 1 (x m) m = x 0 1 n C = (c i,j ) n i,j=1, c i,j = σ2 c exp{b t j t i } Inverse and determinant of C can be computed explicitly. Therefore we can get the maximum likelihood estimators of the parameters. Explicit calculation od the estimators is not practical, therefore we need a numerical optimisation procedure (Cross-Entropy method). Information about parameters from the data E m = cx 0, where m = 1 n n i=1 m c(t i ) E ( 1 n n i=1 (m c(t i ) m) 2) = cσ 2 Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
10 Estimation for queues from queue length data OU approximation for M/M/c queue OU approximation for M/M/c queue Suppose that λ = O(c) (that is λ αc, where α is a constant). Process is density-dependent with F (x) = α µx and G(x) = α + µx. F c (x) = l lf c(x, l). f c (x, l) is non-zero for l = ±1 (and l = 0) F c (x) = 1 f c (x, 1) 1 f c (x, 1) = λ c F (x) = α µx G c (x) = l l 2 f c (x, l). Similarly µ min{m, c} c G c (x) = 1 2 f c (x, 1) + ( 1) 2 f c (x, 1) = λ c G(x) = α + µx + µ min{m, c} c x 0 = α/µ. Traffic intensity λ/(µc) = ρ/c tends to x 0 as c. Therefore, we may interpret x 0 as the asymptotic traffic intensity. Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
11 Estimation for queues from queue length data OU approximation for M/M/c queue We consider the OU approximation about the stable equilibrium x 0 < 1 Local drift is B = F (x 0 ) = µ Local variance is V = G(x 0 ) = 2α, that is approximately 2λ/c For large c, X c (t) has an approximate normal distribution with EX c (t) x 0 + e µt (X c (0) x 0 ), Var X c (t) λ µc 2 (1 e 2µt ) Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
12 Estimation for queues from queue length data OU approximation for M/M/c queue X c (0) N(ρ/c, ρ/c 2 ), cov (X c (s), X c (s + t)) ρ c 2 e µ t Likelihood for the Gaussian vector (X c (t 1 ),..., X c (t n )) is f (x) = { 1 (2π) n C exp 1 } 2 (x m)c 1 (x m), where m = ρ/c1 n and C = (c i,j ) n i,j=1, c i,j = ρ/c 2 exp{ µ t i+s t i } We minimize log-likelihood l(x) = n 2 log(2π) 1 2 log( C ) 1 2 (x m)c 1 (x m). For M/M/c queue, the expected values of the mean and variance of our data set are E [ 1 n E m = ρ (= cx 0 ) ] n (m c (t i ) m) 2 = ρ c (= cσ2 ) i=1 Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
13 Estimation for queues from queue length data Simulation results Large telecommunication system (M/M/300) λ = 25, µ = 0.09, c = 300 Hypothetical system with 300 servers, calls arriving at rate 50 per second, each server being able to process 0.18 calls per second and being sampled every 0.5 of second. Traffic intensity minute s worth of observations (600 data points), 5 runs The relative errors in the average of this estimates for λ and µ are small, being 2.6% and 2%. Small telecommunication system (M/M/50) λ = 4250/600, µ = 1/6, c = 50 Hypothetical system with 50 servers, calls arriving at rate 4250/300 per minute, each server being able to process 1/3 of a call per minute and being sampled every 30 seconds. Traffic intensity hours worth of observations (240 data points), 5 runs The relative errors in the average of this estimates for λ and µ are 2.9% and 1.5%. Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
14 Estimation for queues from queue length data Simulation results Shopping queue (M/M/5) λ = 0.75, µ = 0.175, c = 5 Hypothetical system with 5 servers, customers arriving at rate 3 per minute, each server being able to process 0.7 customer per minute and being sampled every 15 seconds. Traffic intensity hour worth of observations (240 data points), 5 runs The relative errors in the average of this estimates for λ is 7.9% and for µ is 35.7%. Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
15 MLE for single server queues from waiting time data The likelihood function based on the waiting times MLE for single server queues from waiting time data Notation: T n... the time between the arrivals of the n-th and n + 1-th customer W n... the remaining work of the system just before n-th customer arrive U n... the service time of the n-th customer X n = U n T n... i.i.d. random variables, independent of W n W n+1 = (W n + U n T n ) + = (W n + X n ) + Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
16 MLE for single server queues from waiting time data The likelihood function based on the waiting times {W n } is a Markov chain on the state space S = [0, ). The transition distribution function: P (W n+1 y W n = w) = P (W n+1 = 0 W n = w) + + P (0 < W n+1 y W n = w) = = P (w + X n 0) + P (0 < w + X n y) = = F X (y w)i [y 0] Jump at y = 0, with jump-size equal to F X ( w). Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
17 MLE for single server queues from waiting time data The likelihood function based on the waiting times Transition density of W n+1 given W n = w 1 α(w) y = 0 p(y w) = f X (y w) y > 0 0 y < 0 where α(w) = 1 F X ( w) and f X ( ) is the density corresponding to F X ( ) (assumed to exists). p(w n+1 W n ) = (1 α(w n )) 1 Z n+1 f X (W n+1 W n ) Z n+1, where Z n+1 = { 0 W n+1 = 0 1 W n+1 > 0. Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
18 MLE for single server queues from waiting time data The likelihood function based on the waiting times The likelihood function based on sample (W 1, W 2,..., W k ) k 1 L(θ) = p(w 1 ) p(w n+1 W n ) = n=1 k 1 = p(w 1 ) (1 α(w n ; θ)) 1 Z n+1 f X (W n+1 W n ; θ) Z n+1 n=1 where p(w 1 ) is the marginal density of W 1 and θ = (θ 1,..., θ r ) is unknown parameter vector. The likelihood score vector: d ln L dθ = d ln p(w 1; θ) dθ k 1 + n=1 (1 Z n+1 ) dα(w t; θ) (1 α(w t ; θ)) 1 + dθ k 1 n=1 Z t+1 d ln f X (W n+1 W n ; θ) dθ Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
19 MLE for single server queues from waiting time data The consistency and the asymptotic normality of the MLE ˆθ n MLE of θ is typically obtained as a solution of the likelihood equations d ln L dθ = 0. α = EX t = EU t ET t. We assume that α < 0 The first term in the likelihood score vector is asymptotically negligible Define the likelihood score function S n (θ) = n 1 t=1 U t(θ), where U t (θ) =(Z t+1 1) dα(w t; θ) (1 α(w t ; θ)) 1 + dθ d ln f X (W n+1 W n ; θ) + Z t+1 dθ {S n (θ)} is a zero-mean martingale with respect to the σ-field, σ(w n 1,..., W 1 ) The conditional information matrix ξ t (θ) = E[U t (θ)u t(θ) W t ] J n (θ) = n 1 t=1 ξ t(θ). From ergodic theorem n 1 J n (θ) a.s. J(θ) = E[ξ t (θ)] Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
20 MLE for single server queues from waiting time data The consistency and the asymptotic normality of the MLE Consistency and asymptotic normality Next two theorems holds under regularity conditions. Theorem There exists a consistent solution ˆθ n of the equations S n (θ) = 0, with probability tending to 1 as n. Moreover, ˆθ n is a local maximum of L n (θ) with probability tending to 1 as n. Theorem If ˆθ n is any consistent solution of the likelihood equations S n (θ) = 0, then 1 n(ˆθn θ) d N r (0, J 1 (θ)) 2 2(ln L n (ˆθ n ) ln L n (θ)) d χ 2 (r) Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
21 MLE for single server queues from waiting time data Application to M/M/1 λ and µ denote the arrival and service rates The stationarity condition ρ = λ/µ < 1 The density function of X t = U t T t is { λµ λ+µ f x (x) = e µx, 0 x <, λµ λ+µ eλx, < x < 0. F x (x) = P(X t x) = B(x + y)a(dy) 0 B(x + y) = P(U t x + y) = (1 e µ(x+y) )I [x + y 0] A(y) = P(T t y) = (1 e λy )I [y 0] f x (x) = F x(x) α(w) = w f x(y)dy = 1 µ λ+µ e λw To evaluate S n (θ) we need f x (x), α(x) and their partial derivatives, which are easy to compute. Katarína Starinská (starinskak@gmail.com) Arrival and Service Rates / 22
22 MLE for single server queues from waiting time data Application to M/M/1 Asmussen S. (2003): Applied Probability and Queues Basawa I.V., Bhat U.N. and Lund R.(1996):Maximum Likelihood Estimation for Single Server Queues from Waiting Time Data. Queueing Systems 24, p Ross J.V., Taimre T. and Pollet P.K.(2007):Estimation of Queues from Queue Length Data. Queueing Systems 55, p Katarína Starinská Arrival and Service Rates / 22
Constructing estimating equations for queue length data
Constructing estimating equations for queue length data R. McVinish, P.K. Pollett School of Mathematics and Physics, University of Queensland, St Lucia, Brisbane, Queensland 4072, Australia Abstract In
More informationI forgot to mention last time: in the Ito formula for two standard processes, putting
I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy
More informationMarkov Chains. X(t) is a Markov Process if, for arbitrary times t 1 < t 2 <... < t k < t k+1. If X(t) is discrete-valued. If X(t) is continuous-valued
Markov Chains X(t) is a Markov Process if, for arbitrary times t 1 < t 2
More 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 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 informationTHE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA
THE ROYAL STATISTICAL SOCIETY 4 EXAINATIONS SOLUTIONS GRADUATE DIPLOA PAPER I STATISTICAL THEORY & ETHODS The Societ provides these solutions to assist candidates preparing for the examinations in future
More information1. Aufgabenblatt zur Vorlesung Probability Theory
24.10.17 1. Aufgabenblatt zur Vorlesung By (Ω, A, P ) we always enote the unerlying probability space, unless state otherwise. 1. Let r > 0, an efine f(x) = 1 [0, [ (x) exp( r x), x R. a) Show that p f
More informationRenewal processes and Queues
Renewal processes and Queues Last updated by Serik Sagitov: May 6, 213 Abstract This Stochastic Processes course is based on the book Probabilities and Random Processes by Geoffrey Grimmett and David Stirzaker.
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 informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationInference in non-linear time series
Intro LS MLE Other Erik Lindström Centre for Mathematical Sciences Lund University LU/LTH & DTU Intro LS MLE Other General Properties Popular estimatiors Overview Introduction General Properties Estimators
More informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY (formerly the Examinations of the Institute of Statisticians) GRADUATE DIPLOMA, 2004
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY (formerly the Examinations of the Institute of Statisticians) GRADUATE DIPLOMA, 004 Statistical Theory and Methods I Time Allowed: Three Hours Candidates should
More informationContinuous Time Markov Chains
Continuous Time Markov Chains Stochastic Processes - Lecture Notes Fatih Cavdur to accompany Introduction to Probability Models by Sheldon M. Ross Fall 2015 Outline Introduction Continuous-Time Markov
More informationMultivariate Simulations
Multivariate Simulations Katarína Starinská starinskak@gmail.com Charles University Faculty of Mathematics and Physics Prague, Czech Republic 21.11.2011 Katarína Starinská Multivariate Simulations 1 /
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 informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationLesson 4: Stationary stochastic processes
Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@univaq.it Stationary stochastic processes Stationarity is a rather intuitive concept, it means
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 informationn! (k 1)!(n k)! = F (X) U(0, 1). (x, y) = n(n 1) ( F (y) F (x) ) n 2
Order statistics Ex. 4.1 (*. Let independent variables X 1,..., X n have U(0, 1 distribution. Show that for every x (0, 1, we have P ( X (1 < x 1 and P ( X (n > x 1 as n. Ex. 4.2 (**. By using induction
More informationFormulas for probability theory and linear models SF2941
Formulas for probability theory and linear models SF2941 These pages + Appendix 2 of Gut) are permitted as assistance at the exam. 11 maj 2008 Selected formulae of probability Bivariate probability Transforms
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 information1 Continuous-time chains, finite state space
Université Paris Diderot 208 Markov chains Exercises 3 Continuous-time chains, finite state space Exercise Consider a continuous-time taking values in {, 2, 3}, with generator 2 2. 2 2 0. Draw the diagramm
More information(Practice Version) Midterm Exam 2
EECS 126 Probability and Random Processes University of California, Berkeley: Fall 2014 Kannan Ramchandran November 7, 2014 (Practice Version) Midterm Exam 2 Last name First name SID Rules. DO NOT open
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 informationMAT SYS 5120 (Winter 2012) Assignment 5 (not to be submitted) There are 4 questions.
MAT 4371 - SYS 5120 (Winter 2012) Assignment 5 (not to be submitted) There are 4 questions. Question 1: Consider the following generator for a continuous time Markov chain. 4 1 3 Q = 2 5 3 5 2 7 (a) Give
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationHypothesis testing: theory and methods
Statistical Methods Warsaw School of Economics November 3, 2017 Statistical hypothesis is the name of any conjecture about unknown parameters of a population distribution. The hypothesis should be verifiable
More informationTwelfth Problem Assignment
EECS 401 Not Graded PROBLEM 1 Let X 1, X 2,... be a sequence of independent random variables that are uniformly distributed between 0 and 1. Consider a sequence defined by (a) Y n = max(x 1, X 2,..., X
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 informationQuick Tour of Basic Probability Theory and Linear Algebra
Quick Tour of and Linear Algebra Quick Tour of and Linear Algebra CS224w: Social and Information Network Analysis Fall 2011 Quick Tour of and Linear Algebra Quick Tour of and Linear Algebra Outline Definitions
More informationVariance reduction techniques
Variance reduction techniques Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/ moltchan/modsim/ http://www.cs.tut.fi/kurssit/tlt-2706/ OUTLINE: Simulation with a given confidence;
More informationBrief Review on Estimation Theory
Brief Review on Estimation Theory K. Abed-Meraim ENST PARIS, Signal and Image Processing Dept. abed@tsi.enst.fr This presentation is essentially based on the course BASTA by E. Moulines Brief review on
More informationJesper Møller ) and Kateřina Helisová )
Jesper Møller ) and ) ) Aalborg University (Denmark) ) Czech Technical University/Charles University in Prague 5 th May 2008 Outline 1. Describing model 2. Simulation 3. Power tessellation of a union of
More informationACM 116: Lectures 3 4
1 ACM 116: Lectures 3 4 Joint distributions The multivariate normal distribution Conditional distributions Independent random variables Conditional distributions and Monte Carlo: Rejection sampling Variance
More informationThis exam is closed book and closed notes. (You will have access to a copy of the Table of Common Distributions given in the back of the text.
TEST #3 STA 536 December, 00 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. You will have access to a copy
More informationPoisson Approximation for Two Scan Statistics with Rates of Convergence
Poisson Approximation for Two Scan Statistics with Rates of Convergence Xiao Fang (Joint work with David Siegmund) National University of Singapore May 28, 2015 Outline The first scan statistic The second
More informationMaster s Written Examination
Master s Written Examination Option: Statistics and Probability Spring 016 Full points may be obtained for correct answers to eight questions. Each numbered question which may have several parts is worth
More informationDeterministic. Deterministic data are those can be described by an explicit mathematical relationship
Random data Deterministic Deterministic data are those can be described by an explicit mathematical relationship Deterministic x(t) =X cos r! k m t Non deterministic There is no way to predict an exact
More informationNon Markovian Queues (contd.)
MODULE 7: RENEWAL PROCESSES 29 Lecture 5 Non Markovian Queues (contd) For the case where the service time is constant, V ar(b) = 0, then the P-K formula for M/D/ queue reduces to L s = ρ + ρ 2 2( ρ) where
More informationExercises with solutions (Set D)
Exercises with solutions Set D. A fair die is rolled at the same time as a fair coin is tossed. Let A be the number on the upper surface of the die and let B describe the outcome of the coin toss, where
More informationn! (k 1)!(n k)! = F (X) U(0, 1). (x, y) = n(n 1) ( F (y) F (x) ) n 2
Order statistics Ex. 4. (*. Let independent variables X,..., X n have U(0, distribution. Show that for every x (0,, we have P ( X ( < x and P ( X (n > x as n. Ex. 4.2 (**. By using induction or otherwise,
More informationP (L d k = n). P (L(t) = n),
4 M/G/1 queue In the M/G/1 queue customers arrive according to a Poisson process with rate λ and they are treated in order of arrival The service times are independent and identically distributed 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 informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More 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 informationTHE QUEEN S UNIVERSITY OF BELFAST
THE QUEEN S UNIVERSITY OF BELFAST 0SOR20 Level 2 Examination Statistics and Operational Research 20 Probability and Distribution Theory Wednesday 4 August 2002 2.30 pm 5.30 pm Examiners { Professor R M
More informationStochastic Processes
Stochastic Processes 8.445 MIT, fall 20 Mid Term Exam Solutions October 27, 20 Your Name: Alberto De Sole Exercise Max Grade Grade 5 5 2 5 5 3 5 5 4 5 5 5 5 5 6 5 5 Total 30 30 Problem :. True / False
More informationConsistency of the maximum likelihood estimator for general hidden Markov models
Consistency of the maximum likelihood estimator for general hidden Markov models Jimmy Olsson Centre for Mathematical Sciences Lund University Nordstat 2012 Umeå, Sweden Collaborators Hidden Markov models
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 and Jiekun Feng Cornell University Workshop on Congestion Games Institute
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationSystem Identification, Lecture 4
System Identification, Lecture 4 Kristiaan Pelckmans (IT/UU, 2338) Course code: 1RT880, Report code: 61800 - Spring 2016 F, FRI Uppsala University, Information Technology 13 April 2016 SI-2016 K. Pelckmans
More informationFinal Examination Statistics 200C. T. Ferguson June 11, 2009
Final Examination Statistics 00C T. Ferguson June, 009. (a) Define: X n converges in probability to X. (b) Define: X m converges in quadratic mean to X. (c) Show that if X n converges in quadratic mean
More informationTesting Hypothesis. Maura Mezzetti. Department of Economics and Finance Università Tor Vergata
Maura Department of Economics and Finance Università Tor Vergata Hypothesis Testing Outline It is a mistake to confound strangeness with mystery Sherlock Holmes A Study in Scarlet Outline 1 The Power Function
More informationPractical conditions on Markov chains for weak convergence of tail empirical processes
Practical conditions on Markov chains for weak convergence of tail empirical processes Olivier Wintenberger University of Copenhagen and Paris VI Joint work with Rafa l Kulik and Philippe Soulier Toronto,
More informationMaster s Written Examination
Master s Written Examination Option: Statistics and Probability Fall 2013 Full points may be obtained for correct answers to eight questions. Each numbered question (which may have several parts) is worth
More informationMarkovian N-Server Queues (Birth & Death Models)
Markovian -Server Queues (Birth & Death Moels) - Busy Perio Arrivals Poisson (λ) ; Services exp(µ) (E(S) = /µ) Servers statistically ientical, serving FCFS. Offere loa R = λ E(S) = λ/µ Erlangs Q(t) = number
More informationStatistical Problems Related to Excitation Threshold and Reset Value of Membrane Potentials
Statistical Problems Related to Excitation Threshold and Reset Value of Membrane Potentials jahn@mathematik.uni-mainz.de Le Mans - March 18, 2009 Contents 1 Biological Background and the Problem Definition
More informationEECS 126 Probability and Random Processes University of California, Berkeley: Fall 2014 Kannan Ramchandran November 13, 2014.
EECS 126 Probability and Random Processes University of California, Berkeley: Fall 2014 Kannan Ramchandran November 13, 2014 Midterm Exam 2 Last name First name SID Rules. DO NOT open the exam until instructed
More information2. Variance and Covariance: We will now derive some classic properties of variance and covariance. Assume real-valued random variables X and Y.
CS450 Final Review Problems Fall 08 Solutions or worked answers provided Problems -6 are based on the midterm review Identical problems are marked recap] Please consult previous recitations and textbook
More informationMAS223 Statistical Inference and Modelling Exercises
MAS223 Statistical Inference and Modelling Exercises The exercises are grouped into sections, corresponding to chapters of the lecture notes Within each section exercises are divided into warm-up questions,
More informationLecture 25: Review. Statistics 104. April 23, Colin Rundel
Lecture 25: Review Statistics 104 Colin Rundel April 23, 2012 Joint CDF F (x, y) = P [X x, Y y] = P [(X, Y ) lies south-west of the point (x, y)] Y (x,y) X Statistics 104 (Colin Rundel) Lecture 25 April
More informationLAN property for ergodic jump-diffusion processes with discrete observations
LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &
More informationT. Liggett Mathematics 171 Final Exam June 8, 2011
T. Liggett Mathematics 171 Final Exam June 8, 2011 1. The continuous time renewal chain X t has state space S = {0, 1, 2,...} and transition rates (i.e., Q matrix) given by q(n, n 1) = δ n and q(0, n)
More informationChapter 3 sections. SKIP: 3.10 Markov Chains. SKIP: pages Chapter 3 - continued
Chapter 3 sections Chapter 3 - continued 3.1 Random Variables and Discrete Distributions 3.2 Continuous Distributions 3.3 The Cumulative Distribution Function 3.4 Bivariate Distributions 3.5 Marginal Distributions
More informationMATH c UNIVERSITY OF LEEDS Examination for the Module MATH2715 (January 2015) STATISTICAL METHODS. Time allowed: 2 hours
MATH2750 This question paper consists of 8 printed pages, each of which is identified by the reference MATH275. All calculators must carry an approval sticker issued by the School of Mathematics. c UNIVERSITY
More informationHomework for 1/13 Due 1/22
Name: ID: Homework for 1/13 Due 1/22 1. [ 5-23] An irregularly shaped object of unknown area A is located in the unit square 0 x 1, 0 y 1. Consider a random point distributed uniformly over the square;
More informationQueueing Review. Christos Alexopoulos and Dave Goldsman 10/6/16. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA
1 / 24 Queueing Review (mostly from BCNN) Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/6/16 2 / 24 Outline 1 Introduction 2 Queueing Notation 3 Transient
More informationf(x θ)dx with respect to θ. Assuming certain smoothness conditions concern differentiating under the integral the integral sign, we first obtain
0.1. INTRODUCTION 1 0.1 Introduction R. A. Fisher, a pioneer in the development of mathematical statistics, introduced a measure of the amount of information contained in an observaton from f(x θ). Fisher
More informationIntroduction to queuing theory
Introduction to queuing theory Claude Rigault ENST claude.rigault@enst.fr Introduction to Queuing theory 1 Outline The problem The number of clients in a system The client process Delay processes Loss
More information7 Variance Reduction Techniques
7 Variance Reduction Techniques In a simulation study, we are interested in one or more performance measures for some stochastic model. For example, we want to determine the long-run average waiting time,
More informationStatistics & Data Sciences: First Year Prelim Exam May 2018
Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book
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 concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationProbability Models in Electrical and Computer Engineering Mathematical models as tools in analysis and design Deterministic models Probability models
Probability Models in Electrical and Computer Engineering Mathematical models as tools in analysis and design Deterministic models Probability models Statistical regularity Properties of relative frequency
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationRecap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks
Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution
More 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 informationSTAT FINAL EXAM
STAT101 2013 FINAL EXAM This exam is 2 hours long. It is closed book but you can use an A-4 size cheat sheet. There are 10 questions. Questions are not of equal weight. You may need a calculator for some
More informationInterest Rate Models:
1/17 Interest Rate Models: from Parametric Statistics to Infinite Dimensional Stochastic Analysis René Carmona Bendheim Center for Finance ORFE & PACM, Princeton University email: rcarmna@princeton.edu
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 informationSpring 2012 Math 541B Exam 1
Spring 2012 Math 541B Exam 1 1. A sample of size n is drawn without replacement from an urn containing N balls, m of which are red and N m are black; the balls are otherwise indistinguishable. Let X denote
More informationEconometría 2: Análisis de series de Tiempo
Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 II. Basic definitions A time series is a set of observations X t, each
More informationSimulation of stationary processes. Timo Tiihonen
Simulation of stationary processes Timo Tiihonen 2014 Tactical aspects of simulation Simulation has always a goal. How to organize simulation to reach the goal sufficiently well and without unneccessary
More informationIntroduction to Estimation Methods for Time Series models Lecture 2
Introduction to Estimation Methods for Time Series models Lecture 2 Fulvio Corsi SNS Pisa Fulvio Corsi Introduction to Estimation () Methods for Time Series models Lecture 2 SNS Pisa 1 / 21 Estimators:
More informationHEAVY-TRAFFIC EXTREME-VALUE LIMITS FOR QUEUES
HEAVY-TRAFFIC EXTREME-VALUE LIMITS FOR QUEUES by Peter W. Glynn Department of Operations Research Stanford University Stanford, CA 94305-4022 and Ward Whitt AT&T Bell Laboratories Murray Hill, NJ 07974-0636
More informationBIRTH DEATH PROCESSES AND QUEUEING SYSTEMS
BIRTH DEATH PROCESSES AND QUEUEING SYSTEMS Andrea Bobbio Anno Accademico 999-2000 Queueing Systems 2 Notation for Queueing Systems /λ mean time between arrivals S = /µ ρ = λ/µ N mean service time traffic
More informationISyE 6650 Test 2 Solutions
1 NAME ISyE 665 Test 2 Solutions Summer 2 This test is open notes, open books. The test is two hours long. 1. Consider an M/M/3/4 queueing system in steady state with arrival rate λ = 3 and individual
More informationChapter 7. Basic Probability Theory
Chapter 7. Basic Probability Theory I-Liang Chern October 20, 2016 1 / 49 What s kind of matrices satisfying RIP Random matrices with iid Gaussian entries iid Bernoulli entries (+/ 1) iid subgaussian entries
More informationLinear Regression. Junhui Qian. October 27, 2014
Linear Regression Junhui Qian October 27, 2014 Outline The Model Estimation Ordinary Least Square Method of Moments Maximum Likelihood Estimation Properties of OLS Estimator Unbiasedness Consistency Efficiency
More informationStatistics Ph.D. Qualifying Exam: Part I October 18, 2003
Statistics Ph.D. Qualifying Exam: Part I October 18, 2003 Student Name: 1. Answer 8 out of 12 problems. Mark the problems you selected in the following table. 1 2 3 4 5 6 7 8 9 10 11 12 2. Write your answer
More informationGaussian Processes. 1. Basic Notions
Gaussian Processes 1. Basic Notions Let T be a set, and X : {X } T a stochastic process, defined on a suitable probability space (Ω P), that is indexed by T. Definition 1.1. We say that X is a Gaussian
More informationSystem Identification, Lecture 4
System Identification, Lecture 4 Kristiaan Pelckmans (IT/UU, 2338) Course code: 1RT880, Report code: 61800 - Spring 2012 F, FRI Uppsala University, Information Technology 30 Januari 2012 SI-2012 K. Pelckmans
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 informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2008 Prof. Gesine Reinert 1 Data x = x 1, x 2,..., x n, realisations of random variables X 1, X 2,..., X n with distribution (model)
More informationLink lecture - Lagrange Multipliers
Link lecture - Lagrange Multipliers Lagrange multipliers provide a method for finding a stationary point of a function, say f(x, y) when the variables are subject to constraints, say of the form g(x, y)
More informationMaster s Written Examination - Solution
Master s Written Examination - Solution Spring 204 Problem Stat 40 Suppose X and X 2 have the joint pdf f X,X 2 (x, x 2 ) = 2e (x +x 2 ), 0 < x < x 2
More informationM.Sc.(Mathematics with Applications in Computer Science) Probability and Statistics
MMT-008 Assignment Booklet M.Sc.(Mathematics with Applications in Computer Science) Probability and Statistics (Valid from 1 st July, 013 to 31 st May, 014) It is compulsory to submit the assignment before
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 informationECE353: Probability and Random Processes. Lecture 18 - Stochastic Processes
ECE353: Probability and Random Processes Lecture 18 - Stochastic Processes Xiao Fu School of Electrical Engineering and Computer Science Oregon State University E-mail: xiao.fu@oregonstate.edu From RV
More information. Find E(V ) and var(v ).
Math 6382/6383: Probability Models and Mathematical Statistics Sample Preliminary Exam Questions 1. A person tosses a fair coin until she obtains 2 heads in a row. She then tosses a fair die the same number
More information