A comparison of time-based maintenance and condition-based maintenance
|
|
- Shauna Wood
- 5 years ago
- Views:
Transcription
1 A comparison of time-based maintenance and condition-based maintenance Bram de Jonge Ruud Teunter Warse Klingenberg Tiedo Tinga University of Groningen The Netherlands COPE congres 2013 November 1, 2013
2 Introduction
3 Introduction We consider the following setting:
4 Introduction We consider the following setting: A single unit.
5 Introduction We consider the following setting: A single unit. Gradual deterioration, condition monitoring is possible.
6 Introduction We consider the following setting: A single unit. Gradual deterioration, condition monitoring is possible. The cost of a breakdown is 1, after which a new unit is putting into use.
7 Introduction We consider the following setting: A single unit. Gradual deterioration, condition monitoring is possible. The cost of a breakdown is 1, after which a new unit is putting into use. Preventive maintenance can be performed at cost c < 1 and makes the unit as good as new.
8 Maintenance policies
9 Maintenance policies Time-based maintenance (TBM)
10 Maintenance policies Time-based maintenance (TBM) Easy to implement, no condition monitoring needed.
11 Maintenance policies Time-based maintenance (TBM) Easy to implement, no condition monitoring needed. Decision variable: T (maintenance age).
12 Maintenance policies Time-based maintenance (TBM) Easy to implement, no condition monitoring needed. Decision variable: T (maintenance age). Maintenance is performed when the unit reaches age T.
13 Maintenance policies Time-based maintenance (TBM) Easy to implement, no condition monitoring needed. Decision variable: T (maintenance age). Maintenance is performed when the unit reaches age T. Condition-based maintenance (CBM)
14 Maintenance policies Time-based maintenance (TBM) Easy to implement, no condition monitoring needed. Decision variable: T (maintenance age). Maintenance is performed when the unit reaches age T. Condition-based maintenance (CBM) Condition monitoring required, more appropriate scheduling.
15 Maintenance policies Time-based maintenance (TBM) Easy to implement, no condition monitoring needed. Decision variable: T (maintenance age). Maintenance is performed when the unit reaches age T. Condition-based maintenance (CBM) Condition monitoring required, more appropriate scheduling. Decision variable: M (threshold deterioration level).
16 Maintenance policies Time-based maintenance (TBM) Easy to implement, no condition monitoring needed. Decision variable: T (maintenance age). Maintenance is performed when the unit reaches age T. Condition-based maintenance (CBM) Condition monitoring required, more appropriate scheduling. Decision variable: M (threshold deterioration level). Maintenance is performed when the deterioration level M is reached.
17 Gamma deterioration process
18 Gamma deterioration process The state of the unit is assumed to deteriorate gradually according to a stationary/homogeneous gamma process X (t).
19 Gamma deterioration process The state of the unit is assumed to deteriorate gradually according to a stationary/homogeneous gamma process X (t). Properties:
20 Gamma deterioration process The state of the unit is assumed to deteriorate gradually according to a stationary/homogeneous gamma process X (t). Properties: Continuous-time process.
21 Gamma deterioration process The state of the unit is assumed to deteriorate gradually according to a stationary/homogeneous gamma process X (t). Properties: Continuous-time process. X (0) = 0 with probability 1.
22 Gamma deterioration process The state of the unit is assumed to deteriorate gradually according to a stationary/homogeneous gamma process X (t). Properties: Continuous-time process. X (0) = 0 with probability 1. X (τ) X (t) f a(τ t),b for τ > t 0 (f is the gamma density function).
23 Gamma deterioration process The state of the unit is assumed to deteriorate gradually according to a stationary/homogeneous gamma process X (t). Properties: Continuous-time process. X (0) = 0 with probability 1. X (τ) X (t) f a(τ t),b for τ > t 0 (f is the gamma density function). X (t) has independent increments.
24 Gamma deterioration process The state of the unit is assumed to deteriorate gradually according to a stationary/homogeneous gamma process X (t). Properties: Continuous-time process. X (0) = 0 with probability 1. X (τ) X (t) f a(τ t),b for τ > t 0 (f is the gamma density function). X (t) has independent increments. Jump process with infinitely many jumps in any time interval.
25 Gamma deterioration process
26 Gamma deterioration process The stationary gamma process has two parameters:
27 Gamma deterioration process The stationary gamma process has two parameters: Shape parameter a.
28 Gamma deterioration process The stationary gamma process has two parameters: Shape parameter a. Scale parameter b.
29 Gamma deterioration process The stationary gamma process has two parameters: Shape parameter a. Scale parameter b. Properties:
30 Gamma deterioration process The stationary gamma process has two parameters: Shape parameter a. Scale parameter b. Properties: E(X (t)) = abt.
31 Gamma deterioration process The stationary gamma process has two parameters: Shape parameter a. Scale parameter b. Properties: E(X (t)) = abt. Var(X (t)) = ab 2 t.
32 Gamma deterioration process The stationary gamma process has two parameters: Shape parameter a. Scale parameter b. Properties: E(X (t)) = abt. Var(X (t)) = ab 2 t. A breakdown is assumed to occur if the level of deterioration exceeds 1. We choose a and b such that the mean time to failure equals 1 (this does not imply ab = 1).
33 Gamma deterioration process X(t) t
34 Gamma deterioration process X(t) 1 Low a, high b t
35 Gamma deterioration process X(t) 1 High a, low b Low a, high b t
36 Gamma deterioration process X(t) 1 In between High a, low b Low a, high b t
37 Research methodology
38 Research methodology Simulation:
39 Research methodology Simulation: Exact calculations are practically impossible.
40 Research methodology Simulation: Exact calculations are practically impossible. If the number of iterations is sufficiently high, the outcomes are very accurate.
41 Research methodology Simulation: Exact calculations are practically impossible. If the number of iterations is sufficiently high, the outcomes are very accurate. Simulation details: Time step: = Number of iterations: n = 100,000.
42 Base case
43 Base case Gamma deterioration process: a = 5, b
44 Base case Gamma deterioration process: a = 5, b Relative cost preventive maintenance: c = 0.2.
45 Cost analysis
46 Cost analysis Cost saving TBM Cost saving CBM Failure cost = 1 Cost due to uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
47 Cost analysis Cost T or M
48 Cost analysis Cost saving TBM Cost saving CBM Failure cost = 1 Cost due to uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
49 Cost analysis Cost T or M
50 Cost analysis Cost 1 TBM T or M
51 Cost analysis Cost saving TBM Cost saving CBM Failure cost = 1 Cost due to uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
52 Cost analysis Cost 1 TBM T or M
53 Cost analysis Cost 1 TBM CBM T or M
54 Cost analysis Cost saving TBM Failure cost = 1 Cost due to uncertain failure level Cost saving CBM Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
55 Cost analysis Cost saving TBM Failure cost = 1 Cost due saving to CBM uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
56 Effect of deterioration process Cost σ
57 Effect of deterioration process Cost 1 TBM σ
58 Effect of deterioration process Cost 1 TBM CBM σ
59 Cost benefit CBM compared with TBM Cost σ
60 Cost analysis Cost saving TBM Failure cost = 1 Cost due saving to CBM uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
61 Planning time
62 Planning time Preventive maintenance needs to be planned s time units in advance.
63 Planning time Preventive maintenance needs to be planned s time units in advance. The cost c of preventive maintenance does not have to be paid if failure occurs between the moment that maintenance is planned and the moment that maintenance would have been performed.
64 Planning time Preventive maintenance needs to be planned s time units in advance. The cost c of preventive maintenance does not have to be paid if failure occurs between the moment that maintenance is planned and the moment that maintenance would have been performed. Under the CBM policy, maintenance is performed at time t(m) + s, with t(m) the time at which deterioration level M is reached.
65 Planning time Preventive maintenance needs to be planned s time units in advance. The cost c of preventive maintenance does not have to be paid if failure occurs between the moment that maintenance is planned and the moment that maintenance would have been performed. Under the CBM policy, maintenance is performed at time t(m) + s, with t(m) the time at which deterioration level M is reached. The planning time s has no influence on the TBM policy.
66 Planning time Cost 1 TBM CBM (s=0) T or M
67 Planning time Cost 1 TBM CBM (s=0.1) CBM (s=0) T or M
68 Planning time Cost s
69 Planning time Cost 0.6 TBM s
70 Planning time Cost 0.6 TBM CBM s
71 Cost analysis Cost saving TBM Failure cost = 1 Cost due saving to CBM uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
72 Cost analysis Cost saving TBM Failure cost = 1 Cost saving CBM Cost due to uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
73 Imperfect condition monitoring
74 Imperfect condition monitoring The difference between the true level of deterioration and the observed level of deterioration is modelled by σ p W (t).
75 Imperfect condition monitoring The difference between the true level of deterioration and the observed level of deterioration is modelled by σ p W (t). W (t) is a standard Brownian motion.
76 Imperfect condition monitoring Prognostic error t Low σ p 0.4
77 Imperfect condition monitoring Prognostic error 0.4 High σ p t 0.4
78 Imperfect condition monitoring Cost σ p
79 Imperfect condition monitoring Cost 0.6 TBM σ p
80 Imperfect condition monitoring Cost 0.6 TBM CBM σ p
81 Cost analysis Cost saving TBM Failure cost = 1 Cost saving CBM Cost due to uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
82 Cost analysis Cost saving TBM Failure cost = 1 Cost saving CBM Cost due to uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
83 Uncertainty deterioration level failure
84 Uncertainty deterioration level failure The deterioration level at which failure occurs is normally distributed with mean 1 and standard deviation σ f.
85 Uncertainty deterioration level failure The deterioration level at which failure occurs is normally distributed with mean 1 and standard deviation σ f. The normal distribution is left truncated at 0 and right truncated at 2.
86 Uncertainty deterioration level failure Cost σ f
87 Uncertainty deterioration level failure Cost 0.8 TBM σ f
88 Uncertainty deterioration level failure Cost 0.8 TBM CBM σ f
89 Cost analysis Cost saving TBM Failure cost = 1 Cost saving CBM Cost due to uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
90 Cost analysis Cost saving TBM Cost saving CBM Failure cost = 1 Cost due to uncertain failure level Cost due to imperfect condition monitoring Cost due to maintenance planning Cost due to noncontinuous deterioration information Cost due to deterioration process with jumps Minimal cost (PM just before failure) = c
91 Questions?
University of Groningen. Maintenance Optimization based on Mathematical Modeling de Jonge, Bram
University of Groningen Maintenance Optimization based on Mathematical Modeling de Jonge, Bram IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from
More informationFinancial Engineering. Craig Pirrong Spring, 2006
Financial Engineering Craig Pirrong Spring, 2006 March 8, 2006 1 Levy Processes Geometric Brownian Motion is very tractible, and captures some salient features of speculative price dynamics, but it is
More informationUniversity of Groningen. Maintenance Optimization based on Mathematical Modeling de Jonge, Bram
University of Groningen Maintenance Optimization based on Mathematical Modeling de Jonge, Bram IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from
More informationVersion A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.
Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x
More informationDefinition Pricing Risk management Second generation barrier options. Barrier Options. Arfima Financial Solutions
Arfima Financial Solutions Contents Definition 1 Definition 2 3 4 Contenido Definition 1 Definition 2 3 4 Definition Definition: A barrier option is an option on the underlying asset that is activated
More informationSemi-Markov model for market microstructure and HFT
Semi-Markov model for market microstructure and HFT LPMA, University Paris Diderot EXQIM 6th General AMaMeF and Banach Center Conference 10-15 June 2013 Joint work with Huyên PHAM LPMA, University Paris
More information1. For a special whole life insurance on (x), payable at the moment of death:
**BEGINNING OF EXAMINATION** 1. For a special whole life insurance on (x), payable at the moment of death: µ () t = 0.05, t > 0 (ii) δ = 0.08 x (iii) (iv) The death benefit at time t is bt 0.06t = e, t
More informationA Generic One-Factor Lévy Model for Pricing Synthetic CDOs
A Generic One-Factor Lévy Model for Pricing Synthetic CDOs Wim Schoutens - joint work with Hansjörg Albrecher and Sophie Ladoucette Maryland 30th of September 2006 www.schoutens.be Abstract The one-factor
More informationI. Time Series and Stochastic Processes
I. Time Series and Stochastic Processes Purpose of this Module Introduce time series analysis as a method for understanding real-world dynamic phenomena Define different types of time series Explain the
More informationUQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.
UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.
More informationMonte Carlo Simulation of Stochastic Processes
Monte Carlo Simulation of Stochastic Processes Last update: January 10th, 2004. In this section is presented the steps to perform the simulation of the main stochastic processes used in real options applications,
More informationHedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo
Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor
More informationOptimal rebalancing of portfolios with transaction costs assuming constant risk aversion
Optimal rebalancing of portfolios with transaction costs assuming constant risk aversion Lars Holden PhD, Managing director t: +47 22852672 Norwegian Computing Center, P. O. Box 114 Blindern, NO 0314 Oslo,
More informationChapter 5. Continuous Random Variables and Probability Distributions. 5.1 Continuous Random Variables
Chapter 5 Continuous Random Variables and Probability Distributions 5.1 Continuous Random Variables 1 2CHAPTER 5. CONTINUOUS RANDOM VARIABLES AND PROBABILITY DISTRIBUTIONS Probability Distributions Probability
More informationBandit Problems with Lévy Payoff Processes
Bandit Problems with Lévy Payoff Processes Eilon Solan Tel Aviv University Joint with Asaf Cohen Multi-Arm Bandits A single player sequential decision making problem. Time is continuous or discrete. The
More informationChapter 4 Continuous Random Variables and Probability Distributions
Chapter 4 Continuous Random Variables and Probability Distributions Part 2: More on Continuous Random Variables Section 4.5 Continuous Uniform Distribution Section 4.6 Normal Distribution 1 / 27 Continuous
More informationFinancial Models with Levy Processes and Volatility Clustering
Financial Models with Levy Processes and Volatility Clustering SVETLOZAR T. RACHEV # YOUNG SHIN ICIM MICHELE LEONARDO BIANCHI* FRANK J. FABOZZI WILEY John Wiley & Sons, Inc. Contents Preface About the
More informationFractional Black - Scholes Equation
Chapter 6 Fractional Black - Scholes Equation 6.1 Introduction The pricing of options is a central problem in quantitative finance. It is both a theoretical and practical problem since the use of options
More informationStochastic Modelling Unit 3: Brownian Motion and Diffusions
Stochastic Modelling Unit 3: Brownian Motion and Diffusions Russell Gerrard and Douglas Wright Cass Business School, City University, London June 2004 Contents of Unit 3 1 Introduction 2 Brownian Motion
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #3 1 Maximum likelihood of the exponential distribution 1. We assume
More informationSampling Distribution
MAT 2379 (Spring 2012) Sampling Distribution Definition : Let X 1,..., X n be a collection of random variables. We say that they are identically distributed if they have a common distribution. Definition
More informationIntroduction to Statistics I
Introduction to Statistics I Keio University, Faculty of Economics Continuous random variables Simon Clinet (Keio University) Intro to Stats November 1, 2018 1 / 18 Definition (Continuous random variable)
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationContinuous Processes. Brownian motion Stochastic calculus Ito calculus
Continuous Processes Brownian motion Stochastic calculus Ito calculus Continuous Processes The binomial models are the building block for our realistic models. Three small-scale principles in continuous
More informationLecture 15: Exotic Options: Barriers
Lecture 15: Exotic Options: Barriers Dr. Hanqing Jin Mathematical Institute University of Oxford Lecture 15: Exotic Options: Barriers p. 1/10 Barrier features For any options with payoff ξ at exercise
More informationNormal Inverse Gaussian (NIG) Process
With Applications in Mathematical Finance The Mathematical and Computational Finance Laboratory - Lunch at the Lab March 26, 2009 1 Limitations of Gaussian Driven Processes Background and Definition IG
More informationNormal Distribution. Notes. Normal Distribution. Standard Normal. Sums of Normal Random Variables. Normal. approximation of Binomial.
Lecture 21,22, 23 Text: A Course in Probability by Weiss 8.5 STAT 225 Introduction to Probability Models March 31, 2014 Standard Sums of Whitney Huang Purdue University 21,22, 23.1 Agenda 1 2 Standard
More informationMONTE CARLO EXTENSIONS
MONTE CARLO EXTENSIONS School of Mathematics 2013 OUTLINE 1 REVIEW OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO 3 SUMMARY MONTE CARLO SO FAR... Simple to program
More informationStatistics for Business and Economics
Statistics for Business and Economics Chapter 5 Continuous Random Variables and Probability Distributions Ch. 5-1 Probability Distributions Probability Distributions Ch. 4 Discrete Continuous Ch. 5 Probability
More informationEVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS. Rick Katz
1 EVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu
More informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More informationChapter 3 - Lecture 4 Moments and Moment Generating Funct
Chapter 3 - Lecture 4 and s October 7th, 2009 Chapter 3 - Lecture 4 and Moment Generating Funct Central Skewness Chapter 3 - Lecture 4 and Moment Generating Funct Central Skewness The expected value of
More informationOptimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models
Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University University of Missouri-Kansas City Department of Mathematics
More informationMAS3904/MAS8904 Stochastic Financial Modelling
MAS3904/MAS8904 Stochastic Financial Modelling Dr Andrew (Andy) Golightly a.golightly@ncl.ac.uk Semester 1, 2018/19 Administrative Arrangements Lectures on Tuesdays at 14:00 (PERCY G13) and Thursdays at
More informationChapter 4 Continuous Random Variables and Probability Distributions
Chapter 4 Continuous Random Variables and Probability Distributions Part 2: More on Continuous Random Variables Section 4.5 Continuous Uniform Distribution Section 4.6 Normal Distribution 1 / 28 One more
More informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
More informationEE266 Homework 5 Solutions
EE, Spring 15-1 Professor S. Lall EE Homework 5 Solutions 1. A refined inventory model. In this problem we consider an inventory model that is more refined than the one you ve seen in the lectures. The
More informationMAFS Computational Methods for Pricing Structured Products
MAFS550 - Computational Methods for Pricing Structured Products Solution to Homework Two Course instructor: Prof YK Kwok 1 Expand f(x 0 ) and f(x 0 x) at x 0 into Taylor series, where f(x 0 ) = f(x 0 )
More informationBrownian Motion. Richard Lockhart. Simon Fraser University. STAT 870 Summer 2011
Brownian Motion Richard Lockhart Simon Fraser University STAT 870 Summer 2011 Richard Lockhart (Simon Fraser University) Brownian Motion STAT 870 Summer 2011 1 / 33 Purposes of Today s Lecture Describe
More informationPart I - Normally-Distributed Asset Values
Part I - Normally-Distributed Asset Values Gary Schurman, MB, CFA September, 2018 In Part I of this sereies we will build an asset value model that assumes that asset values are normally-distributed. To
More informationCS 774 Project: Fall 2009 Version: November 27, 2009
CS 774 Project: Fall 2009 Version: November 27, 2009 Instructors: Peter Forsyth, paforsyt@uwaterloo.ca Office Hours: Tues: 4:00-5:00; Thurs: 11:00-12:00 Lectures:MWF 3:30-4:20 MC2036 Office: DC3631 CS
More informationLecture 23. STAT 225 Introduction to Probability Models April 4, Whitney Huang Purdue University. Normal approximation to Binomial
Lecture 23 STAT 225 Introduction to Probability Models April 4, 2014 approximation Whitney Huang Purdue University 23.1 Agenda 1 approximation 2 approximation 23.2 Characteristics of the random variable:
More informationMATH 3200 Exam 3 Dr. Syring
. Suppose n eligible voters are polled (randomly sampled) from a population of size N. The poll asks voters whether they support or do not support increasing local taxes to fund public parks. Let M be
More informationApplications of Lévy processes
Applications of Lévy processes Graduate lecture 29 January 2004 Matthias Winkel Departmental lecturer (Institute of Actuaries and Aon lecturer in Statistics) 6. Poisson point processes in fluctuation theory
More informationForecasting: an introduction. There are a variety of ad hoc methods as well as a variety of statistically derived methods.
Forecasting: an introduction Given data X 0,..., X T 1. Goal: guess, or forecast, X T or X T+r. There are a variety of ad hoc methods as well as a variety of statistically derived methods. Illustration
More informationChapter 1. Bond Pricing (continued)
Chapter 1 Bond Pricing (continued) How does the bond pricing illustrated here help investors in their investment decisions? This pricing formula can allow the investors to decide for themselves what the
More informationLinear Regression with One Regressor
Linear Regression with One Regressor Michael Ash Lecture 9 Linear Regression with One Regressor Review of Last Time 1. The Linear Regression Model The relationship between independent X and dependent Y
More informationAn Introduction to Stochastic Calculus
An Introduction to Stochastic Calculus Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 2-3 Haijun Li An Introduction to Stochastic Calculus Week 2-3 1 / 24 Outline
More informationOutline Brownian Process Continuity of Sample Paths Differentiability of Sample Paths Simulating Sample Paths Hitting times and Maximum
Normal Distribution and Brownian Process Page 1 Outline Brownian Process Continuity of Sample Paths Differentiability of Sample Paths Simulating Sample Paths Hitting times and Maximum Searching for a Continuous-time
More informationLast Time. Martingale inequalities Martingale convergence theorem Uniformly integrable martingales. Today s lecture: Sections 4.4.1, 5.
MATH136/STAT219 Lecture 21, November 12, 2008 p. 1/11 Last Time Martingale inequalities Martingale convergence theorem Uniformly integrable martingales Today s lecture: Sections 4.4.1, 5.3 MATH136/STAT219
More informationContinous time models and realized variance: Simulations
Continous time models and realized variance: Simulations Asger Lunde Professor Department of Economics and Business Aarhus University September 26, 2016 Continuous-time Stochastic Process: SDEs Building
More informationAPPENDIX C: STRESS-RANGE HISTOGRAM DATA AND REGRESSION
APPENDIX C: STRESS-RANGE HISTOGRAM DATA AND REGRESSION C-1 To determine the appropriate fatigue load for infinite life design, the fatigue-limit-state load with a probability of exceedance of 1:10,000
More informationContinuous Probability Distributions & Normal Distribution
Mathematical Methods Units 3/4 Student Learning Plan Continuous Probability Distributions & Normal Distribution 7 lessons Notes: Students need practice in recognising whether a problem involves a discrete
More informationModeling Credit Exposure for Collateralized Counterparties
Modeling Credit Exposure for Collateralized Counterparties Michael Pykhtin Credit Analytics & Methodology Bank of America Fields Institute Quantitative Finance Seminar Toronto; February 25, 2009 Disclaimer
More informationEconomics 2010c: Lecture 4 Precautionary Savings and Liquidity Constraints
Economics 2010c: Lecture 4 Precautionary Savings and Liquidity Constraints David Laibson 9/11/2014 Outline: 1. Precautionary savings motives 2. Liquidity constraints 3. Application: Numerical solution
More informationLecture 9 - Sampling Distributions and the CLT
Lecture 9 - Sampling Distributions and the CLT Sta102/BME102 Colin Rundel September 23, 2015 1 Variability of Estimates Activity Sampling distributions - via simulation Sampling distributions - via CLT
More informationProb and Stats, Nov 7
Prob and Stats, Nov 7 The Standard Normal Distribution Book Sections: 7.1, 7.2 Essential Questions: What is the standard normal distribution, how is it related to all other normal distributions, and how
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2019 Last Time: Markov Chains We can use Markov chains for density estimation, d p(x) = p(x 1 ) p(x }{{}
More informationOrnstein-Uhlenbeck Theory
Beatrice Byukusenge Department of Technomathematics Lappeenranta University of technology January 31, 2012 Definition of a stochastic process Let (Ω,F,P) be a probability space. A stochastic process is
More informationAMS Spring Western Sectional Meeting April 4, 2004 (USC campus) Two-sided barrier problems with jump-diffusions Alan L. Lewis
AMS Spring Western Sectional Meeting April 4, 2004 (USC campus) Two-sided barrier problems with jump-diffusions Alan L. Lewis No preprint yet, but overheads to be posted at www.optioncity.net (Publications)
More informationSustainable Maintenance Strategy Under Uncertainty in the Lifetime Distribution of Deteriorating Assets
Sustainable Maintenance Strategy Under Uncertainty in the Lifetime Distribution of Deteriorating Assets Alireza Daneshkhah, Amin Hosseinian-Far and Omid Chatrabgoun Abstract In the life-cycle management
More informationReverse Sensitivity Testing: What does it take to break the model? Silvana Pesenti
Reverse Sensitivity Testing: What does it take to break the model? Silvana Pesenti Silvana.Pesenti@cass.city.ac.uk joint work with Pietro Millossovich and Andreas Tsanakas Insurance Data Science Conference,
More informationShort & Long Run impact of volatility on the effect monetary shocks
Short & Long Run impact of volatility on the effect monetary shocks Fernando Alvarez University of Chicago & NBER Inflation: Drivers & Dynamics Conference 218 Cleveland Fed Alvarez Volatility & Monetary
More informationNORMAL APPROXIMATION. In the last chapter we discovered that, when sampling from almost any distribution, e r2 2 rdrdϕ = 2π e u du =2π.
NOMAL APPOXIMATION Standardized Normal Distribution Standardized implies that its mean is eual to and the standard deviation is eual to. We will always use Z as a name of this V, N (, ) will be our symbolic
More informationApplications of Good s Generalized Diversity Index. A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK
Applications of Good s Generalized Diversity Index A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK Internal Report STAT 98/11 September 1998 Applications of Good s Generalized
More informationEffectiveness of CPPI Strategies under Discrete Time Trading
Effectiveness of CPPI Strategies under Discrete Time Trading S. Balder, M. Brandl 1, Antje Mahayni 2 1 Department of Banking and Finance, University of Bonn 2 Department of Accounting and Finance, Mercator
More informationSection 1.3: More Probability and Decisions: Linear Combinations and Continuous Random Variables
Section 1.3: More Probability and Decisions: Linear Combinations and Continuous Random Variables Jared S. Murray The University of Texas at Austin McCombs School of Business OpenIntro Statistics, Chapters
More informationThe Vasicek Distribution
The Vasicek Distribution Dirk Tasche Lloyds TSB Bank Corporate Markets Rating Systems dirk.tasche@gmx.net Bristol / London, August 2008 The opinions expressed in this presentation are those of the author
More informationMarch 30, Preliminary Monte Carlo Investigations. Vivek Bhattacharya. Outline. Mathematical Overview. Monte Carlo. Cross Correlations
March 30, 2011 Motivation (why spend so much time on simulations) What does corr(rj 1, RJ 2 ) really represent? Results and Graphs Future Directions General Questions ( corr RJ (1), RJ (2)) = corr ( µ
More informationProbability Theory and Simulation Methods. April 9th, Lecture 20: Special distributions
April 9th, 2018 Lecture 20: Special distributions Week 1 Chapter 1: Axioms of probability Week 2 Chapter 3: Conditional probability and independence Week 4 Chapters 4, 6: Random variables Week 9 Chapter
More informationLecture 9 - Sampling Distributions and the CLT. Mean. Margin of error. Sta102/BME102. February 6, Sample mean ( X ): x i
Lecture 9 - Sampling Distributions and the CLT Sta102/BME102 Colin Rundel February 6, 2015 http:// pewresearch.org/ pubs/ 2191/ young-adults-workers-labor-market-pay-careers-advancement-recession Sta102/BME102
More informationOn the value of European options on a stock paying a discrete dividend at uncertain date
A Work Project, presented as part of the requirements for the Award of a Master Degree in Finance from the NOVA School of Business and Economics. On the value of European options on a stock paying a discrete
More informationDealing with forecast uncertainty in inventory models
Dealing with forecast uncertainty in inventory models 19th IIF workshop on Supply Chain Forecasting for Operations Lancaster University Dennis Prak Supervisor: Prof. R.H. Teunter June 29, 2016 Dennis Prak
More informationRandom Variables Handout. Xavier Vilà
Random Variables Handout Xavier Vilà Course 2004-2005 1 Discrete Random Variables. 1.1 Introduction 1.1.1 Definition of Random Variable A random variable X is a function that maps each possible outcome
More informationVOLATILITY FORECASTING IN A TICK-DATA MODEL L. C. G. Rogers University of Bath
VOLATILITY FORECASTING IN A TICK-DATA MODEL L. C. G. Rogers University of Bath Summary. In the Black-Scholes paradigm, the variance of the change in log price during a time interval is proportional to
More informationCHAPTER 6 Random Variables
CHAPTER 6 Random Variables 6.1 Discrete and Continuous Random Variables The Practice of Statistics, 5th Edition Starnes, Tabor, Yates, Moore Bedford Freeman Worth Publishers Discrete and Continuous Random
More informationAn Introduction to Stochastic Calculus
An Introduction to Stochastic Calculus Haijun Li lih@math.wsu.edu Department of Mathematics and Statistics Washington State University Lisbon, May 218 Haijun Li An Introduction to Stochastic Calculus Lisbon,
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationStatistics 6 th Edition
Statistics 6 th Edition Chapter 5 Discrete Probability Distributions Chap 5-1 Definitions Random Variables Random Variables Discrete Random Variable Continuous Random Variable Ch. 5 Ch. 6 Chap 5-2 Discrete
More informationImpressum ( 5 TMG) Herausgeber: Fakultät für Wirtschaftswissenschaft Der Dekan. Verantwortlich für diese Ausgabe:
WORKNG PAPER SERES mpressum ( 5 TMG) Herausgeber: Otto-von-Guericke-Universität Magdeburg Fakultät für Wirtschaftswissenschaft er ekan Verantwortlich für diese Ausgabe: Otto-von-Guericke-Universität Magdeburg
More informationIntroduction to Business Statistics QM 120 Chapter 6
DEPARTMENT OF QUANTITATIVE METHODS & INFORMATION SYSTEMS Introduction to Business Statistics QM 120 Chapter 6 Spring 2008 Chapter 6: Continuous Probability Distribution 2 When a RV x is discrete, we can
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationS t d with probability (1 p), where
Stochastic Calculus Week 3 Topics: Towards Black-Scholes Stochastic Processes Brownian Motion Conditional Expectations Continuous-time Martingales Towards Black Scholes Suppose again that S t+δt equals
More informationEconomics 883: The Basic Diffusive Model, Jumps, Variance Measures. George Tauchen. Economics 883FS Spring 2015
Economics 883: The Basic Diffusive Model, Jumps, Variance Measures George Tauchen Economics 883FS Spring 2015 Main Points 1. The Continuous Time Model, Theory and Simulation 2. Observed Data, Plotting
More informationLifetime Portfolio Selection: A Simple Derivation
Lifetime Portfolio Selection: A Simple Derivation Gordon Irlam (gordoni@gordoni.com) July 9, 018 Abstract Merton s portfolio problem involves finding the optimal asset allocation between a risky and a
More informationComputational Finance Improving Monte Carlo
Computational Finance Improving Monte Carlo School of Mathematics 2018 Monte Carlo so far... Simple to program and to understand Convergence is slow, extrapolation impossible. Forward looking method ideal
More informationContinuous-Time Pension-Fund Modelling
. Continuous-Time Pension-Fund Modelling Andrew J.G. Cairns Department of Actuarial Mathematics and Statistics, Heriot-Watt University, Riccarton, Edinburgh, EH4 4AS, United Kingdom Abstract This paper
More informationA Simulation Study of Bipower and Thresholded Realized Variations for High-Frequency Data
Washington University in St. Louis Washington University Open Scholarship Arts & Sciences Electronic Theses and Dissertations Arts & Sciences Spring 5-18-2018 A Simulation Study of Bipower and Thresholded
More informationMean Reverting Asset Trading. Research Topic Presentation CSCI-5551 Grant Meyers
Mean Reverting Asset Trading Research Topic Presentation CSCI-5551 Grant Meyers Table of Contents 1. Introduction + Associated Information 2. Problem Definition 3. Possible Solution 1 4. Problems with
More informationNo-arbitrage theorem for multi-factor uncertain stock model with floating interest rate
Fuzzy Optim Decis Making 217 16:221 234 DOI 117/s17-16-9246-8 No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Xiaoyu Ji 1 Hua Ke 2 Published online: 17 May 216 Springer
More informationOptimal inventory policies with an exact cost function under large demand uncertainty
MPRA Munich Personal RePEc Archive Optimal inventory policies with an exact cost function under large demand uncertainty George Halkos and Ilias Kevork and Chris Tziourtzioumis Department of Economics,
More informationAnalysis of truncated data with application to the operational risk estimation
Analysis of truncated data with application to the operational risk estimation Petr Volf 1 Abstract. Researchers interested in the estimation of operational risk often face problems arising from the structure
More informationMgr. Jakub Petrásek 1. May 4, 2009
Dissertation Report - First Steps Petrásek 1 2 1 Department of Probability and Mathematical Statistics, Charles University email:petrasek@karlin.mff.cuni.cz 2 RSJ Invest a.s., Department of Probability
More informationTheoretical Problems in Credit Portfolio Modeling 2
Theoretical Problems in Credit Portfolio Modeling 2 David X. Li Shanghai Advanced Institute of Finance (SAIF) Shanghai Jiaotong University(SJTU) November 3, 2017 Presented at the University of South California
More informationValue at Risk and Self Similarity
Value at Risk and Self Similarity by Olaf Menkens School of Mathematical Sciences Dublin City University (DCU) St. Andrews, March 17 th, 2009 Value at Risk and Self Similarity 1 1 Introduction The concept
More informationNumerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps
Numerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps, Senior Quantitative Analyst Motivation: Swing Options An electricity or gas SUPPLIER needs to be capable,
More informationNon-informative Priors Multiparameter Models
Non-informative Priors Multiparameter Models Statistics 220 Spring 2005 Copyright c 2005 by Mark E. Irwin Prior Types Informative vs Non-informative There has been a desire for a prior distributions that
More informationApplications of short-time asymptotics to the statistical estimation and option pricing of Lévy-driven models
Applications of short-time asymptotics to the statistical estimation and option pricing of Lévy-driven models José Enrique Figueroa-López 1 1 Department of Statistics Purdue University CIMAT and Universidad
More informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.
UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: STK4540 Non-Life Insurance Mathematics Day of examination: Wednesday, December 4th, 2013 Examination hours: 14.30 17.30 This
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More information