Estimating and Simulating Nonhomogeneous Poisson Processes
|
|
- Eleanore Dulcie Morris
- 6 years ago
- Views:
Transcription
1 Estimating and Simulating Nonhomogeneous Poisson Processes Larry Leemis Department of Mathematics The College of William & Mary Williamsburg, VA USA May 23, 23 Outline. Motivation 2. Probabilistic properties 3. Estimating Λ(t) from k realizations on (, S] 4. Estimating Λ(t) from overlapping realizations 5. Software 6. Conclusions Note: Portions of this work are with Brad Arkin (RST Corporation), Andy Glen (United States Military Academy), John Drew (William & Mary), and Diane Evans (Rose Hulman). Part of this work was supported by an NSF grant supporting the UMSA (Undergraduate Modeling and Simulation Analysis) REU (Research Experience for Undergraduates) at The College of William & Mary.
2 . Motivation Although easy to estimate and simulate, HPPs and renewal processes do not allow for varying rates. The use of an NHPP is often more appropriate for modeling. Example: Customer arrivals to a fast food restaurant rate Breakfast Lunch Dinner time Other examples: Cyclone arrival times in the Arctic Sea (Lee, Wilson, and Crawford, 99) Database transaction times (Lewis and Shedler, 976) Calls for on-line analysis of electrocardiograms at a hospital in Houston (Kao and Chang, 988) Respiratory cancer deaths near a steel complex in Scotland (Lawson, 988) Repairable systems: blood analyzers, fan motors, power supplies, turbines (Nelson, 988)
3 2. Probabilistic properties Notation t time N(t) number of events by time t λ(t) instantaneous arrival rate at time t (intensity function) Λ(t) = t λ(τ)dτ cumulative intensity function Property [ b a λ(τ)dτ ] n e b a λ(τ)dτ Pr (N(b) N(a) = n) = n! Property 2 E[N(t)] = Λ(t) n =,,... Property 3 (Çinlar, 975) If E, E 2,... are event times in a unit HPP then Λ (E ), Λ (E 2 ),... are event times in an NHPP with cumulative intensity function Λ(t). Λ( t ) E 5 E 4 E 3 E 2 E T T 2 T 3 T 4 T 5 t
4 3. Estimating Λ(t) from k realizations on (, S] Data t time (, S] time interval where observations are collected k number of realizations collected n, n 2,..., n k number of observations per realization t (), t (2),..., t (n) superposition of observations n = k i= n i total number of observations collected Intuitive solution (Law and Kelton, 2): partition time axis and let λ(t) be piecewise constant. Problems (a) Determining cell width Small cell width sampling variability Large cell width miss trend (b) Subjective due to arbitrary parameters rate Breakfast Lunch Dinner time
5 O O O O O O O O O Piecewise-linear nonparametric cumulative intensity function estimator in ˆΛ(t) = (n + )k + n ( ) t t (i) (n + ) k ( ) t (i) < t t (i+) t (i+) t (i) for i =,, 2,..., n. Rationale: ˆΛ(S) = n/k n k Λ( t ) 2n (n+)k n (n+)k t () t (2) S t
6 Properties Handles ties as expected Consistency lim ˆΛ(t) = Λ(t) k Confidence interval (asymptotically exact) for Λ(t) ˆΛ(t) ± z α/2 ˆΛ(t) Variate generation straightforward Input: Number of observed arrivals n Number of active realizations k Superpositioned observations t (), t (2),..., t (n+) Output: Event times T, T 2,..., T i on (, S] k i generate U i U(, ) E i log e ( U i ) while E i < n/k do begin m (n+)kei n [initialize variate counter] [generate initial random number] [generate initial exponential variate] [set m ˆΛ(t (m) ) < E i ˆΛ(t (m+) )] ( T i t (m) + [t (m+) t (m) ] (n+)kei n m [generate event time] i i + [increment variate counter] generate U i U(, ) [generate next random number] E i E i log e ( U i ) [generate next HPP event time] end )
7 4. Estimating Λ(t) from overlapping realizations Data t time (, S] time interval where observations are collected r # time intervals where the # realizations is fixed (s j, s j+ ] interval j +, j =,,..., r k j+ # realizations observed on (s j, s j+ ], j =,,..., r n j+ number of observations on (s j, s j+ ] t (), t (),..., t (n+r) superposition of observations, s, s,..., s r Note: s =, s r = S ˆΛ(t) = j q= n q k q + ( i j q= (n q + ) ) ( ) n j+ n j+ + t t(i) ( ) (n j+ + ) k j+ (n j+ + ) k j+ t(i+) t (i) t (i) < t t (i+) ; i =,, 2,..., n + r, Rationale: ˆΛ(s j+ ) = j+ s j < t s j+ ; j =,,..., r, q= n q /k q Single segment of ˆΛ(t) in the (j + )st region:, j + q = n q k q k j+ realizations n j+ observations ˆΛ(t (i+) ) ˆΛ(t (i) ) jq = n q k q s j t (i) t (i+) s j+
8 Example: Lunchwagon arrivals (Klein and Roberts, 984) Depiction of three regions for lunchwagon arrivals from : AM to 2:3 PM for k =, k 2 = 2, and k 3 = ; s =.5, s 2 = 3, and s 3 = Λ(t) t
9 An asymptotically exact ( α)% confidence interval for Λ(t) Λ(t) ˆΛ(t) < z α/2 ˆΛ(t) ˆΛ(s j ) k j+ + j q= ˆΛ(s q ) ˆΛ(s q ) Parent cumulative intensity function, nonparametric estimator, and 95% confidence bands for lunchwagon arrivals from : AM to 2:3 PM for k =, k 2 = 2, and k 3 = ; s =.5, s 2 = 3, and s 3 = 4.5. k q Λ(t) t
10 Variate Generation Input: Number of partitions r Number of active realizations k, k 2,..., k r Number of observed arrivals per partition n, n 2,..., n r Superpositioned values t (), t (),..., t (n+r) Output: Event times T, T 2,..., T i on (, S] i j MAX r q= n q /k q generate U i U(, ) E i log e ( U i ) [initialize variate counter] [initialize region counter] [set MAX to ˆΛ(S)] [generate initial random number] [generate initial exponential variate] while E i < MAX do begin while E i > j+ q= n q /k q do [update j if necessary] begin j j + [increment region counter] end m (nj+ +)k j+(e i j q= n q/k q) n j+ ( T i t (m) +[t (m+) t (m) ] (n j++)k j+ + j q= (n q + ) [set m ˆΛ(t (m) ) < E i ˆΛ(t (m+) )] ) j E i q= nq/kq n j+ ( m j q= (n q + ) ) i i + generate U i U(, ) E i E i log e ( U i ) end [generate event time] [increment variate counter] [generate next random number] [generate next HPP event time]
11 Example : Monte Carlo evaluation of the confidence interval for Λ(t) Coverages in the lunchwagon example (nominal coverage.95;, replications; k =, k 2 = 2, k 3 = ; s =, s =.5, s 2 = 3, s 3 = 4.5). Time Actual Coverage Misses High Misses Low
12 Example 2: Failure times for 2 copy machines (Zaino and Berke 992) Failure times 45 Machine Number Actuations ˆΛ(t) for each machine ˆΛ(t) Actuations
13 ˆΛ(t) for the copy machine failure times 8 ˆΛ(t) Actuations
14 Example 3: Failure times of heat pump compressors (Nelson 99) The compressors are located in five separate buildings, each under repair contract for a time span (a, b], indicated by the boldface values. The data set consists of n = 28 failure times, and yields r = 29 regions. Bldg Num of Comp Entry time, Failure Times, Exit time B , 3.3, 4.62, 4.62, 5.75, 5.75, 7.42, 7.42, 8.77, 9.27, D , 4.47, 4.47, 5.56, 5.57, 5.8, 6.3, 7.2, 7. E 458., 2,85, 4.65, 4.79, 5.85, 6.73, 7.33 H 49.,.7,.7,.34, 5.9 K 95., 2.7, 3.65, 4.4, 4.4 ˆΛ(t) for the heat pump compressor failure times..8 ˆΛ(t) t (years)
15 5. Software Civilization advances by extending the number of important operations which we can perform without thinking about them. Alfred North Whitehead (86 947) APPL (A Probability Programming Language) is a Maplebased language with data structures for discrete and continuous random variables and algorithms for their manipulation. Example : Let X, X 2,..., X be independent and identically distributed U(,) random variables. Find Typical approaches Pr Central limit theorem Simulation 4 < i= n := ; X := UniformRV(, ); Y := Convolution(X, n); CDF(Y, 6) - CDF(Y, 4); X i < 6
16 Example 2: X T riangular(, 2, 4) Y T riangular(, 2, 3) Find the distribution of V = XY. 4 y 3 xy = xy = 8 xy = 6 xy = 4 xy = 3 xy = 2 xy = x X := TriangularRV(, 2, 4); Y := TriangularRV(, 2, 3); V := Product(X, Y); which returns the probability density function of V as f V (v) = 4 3 v ln v + 2v 3 ln v < v ln 2 + 7v 3 ln v 4 ln v 5v 3 ln v 2 < v ln 2 + 7v 3 ln 2 + 2v 2 ln v v ln v 2 ln 3 2v 3 ln < v ln 2 7v 3 ln v 2 ln ln v 2v 3 ln 3 + 4v 3 ln v 4 < v ln 2 4v 3 ln v ln v + v 3 ln v + 4 ln 3 + v 3 ln 3 6 < v ln 2 + 2v 3 ln v + 4 ln 3 4 ln v + v 3 ln 3 v 3 ln v 8 < v < 2
17 Example 3: Kolmogorov Smirnov test statistic (all parameters known) Defining formula: D n = sup x F (x) F n (x), Computational formula: D n = max i=,2,...,n i n x (i), i n x (i) The cdf for the test statistic is (Birnbaum, 952) P D n < 2n + v for v 2n 2n, where for u u 2 u n. = n! 2n +v 2n v 3 2n +v 3 2n +v v... 2n 2n v 2n 2n g(u, u 2,..., u n ) du n... du 2 du g(u, u 2,..., u n ) = CASE I: n = F D (t) = Pr(D t) = t 2 2t 2 < t < t CASE II: n = 2 F D2 (t) = Pr(D 2 t) = t 4 8 ( t ) < t < 2 2( t) 2 2 < t < t
18 CASE I: n =. F(x) x CASE II: n = 2. F(x) x
19 Goal: X := KSRV(n); CASE III: n = 6. F(x) x y < y y y y y2 9 y y < y 6 48 y y y y y y < 4 32 y y y y y y y < 3 28 y F D6 (y) = y y y y y y < y 6 24 y y y y y y < 2 2 y y y y y y < 2 3 y 6 38 y y y3 5 8 y y 2 3 y < y y 5 3 y y 3 3 y y 6 y < y.
20 Example 4:. Let X, X 2,..., X be iid geometric( / 4) random variables (parameterized from ). Find the mean and variance of X (2). Y := OrderStat(GeometricRV( / 4),, 2); Mean(Y); Variance(Y); yielding µ = and σ 2 = Example 5:. Fit a power law process to the odometer failure data over (, ]: 2,942 28,489 65,56 78,254 83,639 85,63 88,43 9,89 92,36 94,78 98,23 99,9. CarFailures := [2942, 28489, 6556, 78254, 83639, 8563, 8843, 989, 9236, 9478, 9823, 999]; X := WeibullRV(lambda, kappa); hat := MLENHPP(X, CarFailures, [lambda, kappa], ); PlotEmpVsFittedCIF(X, Sample, [lambda = hat[], kappa = hat[2]],, ); ˆλ =.2637 ˆκ = 2.568
21 2 8 CIF x 6. Conclusions (a) Nonparametric estimation and simulation for NHPPs is straightforward. (b) Collecting data across overlapping intervals does not pose any significant problems. (c) Once coded, this approach requires less effort than a parametric renewal process in order to simulate the observations. (d) There may be potential for a probability package analogous to statistical packages such as SAS, SPSS, or S-Plus. (e) I am searching for situations where an exact probability calculation is needed (typically not the case in economics, OR, engineering, classical statistics; possibly the case in biology, chemistry, physics).
22 Bibliography Arkin, B., and Leemis, L., Nonparametric Estimation of the Cumulative Intensity Function for a Nonhomogeneous Poisson Process from Overlapping Intervals, Management Science, 46, 7, 2, Drew, J.H., Glen, A.G., Leemis, L.M., Computing the Cumulative Distribution Function of the Kolmogorov-Smirnov Statistic, Computational Statistics and Data Analysis, 34,, 2, 5. Evans, D.L. and Leemis, L.M., Input Modeling Using a Computer Algebra System, in Proceedings of the 2 Winter Simulation Conference, J.A. Joines, R.R. Barton, K. Kang, P.A. Fishwick, eds., Institute of Electrical and Electronics Engineers, Orlando, Florida, 2, Glen, A.G., Leemis, L.M., and Drew, J.H., A Generalized Univariate Change-of-Variable Transformation Technique, INFORMS Journal on Computing, 9, 3, 997, Glen, A.G., Evans, D.L., and Leemis, L.M., APPL: A Probability Programming Language, The American Statistician, 55, 2, 2, Glen, A.G., Leemis, L.M., and Drew, J.H., Computing the Distribution of the Product of Two Continuous Random Variables, forthcoming, Computational Statistics and Data Analysis. Leemis, L.M., Nonparametric Estimation of the Intensity Function for a Nonhomogeneous Poisson Process, Management Science, 37, 7, 99,
Nonparametric estimation and variate generation for a nonhomogeneous Poisson process from event count data
IIE Transactions (2004) 36, 1155 1160 Copyright C IIE ISSN: 0740-817X print / 1545-8830 online DOI: 10.1080/07408170490507693 TECHNICAL NOTE Nonparametric estimation and variate generation for a nonhomogeneous
More informationEstimation for Nonhomogeneous Poisson Processes from Aggregated Data
Estimation for Nonhomogeneous Poisson Processes from Aggregated Data Shane G. Henderson School of Operations Research and Industrial Engineering, Cornell University, Ithaca, NY 14853. November 22, 2002
More informationA TEST OF RANDOMNESS BASED ON THE DISTANCE BETWEEN CONSECUTIVE RANDOM NUMBER PAIRS. Matthew J. Duggan John H. Drew Lawrence M.
Proceedings of the 2005 Winter Simulation Conference M. E. Kuhl, N. M. Steiger, F. B. Armstrong, and J. A. Joines, eds. A TEST OF RANDOMNESS BASED ON THE DISTANCE BETWEEN CONSECUTIVE RANDOM NUMBER PAIRS
More informationSection 7.5: Nonstationary Poisson Processes
Section 7.5: Nonstationary Poisson Processes Discrete-Event Simulation: A First Course c 2006 Pearson Ed., Inc. 0-13-142917-5 Discrete-Event Simulation: A First Course Section 7.5: Nonstationary Poisson
More informationReliability Growth in JMP 10
Reliability Growth in JMP 10 Presented at Discovery Summit 2012 September 13, 2012 Marie Gaudard and Leo Wright Purpose of Talk The goal of this talk is to provide a brief introduction to: The area of
More informationMahdi karbasian* & Zoubi Ibrahim
International Journal of Industrial Engineering & Production Research (010) pp. 105-110 September 010, Volume 1, Number International Journal of Industrial Engineering & Production Research ISSN: 008-4889
More informationAn Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems
An Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems M.P. Kaminskiy and V.V. Krivtsov Abstract This paper introduces a simple index that helps to assess the degree of aging
More informationIntroduction to repairable systems STK4400 Spring 2011
Introduction to repairable systems STK4400 Spring 2011 Bo Lindqvist http://www.math.ntnu.no/ bo/ bo@math.ntnu.no Bo Lindqvist Introduction to repairable systems Definition of repairable system Ascher and
More informationRepairable Systems Reliability Trend Tests and Evaluation
Repairable Systems Reliability Trend Tests and Evaluation Peng Wang, Ph.D., United Technologies Research Center David W. Coit, Ph.D., Rutgers University Keywords: repairable system, reliability trend test,
More informationAN INTEGRAL MEASURE OF AGING/REJUVENATION FOR REPAIRABLE AND NON REPAIRABLE SYSTEMS
R&RAA # 1 (Vol.1) 8, March AN INEGRAL MEASURE OF AGING/REJUVENAION FOR REPAIRABLE AND NON REPAIRABLE SYSEMS M.P. Kaminskiy and V.V. Krivtsov Abstract his paper introduces a simple index that helps to assess
More informationIntroduction to Modeling and Generating Probabilistic Input Processes for Simulation
Introduction to Modeling and Generating Probabilistic Input Processes for Simulation Michael E. Kuhl RIT Natalie M. Steiger University of Maine Emily K. Lada SAS Institute Inc. Mary Ann Wagner SAIC James
More informationISSN (ONLINE): , ISSN (PRINT):
www.ijemr.net ISSN ONLINE): 2250-0758, ISSN PRINT): 2394-6962 Volume-8, Issue-2, April 2018 International Journal of Engineering and Management Research Page Number: 232-236 DOI: doi.org/10.31033/ijemr.v8i02.11652
More informationTime-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation
Time-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation H. Zhang, E. Cutright & T. Giras Center of Rail Safety-Critical Excellence, University of Virginia,
More informationConfidence Intervals for Reliability Growth Models with Small Sample Sizes. Reliability growth models, Order statistics, Confidence intervals
Confidence Intervals for Reliability Growth Models with Small Sample Sizes John Quigley, Lesley Walls University of Strathclyde, Glasgow, Scotland Key Words Reliability growth models, Order statistics,
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 informationLEAST SQUARES ESTIMATION OF NONHOMOGENEOUS POISSON PROCESSES
Proceedings of the 1998 Winter Simulation Conference D.J. Medeiros, E.F. Watson, J.S. Carson and M.S. Manivannan, eds. LEAST SQUARES ESTIMATION OF NONHOMOGENEOUS POISSON PROCESSES Michael E. Kuhl Department
More informationComputer Simulation of Repairable Processes
SEMATECH 1996 Applied Reliability Tools Workshop (ARTWORK IX) Santa Fe Computer Simulation of Repairable Processes Dave Trindade, Ph.D. Senior AMD Fellow Applied Statistics Introduction Computer simulation!
More informationProceedings of the 2014 Winter Simulation Conference A. Tolk, S. D. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds.
Proceedings of the 2014 Winter Simulation Conference A. Tolk, S. D. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds. A Continuous Piecewise-Linear NHPP Intensity Function Estimator David
More informationBAYESIAN MODELING OF DYNAMIC SOFTWARE GROWTH CURVE MODELS
BAYESIAN MODELING OF DYNAMIC SOFTWARE GROWTH CURVE MODELS Zhaohui Liu, Nalini Ravishanker, University of Connecticut Bonnie K. Ray, IBM Watson Research Center Department of Mathematical Sciences, IBM Watson
More informationStat 710: Mathematical Statistics Lecture 31
Stat 710: Mathematical Statistics Lecture 31 Jun Shao Department of Statistics University of Wisconsin Madison, WI 53706, USA Jun Shao (UW-Madison) Stat 710, Lecture 31 April 13, 2009 1 / 13 Lecture 31:
More informationIndependent Events. Two events are independent if knowing that one occurs does not change the probability of the other occurring
Independent Events Two events are independent if knowing that one occurs does not change the probability of the other occurring Conditional probability is denoted P(A B), which is defined to be: P(A and
More informationSolution: The process is a compound Poisson Process with E[N (t)] = λt/p by Wald's equation.
Solutions Stochastic Processes and Simulation II, May 18, 217 Problem 1: Poisson Processes Let {N(t), t } be a homogeneous Poisson Process on (, ) with rate λ. Let {S i, i = 1, 2, } be the points of the
More informationThe exponential distribution and the Poisson process
The exponential distribution and the Poisson process 1-1 Exponential Distribution: Basic Facts PDF f(t) = { λe λt, t 0 0, t < 0 CDF Pr{T t) = 0 t λe λu du = 1 e λt (t 0) Mean E[T] = 1 λ Variance Var[T]
More informationNotes largely based on Statistical Methods for Reliability Data by W.Q. Meeker and L. A. Escobar, Wiley, 1998 and on their class notes.
Unit 2: Models, Censoring, and Likelihood for Failure-Time Data Notes largely based on Statistical Methods for Reliability Data by W.Q. Meeker and L. A. Escobar, Wiley, 1998 and on their class notes. Ramón
More informationSTAT331 Lebesgue-Stieltjes Integrals, Martingales, Counting Processes
STAT331 Lebesgue-Stieltjes Integrals, Martingales, Counting Processes This section introduces Lebesgue-Stieltjes integrals, and defines two important stochastic processes: a martingale process and a counting
More informationNorthwestern University Department of Electrical Engineering and Computer Science
Northwestern University Department of Electrical Engineering and Computer Science EECS 454: Modeling and Analysis of Communication Networks Spring 2008 Probability Review As discussed in Lecture 1, probability
More informationRandom variables. DS GA 1002 Probability and Statistics for Data Science.
Random variables DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Motivation Random variables model numerical quantities
More informationStat 516, Homework 1
Stat 516, Homework 1 Due date: October 7 1. Consider an urn with n distinct balls numbered 1,..., n. We sample balls from the urn with replacement. Let N be the number of draws until we encounter a ball
More informationIEOR 3106: Second Midterm Exam, Chapters 5-6, November 7, 2013
IEOR 316: Second Midterm Exam, Chapters 5-6, November 7, 13 SOLUTIONS Honor Code: Students are expected to behave honorably, following the accepted code of academic honesty. You may keep the exam itself.
More informationMartín Gastón Teresa León Fermín Mallor
Proceedings of the 2008 Winter Simulation Conference S. J. Mason, R. R. Hill, L. Mönch, O. Rose, T. Jefferson, J. W. Fowler eds. FUNCTIONAL DATA ANALYSIS FOR NON HOMOGENEOUS POISSON PROCESSES Martín Gastón
More informationMAS1302 Computational Probability and Statistics
MAS1302 Computational Probability and Statistics April 23, 2008 3. Simulating continuous random behaviour 3.1 The Continuous Uniform U(0,1) Distribution We have already used this random variable a great
More informationMODELING AND ANALYZING TRANSIENT MILITARY AIR TRAFFIC CONTROL
Proceedings of the 2010 Winter Simulation Conference B Johansson, S Jain, J Montoya-Torres, J Hugan, and E Yücesan, eds MODELING AND ANALYZING TRANSIENT MILITARY AIR TRAFFIC CONTROL William Kaczynski United
More informationSimulating Multivariate Nonhomogeneous Poisson Processes Using Projections
Simulating Multivariate Nonhomogeneous Poisson Processes Using Projections EVAN A. SALTZMAN, RAND Corporation JOHN H. DREW and LAWRENCE M. LEEMIS, The College of William & Mary SHANE G. HENDERSON, Cornell
More informationIE 581 Introduction to Stochastic Simulation
1. List criteria for choosing the majorizing density r (x) when creating an acceptance/rejection random-variate generator for a specified density function f (x). 2. Suppose the rate function of a nonhomogeneous
More informationEconomics 583: Econometric Theory I A Primer on Asymptotics
Economics 583: Econometric Theory I A Primer on Asymptotics Eric Zivot January 14, 2013 The two main concepts in asymptotic theory that we will use are Consistency Asymptotic Normality Intuition consistency:
More informationStochastic Renewal Processes in Structural Reliability Analysis:
Stochastic Renewal Processes in Structural Reliability Analysis: An Overview of Models and Applications Professor and Industrial Research Chair Department of Civil and Environmental Engineering University
More informationPoisson Processes. Particles arriving over time at a particle detector. Several ways to describe most common model.
Poisson Processes Particles arriving over time at a particle detector. Several ways to describe most common model. Approach 1: a) numbers of particles arriving in an interval has Poisson distribution,
More information3 Continuous Random Variables
Jinguo Lian Math437 Notes January 15, 016 3 Continuous Random Variables Remember that discrete random variables can take only a countable number of possible values. On the other hand, a continuous random
More informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More informationVariability within multi-component systems. Bayesian inference in probabilistic risk assessment The current state of the art
PhD seminar series Probabilistics in Engineering : g Bayesian networks and Bayesian hierarchical analysis in engeering g Conducted by Prof. Dr. Maes, Prof. Dr. Faber and Dr. Nishijima Variability within
More informationS1. Proofs of the Key Properties of CIATA-Ph
Article submitted to INFORMS Journal on Computing; manuscript no. (Please, provide the mansucript number!) S-1 Online Supplement to Modeling and Simulation of Nonstationary and Non-Poisson Processes by
More informationEstimation for generalized half logistic distribution based on records
Journal of the Korean Data & Information Science Society 202, 236, 249 257 http://dx.doi.org/0.7465/jkdi.202.23.6.249 한국데이터정보과학회지 Estimation for generalized half logistic distribution based on records
More informationSimulation and Calibration of a Fully Bayesian Marked Multidimensional Hawkes Process with Dissimilar Decays
Simulation and Calibration of a Fully Bayesian Marked Multidimensional Hawkes Process with Dissimilar Decays Kar Wai Lim, Young Lee, Leif Hanlen, Hongbiao Zhao Australian National University Data61/CSIRO
More informationModeling and Simulation of Nonstationary Non-Poisson Arrival Processes
Submitted to INFORMS Journal on Computing manuscript (Please, provide the mansucript number!) Authors are encouraged to submit new papers to INFORMS journals by means of a style file template, which includes
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 informationLecture 2: CDF and EDF
STAT 425: Introduction to Nonparametric Statistics Winter 2018 Instructor: Yen-Chi Chen Lecture 2: CDF and EDF 2.1 CDF: Cumulative Distribution Function For a random variable X, its CDF F () contains all
More informationOptimal maintenance for minimal and imperfect repair models
Optimal maintenance for minimal and imperfect repair models Gustavo L. Gilardoni Universidade de Brasília AMMSI Workshop Grenoble, January 2016 Belo Horizonte Enrico Colosimo Marta Freitas Rodrigo Padilha
More informationChoosing Arrival Process Models for Service Systems: Tests of a Nonhomogeneous Poisson Process
Choosing Arrival Process Models for Service Systems: Tests of a Nonhomogeneous Poisson Process Song-Hee Kim and Ward Whitt Industrial Engineering and Operations Research Columbia University New York, NY,
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 informationISE/OR 762 Stochastic Simulation Techniques
ISE/OR 762 Stochastic Simulation Techniques Topic 0: Introduction to Discrete Event Simulation Yunan Liu Department of Industrial and Systems Engineering NC State University January 9, 2018 Yunan Liu (NC
More informationComputer simulation on homogeneity testing for weighted data sets used in HEP
Computer simulation on homogeneity testing for weighted data sets used in HEP Petr Bouř and Václav Kůs Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University
More informationIntroduction to Reliability Theory (part 2)
Introduction to Reliability Theory (part 2) Frank Coolen UTOPIAE Training School II, Durham University 3 July 2018 (UTOPIAE) Introduction to Reliability Theory 1 / 21 Outline Statistical issues Software
More informationI I FINAL, 01 Jun 8.4 to 31 May TITLE AND SUBTITLE 5 * _- N, '. ', -;
R AD-A237 850 E........ I N 11111IIIII U 1 1I!til II II... 1. AGENCY USE ONLY Leave 'VanK) I2. REPORT DATE 3 REPORT TYPE AND " - - I I FINAL, 01 Jun 8.4 to 31 May 88 4. TITLE AND SUBTITLE 5 * _- N, '.
More informationContinuous-time Markov Chains
Continuous-time Markov Chains Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ October 23, 2017
More informationA STATISTICAL TEST FOR MONOTONIC AND NON-MONOTONIC TREND IN REPAIRABLE SYSTEMS
A STATISTICAL TEST FOR MONOTONIC AND NON-MONOTONIC TREND IN REPAIRABLE SYSTEMS Jan Terje Kvaløy Department of Mathematics and Science, Stavanger University College, P.O. Box 2557 Ullandhaug, N-491 Stavanger,
More information1 Poisson processes, and Compound (batch) Poisson processes
Copyright c 2007 by Karl Sigman 1 Poisson processes, and Compound (batch) Poisson processes 1.1 Point Processes Definition 1.1 A simple point process ψ = {t n : n 1} is a sequence of strictly increasing
More informationOnline Supplement to Are Call Center and Hospital Arrivals Well Modeled by Nonhomogeneous Poisson Processes?
Online Supplement to Are Call Center and Hospital Arrivals Well Modeled by Nonhomogeneous Poisson Processes? Song-Hee Kim and Ward Whitt Industrial Engineering and Operations Research Columbia University
More informationRecognition Performance from SAR Imagery Subject to System Resource Constraints
Recognition Performance from SAR Imagery Subject to System Resource Constraints Michael D. DeVore Advisor: Joseph A. O SullivanO Washington University in St. Louis Electronic Systems and Signals Research
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 informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationABSTRACT. LIU, RAN. Modeling and Simulation of Nonstationary Non-Poisson Processes. (Under the direction of Dr. James R. Wilson.)
ABSTRACT LIU, RAN. Modeling and Simulation of Nonstationary Non-Poisson Processes. (Under the direction of Dr. James R. Wilson.) We develop and experimentally evaluate various methods for modeling and
More informationST745: Survival Analysis: Nonparametric methods
ST745: Survival Analysis: Nonparametric methods Eric B. Laber Department of Statistics, North Carolina State University February 5, 2015 The KM estimator is used ubiquitously in medical studies to estimate
More informationXLVII SIMPÓSIO BRASILEIRO DE PESQUISA OPERACIONAL
Static and Dynamics Maintenance Policies for Repairable Systems under Imperfect Repair Models Maria Luíza Guerra de Toledo Escola Nacional de Ciências Estatísticas - IBGE Rua André Cavalcanti, 106 - Santa
More information3. Poisson Processes (12/09/12, see Adult and Baby Ross)
3. Poisson Processes (12/09/12, see Adult and Baby Ross) Exponential Distribution Poisson Processes Poisson and Exponential Relationship Generalizations 1 Exponential Distribution Definition: The continuous
More informationStatistical Methods for Astronomy
Statistical Methods for Astronomy If your experiment needs statistics, you ought to have done a better experiment. -Ernest Rutherford Lecture 1 Lecture 2 Why do we need statistics? Definitions Statistical
More informationMultivariate Distributions
Copyright Cosma Rohilla Shalizi; do not distribute without permission updates at http://www.stat.cmu.edu/~cshalizi/adafaepov/ Appendix E Multivariate Distributions E.1 Review of Definitions Let s review
More informationMaster of Science in Statistics A Proposal
1 Master of Science in Statistics A Proposal Rationale of the Program In order to cope up with the emerging complexity on the solutions of realistic problems involving several phenomena of nature it is
More informationSlides 5: Random Number Extensions
Slides 5: Random Number Extensions We previously considered a few examples of simulating real processes. In order to mimic real randomness of events such as arrival times we considered the use of random
More informationPrerequisite: STATS 7 or STATS 8 or AP90 or (STATS 120A and STATS 120B and STATS 120C). AP90 with a minimum score of 3
University of California, Irvine 2017-2018 1 Statistics (STATS) Courses STATS 5. Seminar in Data Science. 1 Unit. An introduction to the field of Data Science; intended for entering freshman and transfers.
More informationThis paper investigates the impact of dependence among successive service times on the transient and
MANUFACTURING & SERVICE OPERATIONS MANAGEMENT Vol. 14, No. 2, Spring 212, pp. 262 278 ISSN 1523-4614 (print) ISSN 1526-5498 (online) http://dx.doi.org/1.1287/msom.111.363 212 INFORMS The Impact of Dependent
More informationRENEWAL PROCESSES AND POISSON PROCESSES
1 RENEWAL PROCESSES AND POISSON PROCESSES Andrea Bobbio Anno Accademico 1997-1998 Renewal and Poisson Processes 2 Renewal Processes A renewal process is a point process characterized by the fact that the
More informationChapter 5. Statistical Models in Simulations 5.1. Prof. Dr. Mesut Güneş Ch. 5 Statistical Models in Simulations
Chapter 5 Statistical Models in Simulations 5.1 Contents Basic Probability Theory Concepts Discrete Distributions Continuous Distributions Poisson Process Empirical Distributions Useful Statistical Models
More informationFinite-Horizon Statistics for Markov chains
Analyzing FSDT Markov chains Friday, September 30, 2011 2:03 PM Simulating FSDT Markov chains, as we have said is very straightforward, either by using probability transition matrix or stochastic update
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 informationA Simulation Study on Confidence Interval Procedures of Some Mean Cumulative Function Estimators
Statistics Preprints Statistics -00 A Simulation Study on Confidence Interval Procedures of Some Mean Cumulative Function Estimators Jianying Zuo Iowa State University, jiyizu@iastate.edu William Q. Meeker
More informationModelling and simulating non-stationary arrival processes to facilitate analysis
Journal of Simulation (211) 5, 3 8 r 211 Operational Research Society Ltd. All rights reserved. 1747-7778/11 www.palgrave-journals.com/jos/ Modelling and simulating non-stationary arrival processes to
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
More informationUniform random numbers generators
Uniform random numbers generators Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/tlt-2707/ OUTLINE: The need for random numbers; Basic steps in generation; Uniformly
More informationGradient-based Adaptive Stochastic Search
1 / 41 Gradient-based Adaptive Stochastic Search Enlu Zhou H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology November 5, 2014 Outline 2 / 41 1 Introduction
More informationNotation Precedence Diagram
Notation Precedence Diagram xix xxi CHAPTER 1 Introduction 1 1.1. Systems, Models, and Simulation 1 1.2. Verification, Approximation, and Validation 8 1.2.1. Verifying a Program 9 1.2.2. Approximation
More informationExamination paper for TMA4275 Lifetime Analysis
Department of Mathematical Sciences Examination paper for TMA4275 Lifetime Analysis Academic contact during examination: Ioannis Vardaxis Phone: 95 36 00 26 Examination date: Saturday May 30 2015 Examination
More informationMultistate models and recurrent event models
Multistate models Multistate models and recurrent event models Patrick Breheny December 10 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/22 Introduction Multistate models In this final lecture,
More informationPermutation tests. Patrick Breheny. September 25. Conditioning Nonparametric null hypotheses Permutation testing
Permutation tests Patrick Breheny September 25 Patrick Breheny STA 621: Nonparametric Statistics 1/16 The conditioning idea In many hypothesis testing problems, information can be divided into portions
More informationEstimating a parametric lifetime distribution from superimposed renewal process data
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2013 Estimating a parametric lifetime distribution from superimposed renewal process data Ye Tian Iowa State
More informationDepartment of Mathematics
Department of Mathematics Ma 3/103 KC Border Introduction to Probability and Statistics Winter 2017 Supplement 2: Review Your Distributions Relevant textbook passages: Pitman [10]: pages 476 487. Larsen
More informationerrors every 1 hour unless he falls asleep, in which case he just reports the total errors
I. First Definition of a Poisson Process A. Definition: Poisson process A Poisson Process {X(t), t 0} with intensity λ > 0 is a counting process with the following properties. INDEPENDENT INCREMENTS. For
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationAn algorithm for computing the PDF of order statistics drawn from discrete parent populations is presented,
INFORMS Journal on Computing Vol. 18, No. 1, Winter 2006, pp. 19 30 issn 1091-9856 eissn 1526-5528 06 1801 0019 informs doi 10.1287/ijoc.100.0105 2006 INFORMS The Distribution of Order Statistics for Discrete
More informationDepartment of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim
Tests for trend in more than one repairable system. Jan Terje Kvaly Department of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim ABSTRACT: If failure time data from several
More informationProceedings of the 2012 Winter Simulation Conference C. Laroque, J. Himmelspach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds.
Proceedings of the 2012 Winter Simulation Conference C. Laroque, J. Himmelspach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds. TUTORIAL: TEACHING AN ADVANCED SIMULATION TOPIC Shane G. Henderson School
More informationExpected Values, Exponential and Gamma Distributions
Expected Values, Exponential and Gamma Distributions Sections 5.2 & 5.4 Cathy Poliak, Ph.D. cathy@math.uh.edu Office in Fleming 11c Department of Mathematics University of Houston Lecture 13-3339 Cathy
More informationA nonparametric monotone maximum likelihood estimator of time trend for repairable systems data
A nonparametric monotone maximum likelihood estimator of time trend for repairable systems data Knut Heggland 1 and Bo H Lindqvist Department of Mathematical Sciences Norwegian University of Science and
More informationSTATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL. A Thesis. Presented to the. Faculty of. San Diego State University
STATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL A Thesis Presented to the Faculty of San Diego State University In Partial Fulfillment of the Requirements for the Degree
More informationStock Sampling with Interval-Censored Elapsed Duration: A Monte Carlo Analysis
Stock Sampling with Interval-Censored Elapsed Duration: A Monte Carlo Analysis Michael P. Babington and Javier Cano-Urbina August 31, 2018 Abstract Duration data obtained from a given stock of individuals
More informationAnalytical Bootstrap Methods for Censored Data
JOURNAL OF APPLIED MATHEMATICS AND DECISION SCIENCES, 6(2, 129 141 Copyright c 2002, Lawrence Erlbaum Associates, Inc. Analytical Bootstrap Methods for Censored Data ALAN D. HUTSON Division of Biostatistics,
More information2 Continuous Random Variables and their Distributions
Name: Discussion-5 1 Introduction - Continuous random variables have a range in the form of Interval on the real number line. Union of non-overlapping intervals on real line. - We also know that for any
More informationSlides 8: Statistical Models in Simulation
Slides 8: Statistical Models in Simulation Purpose and Overview The world the model-builder sees is probabilistic rather than deterministic: Some statistical model might well describe the variations. An
More informationConditioning Nonparametric null hypotheses Permutation testing. Permutation tests. Patrick Breheny. October 5. STA 621: Nonparametric Statistics
Permutation tests October 5 The conditioning idea In most hypothesis testing problems, information can be divided into portions that pertain to the hypothesis and portions that do not The usual approach
More informationPoisson Processes. Stochastic Processes. Feb UC3M
Poisson Processes Stochastic Processes UC3M Feb. 2012 Exponential random variables A random variable T has exponential distribution with rate λ > 0 if its probability density function can been written
More informationBirth-Death Processes
Birth-Death Processes Birth-Death Processes: Transient Solution Poisson Process: State Distribution Poisson Process: Inter-arrival Times Dr Conor McArdle EE414 - Birth-Death Processes 1/17 Birth-Death
More information