Chapter 11 Estimation of Absolute Performance
|
|
- Prudence Cannon
- 5 years ago
- Views:
Transcription
1 Chapter 11 Estimation of Absolute Performance
2 Purpose Objective: Estimate system performance via simulation If q is the system performance, the precision of the estimator qˆ can be measured by: The standard error of qˆ. The width of a confidence interval (CI) for q. Purpose of statistical analysis: To estimate the standard error or CI. To figure out the number of observations required to achieve desired error/ci. Potential issues to overcome: Autocorrelation, e.g. inventory cost for subsequent weeks lack statistical independence. Initial conditions, e.g. inventory on hand and # of backorders at time 0 would most likely influence the performance of week 1. 2
3 Outline Distinguish the two types of simulation: transient vs. steady state. Illustrate the inherent variability in a stochastic discreteevent simulation. Cover the statistical estimation of performance measures. Discusses the analysis of transient simulations. Discusses the analysis of steady-state simulations. 3
4 Introduction to simulation From System to Model: Simulation based problem solving DISCRETE EVENT SIMULATION WITH APPLICATION TO COMPUTER COMMUNICATION SYSTEMS PERFORMANCE Introduction to Simulation Helena Szczerbicka, Kishor S. Trivedi and Pawan K. Choudhary
5 Type of Simulations Terminating verses non-terminating simulations Terminating simulation: Runs for some duration of time T E, where E is a specified event that stops the simulation. Starts at time 0 under well-specified initial conditions. Ends at the stopping time T E. T E might be random/unknown or predetermined Bank example: Opens at 8:30 am (time 0) with no customers present and 8 of the 11 teller working (initial conditions), and closes at 4:30 pm (Time T E = 480 minutes). The simulation analyst chooses to consider it a terminating system because the object of interest is one day s operation. 5
6 Type of Simulations Non-terminating simulation: Runs continuously, or at least over a very long period of time. Examples: Earth Simulator, assembly lines that shut down infrequently, telephone systems, networks, hospital emergency rooms. Initial conditions defined by the analyst. Runs for some analyst-specified period of time T E. Study the steady-state (long-run) properties of the system, properties that are not influenced by the initial conditions of the model. Whether a simulation is considered to be terminating or non-terminating depends on both The objectives of the simulation study and The nature of the system. 6
7 Stochastic Nature of Output Data Model output consists of one or more random variables because the model is an input-output transformation and the input variables are RVs M/M/1 queueing example: Poisson arrival rate = 5.0 per second (r i = 0.2 second); service time ~ Exp(6.67 or 0.15 second) System performance: utilization, time-avg queue length, L Q,avg. response time, & avg queueing delay Using OMNeT++, we run 5 statistically independent replications Each replication/run uses a distinct stream of random numbers Facilitated by OMNET++ ``repeat = 5 statement 7
8 Stochastic Nature of Output Data M/M/1 queueing example (cont.): Run/Replication (r) Uilization L_Q w d Average
9 Stochastic Nature of Output Data M/M/1 queueing example (cont.): U L_Q w d
10 Measures of performance Consider the estimation of a performance parameter, q (or f), of a simulated system in a single replication Discrete time data: [Y 1, Y 2,, Y n ], with ordinary mean: q Continuous-time data: {Y(t), 0 t T E } with time-weighted mean: f Point estimation for discrete time data. The point estimator: qˆ = 1 n n i= 1 Y i Is unbiased if its expected value is q, that is if: Is biased if: E( ˆ) q q E( ˆ) q =q Desired 10
11 Point Estimator [Performance Measures] Point estimation for continuous-time data. The point estimator: Is biased in general where: E( ˆ) f f. An unbiased or low-bias estimator is desired. Usually, system performance measures can be put into the common framework of q or f: e.g., the proportion of customers whose queueing delays equal, less, or more than 2 seconds % of time queue > 7 fˆ = ìï Y i = í îï 1 T E Y( t) dt T 0 E 1, if d i = {>,<} 2 0, otherwise ì ï 1, if L Q (t) > 7 Y(t) = í îï 0, otherwise 11
12 Point Estimator [Performance Measures] Performance measure that does not fit: quantile or percentile: P[Y q]= 0.85 E.g., the queueing delay for which 85% of customers will experience it or less than it Estimating quantiles: the inverse of the problem of estimating a proportion or probability. Consider a histogram of the observed values Y: Find qˆ such that 100p% of the histogram is to the left of (smaller than) qˆ. 12
13 Confidence-Interval Estimation [Performance Measures] To understand confidence intervals fully, it is important to distinguish between measures of error, and measures of risk e.g., confidence interval versus prediction interval. Suppose the model is the normal distribution with mean q, variance s 2 (both unknown). Let Y i be the average cycle time for parts produced on the i th replication of the simulation (its mathematical expectation is q). Average cycle time will vary from day to day, but over the long-run the average of the averages will be close to q. Sample variance across R replications: S 2 = 1 R 1 R i= 1 ( Y i Y )... 2 Y.. = 1 R R å i=1 Y i. 13
14 Confidence-Interval Estimation [Performance Measures] Confidence Interval (CI): Y.. is not q; Y.. is an estimate of q based on the sample; & it has error A CI is measure of error between q and the estimator Y.. Where Y i. are normally distributed. S Y.. t / 2, R 1 R We cannot know for certain how far Y.. is from q but CI attempts to bound that error. A CI may be wrong So, a CI, such as 95%, tells us how much we can trust the interval to actually bound the error between Y.. and q. The more replications we make, the less error there is in (converging to 0 as R goes to infinity). Y.. 14
15 Confidence-Interval Estimation [Performance Measures] Prediction Interval (PI): A measure of risk. A good guess for the average cycle time on a particular day is our estimator (Y.. ) but it is unlikely to be exactly right. PI is designed to be wide enough to contain the actual average cycle time on any particular day with high probability. Normal-theory prediction interval: Y.. t / 2, R 1S 1 1 R The length of PI will not go to 0 as R increases because we can never simulate away risk. PI s limit is: q z / 2s 15
16 CI & PI: An Example [Performance Measures] M/M/1 System Run/Replication (r) Uilization L_Q w d Average Standard Deviation t * S /SQRT(R) CI lower CI upper t * S * SQRT(1 + 1/R) PI lower PI upper Theoretical
17 CI & PI: An Example [Performance Measures] M/M/1 System Replications (R) = 5 A 95% CI.t 0.05/2, 5-1 = t 0.025, 4 = 2.78 A CI for Utilization = ± 2.78(0.1242/sqrt(5)) = ± = [0.5505, ] Our estimate of the long-run utilization is , but there could be an error as much as ± in this estimate A PI for Utilization = ± 2.78(0.1242/sqrt(1+1/5)) = ± = [0.3897, ] In a particular run, we are 95% confident that the utilization will be = ±
18 t-distribution Table [Performance Measures] 18
19 CI & PI: An Example [Performance Measures] M/M/1 System: CIs Avg Lo Hi U 1 L_Q 2 w 3 4 d 19
20 Output Analysis for Terminating Simulations A terminating simulation: runs over a simulated time interval [0, T E ]. A common goal is to estimate: 1 q = E n 1 f = E TE n i= 1 0 Better, we use independent replications/runs To eliminate the bias introduced by a particular random number stream Each run using a different random number stream Independently chosen initial conditions Y T E i, Y ( t) dt, for discrete output for continuous output Y ( t),0 t T E 20
21 Statistical Background [Terminating Simulations] Across Replication: For example: the response time averages (discrete time data) R The average: 1 Y.. = Y i. R i= 1 R The sample variance: S = ( Y i. Y.. ) R 1 i= 1 The confidence-interval half-width: S H = t / 2, R 1 R Within replication: For example: the L Q (a continuous time data) The average: 1 T Ei Yi. = Yi ( t) dt T 0 The sample variance: Ei S 2 i = 1 T Ei T Ei 0 Y ( t) Y i i. 2 dt 21
22 Statistical Background [Terminating Simulations] Important to distinguish within-replication data from across-replication data. For example, simulation of a computer system Two performance measures of that system: response time and length of the queue at the CPU. Let Y ij be the response time for the j th process executed in the i th replication. Across-replication data are formed by summarizing withinreplication data. Y i. 22
23 Statistical Background [Terminating Simulations] Within-Replication Data Across-Replication Data Y 11 Y 12 Y 1n1 Y 1. S 1 2 H 1 Y 21 Y 22 Y 2n2 Y 2. S 2 2 H 2 Y R1 Y R2 Y RnR Y R. S 2 R H R Y.. S 2 H Within-Replication Data Across-Replication Data Y 1 (t), 0 < t < T E1 Y 1. S 1 2 H 1 Y 2 (t), 0 < t < T E2 Y 2. S 2 2 H 2 Y R (t), 0 < t < T ER Y R. S R 2 H R Y.. S 2 H 23
24 Statistical Background [Terminating Simulations] Overall sample average, Y.., and the individual replication sample averages,, are always unbiased estimators. Y i. Across-replication data are Independent (different random numbers), and Identically distributed (same model). But, within-replication data do not have these properties. 24
25 Output Analysis for Terminating Simulations Example: the M/M/1 model previously discussed See Slides # 7 9, And, the example(s) in the textbook 25
26 C.I. with Specified Precision [Terminating Simulations] The half-length H of a 100(1 )% confidence interval for a mean q, based on the t distribution, is given by: H = t / 2, R 1 R is the # of replications Suppose that an error criterion e is specified with probability 1 -, a sufficiently large sample size should satisfy: P What is this sufficiently large R? Problem: Both S & H are not-so-straightforward fns of R S R q < e 1 Y.. S 2 is the sample variance 26
27 C.I. with Specified Precision [Terminating Simulations] Make an initial sample of size R 0 independent replications. Based on R 0, obtain an initial estimate S 0 2 of the population variance s 2. Then, choose a sample size R such that R R 0 : Since t /2, R-1 z /2, an initial estimate of R: 2 z / 2S0 R, z / 2 is the standard normal distribution. e t / 2, 1 0 R is the smallest integer satisfying R R 0 and R S R e Collect/conduct R - R 0 additional observations/replications. Form the 100(1- )% C.I. for q based on R: S Y.. t / 2, R 1 R Since S could differ than S 0, H would differ (<, >) than e 2 27
28 C.I. with Specified Precision [Terminating Simulations] Density of the t-distribution for 1, 2, 3, 5, 10 (all in green), and 30 (in red) degrees of freedom vs. standard normal distribution (blue) 28
29 C.I. with Specified Precision: An Example [Terminating Simulations] M/M/1 Example Objective: estimate utilization r to within ±0.04 with probability 0.95 Initial sample of size R 0 = 5 is taken and an initial estimate of the population variance is S 0 = The error criterion is e = 0.04 and confidence coefficient is 1- = 0.95, hence, the final sample size must be at least: æ z S 0 ö 2 æ ö 2 ç = ç = è e ø è 0.04 ø Next, we test the possible candidates (R = 37, 38, 39, 40, ) with the inequality of t /2, R-1 for the final sample size: R t 0.025, R R = 40 is the smallest integer satisfying the error criterion So, R - R 0 = 35 additional replications are needed. After obtaining additional outputs, H should be checked. 29
30 z-distribution Table [Performance Measures] 30
31 Output Analysis for Steady-State Simulation Consider a single run of a simulation model to estimate a steady-state or long-run characteristics of the system. The single run produces observations Y 1, Y 2,... (generally the samples of an autocorrelated time series). Steady-state/long-run Performance measure: q = f = lim 1 n n 1 T lim 0 T E E n å Y i, i=1 T E Y( t) dt, for discrete measure for continuous measure (with probability 1) (with probability 1) Independent of the initial conditions. 31
32 Output Analysis for Steady-State Simulation The sample size is a design choice, with several considerations in mind: Any bias in the point estimator that is due to artificial or arbitrary initial conditions (bias can be severe if run length is too short). Desired precision of the point estimator. Budget constraints on computer resources. Notation: the estimation of q from a discrete-time output process. One replication (or run), the output data: Y 1, Y 2, Y 3, With several replications, the output data for replication r: Y r1, Y r2, Y r3, 32
33 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d) 20 Experiments (i = 1,, 20); each begins with empty & idle state Simulation runs until 2 i jobs departs the system (each 1 replication) Experiment # No. of Jobs Uilization L_Q w d
34 Initialization Bias [Steady-State Simulations] Utilization
35 Initialization Bias [Steady-State Simulations] 4.50 L Q
36 Initialization Bias [Steady-State Simulations] w
37 Initialization Bias [Steady-State Simulations] d
38 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d) Five initial states: X(t = 0) = 0, 1, 2, 3, 4, 5, 10, 20 X(t): # of jobs in system at t 20 Experiments (i = 1,, 20) for each initial state Results only for L Q Green colored cells correspond for values within [Theoretical value 5%, Theoretical value + 5%] 38
39 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d) X(0) = 0 X(0) = 1 X(0) = 2 X(0) = 3 X(0) = 4 X(0) = 5 X(0) = 10 X(0) =
40 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d)
41 Initialization Bias [Steady-State Simulations] Methods to reduce the point-estimator bias caused by using artificial and unrealistic initial conditions: Intelligent initialization. Divide simulation into an initialization phase and data-collection phase. Intelligent initialization Initialize the simulation in a state that is more representative of long-run conditions. If the system exists, collect data on it and use these data to specify more nearly typical initial conditions, e.g., measure/sample a router at some time. If the system can be simplified enough to be mathematically solvable, solve the simplified model to find long-run expected conditions. E.g. queueing models 41
42 Initialization Bias [Steady-State Simulations] Divide each simulation into two phases: An initialization phase, from time 0 to time T 0. A data-collection phase, from T 0 to the stopping time T 0 +T E. The choice of T 0 is important: After T 0, system should be more nearly representative of steady-state behavior. System has reached steady state: the probability distribution of the system state is close to the steady-state probability distribution (bias of response variable is negligible). In OMNeT++ terminology, warmup-period E.g., warmup-period = 20s 42
43 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d) Three initial states: X(t = 0) = 0, 3 & 10 X(t): # of jobs in system at t 5 Experiments for each initial state # of replications in each experiment: 5 & 20 Results only for d 95% confidence interval: Y.. t / 2, R 1 S R 43
44 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d) Initial states: X(t = 0) = 0 # of replications in each experiment: 5 X(0) = 0 R Average Standard Deviation t * S /SQRT(R) CI lower CI upper Theoretical
45 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d) Initial states: X(t = 0) = 0 & 3 # of replications in each experiment:
46 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d) Initial states: X(t = 0) = 10 & 3 # of replications in each experiment:
47 Initialization Bias [Steady-State Simulations] M/M/1 Example (cont d) Initial states: X(t = 0) = 3 # of replications in each experiment: 5 &
48 Replication Method [Steady-State Simulations] Use to estimate point-estimator variability and to construct a confidence interval. Approach: make R replications, initializing and deleting from each one the same way. Important to do a thorough job of investigating the initial-condition bias: Bias is not affected by the number of replications, instead, it is affected only by deleting more data (i.e., increasing T 0 ) or extending the length of each run (i.e. increasing T E ). Basic raw output data {Y rj, r = 1,..., R; j = 1,, n} is derived by: Individual observation from within replication r. 48
49 Replication Method [Steady-State Simulations] Each replication is regarded as a single sample for n 1 estimating q. For replication r: Y ( n, d) = The overall point estimator: Y R.. n, d) = Yr.( n, d) and E[ Y.. ( n, d)] R r= 1 If d and n are chosen sufficiently large: q n,d ~ q. Y ( n, ) is an approximately unbiased estimator of q... d To estimate standard error of Y.., the sample variance and standard error: S 2 Y r. n d rj j= d 1 1 ( = q R R = ( Yr. Y.. ) = Yr. RY.. and s. e.( Y.. ) = R 1 r= 1 R 1 r= 1 n, d S R 49
50 Replication Method [Steady-State Simulations] Length of each replication (n) beyond deletion point (d): (n - d) > 10d Number of replications (R) should be as many as time permits, up to about 25 replications. For a fixed total sample size (n), as fewer data are deleted ( d): C.I. shifts: greater bias. Standard error of Y ( n,.. d ) decreases: decrease variance. Reducing bias Trade off Increasing variance 50
51 Summary Stochastic discrete-event simulation is a statistical experiment. Purpose of statistical experiment: obtain estimates of the performance measures of the system. Purpose of statistical analysis: acquire some assurance that these estimates are sufficiently precise. Distinguish: terminating simulations and steady-state simulations. Steady-state output data are more difficult to analyze Decisions: initial conditions and run length Possible solutions to bias: deletion of data and increasing run length Statistical precision of point estimators are estimated by standarderror or confidence interval Method of independent replications was emphasized. 51
Chapter 11 Output Analysis for a Single Model. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation
Chapter 11 Output Analysis for a Single Model Banks, Carson, Nelson & Nicol Discrete-Event System Simulation Purpose Objective: Estimate system performance via simulation If q is the system performance,
More informationSlides 12: Output Analysis for a Single Model
Slides 12: Output Analysis for a Single Model Objective: Estimate system performance via simulation. If θ is the system performance, the precision of the estimator ˆθ can be measured by: The standard error
More informationComputer Science, Informatik 4 Communication and Distributed Systems. Simulation. Discrete-Event System Simulation. Dr.
Simulation Discrete-Event System Simulation Chapter 0 Output Analysis for a Single Model Purpose Objective: Estimate system performance via simulation If θ is the system performance, the precision of the
More informationChapter 11. Output Analysis for a Single Model Prof. Dr. Mesut Güneş Ch. 11 Output Analysis for a Single Model
Chapter Output Analysis for a Single Model. Contents Types of Simulation Stochastic Nature of Output Data Measures of Performance Output Analysis for Terminating Simulations Output Analysis for Steady-state
More informationOutput Analysis for a Single Model
Output Analysis for a Single Model Output Analysis for a Single Model Output analysis is the examination of data generated by a simulation. Its purpose is to predict the performance of a system or to compare
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 Quote of the Day A person with one watch knows what time it is. A person with two
More informationEE/PEP 345. Modeling and Simulation. Spring Class 11
EE/PEP 345 Modeling and Simulation Class 11 11-1 Output Analysis for a Single Model Performance measures System being simulated Output Output analysis Stochastic character Types of simulations Output analysis
More informationOutput Data Analysis for a Single System
Output Data Analysis for a Single System Chapter 9 Based on the slides provided with the textbook 2 9.1 Introduction Output data analysis is often not conducted appropriately Treating output of a single
More informationOverall Plan of Simulation and Modeling I. Chapters
Overall Plan of Simulation and Modeling I Chapters Introduction to Simulation Discrete Simulation Analytical Modeling Modeling Paradigms Input Modeling Random Number Generation Output Analysis Continuous
More informationB. Maddah INDE 504 Discrete-Event Simulation. Output Analysis (1)
B. Maddah INDE 504 Discrete-Event Simulation Output Analysis (1) Introduction The basic, most serious disadvantage of simulation is that we don t get exact answers. Two different runs of the same model
More informationComputer Science, Informatik 4 Communication and Distributed Systems. Simulation. Discrete-Event System Simulation. Dr.
Simulation Discrete-Event System Simulation Chapter Comparison and Evaluation of Alternative System Designs Purpose Purpose: comparison of alternative system designs. Approach: discuss a few of many statistical
More informationSimulation. Where real stuff starts
1 Simulation Where real stuff starts ToC 1. What is a simulation? 2. Accuracy of output 3. Random Number Generators 4. How to sample 5. Monte Carlo 6. Bootstrap 2 1. What is a simulation? 3 What is a simulation?
More informationConstruction Operation Simulation
Construction Operation Simulation Lecture #8 Output analysis Amin Alvanchi, PhD Construction Engineering and Management Department of Civil Engineering, Sharif University of Technology Outline 2 Introduction
More informationNetwork Simulation Chapter 6: Output Data Analysis
Network Simulation Chapter 6: Output Data Analysis Prof. Dr. Jürgen Jasperneite 1 Contents Introduction Types of simulation output Transient detection When to terminate a simulation 2 Prof. Dr. J ürgen
More informationRedacted for Privacy
AN ABSTRACT OF THE THESIS OF Lori K. Baxter for the degree of Master of Science in Industrial and Manufacturing Engineering presented on June 4, 1990. Title: Truncation Rules in Simulation Analysis: Effect
More informationSimulation. Where real stuff starts
Simulation Where real stuff starts March 2019 1 ToC 1. What is a simulation? 2. Accuracy of output 3. Random Number Generators 4. How to sample 5. Monte Carlo 6. Bootstrap 2 1. What is a simulation? 3
More informationOutput Data Analysis for a Single System
CHAPTER 9 Output Data Analysis for a Single System 9.1 Introduction... 9. Transient and Steady-State Behavior of a Stochastic Process...10 9.3 Types of Simulations with Regard to Output Analysis...1 9.4
More informationPBW 654 Applied Statistics - I Urban Operations Research
PBW 654 Applied Statistics - I Urban Operations Research Lecture 2.I Queuing Systems An Introduction Operations Research Models Deterministic Models Linear Programming Integer Programming Network Optimization
More informationAnalysis of Simulated Data
Analysis of Simulated Data Media Media di una popolazione: somma di tutti i valori delle variabili della popolazione diviso il numero di unità della popolazione (N) µ N i= = 1 N X i Dove: - N = numero
More informationQueueing systems. Renato Lo Cigno. Simulation and Performance Evaluation Queueing systems - Renato Lo Cigno 1
Queueing systems Renato Lo Cigno Simulation and Performance Evaluation 2014-15 Queueing systems - Renato Lo Cigno 1 Queues A Birth-Death process is well modeled by a queue Indeed queues can be used to
More information2WB05 Simulation Lecture 7: Output analysis
2WB05 Simulation Lecture 7: Output analysis Marko Boon http://www.win.tue.nl/courses/2wb05 December 17, 2012 Outline 2/33 Output analysis of a simulation Confidence intervals Warm-up interval Common random
More informationComputer Science, Informatik 4 Communication and Distributed Systems. Simulation. Discrete-Event System Simulation. Dr.
Simulation Discrete-Event System Simulation Chapter 9 Verification and Validation of Simulation Models Purpose & Overview The goal of the validation process is: To produce a model that represents true
More informationProbabilities & Statistics Revision
Probabilities & Statistics Revision Christopher Ting Christopher Ting http://www.mysmu.edu/faculty/christophert/ : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 January 6, 2017 Christopher Ting QF
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 informationDiscrete-event simulations
Discrete-event simulations Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/elt-53606/ OUTLINE: Why do we need simulations? Step-by-step simulations; Classifications;
More informationChapter 6 Queueing Models. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation
Chapter 6 Queueing Models Banks, Carson, Nelson & Nicol Discrete-Event System Simulation Purpose Simulation is often used in the analysis of queueing models. A simple but typical queueing model: Queueing
More informationReview of Queuing Models
Review of Queuing Models Recitation, Apr. 1st Guillaume Roels 15.763J Manufacturing System and Supply Chain Design http://michael.toren.net/slides/ipqueue/slide001.html 2005 Guillaume Roels Outline Overview,
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 informationL = λ W time average number in Line or system arrival rate W average Waiting time per customer
IEOR 4615, Lecture 3, January 27, 2015 L λ Little s Law* L = λ W time average number in Line or system arrival rate W average Waiting time per customer *J. D. C. Little, A proof of the queueing formula:
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 information14 Random Variables and Simulation
14 Random Variables and Simulation In this lecture note we consider the relationship between random variables and simulation models. Random variables play two important roles in simulation models. We assume
More informationSEQUENTIAL ESTIMATION OF THE STEADY-STATE VARIANCE IN DISCRETE EVENT SIMULATION
SEQUENTIAL ESTIMATION OF THE STEADY-STATE VARIANCE IN DISCRETE EVENT SIMULATION Adriaan Schmidt Institute for Theoretical Information Technology RWTH Aachen University D-5056 Aachen, Germany Email: Adriaan.Schmidt@rwth-aachen.de
More informationPerformance Evaluation
Performance Evaluation Chapter 5 Short Probability Primer & Obtaining Results (Acknowledgement: these slides have mostly been compiled from [Kar04, Rob0, BLK0]) otivation! Recall:! Usually, a simulation
More informationIE 581 Introduction to Stochastic Simulation. Three pages of hand-written notes, front and back. Closed book. 120 minutes.
Three pages of hand-written notes, front and back. Closed book. 120 minutes. 1. SIMLIB. True or false. (2 points each. If you wish, write an explanation of your thinking.) (a) T F The time of the next
More informationDesign of discrete-event simulations
Design of discrete-event simulations Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/tlt-2707/ OUTLINE: Discrete event simulation; Event advance design; Unit-time
More informationDiscrete-Event System Simulation
Discrete-Event System Simulation FOURTH EDITION Jerry Banks Independent Consultant John S. Carson II Brooks Automation Barry L. Nelson Northwestern University David M. Nicol University of Illinois, Urbana-Champaign
More informationFormal Statement of Simple Linear Regression Model
Formal Statement of Simple Linear Regression Model Y i = β 0 + β 1 X i + ɛ i Y i value of the response variable in the i th trial β 0 and β 1 are parameters X i is a known constant, the value of the predictor
More informationF O R SOCI AL WORK RESE ARCH
7 TH EUROPE AN CONFERENCE F O R SOCI AL WORK RESE ARCH C h a l l e n g e s i n s o c i a l w o r k r e s e a r c h c o n f l i c t s, b a r r i e r s a n d p o s s i b i l i t i e s i n r e l a t i o n
More informationESTIMATION AND OUTPUT ANALYSIS (L&K Chapters 9, 10) Review performance measures (means, probabilities and quantiles).
ESTIMATION AND OUTPUT ANALYSIS (L&K Chapters 9, 10) Set up standard example and notation. Review performance measures (means, probabilities and quantiles). A framework for conducting simulation experiments
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 informationMore on Input Distributions
More on Input Distributions Importance of Using the Correct Distribution Replacing a distribution with its mean Arrivals Waiting line Processing order System Service mean interarrival time = 1 minute mean
More informationprobability of k samples out of J fall in R.
Nonparametric Techniques for Density Estimation (DHS Ch. 4) n Introduction n Estimation Procedure n Parzen Window Estimation n Parzen Window Example n K n -Nearest Neighbor Estimation Introduction Suppose
More informationStatistical inference
Statistical inference Contents 1. Main definitions 2. Estimation 3. Testing L. Trapani MSc Induction - Statistical inference 1 1 Introduction: definition and preliminary theory In this chapter, we shall
More informationChapter 10 Verification and Validation of Simulation Models. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation
Chapter 10 Verification and Validation of Simulation Models Banks, Carson, Nelson & Nicol Discrete-Event System Simulation Purpose & Overview The goal of the validation process is: To produce a model that
More informationCHAPTER 3 STOCHASTIC MODEL OF A GENERAL FEED BACK QUEUE NETWORK
CHAPTER 3 STOCHASTIC MODEL OF A GENERAL FEED BACK QUEUE NETWORK 3. INTRODUCTION: Considerable work has been turned out by Mathematicians and operation researchers in the development of stochastic and simulation
More informationPower Changes in a Critical Reactor. The Critical Reactor
Chapter 8 Power Changes in a Critical Reactor n For very small reactivity increases n For small reactivity increases n For large reactivity increases/decreases The Critical Reactor < k = hfpel f L t =
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 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 informationVerification and Validation. CS1538: Introduction to Simulations
Verification and Validation CS1538: Introduction to Simulations Steps in a Simulation Study Problem & Objective Formulation Model Conceptualization Data Collection Model translation, Verification, Validation
More informationChapter 2 SIMULATION BASICS. 1. Chapter Overview. 2. Introduction
Chapter 2 SIMULATION BASICS 1. Chapter Overview This chapter has been written to introduce the topic of discreteevent simulation. To comprehend the material presented in this chapter, some background in
More information2.3 Estimating PDFs and PDF Parameters
.3 Estimating PDFs and PDF Parameters estimating means - discrete and continuous estimating variance using a known mean estimating variance with an estimated mean estimating a discrete pdf estimating a
More informationCOSC460 Honours Report: A Sequential Steady-State Detection Method for Quantitative Discrete-Event Simulation
COSC460 Honours Report: A Sequential Steady-State Detection Method for Quantitative Discrete-Event Simulation December 1, 2012 Adam Freeth afr58@uclive.ac.nz Department of Computer Science and Software
More informationSystems Simulation Chapter 6: Queuing Models
Systems Simulation Chapter 6: Queuing Models Fatih Cavdur fatihcavdur@uludag.edu.tr April 2, 2014 Introduction Introduction Simulation is often used in the analysis of queuing models. A simple but typical
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 informationMarkov Processes and Queues
MIT 2.853/2.854 Introduction to Manufacturing Systems Markov Processes and Queues Stanley B. Gershwin Laboratory for Manufacturing and Productivity Massachusetts Institute of Technology Markov Processes
More informationSection 8.1: Interval Estimation
Section 8.1: Interval Estimation Discrete-Event Simulation: A First Course c 2006 Pearson Ed., Inc. 0-13-142917-5 Discrete-Event Simulation: A First Course Section 8.1: Interval Estimation 1/ 35 Section
More informationPractice Problems Section Problems
Practice Problems Section 4-4-3 4-4 4-5 4-6 4-7 4-8 4-10 Supplemental Problems 4-1 to 4-9 4-13, 14, 15, 17, 19, 0 4-3, 34, 36, 38 4-47, 49, 5, 54, 55 4-59, 60, 63 4-66, 68, 69, 70, 74 4-79, 81, 84 4-85,
More information1.225J J (ESD 205) Transportation Flow Systems
1.225J J (ESD 25) Transportation Flow Systems Lecture 9 Simulation Models Prof. Ismail Chabini and Prof. Amedeo R. Odoni Lecture 9 Outline About this lecture: It is based on R16. Only material covered
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 informationARD: AN AUTOMATED REPLICATION-DELETION METHOD FOR SIMULATION ANALYSIS
ARD: AN AUTOMATED REPLICATION-DELETION METHOD FOR SIMULATION ANALYSIS Emily K. Lada Anup C. Mokashi SAS Institute Inc. 100 SAS Campus Drive, R5413 Cary, NC 27513-8617, USA James R. Wilson North Carolina
More informationREAL-TIME DELAY ESTIMATION BASED ON DELAY HISTORY SUPPLEMENTARY MATERIAL
REAL-TIME DELAY ESTIMATION BASED ON DELAY HISTORY SUPPLEMENTARY MATERIAL by Rouba Ibrahim and Ward Whitt IEOR Department Columbia University {rei2101, ww2040}@columbia.edu Abstract Motivated by interest
More informationExercises Solutions. Automation IEA, LTH. Chapter 2 Manufacturing and process systems. Chapter 5 Discrete manufacturing problems
Exercises Solutions Note, that we have not formulated the answers for all the review questions. You will find the answers for many questions by reading and reflecting about the text in the book. Chapter
More informationModeling and Performance Analysis with Discrete-Event Simulation
Simulation Modeling and Performance Analysis with Discrete-Event Simulation Chapter 9 Input Modeling Contents Data Collection Identifying the Distribution with Data Parameter Estimation Goodness-of-Fit
More informationM/G/1 and Priority Queueing
M/G/1 and Priority Queueing Richard T. B. Ma School of Computing National University of Singapore CS 5229: Advanced Compute Networks Outline PASTA M/G/1 Workload and FIFO Delay Pollaczek Khinchine Formula
More informationVariance reduction techniques
Variance reduction techniques Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/elt-53606/ OUTLINE: Simulation with a given accuracy; Variance reduction techniques;
More informationSingle-part-type, multiple stage systems. Lecturer: Stanley B. Gershwin
Single-part-type, multiple stage systems Lecturer: Stanley B. Gershwin Flow Line... also known as a Production or Transfer Line. M 1 B 1 M 2 B 2 M 3 B 3 M 4 B 4 M 5 B 5 M 6 Machine Buffer Machines are
More informationSlides 9: Queuing Models
Slides 9: Queuing Models Purpose Simulation is often used in the analysis of queuing models. A simple but typical queuing model is: Queuing models provide the analyst with a powerful tool for designing
More informationComputer Science, Informatik 4 Communication and Distributed Systems. Simulation. Discrete-Event System Simulation. Dr.
Simulation Discrete-Event System Simulation Chapter 8 Input Modeling Purpose & Overview Input models provide the driving force for a simulation model. The quality of the output is no better than the quality
More information0utline. 1. Tools from Operations Research. 2. Applications
0utline 1. Tools from Operations Research Little s Law (average values) Unreliable Machine(s) (operation dependent) Buffers (zero buffers & infinite buffers) M/M/1 Queue (effects of variation) 2. Applications
More informationInference for the Regression Coefficient
Inference for the Regression Coefficient Recall, b 0 and b 1 are the estimates of the slope β 1 and intercept β 0 of population regression line. We can shows that b 0 and b 1 are the unbiased estimates
More informationSimulation optimization via bootstrapped Kriging: Survey
Simulation optimization via bootstrapped Kriging: Survey Jack P.C. Kleijnen Department of Information Management / Center for Economic Research (CentER) Tilburg School of Economics & Management (TiSEM)
More informationOn the Partitioning of Servers in Queueing Systems during Rush Hour
On the Partitioning of Servers in Queueing Systems during Rush Hour This paper is motivated by two phenomena observed in many queueing systems in practice. The first is the partitioning of server capacity
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 informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
More informationMonte Carlo Methods for Computation and Optimization (048715)
Technion Department of Electrical Engineering Monte Carlo Methods for Computation and Optimization (048715) Lecture Notes Prof. Nahum Shimkin Spring 2015 i PREFACE These lecture notes are intended for
More informationSession-Based Queueing Systems
Session-Based Queueing Systems Modelling, Simulation, and Approximation Jeroen Horters Supervisor VU: Sandjai Bhulai Executive Summary Companies often offer services that require multiple steps on the
More informationDifferentiating Functions & Expressions - Edexcel Past Exam Questions
- Edecel Past Eam Questions. (a) Differentiate with respect to (i) sin + sec, (ii) { + ln ()}. 5-0 + 9 Given that y =, ¹, ( -) 8 (b) show that = ( -). (6) June 05 Q. f() = e ln, > 0. (a) Differentiate
More informationARD: AN AUTOMATED REPLICATION-DELETION METHOD FOR SIMULATION ANALYSIS
Proceedings of the 2013 Winter Simulation Conference R. Pasupathy, S.-H. Kim, A. Tolk, R. Hill, and M. E. Kuhl, eds. ARD: AN AUTOMATED REPLICATION-DELETION METHOD FOR SIMULATION ANALYSIS Emily K. Lada
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 informationCDA5530: Performance Models of Computers and Networks. Chapter 4: Elementary Queuing Theory
CDA5530: Performance Models of Computers and Networks Chapter 4: Elementary Queuing Theory Definition Queuing system: a buffer (waiting room), service facility (one or more servers) a scheduling policy
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.262 Discrete Stochastic Processes Midterm Quiz April 6, 2010 There are 5 questions, each with several parts.
More informationGEOMETRIC -discrete A discrete random variable R counts number of times needed before an event occurs
STATISTICS 4 Summary Notes. Geometric and Exponential Distributions GEOMETRIC -discrete A discrete random variable R counts number of times needed before an event occurs P(X = x) = ( p) x p x =,, 3,...
More informationAnalysis of Call Center Data
University of Pennsylvania ScholarlyCommons Wharton Research Scholars Wharton School 4-1-2004 Analysis of Call Center Data Yu Chu Cheng University of Pennsylvania Follow this and additional works at: http://repository.upenn.edu/wharton_research_scholars
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 informationPart I Stochastic variables and Markov chains
Part I Stochastic variables and Markov chains Random variables describe the behaviour of a phenomenon independent of any specific sample space Distribution function (cdf, cumulative distribution function)
More informationWooldridge, Introductory Econometrics, 4th ed. Appendix C: Fundamentals of mathematical statistics
Wooldridge, Introductory Econometrics, 4th ed. Appendix C: Fundamentals of mathematical statistics A short review of the principles of mathematical statistics (or, what you should have learned in EC 151).
More informationMultilayer Neural Networks
Multilayer Neural Networks Introduction Goal: Classify objects by learning nonlinearity There are many problems for which linear discriminants are insufficient for minimum error In previous methods, the
More informationContents LIST OF TABLES... LIST OF FIGURES... xvii. LIST OF LISTINGS... xxi PREFACE. ...xxiii
LIST OF TABLES... xv LIST OF FIGURES... xvii LIST OF LISTINGS... xxi PREFACE...xxiii CHAPTER 1. PERFORMANCE EVALUATION... 1 1.1. Performance evaluation... 1 1.2. Performance versus resources provisioning...
More informationThe Scale Factor: A New Degree of Freedom in Phase Type Approximation
The Scale Factor: A New Degree of Freedom in Phase Type Approximation Andrea Bobbio DISTA, Università del Piemonte Orientale, Alessandria, Italy,bobbio@unipmn.it András Horváth, Miklós Telek Dept. of Telecommunications,
More informationChapter 10 Verification and Validation of Simulation Models. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation
Chapter 10 Verification and Validation of Simulation Models Banks, Carson, Nelson & Nicol Discrete-Event System Simulation The Black Box [Bank Example: Validate I-O Transformation] A model was developed
More informationExpected Time Delay in Multi-Item Inventory Systems with Correlated Demands
Expected Time Delay in Multi-Item Inventory Systems with Correlated Demands Rachel Q. Zhang Department of Industrial and Operations Engineering, University of Michigan, Ann Arbor, Michigan 48109 Received
More information6 Solving Queueing Models
6 Solving Queueing Models 6.1 Introduction In this note we look at the solution of systems of queues, starting with simple isolated queues. The benefits of using predefined, easily classified queues will
More informationDiscrete Event Simulation
Discrete Event Simulation Henry Z. Lo June 9, 2014 1 Motivation Suppose you were a branch manager of a bank. How many tellers should you staff to keep customers happy? With prior knowledge of how often
More informationA STAFFING ALGORITHM FOR CALL CENTERS WITH SKILL-BASED ROUTING: SUPPLEMENTARY MATERIAL
A STAFFING ALGORITHM FOR CALL CENTERS WITH SKILL-BASED ROUTING: SUPPLEMENTARY MATERIAL by Rodney B. Wallace IBM and The George Washington University rodney.wallace@us.ibm.com Ward Whitt Columbia University
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 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 informationLECTURE #6 BIRTH-DEATH PROCESS
LECTURE #6 BIRTH-DEATH PROCESS 204528 Queueing Theory and Applications in Networks Assoc. Prof., Ph.D. (รศ.ดร. อน นต ผลเพ ม) Computer Engineering Department, Kasetsart University Outline 2 Birth-Death
More informationDr. Maddah ENMG 617 EM Statistics 10/15/12. Nonparametric Statistics (2) (Goodness of fit tests)
Dr. Maddah ENMG 617 EM Statistics 10/15/12 Nonparametric Statistics (2) (Goodness of fit tests) Introduction Probability models used in decision making (Operations Research) and other fields require fitting
More informationBias Variance Trade-off
Bias Variance Trade-off The mean squared error of an estimator MSE(ˆθ) = E([ˆθ θ] 2 ) Can be re-expressed MSE(ˆθ) = Var(ˆθ) + (B(ˆθ) 2 ) MSE = VAR + BIAS 2 Proof MSE(ˆθ) = E((ˆθ θ) 2 ) = E(([ˆθ E(ˆθ)]
More informationStochastic Models in Computer Science A Tutorial
Stochastic Models in Computer Science A Tutorial Dr. Snehanshu Saha Department of Computer Science PESIT BSC, Bengaluru WCI 2015 - August 10 to August 13 1 Introduction 2 Random Variable 3 Introduction
More information