1 Residual life for gamma and Weibull distributions
|
|
- Nathan Murphy
- 6 years ago
- Views:
Transcription
1 Supplement to Tail Estimation for Window Censored Processes Residual life for gamma and Weibull distributions. Gamma distribution Let Γ(k, x = x yk e y dy be the upper incomplete gamma function, and let Γ(k = Γ(k, 0 be the gamma function. The gamma distribution with shape parameter k > 0 and scale parameter σ > 0 has the probability density and tail functions f 0 (x = σ k Γ(k xk e x/σ and F0 (x = Γ(k, x/σ, Γ(k respectively, and has mean µ 0 = kσ. Thus, using the equality Γ(k + = kγ(k, the residual life probability density and cumulative distribution functions are f r 0 (x = σγ(k + Γ(k, x/σ and F r 0 (x = (Γ(k +, x/σ x Γ(k + σ Γ(k, x/σ, where F r 0 was obtained by changing the integration order. All the above holds for x > 0.. Weibull distribution The probability density and the cumulative distribution functions of the Weibull distribution with shape parameter k > 0 and scale parameter σ > 0 are given by f 0 (x = k ( x k e (x/σ k and F 0 (x = e (x/σk, σ σ respectively, and the mean is µ 0 = σγ( + /k. Thus, using a change of variables to obtain the second expression, and with γ(x, k = x 0 yk e y dy, the lower incomplete gamma function, the density and distribution functions of the residual life time are f r 0 (x = σγ( + /k e (x/σk and F r 0 (x = ( ( x k Γ(/k γ k,. σ
2 Efficient algorithms for numerical computation of the incomplete gamma function are available in e.g. MatLab or Wolfram Mathematica software. Asymptotic normality Here we establish that the tail estimators ˆθ = (ˆσ, ˆγ from Subsection. satisfy the conditions of Theorem 5.3 of van der Vaart (998, and hence are asymptotically normal. For this we assume that θ = (σ, γ belongs to a compact set K = [σ l, σ u ] [γ l, γ u ], with 0 < σ l and 0 < γ l < γ u <. Here the last restriction, as before, is to ensure that means exist. We use notation and equation numbers from the paper and assume that 0 < u < w so that the observation window is longer than the threshold. Since σ is a scale parameter, it is only the ratio between σ and w u which determines asymptotic behavior, and not their individual values. Hence we fixed w u = in the calculations. Additionally, we did them for σ l = 0.05 and σ u =.5, since censoring is too mild, or too severe, for the estimation problem to be practically interesting for σ outside of this range. We further used γ l = 0.0 and γ u = 0.99 However, it is straightforward to consider also wider parameter ranges (and we in fact have done this too. Recall the notation X = L u for the excesses of the threshold u by the observed 0-intervals, and S for the length of the first -interval, and that X and S are assumed to be independent. Let m nc,θ (x = log σ γ + log ( + γ σ x, m rc,θ (x = ( γ log + γ σ x, m lc,θ (x = log σ + log( γ γ log ( + γ σ x m dc,θ = γ log ( + γ σ (w u,, and let nc (X = if X is non-censored and nc (X = 0 otherwise, and define rc, lc and dc similarly. Next, set m 0,θ (X = lc (Xm lc,θ (X + dc (Xm dc,θ m,θ (X = nc (Xm nc,θ (X + rc (Xm rc,θ (X m θ (X = m 0,θ (X + m,θ (X so that, with l u (σ, γ defined in (7, and with X,... X n the entire set of n = n nc + n rc + n lc + n dc
3 observed 0-intervals, we have that l u (θ = n m θ (X i. i= Now, inspection of the partial derivatives of m nc,θ (X with respect to σ and γ (computed with the help of Mathematica, shows that they are continuous on the compact set (θ, x K [0, ], and hence are uniformly bounded. Thus, estimating nc by, we obtain that nc (xm nc,θ (x nc (xm nc,θ (x C( σ σ + γ γ, where C denotes a generic constant which does not depend on θ or x. The corresponding inequalities for the rc, lc, and nc cases are obtained in precisely the same way, and, using Cauchy s inequality, it follows that m θ (x m θ (x C θ θ, with denoting Euclidean norm, so that the first condition in Theorem 5.3 of van der Vaart (998 is satisfied. Further, write E for expectation taken with respect to the true parameter value θ 0 = (σ 0, γ 0. We next use numerical computation to show that the symmetric matrix θ E (m 0,θ (X, evaluated at θ 0 is strictly negative definite, by showing that the determinant is strictly positive and that the diagonal elements are negative. Inspection of all the second order derivatives show they also are continuous on the compact set (θ, x K [0, ], and hence are uniformly bounded. It follows that the operations of differentiation and taking expectations may be freely interchanged. Thus we start by showing that the determinant D 0 (σ 0, σ, γ 0, γ := E σ m 0,θ(X E γ m 0,θ(X ( ( E σ γ m 0,θ(X (S- evaluated at θ = θ 0 is strictly positive. Recall the notation p 0 = E ( lc (X + dc (X and p = p 0 for the probabilities that X(0 = 0 and X(0 =, respectively. Considering the first factor in the righthand side of (S-, we have that E σ m 0,θ(X = E ( lc (X σ m lc,θ(x + dc (X σ m dc,θ ( = p 0 0 σ m lc,θ(xg0(x; r σ 0, γ 0 dx + Ḡr 0(; σ 0, γ 0 σ m dc,θ. (S- 3
4 Here g r 0(x is given by (5, and σ m lc,θ (x and σ m dc,θ are straightforwardly obtained from symbolic Mathematica differentiation. Setting σ = σ 0, γ = γ 0, and inserting numerical values for σ 0, γ 0, the integral in the second line of (S- may be obtained by numerical computation, and hence the entire expression in the righthand side may be computed numerically. Similarly the other three expressions in (S- may be computed numerically, and thus also D 0 (σ 0, σ 0, γ 0, γ 0. Using Mathematica to plot D 0 (σ 0, σ 0, γ 0, γ 0 as a function of σ 0 and γ 0 for (σ 0, γ 0 [σ l, σ u ] [γ l, γ u ], shows that it is strictly positive, see online supplement described in Section 6. These computations also show that the diagonal elements of θ E (m 0,θ (X are negative. Hence θ E (m 0,θ (X is strictly negative definite for θ = θ 0 K. Next, writing E S for conditional expectation given S, we as above have that θ E (m,θ(x = E θ m,θ(x = E (E S θ m,θ(x. Numerical computation of the determinant and diagonal elements of E θ,θ (X similarly showed that E θ,θ (X is non-negative definite for all (σ 0, γ 0, S [σ l, σ u ] [γ l, γ u ] [0, ]. For these computations we additionally used the ( following observations. Let D (σ 0, γ 0 S denote the determinant of E θ,θ (X for σ = σ 0, γ = γ 0. Then, D (σ 0, γ 0 S = D (σ 0 /( S, γ 0 /( S and hence computation over (σ 0, γ 0, S can be transferred to a computation with S = 0 and σ u =.5/( S. We did these computations from σ l = 0.05 up to σ u = 5, corresponding to S = 0.9. Since D (σ 0, γ 0 S 0 as σ 0 the values were negligible for S (0.9, 0.. The computations of the diagonal elements used the corresponding scaling properties of the diagonal elements. Next, θ E(m,θ (X is non-negative definite since integrals of non-negative definite matrices are non-negative definite. Finally, a sum of a strictly negative definite matrix and a non-negative definite matrix is strictly negative definite, and thus θ E(m θ (X = θ E (m 0,θ (X + m,θ (X evaluated at θ = θ 0 is strictly negative definite, and the second condition of Theorem 5.3 of van der Vaart (998 holds. To obtain last condition of Theorem 5.3 of van der Vaart (998, consistency of ˆθ, we use Theorem 5.9 of this reference. To check the first condition of that theorem, that n n i= m θ(x i converges to its mean uniformly for θ K is standard, and is left to the reader. The second condition follows by similar, but simpler, arguments as above (and in fact, the local part, that the condition is satisfied in some surrounding of θ 0 is already established above. As a final comment, since integration and differentiation may be freely exchanged, it is clear that one may use the observed information matrix, instead of the expected one, to measure estimation uncertainty. 4
5 3 Estimation algorithms We provide MatLab scripts that estimate the parameters of the distribution for exponential, gamma, and Weibull distributions, under full-parametric model, and for exponential and Generalized Pareto distributions, under semi-parametric model. One may also compute the corresponding 95% confidence intervals: except for exponential distribution in full-parametric case, where the explicit formula is available, confidence intervals are obtained using standard delta method where the inverse of the hessian matrix of the corresponding log-likelihood function is computed numerically. The scripts are located under MatLab Scripts Estimators folder. 4 Simulation study MatLab Scripts Simulation Study folder contains MatLab scripts that were used to produce Table, Table, and Table 3 of the Simulation Study. In each of the cases considered in the study, SimulateGlances*.m is used to simulate censored glances, and RunSimulation*.m produces 0, 000 estimates of the parameters of the distribution from simulated data (for various choices of the original distribution parameters and sample sizes, i.e. number of simulated censored glances Car data analysis MatLab Scripts Analysis of 00CarData folder contains MatLab scripts that were used to produce Figure 3 and Figure 4. Similar figures can be produced for any other choice of the parametric or semi-parametric models considered in the paper. 6 Numerical verification of asymptotic normality Wolfram Mathematica Scripts Asymptotic Normalit.nb is a numerical study where we verify the non-negative definiteness of the expected information matrix. Reference van der Vaart, A. (998. Asymptotic Statistics. Cambridge University Press, New York. 5
**BEGINNING OF EXAMINATION** A random sample of five observations from a population is:
**BEGINNING OF EXAMINATION** 1. You are given: (i) A random sample of five observations from a population is: 0.2 0.7 0.9 1.1 1.3 (ii) You use the Kolmogorov-Smirnov test for testing the null hypothesis,
More informationChapter 8: Sampling distributions of estimators Sections
Chapter 8 continued Chapter 8: Sampling distributions of estimators Sections 8.1 Sampling distribution of a statistic 8.2 The Chi-square distributions 8.3 Joint Distribution of the sample mean and sample
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 informationProbability. An intro for calculus students P= Figure 1: A normal integral
Probability An intro for calculus students.8.6.4.2 P=.87 2 3 4 Figure : A normal integral Suppose we flip a coin 2 times; what is the probability that we get more than 2 heads? Suppose we roll a six-sided
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 information8.1 Estimation of the Mean and Proportion
8.1 Estimation of the Mean and Proportion Statistical inference enables us to make judgments about a population on the basis of sample information. The mean, standard deviation, and proportions of a population
More informationIntroduction to Algorithmic Trading Strategies Lecture 8
Introduction to Algorithmic Trading Strategies Lecture 8 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationAn Improved Skewness Measure
An Improved Skewness Measure Richard A. Groeneveld Professor Emeritus, Department of Statistics Iowa State University ragroeneveld@valley.net Glen Meeden School of Statistics University of Minnesota Minneapolis,
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 informationGPD-POT and GEV block maxima
Chapter 3 GPD-POT and GEV block maxima This chapter is devoted to the relation between POT models and Block Maxima (BM). We only consider the classical frameworks where POT excesses are assumed to be GPD,
More informationChapter 4: Asymptotic Properties of MLE (Part 3)
Chapter 4: Asymptotic Properties of MLE (Part 3) Daniel O. Scharfstein 09/30/13 1 / 1 Breakdown of Assumptions Non-Existence of the MLE Multiple Solutions to Maximization Problem Multiple Solutions to
More informationAsymptotic results discrete time martingales and stochastic algorithms
Asymptotic results discrete time martingales and stochastic algorithms Bernard Bercu Bordeaux University, France IFCAM Summer School Bangalore, India, July 2015 Bernard Bercu Asymptotic results for discrete
More informationLecture 7: Computation of Greeks
Lecture 7: Computation of Greeks Ahmed Kebaier kebaier@math.univ-paris13.fr HEC, Paris Outline 1 The log-likelihood approach Motivation The pathwise method requires some restrictive regularity assumptions
More informationChapter 4: Commonly Used Distributions. Statistics for Engineers and Scientists Fourth Edition William Navidi
Chapter 4: Commonly Used Distributions Statistics for Engineers and Scientists Fourth Edition William Navidi 2014 by Education. This is proprietary material solely for authorized instructor use. Not authorized
More informationSYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data
SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu September 5, 2015
More informationLikelihood Methods of Inference. Toss coin 6 times and get Heads twice.
Methods of Inference Toss coin 6 times and get Heads twice. p is probability of getting H. Probability of getting exactly 2 heads is 15p 2 (1 p) 4 This function of p, is likelihood function. Definition:
More informationGamma Distribution Fitting
Chapter 552 Gamma Distribution Fitting Introduction This module fits the gamma probability distributions to a complete or censored set of individual or grouped data values. It outputs various statistics
More informationExercise. Show the corrected sample variance is an unbiased estimator of population variance. S 2 = n i=1 (X i X ) 2 n 1. Exercise Estimation
Exercise Show the corrected sample variance is an unbiased estimator of population variance. S 2 = n i=1 (X i X ) 2 n 1 Exercise S 2 = = = = n i=1 (X i x) 2 n i=1 = (X i µ + µ X ) 2 = n 1 n 1 n i=1 ((X
More informationUser s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions
User s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions Yacine Aït-Sahalia Department of Economics and Bendheim Center for Finance Princeton University and NBER This Version: July
More informationA New Hybrid Estimation Method for the Generalized Pareto Distribution
A New Hybrid Estimation Method for the Generalized Pareto Distribution Chunlin Wang Department of Mathematics and Statistics University of Calgary May 18, 2011 A New Hybrid Estimation Method for the GPD
More informationEcon 582 Nonlinear Regression
Econ 582 Nonlinear Regression Eric Zivot June 3, 2013 Nonlinear Regression In linear regression models = x 0 β (1 )( 1) + [ x ]=0 [ x = x] =x 0 β = [ x = x] [ x = x] x = β it is assumed that the regression
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationChapter 8 Statistical Intervals for a Single Sample
Chapter 8 Statistical Intervals for a Single Sample Part 1: Confidence intervals (CI) for population mean µ Section 8-1: CI for µ when σ 2 known & drawing from normal distribution Section 8-1.2: Sample
More informationMeasuring Financial Risk using Extreme Value Theory: evidence from Pakistan
Measuring Financial Risk using Extreme Value Theory: evidence from Pakistan Dr. Abdul Qayyum and Faisal Nawaz Abstract The purpose of the paper is to show some methods of extreme value theory through analysis
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 informationWeek 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals
Week 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg :
More informationSTRESS-STRENGTH RELIABILITY ESTIMATION
CHAPTER 5 STRESS-STRENGTH RELIABILITY ESTIMATION 5. Introduction There are appliances (every physical component possess an inherent strength) which survive due to their strength. These appliances receive
More informationAN EXTREME VALUE APPROACH TO PRICING CREDIT RISK
AN EXTREME VALUE APPROACH TO PRICING CREDIT RISK SOFIA LANDIN Master s thesis 2018:E69 Faculty of Engineering Centre for Mathematical Sciences Mathematical Statistics CENTRUM SCIENTIARUM MATHEMATICARUM
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER
Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.
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 informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
More informationPosterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties
Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where
More informationPoint Estimation. Principle of Unbiased Estimation. When choosing among several different estimators of θ, select one that is unbiased.
Point Estimation Point Estimation Definition A point estimate of a parameter θ is a single number that can be regarded as a sensible value for θ. A point estimate is obtained by selecting a suitable statistic
More information12 The Bootstrap and why it works
12 he Bootstrap and why it works For a review of many applications of bootstrap see Efron and ibshirani (1994). For the theory behind the bootstrap see the books by Hall (1992), van der Waart (2000), Lahiri
More informationEstimation Procedure for Parametric Survival Distribution Without Covariates
Estimation Procedure for Parametric Survival Distribution Without Covariates The maximum likelihood estimates of the parameters of commonly used survival distribution can be found by SAS. The following
More informationHomework Problems Stat 479
Chapter 10 91. * A random sample, X1, X2,, Xn, is drawn from a distribution with a mean of 2/3 and a variance of 1/18. ˆ = (X1 + X2 + + Xn)/(n-1) is the estimator of the distribution mean θ. Find MSE(
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Moments of a distribubon Measures of
More informationLikelihood-based Optimization of Threat Operation Timeline Estimation
12th International Conference on Information Fusion Seattle, WA, USA, July 6-9, 2009 Likelihood-based Optimization of Threat Operation Timeline Estimation Gregory A. Godfrey Advanced Mathematics Applications
More informationASSIGNMENT - 1, MAY M.Sc. (PREVIOUS) FIRST YEAR DEGREE STATISTICS. Maximum : 20 MARKS Answer ALL questions.
(DMSTT 0 NR) ASSIGNMENT -, MAY-04. PAPER- I : PROBABILITY AND DISTRIBUTION THEORY ) a) State and prove Borel-cantelli lemma b) Let (x, y) be jointly distributed with density 4 y(+ x) f( x, y) = y(+ x)
More informationCambridge University Press Risk Modelling in General Insurance: From Principles to Practice Roger J. Gray and Susan M.
adjustment coefficient, 272 and Cramér Lundberg approximation, 302 existence, 279 and Lundberg s inequality, 272 numerical methods for, 303 properties, 272 and reinsurance (case study), 348 statistical
More informationChapter 6. Importance sampling. 6.1 The basics
Chapter 6 Importance sampling 6.1 The basics To movtivate our discussion consider the following situation. We want to use Monte Carlo to compute µ E[X]. There is an event E such that P(E) is small but
More informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
More informationOn the Distribution and Its Properties of the Sum of a Normal and a Doubly Truncated Normal
The Korean Communications in Statistics Vol. 13 No. 2, 2006, pp. 255-266 On the Distribution and Its Properties of the Sum of a Normal and a Doubly Truncated Normal Hea-Jung Kim 1) Abstract This paper
More informationCase Study: Heavy-Tailed Distribution and Reinsurance Rate-making
Case Study: Heavy-Tailed Distribution and Reinsurance Rate-making May 30, 2016 The purpose of this case study is to give a brief introduction to a heavy-tailed distribution and its distinct behaviors in
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 informationME3620. Theory of Engineering Experimentation. Spring Chapter III. Random Variables and Probability Distributions.
ME3620 Theory of Engineering Experimentation Chapter III. Random Variables and Probability Distributions Chapter III 1 3.2 Random Variables In an experiment, a measurement is usually denoted by a variable
More informationRandom Variables and Probability Distributions
Chapter 3 Random Variables and Probability Distributions Chapter Three Random Variables and Probability Distributions 3. Introduction An event is defined as the possible outcome of an experiment. In engineering
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationContinuous random variables
Continuous random variables probability density function (f(x)) the probability distribution function of a continuous random variable (analogous to the probability mass function for a discrete random variable),
More informationChapter 5. Statistical inference for Parametric Models
Chapter 5. Statistical inference for Parametric Models Outline Overview Parameter estimation Method of moments How good are method of moments estimates? Interval estimation Statistical Inference for Parametric
More informationmay be of interest. That is, the average difference between the estimator and the truth. Estimators with Bias(ˆθ) = 0 are called unbiased.
1 Evaluating estimators Suppose you observe data X 1,..., X n that are iid observations with distribution F θ indexed by some parameter θ. When trying to estimate θ, one may be interested in determining
More informationFinite Memory and Imperfect Monitoring
Federal Reserve Bank of Minneapolis Research Department Finite Memory and Imperfect Monitoring Harold L. Cole and Narayana Kocherlakota Working Paper 604 September 2000 Cole: U.C.L.A. and Federal Reserve
More informationChapter 3 Common Families of Distributions. Definition 3.4.1: A family of pmfs or pdfs is called exponential family if it can be expressed as
Lecture 0 on BST 63: Statistical Theory I Kui Zhang, 09/9/008 Review for the previous lecture Definition: Several continuous distributions, including uniform, gamma, normal, Beta, Cauchy, double exponential
More informationPractice Exam 1. Loss Amount Number of Losses
Practice Exam 1 1. You are given the following data on loss sizes: An ogive is used as a model for loss sizes. Determine the fitted median. Loss Amount Number of Losses 0 1000 5 1000 5000 4 5000 10000
More informationQQ PLOT Yunsi Wang, Tyler Steele, Eva Zhang Spring 2016
QQ PLOT INTERPRETATION: Quantiles: QQ PLOT Yunsi Wang, Tyler Steele, Eva Zhang Spring 2016 The quantiles are values dividing a probability distribution into equal intervals, with every interval having
More information2.1 Random variable, density function, enumerative density function and distribution function
Risk Theory I Prof. Dr. Christian Hipp Chair for Science of Insurance, University of Karlsruhe (TH Karlsruhe) Contents 1 Introduction 1.1 Overview on the insurance industry 1.1.1 Insurance in Benin 1.1.2
More informationP VaR0.01 (X) > 2 VaR 0.01 (X). (10 p) Problem 4
KTH Mathematics Examination in SF2980 Risk Management, December 13, 2012, 8:00 13:00. Examiner : Filip indskog, tel. 790 7217, e-mail: lindskog@kth.se Allowed technical aids and literature : a calculator,
More informationWeb-based Supplementary Materials for. A space-time conditional intensity model. for invasive meningococcal disease occurence
Web-based Supplementary Materials for A space-time conditional intensity model for invasive meningococcal disease occurence by Sebastian Meyer 1,2, Johannes Elias 3, and Michael Höhle 4,2 1 Department
More informationLecture 10: Point Estimation
Lecture 10: Point Estimation MSU-STT-351-Sum-17B (P. Vellaisamy: MSU-STT-351-Sum-17B) Probability & Statistics for Engineers 1 / 31 Basic Concepts of Point Estimation A point estimate of a parameter θ,
More informationChapter 8. Introduction to Statistical Inference
Chapter 8. Introduction to Statistical Inference Point Estimation Statistical inference is to draw some type of conclusion about one or more parameters(population characteristics). Now you know that a
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More informationEE641 Digital Image Processing II: Purdue University VISE - October 29,
EE64 Digital Image Processing II: Purdue University VISE - October 9, 004 The EM Algorithm. Suffient Statistics and Exponential Distributions Let p(y θ) be a family of density functions parameterized by
More informationTELECOMMUNICATIONS ENGINEERING
TELECOMMUNICATIONS ENGINEERING STATISTICS 2012-2013 COMPUTER LAB SESSION # 3. PROBABILITY MODELS AIM: Introduction to most common discrete and continuous probability models. Characterization, graphical
More informationModel Uncertainty in Operational Risk Modeling
Model Uncertainty in Operational Risk Modeling Daoping Yu 1 University of Wisconsin-Milwaukee Vytaras Brazauskas 2 University of Wisconsin-Milwaukee Version #1 (March 23, 2015: Submitted to 2015 ERM Symposium
More informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS & STATISTICS SEMESTER /2013 MAS8304. Environmental Extremes: Mid semester test
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS & STATISTICS SEMESTER 2 2012/2013 Environmental Extremes: Mid semester test Time allowed: 50 minutes Candidates should attempt all questions. Marks for each question
More informationDRAFT. 1 exercise in state (S, t), π(s, t) = 0 do not exercise in state (S, t) Review of the Risk Neutral Stock Dynamics
Chapter 12 American Put Option Recall that the American option has strike K and maturity T and gives the holder the right to exercise at any time in [0, T ]. The American option is not straightforward
More information4-1. Chapter 4. Commonly Used Distributions by The McGraw-Hill Companies, Inc. All rights reserved.
4-1 Chapter 4 Commonly Used Distributions 2014 by The Companies, Inc. All rights reserved. Section 4.1: The Bernoulli Distribution 4-2 We use the Bernoulli distribution when we have an experiment which
More informationOPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF FINITE
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 005 Seville, Spain, December 1-15, 005 WeA11.6 OPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF
More informationOn Existence of Equilibria. Bayesian Allocation-Mechanisms
On Existence of Equilibria in Bayesian Allocation Mechanisms Northwestern University April 23, 2014 Bayesian Allocation Mechanisms In allocation mechanisms, agents choose messages. The messages determine
More informationGenerating Random Numbers
Generating Random Numbers Aim: produce random variables for given distribution Inverse Method Let F be the distribution function of an univariate distribution and let F 1 (y) = inf{x F (x) y} (generalized
More informationA Convenient Way of Generating Normal Random Variables Using Generalized Exponential Distribution
A Convenient Way of Generating Normal Random Variables Using Generalized Exponential Distribution Debasis Kundu 1, Rameshwar D. Gupta 2 & Anubhav Manglick 1 Abstract In this paper we propose a very convenient
More informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question
More information[D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright
Faculty and Institute of Actuaries Claims Reserving Manual v.2 (09/1997) Section D7 [D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright 1. Introduction
More informationAn Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process
Computational Statistics 17 (March 2002), 17 28. An Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process Gordon K. Smyth and Heather M. Podlich Department
More informationProbability and Statistics
Kristel Van Steen, PhD 2 Montefiore Institute - Systems and Modeling GIGA - Bioinformatics ULg kristel.vansteen@ulg.ac.be CHAPTER 3: PARAMETRIC FAMILIES OF UNIVARIATE DISTRIBUTIONS 1 Why do we need distributions?
More informationExam M Fall 2005 PRELIMINARY ANSWER KEY
Exam M Fall 005 PRELIMINARY ANSWER KEY Question # Answer Question # Answer 1 C 1 E C B 3 C 3 E 4 D 4 E 5 C 5 C 6 B 6 E 7 A 7 E 8 D 8 D 9 B 9 A 10 A 30 D 11 A 31 A 1 A 3 A 13 D 33 B 14 C 34 C 15 A 35 A
More informationFinancial Risk Forecasting Chapter 9 Extreme Value Theory
Financial Risk Forecasting Chapter 9 Extreme Value Theory Jon Danielsson 2017 London School of Economics To accompany Financial Risk Forecasting www.financialriskforecasting.com Published by Wiley 2011
More informationModeling the dependence between a Poisson process and a continuous semimartingale
1 / 28 Modeling the dependence between a Poisson process and a continuous semimartingale Application to electricity spot prices and wind production modeling Thomas Deschatre 1,2 1 CEREMADE, University
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 informationContinuous Distributions
Quantitative Methods 2013 Continuous Distributions 1 The most important probability distribution in statistics is the normal distribution. Carl Friedrich Gauss (1777 1855) Normal curve A normal distribution
More informationMLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models
MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models Matthew Dixon and Tao Wu 1 Illinois Institute of Technology May 19th 2017 1 https://papers.ssrn.com/sol3/papers.cfm?abstract
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More informationMonitoring Processes with Highly Censored Data
Monitoring Processes with Highly Censored Data Stefan H. Steiner and R. Jock MacKay Dept. of Statistics and Actuarial Sciences University of Waterloo Waterloo, N2L 3G1 Canada The need for process monitoring
More information1 Explicit Euler Scheme (or Euler Forward Scheme )
Numerical methods for PDE in Finance - M2MO - Paris Diderot American options January 2018 Files: https://ljll.math.upmc.fr/bokanowski/enseignement/2017/m2mo/m2mo.html We look for a numerical approximation
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 informationInterpolation. 1 What is interpolation? 2 Why are we interested in this?
Interpolation 1 What is interpolation? For a certain function f (x we know only the values y 1 = f (x 1,,y n = f (x n For a point x different from x 1,,x n we would then like to approximate f ( x using
More informationModelling Environmental Extremes
19th TIES Conference, Kelowna, British Columbia 8th June 2008 Topics for the day 1. Classical models and threshold models 2. Dependence and non stationarity 3. R session: weather extremes 4. Multivariate
More informationAmath 546/Econ 589 Univariate GARCH Models
Amath 546/Econ 589 Univariate GARCH Models Eric Zivot April 24, 2013 Lecture Outline Conditional vs. Unconditional Risk Measures Empirical regularities of asset returns Engle s ARCH model Testing for ARCH
More informationConstruction and behavior of Multinomial Markov random field models
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2010 Construction and behavior of Multinomial Markov random field models Kim Mueller Iowa State University Follow
More informationA Bayesian Control Chart for the Coecient of Variation in the Case of Pooled Samples
A Bayesian Control Chart for the Coecient of Variation in the Case of Pooled Samples R van Zyl a,, AJ van der Merwe b a PAREXEL International, Bloemfontein, South Africa b University of the Free State,
More informationCourse information FN3142 Quantitative finance
Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken
More informationModelling Environmental Extremes
19th TIES Conference, Kelowna, British Columbia 8th June 2008 Topics for the day 1. Classical models and threshold models 2. Dependence and non stationarity 3. R session: weather extremes 4. Multivariate
More informationAsymmetric Type II Compound Laplace Distributions and its Properties
CHAPTER 4 Asymmetric Type II Compound Laplace Distributions and its Properties 4. Introduction Recently there is a growing trend in the literature on parametric families of asymmetric distributions which
More informationGlobal convergence rate analysis of unconstrained optimization methods based on probabilistic models
Math. Program., Ser. A DOI 10.1007/s10107-017-1137-4 FULL LENGTH PAPER Global convergence rate analysis of unconstrained optimization methods based on probabilistic models C. Cartis 1 K. Scheinberg 2 Received:
More informationFolded- and Log-Folded-t Distributions as Models for Insurance Loss Data
Folded- and Log-Folded-t Distributions as Models for Insurance Loss Data Vytaras Brazauskas University of Wisconsin-Milwaukee Andreas Kleefeld University of Wisconsin-Milwaukee Revised: September 009 (Submitted:
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 informationPakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks
Pakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks Spring 2009 Main question: How much are patents worth? Answering this question is important, because it helps
More informationSTART HERE: Instructions. 1 Exponential Family [Zhou, Manzil]
START HERE: Instructions Thanks a lot to John A.W.B. Constanzo and Shi Zong for providing and allowing to use the latex source files for quick preparation of the HW solution. The homework was due at 9:00am
More informationLecture Notes 6. Assume F belongs to a family of distributions, (e.g. F is Normal), indexed by some parameter θ.
Sufficient Statistics Lecture Notes 6 Sufficiency Data reduction in terms of a particular statistic can be thought of as a partition of the sample space X. Definition T is sufficient for θ if the conditional
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More information