DECISION SCIENCES INSTITUTE. Exponentially Derived Antithetic Random Numbers. (Full paper submission)
|
|
- Ellen Welch
- 5 years ago
- Views:
Transcription
1 DECISION SCIENCES INSTITUTE (Full paper submission) Dennis Ridley, Ph.D. SBI, Florida A&M University and Scientific Computing, Florida State University Pierre Ngnepieba, Ph.D. Department of Mathematics, Florida A&M University, Tallahassee, Florida. ABSTRACT Serial correlation in pseudo random numbers that are used in Monte Carlo simulations can create hidden errors in the response variable. The estimate of the response variable becomes biased, does not converge to the true value, and its variance is larger than the true value. A pair of antithetic exponential random numbers is created by power transformation of uniformly distributed random [0,1] numbers. This pair of antithetic numbers is used in two separate computer simulation experiments and response variables from the simulations combined. The combined response has bias and variance that are smaller than those of the individual responses. KEYWORDS: Exponential antithetic random variates, Simulation model bias, Variance reduction, Bias reduction, Inverse correlation.
2 DECISION SCIENCES INSTITUTE (Full paper submission) Dennis Ridley, Ph.D. SBI, Florida A&M University and Scientific Computing, Florida State University Pierre Ngnepieba, Ph.D. Department of Mathematics, Florida A&M University, Tallahassee, Florida. ABSTRACT Serial correlation in pseudo random numbers that are used in Monte Carlo simulations can create hidden errors in the response variable. The estimate of the response variable becomes biased, does not converge to the true value, and its variance is larger than the true value. A pair of antithetic exponential random numbers is created by power transformation of uniformly distributed random [0,1] numbers. This pair of antithetic numbers is used in two separate computer simulation experiments and response variables from the simulations combined. The combined response has bias and variance that are smaller than those of the individual responses. KEYWORDS: Exponential antithetic random variates, Simulation model bias, Variance reduction, Bias reduction, Inverse correlation. INTRODUCTION Background Hammersley & Morton (1956) suggested a very simple method for obtaining inversely correlated random numbers for use in Monte Carlo computer simulation experiments. A computer simulation experiment is conducted with a uniformly distributed random number r ~ U(0,1). The simulation is repeated with the inversely correlated complementary random number 1-r. Finally, the response variables from the two simulations are averaged. This can result in a variance in the response variable that is smaller than that for either simulation. However, there is no guarantee that the variance will decrease. In practice the variance can actually increase. See also Kleijnen (1975) Antithetic random variables by power transformation was introduced by Ridley(1999). The Ridley (1999) antithetic time series theorem states that if XX tt > 0, tt = 1,2,3, is a discrete realization of a lognormal stochastic process, such that ln XX tt ~NN(µ, σσ), then if the correlation between XX tt and XX tt pp is ρρ XXXX pp, then lim pp 0,σσ 0 ρρ XXXX pp = 1. This resulted in large reduction in time series model fitted and forecast values. See also Ridley (1995, 1997, 2013 & 2014) for applications in time series analysis and forecasting. In this paper we are interested in applications in computer simulation where exponential random variates are called for.
3 Proposed research. The common practice in computer simulation is to generate uniformly distributed pseudo random numbers on the interval from 0 to 1. These rr~u(0,1) numbers are then used to generate random variates from the cumulative distribution that correspond to the actual random variable under study. The original random numbers and therefore the random variates are assumed to be identically and independently distributed. In practice however, the random numbers are generated from a mathematical formula. One simple method is the linear congruential generator. Random numbers generated from a mathematical formula have the advantage that very large sets can be reproduced for repetition and comparison of the outcomes from different computer simulations under different conditions. However, if the random numbers contain any traces of serial correlations, they can introduce hidden errors into the response variable. For an example see Ferrenberg, Landau & Wong (1992) where so called good random numbers were used and inaccurate results were obtained, but where better results were obtained using simple methods. The objectives of this research are two-fold. One objective is to correct for the bad effects whether they are caused by good or bad random numbers. The other objective is to reduce the variance of the estimate of the response function. This is particularly important when very long simulations are required. Some physics experiments utilize as much as in length. Variance reduction can reduce the length that is necessary for the same level of accuracy. The paper is organized as follows. In section 2 we investigate the correlation between rr and rr pp as p 0. In section 3 we investigate the correlation between rr pp and rr pp as p 0. In section 4 we investigate the distributions of rr, rr pp and rr pp. In section 5 we carry out computer simulation experiments using exponential antithetic random variates and compare the results with those obtained by traditional random numbers. Section 6 contains conclusions and suggestions for future research. CORRELATION BETWEEN rr AND rr pp vs p. First we consider the correlation between rr and rr pp. Correlations for different values of p as p approaches zero from the left are listed in Table 1 and drawn in Figure 1. The numbers were generated with Matlab (2008) and calculated as the average of 1000 samples of 500 numbers each. The correlation is seen to approach recurring, or - 3/2. Therefore rr and rr pp are antithetic in the sense that they are negatively correlated. Table 1: Correlation between rr and rr pp vs p p Corr(r,rr pp ) Figure 1: Correlation between rr and rr pp vs p
4 CORRELATION BETWEEN rr pp AND rr pp vs p. Next, we consider the correlation between rr pp and rr pp. Correlations for different values of p as p approaches zero from the left are listed in Table 2 and drawn in Figure 2. In this case, the sample correlations converge very rapidly in just one sample of 100 random numbers. The correlation is seen to approach -1. Therefore, rr pp and rr pp, are good candidates for use as antithetic random variates. Table 2: Correlation between rr pp and rr pp vs p p Corr(rr pp,rr pp ) Figure 2: Correlation between rr pp and rr pp vs p DISTRIBUTION r, rr pp AND rr pp. Before we can proceed, we need to understand the distributional characteristics of rr, rr pp and rr pp. The histograms for rr, rr pp and rr pp are given in Figure 3. The histograms are based on 300 random numbers rr~u(0,1). As illustrated, in the limit as pp 0, the distributions of rr pp and rr pp are similar to negative and positive exponential distributions, respectively. Figure 3: Distributional characteristics pp 0 a. Uniform rr~u(0,1) b. rr pp c. rr pp
5 Source of serial correlation Any sequence of random numbers that in its totality is serially perfectly uncorrelated, can contain what appears to be subsets of correlated sequences. The smallest trace of serial correlation in any such subset can create a trend that biases the random numbers to the high or low side of the true value. Even if an opposite trend develops subsequently, a very large number of random numbers can occur before the net bias is corrected. Either way, apparent subsets of traces of serial correlation introduce unreliability into the simulation. Positively Correlated Random Numbers Approach Consider the regression of rr pp on rr pp rr pp = ff(rr pp ) + εε = αα + ββrr pp + εε. Since lim pp 0 Corr(rrpp, rr pp ) = 1, that is rr pp and rr pp are perfectly correlated, it follows that lim pp 0 rrpp = αα + ββrr pp + εε wwheeeeee εε 0. That is, rr pp αα + ββrr pp. Let the fitted values of rr pp be denoted by rr pp, let αα and ββ denote least squares estimates and let s denote standard deviation. Then, rr pp = αα + ββ rr pp = rr pp + Corr(rr pp, rr pp )(ss rr pp/ss rr pp)(rr pp ). rr pp This random variable is negatively correlated with rr. However, unlike rr, the probability pp pp distribution of rr is not uniform. For purpose of comparison and uniformity, rr can be converted back to the original scale of and distribution of rr by, rr = (rr pp ) pp. Based on this formulation of the pair of random numbers rr and rr, p=-0.01, the combined response from two separate computer simulations based on rr and rr was calculated. As it turns out, with rr and rr positively correlated, the result was two identical values for the response variables and no gain from combining. Antithetic Random Numbers Random Variates Approach Consider the negatively correlated random number random variate pair rr and rr pp. This suggest the combination of responses from two computer simulations based on rr and rr pp. Consider the simple queueing system shown in Figure 4. The time between arrivals are exponentially distributed with mean (1/ λ). The mean arrival rate= λ. The exponential probability density function is given by f(x)= λλee λλλλ for XX 0 0 eeeeeeeeeeheeeeee The corresponding cumulative distribution is given by 1 F(x)=P(X<x)= 1 ee λλλλ for xx 0 0 eeeeeeeeeeheeeeee
6 The inter arrival times are calculated by setting rr equal to the cumulative distribution rr = 1 ee λλλλ The inter arrival times are given by x=-(1/ λ)ln(1-r) For the antithetic variates rr pp, the inter arrival times are found by converting the mean and standard deviation. This is accomplished by subtracting its sample average rr pp and dividing by the standard deviation (ss rr pp), then multiply by the population standard deviation (1/ λ) and adding the population mean (1/ λ) time of interest as follows. Antithetic Random Variates Approach xx = rrpp rr pp ss rr ppλλ + 1 λλ Consider the perfectly negatively correlated antithetic random variates rr pp and rr pp. Based on this formulation of 300 pairs of random variates rr pp, rr pp, p=-0.01, the following were calculated: rr pp = , rr pp = , Corr rr pp, rr pp ii = , ss rr pp = , ss rr pp = , rr pp = The following averages were calculated: rr = , rr = , 0.5rr +0.5rr = Therefore, the average of the equally weighted combined random number is closer to the true value of 0.5 for the uniform distribution than either individual average. The corresponding standard deviation deviations were calculated: ss rr = , ss rr = , (0.5ss 2 rr + 0.5ss 2 rr ) 1/2 = Therefore, the standard deviation of the equally weighted combined random number is closer to the true value of 1/ 12= for the uniform distribution than either individual standard deviation. This suggests the combination of responses from two computer simulations based rr pp and rr pp. However, this experiment is beyond the scope of this paper. COMPUTER SIMULATION EXPERIMENT To explore the possibility of reducing bias and variance by combining the estimates of the responses from two Monte Carlo experiments, one based on xx and one based on xx, a simple queue will be analyzed as shown in Figure 4. In the simple queue, the source population for arriving entities is infinite. The queue is not capacitated. There is one server. The entities are processed on a first come served basis. The time between arrivals and the service times are exponentially distributed. The entity arrival rate is λ. The service rate is μ.
7 Figure 4. Simple queueing system. Arriving entities Queue Server Departing entities Arrival rate λ xx =-(1/ λ)ln(1-r) Service rate μ TT Arriving entities Queue Server Departing entities Combined response ωωωω + (1 ωω)tt Arrival rate λ xx = rrpp rr pp ss rr ppλλ + 1 λλ Service rate μ TT The response variable chosen is the average time that an entity spends in the system. The average time that entities spend in the system is well known to be precisely TT ttheeeeeeeeeeeeee =1/( μ - λ). Consider two Monte Carlo simulations. One simulation uses the exponential random numbers xx =-(1/ λ)ln(1-r), rr~uu(0,1) with response TT, and the other simulation uses the exponential random numbers xx = (rr pp rr )/ss pp rr ppλλ + 1 λλ, rr~uu(0,1), pp = 0.01 with response TT. The combined response is ωωωω + (1 ωω)tt where ωω and (1 ωω) are combining weights. The results in Table 3 are based on ωω = 0.5. The length of the simulation is 40,000 arrivals. Before antithetic combining, the average time in the system shows a bias of 1.8 minutes. After combining, the average time in the system has a smaller bias of 1.1. Before antithetic combining, the standard deviation of the time in the system is 14.4 minutes. After combining, the standard deviation of time in the system is 13.5 minutes. Table 3: Computer simulation results Simulation step number 16k 20k 24k 28k 32k 36k 40k TT Bias Standard Deviation of TT TT Standard Deviation of TT ωωωω + (1 ωω)tt Bias Standard Deviation of ωωωω + (1 ωω)tt
8 CONCLUSIONS Antithetic random number random variate pairs were created. These antithetic variables were used to perform two separate Monte Carlo computer simulations. The two estimates of the response variable were combined. The combined estimate exhibited better statistical properties than those of the uncombined estimates. The combined average of the estimated response variable had smaller bias. REFERENCES Ferrenberg, A.M., Landau, D.P., & Wong, Y. J. (1992). Monte Carlo Simulations: Hidden errors from good random number generators, Physical Review Letters, 69(23), Hammersley, J.M., & Morton, K.W. (1956). A new Monte Carlo technique: antithetic variates, Mathematical Proceedings of the Cambridge Philosophical Society, 52(3), Kleijnen, J. P. C. (1975). Antithetic variates, common random numbers and optimal computer time allocation in simulation, Management Science, 21(10), MATLAB (2008). Application Program Interface Reference, Version 8, The MathWorks, Inc. Ridley, A.D. (1995). Combining Global Antithetic Forecasts, International Transactions in Operational Research, 2(4), Ridley, A.D. (1997). Optimal Weights for Combining Antithetic Forecasts, Computers & Industrial Engineering, 32(2), Ridley, A.D. (1999). Optimal Antithetic Weights for Lognormal Time Series Forecasting, Computers and Operations Research, 26(3), Ridley, A.D., Ngnepieba, P., & Duke, D. (2013). Parameter Optimization for Combining Lognormal Antithetic Time Series, European Journal of Mathematical Sciences, 2(2), Ridley, A.D., & Ngnepieba, P. (2014). Antithetic time series analysis and the CompanyX data, Journal of the Royal Statistical Society, A, 177(1),
Computational Methods. Randomness and Monte Carlo Methods
Computational Methods Randomness and Monte Carlo Methods Manfred Huber 2010 1 Randomness and Monte Carlo Methods Introducing randomness in an algorithm can lead to improved efficiencies Random sampling
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 13 Random Numbers and Stochastic Simulation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright
More informationWhat We ll Do... Random
What We ll Do... Random- number generation Random Number Generation Generating random variates Nonstationary Poisson processes Variance reduction Sequential sampling Designing and executing simulation
More informationGenerating random samples from user-defined distributions
The Stata Journal (2011) 11, Number 2, pp. 299 304 Generating random samples from user-defined distributions Katarína Lukácsy Central European University Budapest, Hungary lukacsy katarina@phd.ceu.hu Abstract.
More informationMULTI-DIMENSIONAL MONTE CARLO INTEGRATION
CS580: Computer Graphics KAIST School of Computing Chapter 3 MULTI-DIMENSIONAL MONTE CARLO INTEGRATION 2 1 Monte Carlo Integration This describes a simple technique for the numerical evaluation of integrals
More informationLecture: Simulation. of Manufacturing Systems. Sivakumar AI. Simulation. SMA6304 M2 ---Factory Planning and scheduling. Simulation - A Predictive Tool
SMA6304 M2 ---Factory Planning and scheduling Lecture Discrete Event of Manufacturing Systems Simulation Sivakumar AI Lecture: 12 copyright 2002 Sivakumar 1 Simulation Simulation - A Predictive Tool Next
More informationA simple OMNeT++ queuing experiment using different random number generators
A simple OMNeT++ queuing experiment using different random number generators Bernhard Hechenleitner and Karl Entacher December 5, 2002 Abstract We apply a simple queuing-experiment using parallel streams
More informationCombinatorial Search; Monte Carlo Methods
Combinatorial Search; Monte Carlo Methods Parallel and Distributed Computing Department of Computer Science and Engineering (DEI) Instituto Superior Técnico May 02, 2016 CPD (DEI / IST) Parallel and Distributed
More informationProbability and Statistics for Final Year Engineering Students
Probability and Statistics for Final Year Engineering Students By Yoni Nazarathy, Last Updated: April 11, 2011. Lecture 1: Introduction and Basic Terms Welcome to the course, time table, assessment, etc..
More informationVARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung
POLYTECHNIC UNIVERSITY Department of Computer and Information Science VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung Abstract: Techniques for reducing the variance in Monte Carlo
More informationMonte Carlo Methods; Combinatorial Search
Monte Carlo Methods; Combinatorial Search Parallel and Distributed Computing Department of Computer Science and Engineering (DEI) Instituto Superior Técnico November 22, 2012 CPD (DEI / IST) Parallel and
More informationLearner Expectations UNIT 1: GRAPICAL AND NUMERIC REPRESENTATIONS OF DATA. Sept. Fathom Lab: Distributions and Best Methods of Display
CURRICULUM MAP TEMPLATE Priority Standards = Approximately 70% Supporting Standards = Approximately 20% Additional Standards = Approximately 10% HONORS PROBABILITY AND STATISTICS Essential Questions &
More informationDiscrete Event Simulation & VHDL. Prof. K. J. Hintz Dept. of Electrical and Computer Engineering George Mason University
Discrete Event Simulation & VHDL Prof. K. J. Hintz Dept. of Electrical and Computer Engineering George Mason University Discrete Event Simulation Material from VHDL Programming with Advanced Topics by
More informationQuasi-Monte Carlo Methods Combating Complexity in Cost Risk Analysis
Quasi-Monte Carlo Methods Combating Complexity in Cost Risk Analysis Blake Boswell Booz Allen Hamilton ISPA / SCEA Conference Albuquerque, NM June 2011 1 Table Of Contents Introduction Monte Carlo Methods
More informationRandom Number Generation and Monte Carlo Methods
James E. Gentle Random Number Generation and Monte Carlo Methods With 30 Illustrations Springer Contents Preface vii 1 Simulating Random Numbers from a Uniform Distribution 1 1.1 Linear Congruential Generators
More informationA Random Number Based Method for Monte Carlo Integration
A Random Number Based Method for Monte Carlo Integration J Wang and G Harrell Department Math and CS, Valdosta State University, Valdosta, Georgia, USA Abstract - A new method is proposed for Monte Carlo
More informationYou ve already read basics of simulation now I will be taking up method of simulation, that is Random Number Generation
Unit 5 SIMULATION THEORY Lesson 39 Learning objective: To learn random number generation. Methods of simulation. Monte Carlo method of simulation You ve already read basics of simulation now I will be
More information* Hyun Suk Park. Korea Institute of Civil Engineering and Building, 283 Goyangdae-Ro Goyang-Si, Korea. Corresponding Author: Hyun Suk Park
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 11 (November 2017), PP.47-59 Determination of The optimal Aggregation
More informationThe first few questions on this worksheet will deal with measures of central tendency. These data types tell us where the center of the data set lies.
Instructions: You are given the following data below these instructions. Your client (Courtney) wants you to statistically analyze the data to help her reach conclusions about how well she is teaching.
More informationMCMC Methods for data modeling
MCMC Methods for data modeling Kenneth Scerri Department of Automatic Control and Systems Engineering Introduction 1. Symposium on Data Modelling 2. Outline: a. Definition and uses of MCMC b. MCMC algorithms
More informationIntegers & Absolute Value Properties of Addition Add Integers Subtract Integers. Add & Subtract Like Fractions Add & Subtract Unlike Fractions
Unit 1: Rational Numbers & Exponents M07.A-N & M08.A-N, M08.B-E Essential Questions Standards Content Skills Vocabulary What happens when you add, subtract, multiply and divide integers? What happens when
More informationMonte Carlo Ray Tracing. Computer Graphics CMU /15-662
Monte Carlo Ray Tracing Computer Graphics CMU 15-462/15-662 TODAY: Monte Carlo Ray Tracing How do we render a photorealistic image? Put together many of the ideas we ve studied: - color - materials - radiometry
More informationSimulation: Solving Dynamic Models ABE 5646 Week 12, Spring 2009
Simulation: Solving Dynamic Models ABE 5646 Week 12, Spring 2009 Week Description Reading Material 12 Mar 23- Mar 27 Uncertainty and Sensitivity Analysis Two forms of crop models Random sampling for stochastic
More informationThemes in the Texas CCRS - Mathematics
1. Compare real numbers. a. Classify numbers as natural, whole, integers, rational, irrational, real, imaginary, &/or complex. b. Use and apply the relative magnitude of real numbers by using inequality
More informationMonte Carlo Integration and Random Numbers
Monte Carlo Integration and Random Numbers Higher dimensional integration u Simpson rule with M evaluations in u one dimension the error is order M -4! u d dimensions the error is order M -4/d u In general
More informationChapter 1. Introduction
Chapter 1 Introduction A Monte Carlo method is a compuational method that uses random numbers to compute (estimate) some quantity of interest. Very often the quantity we want to compute is the mean of
More informationTowards Interactive Global Illumination Effects via Sequential Monte Carlo Adaptation. Carson Brownlee Peter S. Shirley Steven G.
Towards Interactive Global Illumination Effects via Sequential Monte Carlo Adaptation Vincent Pegoraro Carson Brownlee Peter S. Shirley Steven G. Parker Outline Motivation & Applications Monte Carlo Integration
More informationComputer Simulation on Structural Plane of Ore-rock By Monte-Carlo Method Pan Dong 1,2, a, Li XiBing 1,b, WangYu 1, c
2nd International Conference on Electronic & Mechanical Engineering and Ination Technology (EMEIT-2012) Computer Simulation on Structural Plane of Ore-rock By Monte-Carlo Method Pan Dong 1,2, a, Li XiBing
More informationUnit WorkBook 2 Level 4 ENG U2 Engineering Maths LO2 Statistical Techniques 2018 UniCourse Ltd. All Rights Reserved. Sample
Pearson BTEC Levels 4 and 5 Higher Nationals in Engineering (RQF) Unit 2: Engineering Maths (core) Unit Workbook 2 in a series of 4 for this unit Learning Outcome 2 Statistical Techniques Page 1 of 37
More informationApplication of Importance Sampling in Simulation of Buffer Policies in ATM networks
Application of Importance Sampling in Simulation of Buffer Policies in ATM networks SAMAD S. KOLAHI School of Computing and Information Systems Unitec New Zealand Carrington Road, Mt Albert, Auckland NEW
More informationTEACHING VARIANCE REDUCTION IN MATLAB
TEACHING VARIANCE REDUCTION IN MATLAB Jaroslav Sklenar Department of Statistics and Operations Research University of Malta Msida MSD 28, Malta jaroslav.sklenar@um.edu.mt ABSTRACT The paper describes an
More informationEstimating survival from Gray s flexible model. Outline. I. Introduction. I. Introduction. I. Introduction
Estimating survival from s flexible model Zdenek Valenta Department of Medical Informatics Institute of Computer Science Academy of Sciences of the Czech Republic I. Introduction Outline II. Semi parametric
More informationRandom Number Generation. Biostatistics 615/815 Lecture 16
Random Number Generation Biostatistics 615/815 Lecture 16 Some Uses of Random Numbers Simulating data Evaluate statistical procedures Evaluate study designs Evaluate program implementations Controlling
More informationSTATISTICS (STAT) Statistics (STAT) 1
Statistics (STAT) 1 STATISTICS (STAT) STAT 2013 Elementary Statistics (A) Prerequisites: MATH 1483 or MATH 1513, each with a grade of "C" or better; or an acceptable placement score (see placement.okstate.edu).
More informationLecture 7: Monte Carlo Rendering. MC Advantages
Lecture 7: Monte Carlo Rendering CS 6620, Spring 2009 Kavita Bala Computer Science Cornell University MC Advantages Convergence rate of O( ) Simple Sampling Point evaluation Can use black boxes General
More informationLeveling Up as a Data Scientist. ds/2014/10/level-up-ds.jpg
Model Optimization Leveling Up as a Data Scientist http://shorelinechurch.org/wp-content/uploa ds/2014/10/level-up-ds.jpg Bias and Variance Error = (expected loss of accuracy) 2 + flexibility of model
More informationMonte Carlo Hospital
Monte Carlo Hospital Let s take a break from springs and things. This model will be inspired by a paper by Schmitz and Kwak [1972]. When a hospital decides to increase the number of beds administrators
More informationVerification and Validation of X-Sim: A Trace-Based Simulator
http://www.cse.wustl.edu/~jain/cse567-06/ftp/xsim/index.html 1 of 11 Verification and Validation of X-Sim: A Trace-Based Simulator Saurabh Gayen, sg3@wustl.edu Abstract X-Sim is a trace-based simulator
More informationNumerical approach estimate
Simulation Nature of simulation Numericalapproachfor investigating models of systems. Data are gathered to estimatethe true characteristics of the model. Garbage in garbage out! One of the techniques of
More informationBehavior of the sample mean. varx i = σ 2
Behavior of the sample mean We observe n independent and identically distributed (iid) draws from a random variable X. Denote the observed values by X 1, X 2,..., X n. Assume the X i come from a population
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. - Statistical Methods -
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics - Statistical Methods - Karsten Heeger heeger@wisc.edu Course Schedule and Reading course website http://neutrino.physics.wisc.edu/teaching/phys736/
More informationDiscrete Event (time) Simulation
Discrete Event (time) Simulation What is a simulation? Simulation is the process of designing a model of a real system and conducting experiments with this model for the purpose either of understanding
More informationStochastic Simulation: Algorithms and Analysis
Soren Asmussen Peter W. Glynn Stochastic Simulation: Algorithms and Analysis et Springer Contents Preface Notation v xii I What This Book Is About 1 1 An Illustrative Example: The Single-Server Queue 1
More informationFast Automated Estimation of Variance in Discrete Quantitative Stochastic Simulation
Fast Automated Estimation of Variance in Discrete Quantitative Stochastic Simulation November 2010 Nelson Shaw njd50@uclive.ac.nz Department of Computer Science and Software Engineering University of Canterbury,
More information1. Estimation equations for strip transect sampling, using notation consistent with that used to
Web-based Supplementary Materials for Line Transect Methods for Plant Surveys by S.T. Buckland, D.L. Borchers, A. Johnston, P.A. Henrys and T.A. Marques Web Appendix A. Introduction In this on-line appendix,
More informationChapter 6 Normal Probability Distributions
Chapter 6 Normal Probability Distributions 6-1 Review and Preview 6-2 The Standard Normal Distribution 6-3 Applications of Normal Distributions 6-4 Sampling Distributions and Estimators 6-5 The Central
More informationLab 5 Monte Carlo integration
Lab 5 Monte Carlo integration Edvin Listo Zec 9065-976 edvinli@student.chalmers.se October 0, 014 Co-worker: Jessica Fredby Introduction In this computer assignment we will discuss a technique for solving
More information(X 1:n η) 1 θ e 1. i=1. Using the traditional MLE derivation technique, the penalized MLEs for η and θ are: = n. (X i η) = 0. i=1 = 1.
EXAMINING THE PERFORMANCE OF A CONTROL CHART FOR THE SHIFTED EXPONENTIAL DISTRIBUTION USING PENALIZED MAXIMUM LIKELIHOOD ESTIMATORS: A SIMULATION STUDY USING SAS Austin Brown, M.S., University of Northern
More informationLearning Objectives. Continuous Random Variables & The Normal Probability Distribution. Continuous Random Variable
Learning Objectives Continuous Random Variables & The Normal Probability Distribution 1. Understand characteristics about continuous random variables and probability distributions 2. Understand the uniform
More informationProbability Models.S4 Simulating Random Variables
Operations Research Models and Methods Paul A. Jensen and Jonathan F. Bard Probability Models.S4 Simulating Random Variables In the fashion of the last several sections, we will often create probability
More informationRadial Basis Function (RBF) Neural Networks Based on the Triple Modular Redundancy Technology (TMR)
Radial Basis Function (RBF) Neural Networks Based on the Triple Modular Redundancy Technology (TMR) Yaobin Qin qinxx143@umn.edu Supervisor: Pro.lilja Department of Electrical and Computer Engineering Abstract
More informationMERBEIN P-10 COLLEGE MATHS SCOPE & SEQUENCE
MERBEIN P-10 COLLEGE MATHS SCOPE & SEQUENCE Year Number & Algebra Measurement & Geometry Statistics & Probability P Numbers to 20 Location Data collection Counting and comparing groups Length Ordinal numbers
More informationCentral Limit Theorem Sample Means
Date Central Limit Theorem Sample Means Group Member Names: Part One Review of Types of Distributions Consider the three graphs below. Match the histograms with the distribution description. Write the
More informationThe Plan: Basic statistics: Random and pseudorandom numbers and their generation: Chapter 16.
Scientific Computing with Case Studies SIAM Press, 29 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit IV Monte Carlo Computations Dianne P. O Leary c 28 What is a Monte-Carlo method?
More informationBland-Altman Plot and Analysis
Chapter 04 Bland-Altman Plot and Analysis Introduction The Bland-Altman (mean-difference or limits of agreement) plot and analysis is used to compare two measurements of the same variable. That is, it
More informationWhy is Statistics important in Bioinformatics?
Why is Statistics important in Bioinformatics? Random processes are inherent in evolution and in sampling (data collection). Errors are often unavoidable in the data collection process. Statistics helps
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit IV Monte Carlo
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit IV Monte Carlo Computations Dianne P. O Leary c 2008 1 What is a Monte-Carlo
More informationData Analysis and Solver Plugins for KSpread USER S MANUAL. Tomasz Maliszewski
Data Analysis and Solver Plugins for KSpread USER S MANUAL Tomasz Maliszewski tmaliszewski@wp.pl Table of Content CHAPTER 1: INTRODUCTION... 3 1.1. ABOUT DATA ANALYSIS PLUGIN... 3 1.3. ABOUT SOLVER PLUGIN...
More informationPerformance Estimation and Regularization. Kasthuri Kannan, PhD. Machine Learning, Spring 2018
Performance Estimation and Regularization Kasthuri Kannan, PhD. Machine Learning, Spring 2018 Bias- Variance Tradeoff Fundamental to machine learning approaches Bias- Variance Tradeoff Error due to Bias:
More informationExcel 2010 with XLSTAT
Excel 2010 with XLSTAT J E N N I F E R LE W I S PR I E S T L E Y, PH.D. Introduction to Excel 2010 with XLSTAT The layout for Excel 2010 is slightly different from the layout for Excel 2007. However, with
More information1) Give decision trees to represent the following Boolean functions:
1) Give decision trees to represent the following Boolean functions: 1) A B 2) A [B C] 3) A XOR B 4) [A B] [C Dl Answer: 1) A B 2) A [B C] 1 3) A XOR B = (A B) ( A B) 4) [A B] [C D] 2 2) Consider the following
More informationROLE OF RANDOM NUMBERS IN SIMULATIONS OF ECONOMIC PROCESSES
562 ROLE OF RANDOM NUMBERS IN SIMULATIONS OF ECONOMIC PROCESSES Dominika Crnjac Milić, Ph.D. J.J. Strossmayer University in Osijek Faculty of Electrical Engineering Ljiljanka Kvesić, M.Sc. University of
More informationBox-Cox Transformation
Chapter 190 Box-Cox Transformation Introduction This procedure finds the appropriate Box-Cox power transformation (1964) for a single batch of data. It is used to modify the distributional shape of a set
More information10.4 Linear interpolation method Newton s method
10.4 Linear interpolation method The next best thing one can do is the linear interpolation method, also known as the double false position method. This method works similarly to the bisection method by
More information>>> SOLUTIONS <<< a) What is a system? Give a short formal definition. (Hint: what you draw a box around is not a formal definition)
Mid-Term Exam for Simulation (CIS 4930) Summer 2009 >>> SOLUTIONS
More informationpacket-switched networks. For example, multimedia applications which process
Chapter 1 Introduction There are applications which require distributed clock synchronization over packet-switched networks. For example, multimedia applications which process time-sensitive information
More informationACCURACY AND EFFICIENCY OF MONTE CARLO METHOD. Julius Goodman. Bechtel Power Corporation E. Imperial Hwy. Norwalk, CA 90650, U.S.A.
- 430 - ACCURACY AND EFFICIENCY OF MONTE CARLO METHOD Julius Goodman Bechtel Power Corporation 12400 E. Imperial Hwy. Norwalk, CA 90650, U.S.A. ABSTRACT The accuracy of Monte Carlo method of simulating
More informationISO INTERNATIONAL STANDARD. Random variate generation methods. Méthodes de génération de nombres pseudo-aléatoires. First edition
INTERNATIONAL STANDARD ISO 28640 First edition 2010-03-15 Random variate generation methods Méthodes de génération de nombres pseudo-aléatoires Reference number ISO 28640:2010(E) ISO 2010 PDF disclaimer
More informationAn Excel Add-In for Capturing Simulation Statistics
2001 Joint Statistical Meetings Atlanta, GA An Excel Add-In for Capturing Simulation Statistics David C. Trindade david.trindade@sun.com Sun Microsystems, Inc. Cupertino, CA David Meade david.meade@amd.com
More informationMathematics Year 9-11 Skills and Knowledge Checklist. Name: Class: Set : 1 Date Year 9 MEG :
Personal targets to help me achieve my grade : AFL Sheet Number 1 : Standard Form, Decimals, Fractions and Percentages Standard Form I can write a number as a product of it s prime factors I can use the
More informationSimulation. Chapter Copyright 2010 Pearson Education, Inc. Publishing as Prentice Hall
Simulation Chapter 14 14-1 Chapter Topics The Monte Carlo Process Computer Simulation with Excel Spreadsheets Simulation of a Queuing System Continuous Probability Distributions Statistical Analysis of
More informationGauss for Econometrics: Simulation
Gauss for Econometrics: Simulation R.G. Pierse 1. Introduction Simulation is a very useful tool in econometric modelling. It allows the economist to examine the properties of models and estimators when
More informationQueuing Systems. 1 Lecturer: Hawraa Sh. Modeling & Simulation- Lecture -4-21/10/2012
Queuing Systems Queuing theory establishes a powerful tool in modeling and performance analysis of many complex systems, such as computer networks, telecommunication systems, call centers, manufacturing
More informationUNIT 1: NUMBER LINES, INTERVALS, AND SETS
ALGEBRA II CURRICULUM OUTLINE 2011-2012 OVERVIEW: 1. Numbers, Lines, Intervals and Sets 2. Algebraic Manipulation: Rational Expressions and Exponents 3. Radicals and Radical Equations 4. Function Basics
More informationDesign of Experiments
Seite 1 von 1 Design of Experiments Module Overview In this module, you learn how to create design matrices, screen factors, and perform regression analysis and Monte Carlo simulation using Mathcad. Objectives
More informationSimulation studies. Patrick Breheny. September 8. Monte Carlo simulation Example: Ridge vs. Lasso vs. Subset
Simulation studies Patrick Breheny September 8 Patrick Breheny BST 764: Applied Statistical Modeling 1/17 Introduction In statistics, we are often interested in properties of various estimation and model
More informationNatasha S. Sharma, PhD
Revisiting the function evaluation problem Most functions cannot be evaluated exactly: 2 x, e x, ln x, trigonometric functions since by using a computer we are limited to the use of elementary arithmetic
More informationThe Bootstrap and Jackknife
The Bootstrap and Jackknife Summer 2017 Summer Institutes 249 Bootstrap & Jackknife Motivation In scientific research Interest often focuses upon the estimation of some unknown parameter, θ. The parameter
More informationIntroduction to Queueing Theory for Computer Scientists
Introduction to Queueing Theory for Computer Scientists Raj Jain Washington University in Saint Louis Jain@eecs.berkeley.edu or Jain@wustl.edu A Mini-Course offered at UC Berkeley, Sept-Oct 2012 These
More informationBASIC SIMULATION CONCEPTS
BASIC SIMULATION CONCEPTS INTRODUCTION Simulation is a technique that involves modeling a situation and performing experiments on that model. A model is a program that imitates a physical or business process
More informationSimulation Input Data Modeling
Introduction to Modeling and Simulation Simulation Input Data Modeling OSMAN BALCI Professor Department of Computer Science Virginia Polytechnic Institute and State University (Virginia Tech) Blacksburg,
More informationBurstiness of Intrusion Detection Process: Empirical Evidence and a Modeling Approach
Burstiness of Intrusion Detection Process: Empirical Evidence and a Modeling Approach Richard Harang, Alexander Kott Abstract We analyze sets of intrusion detection records observed on the networks of
More informationBox-Cox Transformation for Simple Linear Regression
Chapter 192 Box-Cox Transformation for Simple Linear Regression Introduction This procedure finds the appropriate Box-Cox power transformation (1964) for a dataset containing a pair of variables that are
More informationCS6375: Machine Learning Gautam Kunapuli. Mid-Term Review
Gautam Kunapuli Machine Learning Data is identically and independently distributed Goal is to learn a function that maps to Data is generated using an unknown function Learn a hypothesis that minimizes
More information3. Data Analysis and Statistics
3. Data Analysis and Statistics 3.1 Visual Analysis of Data 3.2.1 Basic Statistics Examples 3.2.2 Basic Statistical Theory 3.3 Normal Distributions 3.4 Bivariate Data 3.1 Visual Analysis of Data Visual
More informationSIMULATED ANNEALING TECHNIQUES AND OVERVIEW. Daniel Kitchener Young Scholars Program Florida State University Tallahassee, Florida, USA
SIMULATED ANNEALING TECHNIQUES AND OVERVIEW Daniel Kitchener Young Scholars Program Florida State University Tallahassee, Florida, USA 1. INTRODUCTION Simulated annealing is a global optimization algorithm
More informationExpectation-Maximization Methods in Population Analysis. Robert J. Bauer, Ph.D. ICON plc.
Expectation-Maximization Methods in Population Analysis Robert J. Bauer, Ph.D. ICON plc. 1 Objective The objective of this tutorial is to briefly describe the statistical basis of Expectation-Maximization
More informationAnalysis of Directional Beam Patterns from Firefly Optimization
Analysis of Directional Beam Patterns from Firefly Optimization Nicholas Misiunas, Charles Thompson and Kavitha Chandra Center for Advanced Computation and Telecommunications Department of Electrical and
More informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 2 December 8, 2014 Outline 1 Random Number Generation 2 Uniform Random Variables 3 Normal Random Variables 4 Correlated Random Variables Random Number Generation Monte Carlo
More informationAPPENDIX K MONTE CARLO SIMULATION
APPENDIX K MONTE CARLO SIMULATION K-1. Introduction. Monte Carlo simulation is a method of reliability analysis that should be used only when the system to be analyzed becomes too complex for use of simpler
More informationPackage simed. November 27, 2017
Version 1.0.3 Title Simulation Education Author Barry Lawson, Larry Leemis Package simed November 27, 2017 Maintainer Barry Lawson Imports graphics, grdevices, methods, stats, utils
More informationSimulation Modeling and Analysis
Simulation Modeling and Analysis FOURTH EDITION Averill M. Law President Averill M. Law & Associates, Inc. Tucson, Arizona, USA www. averill-law. com Boston Burr Ridge, IL Dubuque, IA New York San Francisco
More informationInclusion of Aleatory and Epistemic Uncertainty in Design Optimization
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Inclusion of Aleatory and Epistemic Uncertainty in Design Optimization Sirisha Rangavajhala
More informationPart I, Chapters 4 & 5. Data Tables and Data Analysis Statistics and Figures
Part I, Chapters 4 & 5 Data Tables and Data Analysis Statistics and Figures Descriptive Statistics 1 Are data points clumped? (order variable / exp. variable) Concentrated around one value? Concentrated
More information1 More configuration model
1 More configuration model In the last lecture, we explored the definition of the configuration model, a simple method for drawing networks from the ensemble, and derived some of its mathematical properties.
More informationCHAPTER VI BACK PROPAGATION ALGORITHM
6.1 Introduction CHAPTER VI BACK PROPAGATION ALGORITHM In the previous chapter, we analysed that multiple layer perceptrons are effectively applied to handle tricky problems if trained with a vastly accepted
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 WRI C225 Lecture 04 130131 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Histogram Equalization Image Filtering Linear
More informationDealing with Categorical Data Types in a Designed Experiment
Dealing with Categorical Data Types in a Designed Experiment Part II: Sizing a Designed Experiment When Using a Binary Response Best Practice Authored by: Francisco Ortiz, PhD STAT T&E COE The goal of
More informationBiostatistics & SAS programming. Kevin Zhang
Biostatistics & SAS programming Kevin Zhang February 27, 2017 Random variables and distributions 1 Data analysis Simulation study Apply existing methodologies to your collected samples, with the hope to
More informationCommon Core State Standards. August 2010
August 2010 Grade Six 6.RP: Ratios and Proportional Relationships Understand ratio concepts and use ratio reasoning to solve problems. 1. Understand the concept of a ratio and use ratio language to describe
More information