18-660: Numerical Methods for Engineering Design and Optimization
|
|
- Maximillian Shaw
- 6 years ago
- Views:
Transcription
1 18-660: Numerical Methods for Engineering Design and Optimization Xin Li Department of ECE Carnegie Mellon Universit Pittsburgh, PA Slide 1
2 Overview Monte Carlo Analsis Latin hpercube sampling Importance sampling Slide 2
3 Latin Hpercube Sampling (LHS) A great number of samples are tpicall required in traditional Monte Carlo to achieve good accurac Various techniques eist to improve Monte Carlo accurac Controlling sampling points is the ke Latin hpercube sampling is a widel-used method to generate controlled random samples The basic idea is to make sampling point distribution close to probabilit densit function (PDF) M. Mcka, R. Beckman and W. Conover, A comparison of three methods for selecting values of input variables in the analsis of output from a computer code, Technometrics, vol. 21, no. 2, pp , Ma Slide 3
4 Latin Hpercube Sampling (LHS) One dimensional Latin hpercube sampling Evenl partition CDF into N regions Randoml pick up one sampling point in each region One random sample in one region 1 CDF Evenl partitioned regions 0 Avoid the probabilit that all sampling points come from the same local region Slide 4
5 Latin Hpercube Sampling (LHS) Two dimensional Latin hpercube sampling 1 and 2 must be independent Generate one-dimensional LHS samples for 1 Generate one-dimensional LHS samples for 2 Randoml combine the LHS samples to two-dimensional pairs 1 Probabilit of Probabilit of 1 One sample in each row and each column Sampling is random in each grid Higher-dimensional LHS samples can be similarl generated Slide 5
6 Latin Hpercube Sampling (LHS) Matlab code for LHS sampling of independent standard Normal distributions data = rand(nsample,nvar); for i = 1:NVar inde = randperm(nsample); prob = (inde'-data(:,i))/nsample; data(:,i) = sqrt(2)*erfinv(2*prob-1); end; NVar: NSample: data: # of random variables # of samples LHS sampling points Slide 6
7 Latin Hpercube Sampling (LHS) Compare Monte Carlo accurac for a simple eample ~ N(0,1) (standard Normal distribution) Repeatedl estimate the mean value b Monte Carlo analsis Accurac is improved b increasing sampling # Estimated Mean Value Estimated Mean Value # of Samples Random sampling Monte Carlo is not deterministic # of Samples Latin hpercube sampling Slide 7
8 Even with Latin hpercube sampling, Monte Carlo analsis requires a HUGE number of sampling points Eample: rare event estimation ~ N( 0,1) Standard Normal distribution Estimate P( 5 ) =??? The theoretical answer for P( -5) is equal to ~100M sampling points are required if we attempt to estimate this probabilit b random sampling or LHS Slide 8
9 Ke idea: Do not generate random samples from (t) Instead, find a good distorted (t) to improve Monte Carlo sampling accurac Eample: if ~ N(0,1), what is P( -5)? Intuitivel, if we draw sampling points based on (t), more samples will fall into the gre area (t) PDF (t) How do we calculate P( -5) when sampling (t)? P( -5) 0 t Slide 9
10 Assume that we want to estimate the following epected value E [ f ( ) ] f Eample: if we want to estimate P( -5), then + = dt f ( ) = 1 0 ( 5) ( > 5) E 5 + [ f ( ) ] 1 dt + 0 dt = P( 5) = 5 Slide 10
11 Estimate E[f()] where ~ (t) b importance sampling E [ f ( ) ] f dt = f + + = dt E [ f ( ) ] f + = + = dt = g dt E[ g( ) ] Slide 11
12 Estimate E[f()] where ~ (t) b traditional sampling Step 1: draw M random samples {t 1,t 2,...,t M } based on (t) Step 2: calculate f m = f(t m ) at each sampling point m = 1,2,...,M Step 3: calculate E[f] (f 1 +f f M )/M Estimate E[f()] where ~ (t) b importance sampling Step 1: draw M random samples {t 1,t 2,...,t M } based on (t) Step 2: calculate g m = f(t m ) (t m )/ (t m ) at each sampling point m = 1,2,...,M Step 3: calculate E[f] (g 1 +g g M )/M How do we decide the optimal (t) to achieve minimal Monte Carlo analsis error? Slide 12
13 Determine optimal (t) for importance sampling E [ f ] µ = f ( t ) Estimator f 1 M ( tm ) ( t ) The accurac of an estimator can be quantitativel measured b its variance To improve Monte Carlo analsis accurac, we should minimize VAR[μ f ] M m= 1 m [ ] Error ~ VAR µ f m Slide 13
14 Determine optimal (t) for importance sampling E [ f ] µ = f ( t ) m= 1 We achieve the minimal VAR[μ f ] = 0 if f = k f 1 M (Constant) f = (Optimal PDF) k M m ( tm ) ( t ) m f(t m ) (t m )/ (t m ) is equal to the same constant for all sampling points Slide 14
15 How do we decide the value k? = f k K cannot be arbitraril selected (t) must be a valid PDF that satisfies the following condition + dt = + f + k = f = k dt dt E[ f ] =1 Finding the optimal (t) requires to know E[f], i.e., the answer of our Monte Carlo analsis!!! Slide 15
16 In practice, such an optimal (t) cannot be easil applied Instead, we tpicall look for a sub-optimal solution that satisfies the following constraints Eas to construct we do not have to know k = E[f] Eas to sample not all random distributions can be easil sampled b a random number generator Minimal estimator variance the sub-optimal (t) is close to the optimal case as much as possible Slide 16
17 Finding the right (t) is nontrivial for practical problems No magic equation eists in general Engineering approach is based on heuristics Sometimes require a lot of human eperience and numerical optimization The criterion to choose (t) is also application-dependent Slide 17
18 Eample: if ~ N(0,1), what is P( -5)? Shifted Normal PDF (t) Uniform P( -5) 0 Several possible choices for (t) t Slide 18
19 Summar Monte Carlo analsis Latin hpercube sampling Importance sampling Slide 19
Some Formulas neglected in Anderson, Sweeny, and Williams, with a Digression on Statistics and Finance
Some Formulas neglected in Anderson, Sween, and Williams, with a Digression on Statistics and Finance Transformations of a Single Random variable: If ou have a case where a new random variable is defined
More informationA Recommended Financial Model for the Selection of Safest portfolio by using Simulation and Optimization Techniques
Journal of Applied Finance & Banking, vol., no., 20, 3-42 ISSN: 792-6580 (print version), 792-6599 (online) International Scientific Press, 20 A Recommended Financial Model for the Selection of Safest
More informationJoint Distribution of Stock Market Returns and Trading Volume
Rev. Integr. Bus. Econ. Res. Vol 5(3) 0 Joint Distribution of Stock Market Returns and Trading Volume Muhammad Idrees Ahmad * Department of Mathematics and Statistics, Sultan Qaboos Universit, Muscat,
More informationChapter 2 Uncertainty Analysis and Sampling Techniques
Chapter 2 Uncertainty Analysis and Sampling Techniques The probabilistic or stochastic modeling (Fig. 2.) iterative loop in the stochastic optimization procedure (Fig..4 in Chap. ) involves:. Specifying
More informationCentral Limit Theorem, Joint Distributions Spring 2018
Central Limit Theorem, Joint Distributions 18.5 Spring 218.5.4.3.2.1-4 -3-2 -1 1 2 3 4 Exam next Wednesday Exam 1 on Wednesday March 7, regular room and time. Designed for 1 hour. You will have the full
More information2.1 Properties of PDFs
2.1 Properties of PDFs mode median epectation values moments mean variance skewness kurtosis 2.1: 1/13 Mode The mode is the most probable outcome. It is often given the symbol, µ ma. For a continuous random
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 informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationProbability Theory. Mohamed I. Riffi. Islamic University of Gaza
Probability Theory Mohamed I. Riffi Islamic University of Gaza Table of contents 1. Chapter 2 Discrete Distributions The binomial distribution 1 Chapter 2 Discrete Distributions Bernoulli trials and the
More informationCh4. Variance Reduction Techniques
Ch4. Zhang Jin-Ting Department of Statistics and Applied Probability July 17, 2012 Ch4. Outline Ch4. This chapter aims to improve the Monte Carlo Integration estimator via reducing its variance using some
More informationBasic Principles of Probability and Statistics. Lecture notes for PET 472 Spring 2010 Prepared by: Thomas W. Engler, Ph.D., P.E
Basic Principles of Probability and Statistics Lecture notes for PET 472 Spring 2010 Prepared by: Thomas W. Engler, Ph.D., P.E Definitions Risk Analysis Assessing probabilities of occurrence for each possible
More informationMarquette University MATH 1700 Class 8 Copyright 2018 by D.B. Rowe
Class 8 Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science Copyright 208 by D.B. Rowe Agenda: Recap Chapter 4.3-4.5 Lecture Chapter 5. - 5.3 2 Recap Chapter 4.3-4.5 3 4:
More informationBasic Principles of Probability and Statistics. Lecture notes for PET 472 Spring 2012 Prepared by: Thomas W. Engler, Ph.D., P.E
Basic Principles of Probability and Statistics Lecture notes for PET 472 Spring 2012 Prepared by: Thomas W. Engler, Ph.D., P.E Definitions Risk Analysis Assessing probabilities of occurrence for each possible
More information36106 Managerial Decision Modeling Monte Carlo Simulation in Excel: Part IV
36106 Managerial Decision Modeling Monte Carlo Simulation in Excel: Part IV Kipp Martin University of Chicago Booth School of Business November 29, 2017 Reading and Excel Files 2 Reading: Handout: Optimal
More informationEcon 300: Quantitative Methods in Economics. 11th Class 10/19/09
Econ 300: Quantitative Methods in Economics 11th Class 10/19/09 Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write. --H.G. Wells discuss test [do
More informationCIVL Learning Objectives. Definitions. Discrete Distributions
CIVL 3103 Discrete Distributions Learning Objectives Define discrete distributions, and identify common distributions applicable to engineering problems. Identify the appropriate distribution (i.e. binomial,
More informationClass 13. Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science. Marquette University MATH 1700
Class 13 Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science Copyright 017 by D.B. Rowe 1 Agenda: Recap Chapter 6.3 6.5 Lecture Chapter 7.1 7. Review Chapter 5 for Eam 3.
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 informationStats CH 6 Intro Activity 1
Stats CH 6 Intro Activit 1 1. Purpose can ou tell the difference between bottled water and tap water? You will drink water from 3 samples. 1 of these is bottled water.. You must test them in the following
More informationLecture outline. Monte Carlo Methods for Uncertainty Quantification. Importance Sampling. Importance Sampling
Lecture outline Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford KU Leuven Summer School on Uncertainty Quantification Lecture 2: Variance reduction
More informationGenerating Non-Gaussian Vibration for Testing Purposes
Generating Vibration for Testing Purposes David O. Smallwood, Albuquerque, New Meico There is increasing interest in simulating vibration environments that are non-, particularl in the transportation arena.
More informationAnalytic Models of the ROC Curve: Applications to Credit Rating Model Validation
Analtic Models of the ROC Curve: Applications to Credit Rating Model Validation Steve Satchell Facult of Economics and Politics Universit of Cambridge Email: ses@econ.cam.ac.uk Wei Xia Birkbeck College
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2018 Last Time: Markov Chains We can use Markov chains for density estimation, p(x) = p(x 1 ) }{{} d p(x
More informationLean Six Sigma: Training/Certification Books and Resources
Lean Si Sigma Training/Certification Books and Resources Samples from MINITAB BOOK Quality and Si Sigma Tools using MINITAB Statistical Software A complete Guide to Si Sigma DMAIC Tools using MINITAB Prof.
More informationProbability Distributions
Probability Distributions Probability Distributions CEE 2L. Uncertainty, Design, and Optimization Department of Civil and Environmental Engineering Duke University Philip Scott Harvey, Henri P. Gavin and
More informationStratified Sampling in Monte Carlo Simulation: Motivation, Design, and Sampling Error
South Texas Project Risk- Informed GSI- 191 Evaluation Stratified Sampling in Monte Carlo Simulation: Motivation, Design, and Sampling Error Document: STP- RIGSI191- ARAI.03 Revision: 1 Date: September
More informationMonte Carlo Methods in Option Pricing. UiO-STK4510 Autumn 2015
Monte Carlo Methods in Option Pricing UiO-STK4510 Autumn 015 The Basics of Monte Carlo Method Goal: Estimate the expectation θ = E[g(X)], where g is a measurable function and X is a random variable such
More informationImportance Sampling. Sargur N. Srihari
Importance Sampling Sargur N. srihari@cedar.buffalo.edu 1 Topics in Monte Carlo Methods 1. Sampling and Monte Carlo Methods 2. Importance Sampling 3. Markov Chain Monte Carlo Methods 4. Gibbs Sampling
More informationMA Lesson 13 Notes Section 4.1 (calculus part of textbook, page 196) Techniques for Finding Derivatives
Notation for the Derivative: MA 15910 Lesson 13 Notes Section 4.1 (calculus part of tetbook, page 196) Techniques for Finding Derivatives The derivative of a function y f ( ) may be written in any of these
More information"Vibrato" Monte Carlo evaluation of Greeks
"Vibrato" Monte Carlo evaluation of Greeks (Smoking Adjoints: part 3) Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance MCQMC 2008,
More informationGamma. The finite-difference formula for gamma is
Gamma The finite-difference formula for gamma is [ P (S + ɛ) 2 P (S) + P (S ɛ) e rτ E ɛ 2 ]. For a correlation option with multiple underlying assets, the finite-difference formula for the cross gammas
More informationDuality & The Dual Simplex Method & Sensitivity Analysis for Linear Programming. Metodos Cuantitativos M. En C. Eduardo Bustos Farias 1
Dualit & The Dual Simple Method & Sensitivit Analsis for Linear Programming Metodos Cuantitativos M. En C. Eduardo Bustos Farias Dualit EverLP problem has a twin problem associated with it. One problem
More informationInformation Revelation and Market Crashes
Information Revelation and Market Crashes Jan Werner Department of Economics Universit of Minnesota Minneapolis, MN 55455 September 2004 Revised: Ma 2005 Abstract: We show the possibilit of market crash
More informationAccelerated Option Pricing Multiple Scenarios
Accelerated Option Pricing in Multiple Scenarios 04.07.2008 Stefan Dirnstorfer (stefan@thetaris.com) Andreas J. Grau (grau@thetaris.com) 1 Abstract This paper covers a massive acceleration of Monte-Carlo
More information15 : Approximate Inference: Monte Carlo Methods
10-708: Probabilistic Graphical Models 10-708, Spring 2016 15 : Approximate Inference: Monte Carlo Methods Lecturer: Eric P. Xing Scribes: Binxuan Huang, Yotam Hechtlinger, Fuchen Liu 1 Introduction to
More informationChapter 7: Estimation Sections
1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationON PROFIT DENSITY BASED GREEDY ALGORITHM FOR A RESOURCE ALLOCATION PROBLEM IN WEB SERVICES
International Journal of Computers and Applications, Vol. 9, No., 007 ON PROFIT DENSITY BASED GREEDY ALGORITHM FOR A RESOURCE ALLOCATION PROBLEM IN WEB SERVICES J. Sum, J. Wu, and C.-S. Leung Abstract
More informationM1 M1 A1 M1 A1 M1 A1 A1 A1 11 A1 2 B1 B1. B1 M1 Relative efficiency (y) = M1 A1 BEWARE PRINTED ANSWER. 5
Q L e σ π ( W μ e σ π ( W μ M M A Product form. Two Normal terms. Fully correct. (ii ln L const ( W ( W d ln L ( W + ( W dμ 0 σ W σ μ W σ W W ˆ μ σ Chec this is a maximum. d ln L E.g. < 0 dμ σ σ σ μ σ
More informationComputational Finance Least Squares Monte Carlo
Computational Finance Least Squares Monte Carlo School of Mathematics 2019 Monte Carlo and Binomial Methods In the last two lectures we discussed the binomial tree method and convergence problems. One
More informationRisk Neutral Valuation, the Black-
Risk Neutral Valuation, the Black- Scholes Model and Monte Carlo Stephen M Schaefer London Business School Credit Risk Elective Summer 01 C = SN( d )-PV( X ) N( ) N he Black-Scholes formula 1 d (.) : cumulative
More informationMathematics in Finance
Mathematics in Finance Steven E. Shreve Department of Mathematical Sciences Carnegie Mellon University Pittsburgh, PA 15213 USA shreve@andrew.cmu.edu A Talk in the Series Probability in Science and Industry
More informationHow to Consider Risk Demystifying Monte Carlo Risk Analysis
How to Consider Risk Demystifying Monte Carlo Risk Analysis James W. Richardson Regents Professor Senior Faculty Fellow Co-Director, Agricultural and Food Policy Center Department of Agricultural Economics
More informationMEASURES OF RISK-AVERSION
Department of Economics, Universit of California, Davis Professor Giacomo Bonanno Ecn 03 Economics of Uncertaint and Information Lecture MESURES OF RISK-VERSION Can we answer the question: Which of two
More informationMarkov Decision Processes
Markov Decision Processes Robert Platt Northeastern University Some images and slides are used from: 1. CS188 UC Berkeley 2. AIMA 3. Chris Amato Stochastic domains So far, we have studied search Can use
More informationYao s Minimax Principle
Complexity of algorithms The complexity of an algorithm is usually measured with respect to the size of the input, where size may for example refer to the length of a binary word describing the input,
More informationGRAPHS IN ECONOMICS. Appendix. Key Concepts. A Positive Relationship
Appendi GRAPHS IN ECONOMICS Ke Concepts Graphing Data Graphs represent quantit as a distance on a line. On a graph, the horizontal scale line is the -ais, the vertical scale line is the -ais, and the intersection
More informationApplications of Local Gaussian Correlation in Finance
Applications of Local Gaussian Correlation in Finance and towards a parametric version of the LGC Bård Støve Universit of Bergen, Norwa Department of Mathematics Joint work in progress with Dag Tjøstheim,
More informationClass 11. Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science. Marquette University MATH 1700
Class 11 Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science Copyright 2017 by D.B. Rowe 1 Agenda: Recap Chapter 5.3 continued Lecture 6.1-6.2 Go over Eam 2. 2 5: Probability
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 informationNumerical Methods II
Numerical Methods II Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute MC Lecture 3 p. 1 Variance Reduction Monte Carlo starts as a very simple method; much of the complexity
More informationECE 295: Lecture 03 Estimation and Confidence Interval
ECE 295: Lecture 03 Estimation and Confidence Interval Spring 2018 Prof Stanley Chan School of Electrical and Computer Engineering Purdue University 1 / 23 Theme of this Lecture What is Estimation? You
More informationSummary Sampling Techniques
Summary Sampling Techniques MS&E 348 Prof. Gerd Infanger 2005/2006 Using Monte Carlo sampling for solving the problem Monte Carlo sampling works very well for estimating multiple integrals or multiple
More informationStochastic Optimization with cvxpy EE364b Project Final Report
Stochastic Optimization with cvpy EE364b Project Final Report Alnur Ali alnurali@cmu.edu June 5, 2015 1 Introduction A stochastic program is a conve optimization problem that includes random variables,
More informationIntroduction to Sequential Monte Carlo Methods
Introduction to Sequential Monte Carlo Methods Arnaud Doucet NCSU, October 2008 Arnaud Doucet () Introduction to SMC NCSU, October 2008 1 / 36 Preliminary Remarks Sequential Monte Carlo (SMC) are a set
More informationImportance Sampling for Fair Policy Selection
Importance Sampling for Fair Policy Selection Shayan Doroudi Carnegie Mellon University Pittsburgh, PA 15213 shayand@cs.cmu.edu Philip S. Thomas Carnegie Mellon University Pittsburgh, PA 15213 philipt@cs.cmu.edu
More informationThe Course So Far. Atomic agent: uninformed, informed, local Specific KR languages
The Course So Far Traditional AI: Deterministic single agent domains Atomic agent: uninformed, informed, local Specific KR languages Constraint Satisfaction Logic and Satisfiability STRIPS for Classical
More informationA Numerical Approach to the Estimation of Search Effort in a Search for a Moving Object
Proceedings of the 1. Conference on Applied Mathematics and Computation Dubrovnik, Croatia, September 13 18, 1999 pp. 129 136 A Numerical Approach to the Estimation of Search Effort in a Search for a Moving
More informationReal Options Valuation, Inc. Software Technical Support
Real Options Valuation, Inc. Software Technical Support HELPFUL TIPS AND TECHNIQUES Johnathan Mun, Ph.D., MBA, MS, CFC, CRM, FRM, MIFC 1 P a g e Helpful Tips and Techniques The following are some quick
More informationInternational Finance. Estimation Error. Campbell R. Harvey Duke University, NBER and Investment Strategy Advisor, Man Group, plc.
International Finance Estimation Error Campbell R. Harvey Duke University, NBER and Investment Strategy Advisor, Man Group, plc February 17, 2017 Motivation The Markowitz Mean Variance Efficiency is the
More informationMEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL
MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL Isariya Suttakulpiboon MSc in Risk Management and Insurance Georgia State University, 30303 Atlanta, Georgia Email: suttakul.i@gmail.com,
More informationMATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, Student Name (print):
MATH4143 Page 1 of 17 Winter 2007 MATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, 2007 Student Name (print): Student Signature: Student ID: Question
More information1 The principal-agent problems
1 The principal-agent problems The principal-agent problems are at the heart of modern economic theory. One of the reasons for this is that it has widespread applicability. We start with some eamples.
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 informationESD.70J Engineering Economy
ESD.70J Engineering Economy Fall 2010 Session One Xin Zhang xinzhang@mit.edu Prof. Richard de Neufville ardent@mit.edu http://ardent.mit.edu/real_options/rocse_excel_latest/excel_class.html ESD.70J Engineering
More informationand Pricing Problems
Mechanism Design, Machine Learning, and Pricing Problems Maria-Florina Balcan Carnegie Mellon University Overview Pricing and Revenue Maimization Software Pricing Digital Music Pricing Problems One Seller,
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMS091)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMS091) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 3 Importance sampling January 27, 2015 M. Wiktorsson
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulation Efficiency and an Introduction to Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University
More informationAnnuities. Lecture: Weeks 8-9. Lecture: Weeks 8-9 (Math 3630) Annuities Fall Valdez 1 / 41
Annuities Lecture: Weeks 8-9 Lecture: Weeks 8-9 (Math 3630) Annuities Fall 2017 - Valdez 1 / 41 What are annuities? What are annuities? An annuity is a series of payments that could vary according to:
More informationA Macro-Finance Model of the Term Structure: the Case for a Quadratic Yield Model
Title page Outline A Macro-Finance Model of the Term Structure: the Case for a 21, June Czech National Bank Structure of the presentation Title page Outline Structure of the presentation: Model Formulation
More information1.1 Interest rates Time value of money
Lecture 1 Pre- Derivatives Basics Stocks and bonds are referred to as underlying basic assets in financial markets. Nowadays, more and more derivatives are constructed and traded whose payoffs depend on
More informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to 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 information2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises
96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with
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 informationWrite legibly. Unreadable answers are worthless.
MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are
More informationA METHOD FOR STOCHASTIC ESTIMATION OF COST AND COMPLETION TIME OF A MINING PROJECT
A METHOD FOR STOCHASTIC ESTIMATION OF COST AND COMPLETION TIME OF A MINING PROJECT E. Newby, F. D. Fomeni, M. M. Ali and C. Musingwini Abstract The standard methodology used for estimating the cost and
More informationReinforcement Learning 04 - Monte Carlo. Elena, Xi
Reinforcement Learning 04 - Monte Carlo Elena, Xi Previous lecture 2 Markov Decision Processes Markov decision processes formally describe an environment for reinforcement learning where the environment
More informationAlternative VaR Models
Alternative VaR Models Neil Roeth, Senior Risk Developer, TFG Financial Systems. 15 th July 2015 Abstract We describe a variety of VaR models in terms of their key attributes and differences, e.g., parametric
More informationAnnuities. Lecture: Weeks Lecture: Weeks 9-11 (Math 3630) Annuities Fall Valdez 1 / 44
Annuities Lecture: Weeks 9-11 Lecture: Weeks 9-11 (Math 3630) Annuities Fall 2017 - Valdez 1 / 44 What are annuities? What are annuities? An annuity is a series of payments that could vary according to:
More informationMS&E 348 Winter 2011 BOND PORTFOLIO MANAGEMENT: INCORPORATING CORPORATE BOND DEFAULT
MS&E 348 Winter 2011 BOND PORTFOLIO MANAGEMENT: INCORPORATING CORPORATE BOND DEFAULT March 19, 2011 Assignment Overview In this project, we sought to design a system for optimal bond management. Within
More informationChapter 7: Portfolio Theory
Chapter 7: Portfolio Theory 1. Introduction 2. Portfolio Basics 3. The Feasible Set 4. Portfolio Selection Rules 5. The Efficient Frontier 6. Indifference Curves 7. The Two-Asset Portfolio 8. Unrestriceted
More informationValue at Risk Ch.12. PAK Study Manual
Value at Risk Ch.12 Related Learning Objectives 3a) Apply and construct risk metrics to quantify major types of risk exposure such as market risk, credit risk, liquidity risk, regulatory risk etc., and
More informationResults for option pricing
Results for option pricing [o,v,b]=optimal(rand(1,100000 Estimators = 0.4619 0.4617 0.4618 0.4613 0.4619 o = 0.46151 % best linear combination (true value=0.46150 v = 1.1183e-005 %variance per uniform
More informationConsider the following examples: ex: let X = tossing a coin three times and counting the number of heads
Overview Both chapters and 6 deal with a similar concept probability distributions. The difference is that chapter concerns itself with discrete probability distribution while chapter 6 covers continuous
More informationLecture 3: Return vs Risk: Mean-Variance Analysis
Lecture 3: Return vs Risk: Mean-Variance Analysis 3.1 Basics We will discuss an important trade-off between return (or reward) as measured by expected return or mean of the return and risk as measured
More informationA Heuristic Method for Statistical Digital Circuit Sizing
A Heuristic Method for Statistical Digital Circuit Sizing Stephen Boyd Seung-Jean Kim Dinesh Patil Mark Horowitz Microlithography 06 2/23/06 Statistical variation in digital circuits growing in importance
More information5 Profit maximization, Supply
Microeconomics I - Lecture #5, March 17, 2009 5 Profit maximization, Suppl We alread described the technological possibilities now we analze how the firm chooses the amount to produce so as to maximize
More informationCHAPTER 16. SECTION 16.1 (page 1168) SECTION 16.3 (page 1192) SECTION 16.2 (page 1179) Skills Review (page 1168) Skills Review (page 1192)
Answers to Selected Eercises A CHAPTER SECTION. (page ) Skills Review (page )..,..,... nn n n. nn n. n!. n!....... Permutations of seating positions.,,.. (a) (b) (c). A: ; B: ; A.,.. (a), (b).....,.. nn
More informationLESSON 7 INTERVAL ESTIMATION SAMIE L.S. LY
LESSON 7 INTERVAL ESTIMATION SAMIE L.S. LY 1 THIS WEEK S PLAN Part I: Theory + Practice ( Interval Estimation ) Part II: Theory + Practice ( Interval Estimation ) z-based Confidence Intervals for a Population
More informationA Probabilistic Approach to Determining the Number of Widgets to Build in a Yield-Constrained Process
A Probabilistic Approach to Determining the Number of Widgets to Build in a Yield-Constrained Process Introduction Timothy P. Anderson The Aerospace Corporation Many cost estimating problems involve determining
More informationSample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models
Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models Scott Robertson Carnegie Mellon University scottrob@andrew.cmu.edu http://www.math.cmu.edu/users/scottrob June
More informationII. Random Variables
II. Random Variables Random variables operate in much the same way as the outcomes or events in some arbitrary sample space the distinction is that random variables are simply outcomes that are represented
More information9 Some Demographic Topics
9 Some Demographic Topics 9.1 Single Figure Indices How do we compare the mortality eperience of two or more populations? A single figure that summarises the mortality eperience of a population enables
More informationRecovering portfolio default intensities implied by CDO quotes. Rama CONT & Andreea MINCA. March 1, Premia 14
Recovering portfolio default intensities implied by CDO quotes Rama CONT & Andreea MINCA March 1, 2012 1 Introduction Premia 14 Top-down" models for portfolio credit derivatives have been introduced as
More informationContinuous probability distribution
Microarray Center BIOSTATISTICS Lecture 6 Continuous Probability Distributions 16-4-1 Lecture 6. Continuous probability distributions Dr. Petr Nazarov petr.nazarov@crp-sante.lu OUTLINE Lecture 1 Continuous
More informationStatistical analysis and bootstrapping
Statistical analysis and bootstrapping p. 1/15 Statistical analysis and bootstrapping Michel Bierlaire michel.bierlaire@epfl.ch Transport and Mobility Laboratory Statistical analysis and bootstrapping
More informationAnnuities. Lecture: Weeks 8-9. Lecture: Weeks 8-9 (Math 3630) Annuities Fall Valdez 1 / 41
Annuities Lecture: Weeks 8-9 Lecture: Weeks 8-9 (Math 3630) Annuities Fall 2017 - Valdez 1 / 41 What are annuities? What are annuities? An annuity is a series of payments that could vary according to:
More informationJDEP 384H: Numerical Methods in Business
Chapter 4: Numerical Integration: Deterministic and Monte Carlo Methods Chapter 8: Option Pricing by Monte Carlo Methods JDEP 384H: Numerical Methods in Business Instructor: Thomas Shores Department of
More information