Modelling the risk process
|
|
- Della Hardy
- 5 years ago
- Views:
Transcription
1 Modelling the risk process Krzysztof Burnecki Hugo Steinhaus Center Wroc law University of Technology hugo
2 Modelling the risk process 1 Risk process If (Ω, F, P) is a probability space carrying: 1. a point process {N t } t 0, i.e. an integer valued stochastic process with N 0 = 0 a.s., N t < for each t < and nondecreasing realizations, and 2. an independent sequence {X k } k=1 variables, of positive i.i.d. random then the risk process {R t } t 0 is given by: R t = u + c(t) N t i=1 X i. (1)
3 Modelling the risk process 2 Risk process cont. where: R t = u + c(t) N t i=1 X i, {N t } t 0 is the claim arrival point process, {X k } k=1 is an independent claim sequence of positive i.i.d. random variables with common mean µ, u is a nonnegative constant representing the initial capital of the company, c(t) is the premium function.
4 Modelling the risk process 3 Simulation of the risk process The simulation of the risk process R t or the aggregated claim process { N t i=1 X i} reduces to modeling: the claim arrival point process {N t }, the claim size sequence {X k }, Both processes are assumed to be independent, hence can be simulated independently of each other.
5 Modelling the risk process 4 Generating continuous random variables Problem: Generate sample of a random variable X with a given density f. (The sample is called a random variate) The simulation of the claim sizes can be done by: inverse transform method, convolution method, composition method, acceptance-rejection method.
6 Modelling the risk process 5 Inverse transform method Assumption: We can generate U, i.e., uniform (0, 1) random variable. Consider a continuous r.v. having distribution function F. Theorem: For any continuous distribution F the r.v. X defined by X = F 1 (U) has distribution F, where F 1 (u) = inf x : F (x) = u. Thus, the algorithm is: Step 1: Generate a uniform random variable U. Step 2: Set X = F 1 (U).
7 Modelling the risk process 6 Inverse transform method. Example Exponential r.v.: F (x) = 1 exp( βx). If 1 exp( βx) = u, then x = 1 β log(1 u). X = 1 β log(1 U). Algorithm: Step 1: Generate a uniform random variable U. Step 2: Set X = 1 β log(u).
8 Modelling the risk process 7 Rejection Method Suppose we can simulate a r.v. with density function g(x). We would like to simulate a r.v. with density function f(y). Let c be a constant such that f(y)/g(y) c, for all y. Then the rejection algorithm is: Step 1: Generate Y having density g and generate U. Step 2: If U f(y) cg(y), then X = Y ; else go to step 1. Theorem. The r.v. X generated by the rejection method has density function f. Moreover, the number of iterations needed to obtain X is a geometric r.v. with mean c.
9 Modelling the risk process 8 Convolution method Suppose X is a sum of independent random variables Z 1, Z 2,... Z m, i.e. X = Z 1 + Z Z m, where Z i F i and are all independent. Algorithm: Step 1: Generate m random numbers U 1, U 2,... U m. Step 2: Inverse transform method: Z i = F 1 (U i ). Step 3: Set X = m i=1 Z i.
10 Modelling the risk process 9 Convolution method. Example Generate a sample from Erlang(β; m) distribution. Algorithm: Step 1: Generate m random numbers U 1, U 2,... U m. Step 2: Inverse transform method: Z i = 1 β log(u i). Step 3: Set X = m i=1 Z i.
11 Modelling the risk process 10 Composition method Suppose that either the distribution F X or the probability density f X can be represented either of the following two forms: (1) F X (x) = p 1 F Y1 (x) + p 2 F Y2 (x) +... p m F Ym (x) (2) f X (x) = p 1 f Y1 (x) + p 2 f Y2 (x) +... p m f Ym (x) where p 1,..., p m are non-negative and sum to one (so that they form a probability mass function).
12 Modelling the risk process 11 Composition method cont. Then, assuming that the Y i s are relatively easily to generate, we can generate X as follows: Step 1: Generate a discrete random variable I on {1,..., m}, where P (I = j) = p j for 1 j m. Step 2: Generate Y I from F YI (or f YI ). Step 3: Return X = Y I.
13 Modelling the risk process 12 Simulation of the claim arrival point process {N t } is simulated either via: the arrival times {T i }, i.e. moments when the ith claim occurs, or the inter-arrival times (or waiting times) W i = T i T i 1, i.e. the time periods between successive claims. The prominent scenarios for {N t }, are given by: the homogeneous Poisson process (HPP), the non-homogeneous Poisson process (NHPP), the mixed Poisson process, the Cox process (or doubly stochastic Poisson process), the renewal process.
14 Modelling the risk process 13 Homogeneous Poisson process A continuous-time stochastic process {N t : t 0} is a (homogeneous) Poisson process with intensity (or rate) λ > 0 if: 1. {N t } is a point process, and 2. the times between events are independent and identically distributed with an exponential(λ) distribution, i.e. exponential with mean 1/λ.
15 Modelling the risk process 14 Simulation of the HPP Successive arrival times T 1, T 2,..., T n of the (homogeneous) Poisson process can be generated by the following algorithm: Step 1: set T 0 = 0 Step 2: for i = 1, 2,..., n do Step 2a: generate an exponential random variable E with intensity λ Step 2b: set T i = T i 1 + E
16 Modelling the risk process 15 Simulation of the HPP cont. Given that N(t) = n, the n occurrence times T 1, T 2,..., T n have the same distributions at the order statistics corresponding to n i.i.d. random variables uniformly distributed on the interval (0, t]. Hence the arrival times T 1, T 2,..., T n of the HPP on the interval (0, t] can be generated as follows: Step 1: Generate a Poisson random variable N with intensity λ. Let N = n. Step 2: Generate n random variables U i distributed uniformly on (0, 1), i.e. U i U(0, 1), i = 1, 2,..., n. Step 3: (T 1, T 2,..., T n ) = t sort{u 1, U 2,..., U n }.
17 Modelling the risk process 16 Simulation of the HPP cont. Since EN t = λt, it is natural to define the premium function as: c(t) = ct = (1 + θ)µλt, where µ = EX k and θ > 0 is the relative safety loading which guarantees survival of the insurance company. With such a choice of the premium function we obtain the classical form of the risk process: R t = u + ct N t i=1 X i.
18 Modelling the risk process 17 Non-homogeneous Poisson process The non-homogeneous Poisson process (NHPP) can be thought of as a Poisson process with a variable intensity defined by the deterministic intensity (rate) function λ(t). A NHPP can model situations, where claim occurrence epochs are likely to depend on the time of the year or of the week. The increments of a NHPP do not have to be stationary. When λ(t) = λ, the NHPP reduces to the HPP with intensity λ. The simulation of the non-homogeneous Poisson process is slightly more complicated than the homogeneous one.
19 Modelling the risk process 18 Simulation of the NHPP integration method The increment of a NHPP with rate function λ(t) is distributed as a Poisson random variable with intensity λ = t s λ(u)du. Hence, the distribution function F s of the waiting time W s satisfies: F s (t) = P (W s t) = 1 P (W s > t) = 1 P (N s+t N s = 0) = ( s+t ) ( t ) = 1 exp λ(u)du = 1 exp λ(s + v)dv. s If we can find a formula for the inverse Fs 1 then for each s we can easily generate W s using the inverse transform method. 0
20 Modelling the risk process 19 Simulation of the NHPP integration method cont. The resulting algorithm can be summarized as follows: Step 1: set T 0 = 0 Step 2: for i = 1, 2,..., n do Step 2a: generate a random variable U distributed uniformly on (0, 1), i.e. U U(0, 1) Step 2b: set T i = T i 1 + Fs 1 (U)
21 Modelling the risk process 20 Simulation of the NHPP thinning Suppose that there exists a constant λ such that λ(t) λ for all t. Let T 1, T 2,... be the arrival times of a HPP with intensity λ. Accept the ith arrival time with probability λ(ti )/λ, independently of all other arrivals, as part of the thinned process (hence the name of the method). The sequence T 1, T 2,... of the accepted arrival times forms a sequence of the arrival times of a NHPP with rate function λ(t). The algorithm amounts to rejecting (hence the alternative name rejection method) or accepting a particular arrival as part of the thinned process.
22 Modelling the risk process 21 Simulation of the NHPP thinning cont. The resulting algorithm reads as follows: Step 1: set T 0 = 0 and T = 0 Step 2: for i = 1, 2,..., n do Step 2a: generate an exponential random variable E with intensity λ Step 2b: set T = T + E Step 2c: generate a random variable U U(0, 1) Step 2d: if U > λ(t )/λ then return to step 2a ( reject the arrival time) else set T i = T ( accept the arrival time)
23 Modelling the risk process 22 Simulation of the NHPP cont. Given that N(t) = n, the n occurence times T 1, T 2,..., T n have the same distributions at the order statistics corresponding to n independent random variables distributed on the interval (0, t], each with the common density function f(v) = λ(v)/ t λ(u)du, v (0, t]. 0 Hence the arrival times T 1, T 2,..., T n of the NHPP on the interval (0, t] can be generated as follows: Step 1: Generate a Poisson random variable N with intensity t λ(u)du. Let N = n. 0 Step 2: Generate n random variables V i, i = 1, 2,... n given by the densityf(v) = λ(v)/ t 0 λ(u)du. Step 3: (T 1, T 2,..., T n ) = sort{v 1, V 2,..., V n }.
24 Modelling the risk process 23 Simulation of the NHPP cont. Since EN t = t λ(s)ds, it is natural to define the premium function 0 in the non-homogeneous case as: c(t) = (1 + θ)µ t 0 λ(s)ds, where µ = EX k and θ > 0 is the relative safety loading. Then the risk process takes the form: R t = u + (1 + θ)µ t 0 λ(s)ds N t i=1 X i.
25 Modelling the risk process 24 Simulation of the NHPP cont. a = 1, b = 0.01 (blue line), b = 0.1 (red), b = 1 (green) Linear intensity (a+b*t) Seasonal intensity (a+b*sin(2*pi*t)) N(t) N(t) t t STFrisk01.xpl STFrisk02.xpl
26 Modelling the risk process 25 Mixed Poisson process In many situations the portfolio of an insurance company is diversified in the sense that the risks associated with different groups of policy holders are significantly different. For example, in motor insurance we might want to make a difference between male and female drivers or between drivers of different age. We would then assume that the claims come from a heterogeneous group of clients, each one of them generating claims according to a Poisson distribution with the intensity varying from one group to another.
27 Modelling the risk process 26 Mixed Poisson process cont. If Ñ is a HPP with intensity 1 and Λ is a positive random variable independent of Ñ, then the process N = Ñ Λ = (Ñ(Λt)) t is called a mixed Poisson process. The random variable Λ is called a structure variable. A mixed Poisson process has stationary increments, however the independent increments condition is violated. The most common choice for the distribution of the structure variable Λ is the gamma distribution. In such a case the mixed Poisson proces is called a negative binomial process or Pólya process.
28 Modelling the risk process 27 Mixed Poisson process cont. In the mixed Poisson process the distribution of {N t } is given by a mixture of Poisson processes. Conditioning on the extrinsic random variable Λ, the process {N t } behaves like a HPP. Hence, the process can be generated in the following way: first a realization of a non-negative random variable Λ is generated, conditioned upon its realization, {N t } as a HPP with that realization as its intensity is constructed.
29 Modelling the risk process 28 Simulation of the mixed Poisson process Making the algorithm more formal we can write: Step 1: generate a realization λ of the random intensity Λ Step 2: set T 0 = 0 Step 3: for i = 1, 2,..., n do Step 3a: generate an exponential random variable E with intensity λ Step 3b: set T i = T i 1 + E
30 Modelling the risk process 29 Simulation of the mixed Poisson process cont. Since for each t the claim numbers {N t } up to time t are Poisson with intensity Λt, in the mixed case it is reasonable to consider the premium function of the form: c(t) = (1 + θ)µλt, where µ = EX k and θ > 0 is the relative safety loading. Then the risk process takes the form: R t = u + (1 + θ)µλt N t i=1 X i.
31 Modelling the risk process 30 Cox process The Cox process, or doubly stochastic Poisson process, provides flexibility by letting the intensity not only depend on time but also by allowing it to be a stochastic process. Cox processes seem to form a natural class for modeling risk and size fluctuations. IF Ñ is a HPP with intensity 1 and {Λ(t)} is a stochastic process with Λ = 0, non-decreasing sample paths and idependent of Ñ, then the process N = Ñ Λ = (Ñ(Λ)) is called a Cox process or doubly stochastic Poisson process.
32 Modelling the risk process 31 Cox process cont. The intensity process {Λ(t)} is used to generate another process {N t } by acting as its intensity. That is, {N t } is a Poisson process conditional on {Λ(t)} which itself is a stochastic process. If {Λ(t)} is deterministic, then {N t } is a NHPP. This property suggests a simulation method.
33 Modelling the risk process 32 Simulation of the Cox process Step 1: generate a realization λ(t) of the intensity process {Λ(t)} for a sufficiently large time period Step 2: set λ = max {λ(t)} Step 3: set T 0 = 0 and T = 0 Step 4: for i = 1, 2,..., n do Step 4a: generate an exponential random variable E with intensity λ Step 4b: set T = T + E Step 4c: generate a random variable U U(0, 1) Step 4d: if U > λ(t )/λ then return to step 4a ( reject the arrival time) else set T i = T ( accept the arrival time)
34 Modelling the risk process 33 Simulation of the Cox process cont. The premium function is a generalization of the former premium functions: c(t) = (1 + θ)µ t 0 Λ(s)ds, where µ = EX k and θ > 0 is the relative safety loading. Then the risk process takes the form: R t = u + (1 + θ)µ t 0 Λ(s)ds N t i=1 X i.
35 Modelling the risk process 34 Renewal process If the waiting times W i are i.i.d. and nonnegative then the resulting sequence is a renewal process. Note, that the HPP is a renewal process with exponentially distributed inter-arrival times. Hence, we can generate the arrival times of a renewal process by: Step 1: set T 0 = 0 Step 2: for i = 1, 2,..., n do Step 2a: generate a random variable X with an assumed distribution function F Step 2b: set T i = T i 1 + X
36 Modelling the risk process 35 Renewal process cont. For renewal claim arrival processes a constant premium rate allows for a constant safety loading. Let {N t } be a renewal process and assume that EW 1 = 1/λ <. Then the premium function is defined in a natural way as: c(t) = (1 + θ)µλt, like in the homogeneous Poisson process case.
37 Modelling the risk process 36 Empirical analysis Danish fire losses recorded by Copenhagen Re. Losses in profits connected with fires Loss sizes Lognormal, Pareto, Burr, gamma, Weibull and mixture of two exponentials distributions Claim counting process Homogeneous and nonhomogeneous process
38 Modelling the risk process 37 Losses (DKK million) Log(1-F(x)) Time Losses (DKK million) Figure 1: Left panel: Illustration of the major Danish fire losses adjusted for inflation. Right panel: Logarithm of the right tails of the empirical claim size distribution function (thick blue solid line) together with lognormal (red dotted line) and Burr (thin black solid line) fits.
39 Modelling the risk process 38 Danish fire losses. Loss sizes Losses (DKK million) Log(1-F(x)) Time Losses (DKK million) STFrisk03.xpl STFrisk04.xpl
40 Modelling the risk process 39 d.f.: Lognormal Pareto Burr Weibull Para- µ = α = α = α = meters: σ = λ = e6 λ = e7 λ = e-5 τ = χ KS CM AD
41 Modelling the risk process 40 Danish fire losses cont. Claim counting process Mean-value function ACF Time (years) Time lag (qtr) STFrisk05.xpl STFrisk06.xpl
42 Modelling the risk process 41 Danish fire losses cont. Claim counting process The data reveals no seasonality A clear increasing trend can be observed in the number of quarterly losses We tested different exponential and polynomial functional forms A simple linear intensity function λ(s) = a + bs yielded the best fit Applying a least squares procedure we arrived at the values: a = and b = 7.57.
43 Modelling the risk process 42 Danish fire losses cont. The fitted risk process We consider a hypothetical scenario where the insurance company insures losses resulting from fire damage The company s initial capital is assumed to be u = 100 million kr The relative safety loading used is θ = 0.5 We chose two models of the risk process: a non-homogeneous Poisson process with lognormal claim sizes and a non-homogeneous Poisson process with Burr claim sizes.
44 Modelling the risk process 43 Danish fire losses cont. The fitted risk process Capital (DKK million) Capital (DKK million) Time (years) Time (years) STFrisk07.xpl STFrisk08.xpl
HSC Research Report. Ruin Probability in Finite Time HSC/10/04. Krzysztof Burnecki* Marek Teuerle*
HSC Research Report HSC/1/4 Ruin Probability in Finite Time Krzysztof Burnecki* Marek Teuerle* * Hugo Steinhaus Center, Wrocław University of Technology, Poland Hugo Steinhaus Center Wrocław University
More informationPoisson Processes. Stochastic Processes. Feb UC3M
Poisson Processes Stochastic Processes UC3M Feb. 2012 Exponential random variables A random variable T has exponential distribution with rate λ > 0 if its probability density function can been written
More informationCAS Exam MAS-1. Howard Mahler. Stochastic Models
CAS Exam MAS-1 Howard Mahler Stochastic Models 2019 Stochastic Models, HCM 12/31/18, Page 1 The slides are in the same order as the sections of my study guide. Section # Section Name 1 Introduction 2 Exponential
More informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More informationf X (x) = λe λx, , x 0, k 0, λ > 0 Γ (k) f X (u)f X (z u)du
11 COLLECTED PROBLEMS Do the following problems for coursework 1. Problems 11.4 and 11.5 constitute one exercise leading you through the basic ruin arguments. 2. Problems 11.1 through to 11.13 but excluding
More information1 Poisson processes, and Compound (batch) Poisson processes
Copyright c 2007 by Karl Sigman 1 Poisson processes, and Compound (batch) Poisson processes 1.1 Point Processes Definition 1.1 A simple point process ψ = {t n : n 1} is a sequence of strictly increasing
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
More informationIE 303 Discrete-Event Simulation
IE 303 Discrete-Event Simulation 1 L E C T U R E 5 : P R O B A B I L I T Y R E V I E W Review of the Last Lecture Random Variables Probability Density (Mass) Functions Cumulative Density Function Discrete
More informationPoint Process Control
Point Process Control The following note is based on Chapters I, II and VII in Brémaud s book Point Processes and Queues (1981). 1 Basic Definitions Consider some probability space (Ω, F, P). A real-valued
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationStochastic Renewal Processes in Structural Reliability Analysis:
Stochastic Renewal Processes in Structural Reliability Analysis: An Overview of Models and Applications Professor and Industrial Research Chair Department of Civil and Environmental Engineering University
More informationSlides 8: Statistical Models in Simulation
Slides 8: Statistical Models in Simulation Purpose and Overview The world the model-builder sees is probabilistic rather than deterministic: Some statistical model might well describe the variations. An
More informationPractice Exam 1. (A) (B) (C) (D) (E) You are given the following data on loss sizes:
Practice Exam 1 1. Losses for an insurance coverage have the following cumulative distribution function: F(0) = 0 F(1,000) = 0.2 F(5,000) = 0.4 F(10,000) = 0.9 F(100,000) = 1 with linear interpolation
More information2905 Queueing Theory and Simulation PART IV: SIMULATION
2905 Queueing Theory and Simulation PART IV: SIMULATION 22 Random Numbers A fundamental step in a simulation study is the generation of random numbers, where a random number represents the value of a random
More information3. Poisson Processes (12/09/12, see Adult and Baby Ross)
3. Poisson Processes (12/09/12, see Adult and Baby Ross) Exponential Distribution Poisson Processes Poisson and Exponential Relationship Generalizations 1 Exponential Distribution Definition: The continuous
More informationFinite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims
Proceedings of the World Congress on Engineering 29 Vol II WCE 29, July 1-3, 29, London, U.K. Finite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims Tao Jiang Abstract
More informationExam C Solutions Spring 2005
Exam C Solutions Spring 005 Question # The CDF is F( x) = 4 ( + x) Observation (x) F(x) compare to: Maximum difference 0. 0.58 0, 0. 0.58 0.7 0.880 0., 0.4 0.680 0.9 0.93 0.4, 0.6 0.53. 0.949 0.6, 0.8
More informationSPRING 2007 EXAM C SOLUTIONS
SPRING 007 EXAM C SOLUTIONS Question #1 The data are already shifted (have had the policy limit and the deductible of 50 applied). The two 350 payments are censored. Thus the likelihood function is L =
More informationRandom variables. DS GA 1002 Probability and Statistics for Data Science.
Random variables DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Motivation Random variables model numerical quantities
More informationProbability Distributions Columns (a) through (d)
Discrete Probability Distributions Columns (a) through (d) Probability Mass Distribution Description Notes Notation or Density Function --------------------(PMF or PDF)-------------------- (a) (b) (c)
More informationRMSC 2001 Introduction to Risk Management
RMSC 2001 Introduction to Risk Management Tutorial 4 (2011/12) 1 February 20, 2012 Outline: 1. Failure Time 2. Loss Frequency 3. Loss Severity 4. Aggregate Claim ====================================================
More informationStability of the Defect Renewal Volterra Integral Equations
19th International Congress on Modelling and Simulation, Perth, Australia, 12 16 December 211 http://mssanz.org.au/modsim211 Stability of the Defect Renewal Volterra Integral Equations R. S. Anderssen,
More information3 Continuous Random Variables
Jinguo Lian Math437 Notes January 15, 016 3 Continuous Random Variables Remember that discrete random variables can take only a countable number of possible values. On the other hand, a continuous random
More informationChapter 5: Generating Random Numbers from Distributions
Chapter 5: Generating Random Numbers from Distributions See Reading Assignment OR441-DrKhalid Nowibet 1 Review 1 Inverse Transform Generate a number u i between 0 and 1 (one U-axis) and then find the corresponding
More informationLecturer: Olga Galinina
Renewal models Lecturer: Olga Galinina E-mail: olga.galinina@tut.fi Outline Reminder. Exponential models definition of renewal processes exponential interval distribution Erlang distribution hyperexponential
More informationOptional Stopping Theorem Let X be a martingale and T be a stopping time such
Plan Counting, Renewal, and Point Processes 0. Finish FDR Example 1. The Basic Renewal Process 2. The Poisson Process Revisited 3. Variants and Extensions 4. Point Processes Reading: G&S: 7.1 7.3, 7.10
More informationDr. Maddah ENMG 617 EM Statistics 10/15/12. Nonparametric Statistics (2) (Goodness of fit tests)
Dr. Maddah ENMG 617 EM Statistics 10/15/12 Nonparametric Statistics (2) (Goodness of fit tests) Introduction Probability models used in decision making (Operations Research) and other fields require fitting
More informationContinuous-time Markov Chains
Continuous-time Markov Chains Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ October 23, 2017
More informationerrors every 1 hour unless he falls asleep, in which case he just reports the total errors
I. First Definition of a Poisson Process A. Definition: Poisson process A Poisson Process {X(t), t 0} with intensity λ > 0 is a counting process with the following properties. INDEPENDENT INCREMENTS. For
More informationManual for SOA Exam MLC.
Chapter 10. Poisson processes. Section 10.5. Nonhomogenous Poisson processes Extract from: Arcones Fall 2009 Edition, available at http://www.actexmadriver.com/ 1/14 Nonhomogenous Poisson processes Definition
More informationRMSC 2001 Introduction to Risk Management
RMSC 2001 Introduction to Risk Management Tutorial 4 (2011/12) 1 February 20, 2012 Outline: 1. Failure Time 2. Loss Frequency 3. Loss Severity 4. Aggregate Claim ====================================================
More informationModeling Recurrent Events in Panel Data Using Mixed Poisson Models
Modeling Recurrent Events in Panel Data Using Mixed Poisson Models V. Savani and A. Zhigljavsky Abstract This paper reviews the applicability of the mixed Poisson process as a model for recurrent events
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationExponential Distribution and Poisson Process
Exponential Distribution and Poisson Process Stochastic Processes - Lecture Notes Fatih Cavdur to accompany Introduction to Probability Models by Sheldon M. Ross Fall 215 Outline Introduction Exponential
More informationParameter Estimation
Parameter Estimation Chapters 13-15 Stat 477 - Loss Models Chapters 13-15 (Stat 477) Parameter Estimation Brian Hartman - BYU 1 / 23 Methods for parameter estimation Methods for parameter estimation Methods
More informationGenerating Random Variates 2 (Chapter 8, Law)
B. Maddah ENMG 6 Simulation /5/08 Generating Random Variates (Chapter 8, Law) Generating random variates from U(a, b) Recall that a random X which is uniformly distributed on interval [a, b], X ~ U(a,
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 29, 2014 Introduction Introduction The world of the model-builder
More informationCreating New Distributions
Creating New Distributions Section 5.2 Stat 477 - Loss Models Section 5.2 (Stat 477) Creating New Distributions Brian Hartman - BYU 1 / 18 Generating new distributions Some methods to generate new distributions
More informationDiscrete Distributions Chapter 6
Discrete Distributions Chapter 6 Negative Binomial Distribution section 6.3 Consider k r, r +,... independent Bernoulli trials with probability of success in one trial being p. Let the random variable
More information(b) What is the variance of the time until the second customer arrives, starting empty, assuming that we measure time in minutes?
IEOR 3106: Introduction to Operations Research: Stochastic Models Fall 2006, Professor Whitt SOLUTIONS to Final Exam Chapters 4-7 and 10 in Ross, Tuesday, December 19, 4:10pm-7:00pm Open Book: but only
More informationRandom Processes. DS GA 1002 Probability and Statistics for Data Science.
Random Processes DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Aim Modeling quantities that evolve in time (or space)
More informationGB2 Regression with Insurance Claim Severities
GB2 Regression with Insurance Claim Severities Mitchell Wills, University of New South Wales Emiliano A. Valdez, University of New South Wales Edward W. (Jed) Frees, University of Wisconsin - Madison UNSW
More informationUpper and lower bounds for ruin probability
Upper and lower bounds for ruin probability E. Pancheva,Z.Volkovich and L.Morozensky 3 Institute of Mathematics and Informatics, the Bulgarian Academy of Sciences, 3 Sofia, Bulgaria pancheva@math.bas.bg
More informationMethod of Moments. which we usually denote by X or sometimes by X n to emphasize that there are n observations.
Method of Moments Definition. If {X 1,..., X n } is a sample from a population, then the empirical k-th moment of this sample is defined to be X k 1 + + Xk n n Example. For a sample {X 1, X, X 3 } the
More informationEDRP lecture 7. Poisson process. Pawe J. Szab owski
EDRP lecture 7. Poisson process. Pawe J. Szab owski 2007 Counting process Random process fn t ; t 0g is called a counting process, if N t is equal total number of events that have happened up to moment
More informationPROBABILITY DISTRIBUTIONS
Review of PROBABILITY DISTRIBUTIONS Hideaki Shimazaki, Ph.D. http://goo.gl/visng Poisson process 1 Probability distribution Probability that a (continuous) random variable X is in (x,x+dx). ( ) P x < X
More informationA STATISTICAL TEST FOR MONOTONIC AND NON-MONOTONIC TREND IN REPAIRABLE SYSTEMS
A STATISTICAL TEST FOR MONOTONIC AND NON-MONOTONIC TREND IN REPAIRABLE SYSTEMS Jan Terje Kvaløy Department of Mathematics and Science, Stavanger University College, P.O. Box 2557 Ullandhaug, N-491 Stavanger,
More informationLecture Notes on Risk Theory
Lecture Notes on Risk Theory February 2, 21 Contents 1 Introduction and basic definitions 1 2 Accumulated claims in a fixed time interval 3 3 Reinsurance 7 4 Risk processes in discrete time 1 5 The Adjustment
More informationLecture 7. Poisson and lifetime processes in risk analysis
Lecture 7. Poisson and lifetime processes in risk analysis Jesper Rydén Department of Mathematics, Uppsala University jesper.ryden@math.uu.se Statistical Risk Analysis Spring 2014 Example: Life times of
More informationContinuous Probability Spaces
Continuous Probability Spaces Ω is not countable. Outcomes can be any real number or part of an interval of R, e.g. heights, weights and lifetimes. Can not assign probabilities to each outcome and add
More informationIntroduction to Rare Event Simulation
Introduction to Rare Event Simulation Brown University: Summer School on Rare Event Simulation Jose Blanchet Columbia University. Department of Statistics, Department of IEOR. Blanchet (Columbia) 1 / 31
More informationVaR vs. Expected Shortfall
VaR vs. Expected Shortfall Risk Measures under Solvency II Dietmar Pfeifer (2004) Risk measures and premium principles a comparison VaR vs. Expected Shortfall Dependence and its implications for risk measures
More informationEE/CpE 345. Modeling and Simulation. Fall Class 5 September 30, 2002
EE/CpE 345 Modeling and Simulation Class 5 September 30, 2002 Statistical Models in Simulation Real World phenomena of interest Sample phenomena select distribution Probabilistic, not deterministic Model
More informationMixed Poisson Models for Claim Numbers. 1.1 Introduction Poisson Modelling for the Number of Claims
1 Mixed Poisson Models for Claim Numbers 1.1 Introduction 1.1.1 Poisson Modelling for the Number of Claims In view of the economic importance of motor third party liability insurance in industrialized
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationThe exponential distribution and the Poisson process
The exponential distribution and the Poisson process 1-1 Exponential Distribution: Basic Facts PDF f(t) = { λe λt, t 0 0, t < 0 CDF Pr{T t) = 0 t λe λu du = 1 e λt (t 0) Mean E[T] = 1 λ Variance Var[T]
More informationContinuous Distributions
A normal distribution and other density functions involving exponential forms play the most important role in probability and statistics. They are related in a certain way, as summarized in a diagram later
More informationApplying the proportional hazard premium calculation principle
Applying the proportional hazard premium calculation principle Maria de Lourdes Centeno and João Andrade e Silva CEMAPRE, ISEG, Technical University of Lisbon, Rua do Quelhas, 2, 12 781 Lisbon, Portugal
More informationUniversity Of Calgary Department of Mathematics and Statistics
University Of Calgary Department of Mathematics and Statistics Hawkes Seminar May 23, 2018 1 / 46 Some Problems in Insurance and Reinsurance Mohamed Badaoui Department of Electrical Engineering National
More informationFeuille d exercices n 1 : Poisson processes. 1. Give a necessary and su cient condition for having
Processus de Poisson et méthodes actuarielles (25-26) Feuille d exercices n : Poisson processes Exercise. Let ( n,n ) be a sequence of IID (independent and identically distributed) non-negative random
More informationCDA6530: Performance Models of Computers and Networks. Chapter 3: Review of Practical Stochastic Processes
CDA6530: Performance Models of Computers and Networks Chapter 3: Review of Practical Stochastic Processes Definition Stochastic process X = {X(t), t2 T} is a collection of random variables (rvs); one rv
More information1 Delayed Renewal Processes: Exploiting Laplace Transforms
IEOR 6711: Stochastic Models I Professor Whitt, Tuesday, October 22, 213 Renewal Theory: Proof of Blackwell s theorem 1 Delayed Renewal Processes: Exploiting Laplace Transforms The proof of Blackwell s
More informationRuin probabilities in multivariate risk models with periodic common shock
Ruin probabilities in multivariate risk models with periodic common shock January 31, 2014 3 rd Workshop on Insurance Mathematics Universite Laval, Quebec, January 31, 2014 Ruin probabilities in multivariate
More informationCDA5530: Performance Models of Computers and Networks. Chapter 3: Review of Practical
CDA5530: Performance Models of Computers and Networks Chapter 3: Review of Practical Stochastic Processes Definition Stochastic ti process X = {X(t), t T} is a collection of random variables (rvs); one
More informationMixture distributions in Exams MLC/3L and C/4
Making sense of... Mixture distributions in Exams MLC/3L and C/4 James W. Daniel Jim Daniel s Actuarial Seminars www.actuarialseminars.com February 1, 2012 c Copyright 2012 by James W. Daniel; reproduction
More informationS6880 #7. Generate Non-uniform Random Number #1
S6880 #7 Generate Non-uniform Random Number #1 Outline 1 Inversion Method Inversion Method Examples Application to Discrete Distributions Using Inversion Method 2 Composition Method Composition Method
More informationIntroduction to repairable systems STK4400 Spring 2011
Introduction to repairable systems STK4400 Spring 2011 Bo Lindqvist http://www.math.ntnu.no/ bo/ bo@math.ntnu.no Bo Lindqvist Introduction to repairable systems Definition of repairable system Ascher and
More informationA Note On The Erlang(λ, n) Risk Process
A Note On The Erlangλ, n) Risk Process Michael Bamberger and Mario V. Wüthrich Version from February 4, 2009 Abstract We consider the Erlangλ, n) risk process with i.i.d. exponentially distributed claims
More informationChapter 4: CONTINUOUS RANDOM VARIABLES AND PROBABILITY DISTRIBUTIONS
Chapter 4: CONTINUOUS RANDOM VARIABLES AND PROBABILITY DISTRIBUTIONS Part 4: Gamma Distribution Weibull Distribution Lognormal Distribution Sections 4-9 through 4-11 Another exponential distribution example
More informationNon-Life Insurance: Mathematics and Statistics
Exercise sheet 1 Exercise 1.1 Discrete Distribution Suppose the random variable N follows a geometric distribution with parameter p œ (0, 1), i.e. ; (1 p) P[N = k] = k 1 p if k œ N \{0}, 0 else. (a) Show
More informationChapter 2. Poisson Processes. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan
Chapter 2. Poisson Processes Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan Outline Introduction to Poisson Processes Definition of arrival process Definition
More informationSeverity Models - Special Families of Distributions
Severity Models - Special Families of Distributions Sections 5.3-5.4 Stat 477 - Loss Models Sections 5.3-5.4 (Stat 477) Claim Severity Models Brian Hartman - BYU 1 / 1 Introduction Introduction Given that
More informationModèles de dépendance entre temps inter-sinistres et montants de sinistre en théorie de la ruine
Séminaire de Statistiques de l'irma Modèles de dépendance entre temps inter-sinistres et montants de sinistre en théorie de la ruine Romain Biard LMB, Université de Franche-Comté en collaboration avec
More informationThe story of the film so far... Mathematics for Informatics 4a. Continuous-time Markov processes. Counting processes
The story of the film so far... Mathematics for Informatics 4a José Figueroa-O Farrill Lecture 19 28 March 2012 We have been studying stochastic processes; i.e., systems whose time evolution has an element
More informationMahler s Guide to Stochastic Models
HCMSA-S07-3-1A, Stochastic Models, 11/22/06, Page 1 Mahler s Guide to Stochastic Models Copyright 2007 by Howard C. Mahler. Information in bold or sections whose title is in bold are more important for
More informationOrder book modeling and market making under uncertainty.
Order book modeling and market making under uncertainty. Sidi Mohamed ALY Lund University sidi@maths.lth.se (Joint work with K. Nyström and C. Zhang, Uppsala University) Le Mans, June 29, 2016 1 / 22 Outline
More informationA hidden semi-markov model for the occurrences of water pipe bursts
A hidden semi-markov model for the occurrences of water pipe bursts T. Economou 1, T.C. Bailey 1 and Z. Kapelan 1 1 School of Engineering, Computer Science and Mathematics, University of Exeter, Harrison
More informationAn Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems
An Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems M.P. Kaminskiy and V.V. Krivtsov Abstract This paper introduces a simple index that helps to assess the degree of aging
More informationNon-Life Insurance: Mathematics and Statistics
ETH Zürich, D-MATH HS 08 Prof. Dr. Mario V. Wüthrich Coordinator Andrea Gabrielli Non-Life Insurance: Mathematics and Statistics Solution sheet 9 Solution 9. Utility Indifference Price (a) Suppose that
More informationStochastic Simulation
Stochastic Simulation Jan-Pieter Dorsman 1 & Michel Mandjes 1,2,3 1 Korteweg-de Vries Institute for Mathematics, University of Amsterdam 2 CWI, Amsterdam 3 Eurandom, Eindhoven University of Amsterdam,
More informationSince D has an exponential distribution, E[D] = 0.09 years. Since {A(t) : t 0} is a Poisson process with rate λ = 10, 000, A(0.
IEOR 46: Introduction to Operations Research: Stochastic Models Chapters 5-6 in Ross, Thursday, April, 4:5-5:35pm SOLUTIONS to Second Midterm Exam, Spring 9, Open Book: but only the Ross textbook, the
More informationRuin Probabilities of a Discrete-time Multi-risk Model
Ruin Probabilities of a Discrete-time Multi-risk Model Andrius Grigutis, Agneška Korvel, Jonas Šiaulys Faculty of Mathematics and Informatics, Vilnius University, Naugarduko 4, Vilnius LT-035, Lithuania
More informationStochastic Processes
Stochastic Processes Stochastic Process Non Formal Definition: Non formal: A stochastic process (random process) is the opposite of a deterministic process such as one defined by a differential equation.
More informationTest Problems for Probability Theory ,
1 Test Problems for Probability Theory 01-06-16, 010-1-14 1. Write down the following probability density functions and compute their moment generating functions. (a) Binomial distribution with mean 30
More informationIEOR 3106: Second Midterm Exam, Chapters 5-6, November 7, 2013
IEOR 316: Second Midterm Exam, Chapters 5-6, November 7, 13 SOLUTIONS Honor Code: Students are expected to behave honorably, following the accepted code of academic honesty. You may keep the exam itself.
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 6
MATH 56A: STOCHASTIC PROCESSES CHAPTER 6 6. Renewal Mathematically, renewal refers to a continuous time stochastic process with states,, 2,. N t {,, 2, 3, } so that you only have jumps from x to x + and
More informationSolution: The process is a compound Poisson Process with E[N (t)] = λt/p by Wald's equation.
Solutions Stochastic Processes and Simulation II, May 18, 217 Problem 1: Poisson Processes Let {N(t), t } be a homogeneous Poisson Process on (, ) with rate λ. Let {S i, i = 1, 2, } be the points of the
More informationLecture 4a: Continuous-Time Markov Chain Models
Lecture 4a: Continuous-Time Markov Chain Models Continuous-time Markov chains are stochastic processes whose time is continuous, t [0, ), but the random variables are discrete. Prominent examples of continuous-time
More informationSTAT 414: Introduction to Probability Theory
STAT 414: Introduction to Probability Theory Spring 2016; Homework Assignments Latest updated on April 29, 2016 HW1 (Due on Jan. 21) Chapter 1 Problems 1, 8, 9, 10, 11, 18, 19, 26, 28, 30 Theoretical Exercises
More informationA New Fractional Process: A Fractional Non-homogeneous Poisson Process Enrico Scalas
A New Fractional Process: A Fractional Non-homogeneous Poisson Process Enrico Scalas University of Sussex Joint work with Nikolai Leonenko (Cardiff University) and Mailan Trinh Fractional PDEs: Theory,
More informationSolutions to the Spring 2015 CAS Exam ST
Solutions to the Spring 2015 CAS Exam ST (updated to include the CAS Final Answer Key of July 15) There were 25 questions in total, of equal value, on this 2.5 hour exam. There was a 10 minute reading
More informationSTAT 418: Probability and Stochastic Processes
STAT 418: Probability and Stochastic Processes Spring 2016; Homework Assignments Latest updated on April 29, 2016 HW1 (Due on Jan. 21) Chapter 1 Problems 1, 8, 9, 10, 11, 18, 19, 26, 28, 30 Theoretical
More informationIE 581 Introduction to Stochastic Simulation
1. List criteria for choosing the majorizing density r (x) when creating an acceptance/rejection random-variate generator for a specified density function f (x). 2. Suppose the rate function of a nonhomogeneous
More informationMFM Practitioner Module: Quantitative Risk Management. John Dodson. September 23, 2015
MFM Practitioner Module: Quantitative Risk Management September 23, 2015 Mixtures Mixtures Mixtures Definitions For our purposes, A random variable is a quantity whose value is not known to us right now
More informationGraduate Econometrics I: What is econometrics?
Graduate Econometrics I: What is econometrics? Yves Dominicy Université libre de Bruxelles Solvay Brussels School of Economics and Management ECARES Yves Dominicy Graduate Econometrics I: What is econometrics?
More informationNorthwestern University Department of Electrical Engineering and Computer Science
Northwestern University Department of Electrical Engineering and Computer Science EECS 454: Modeling and Analysis of Communication Networks Spring 2008 Probability Review As discussed in Lecture 1, probability
More informationChapter 5. Statistical Models in Simulations 5.1. Prof. Dr. Mesut Güneş Ch. 5 Statistical Models in Simulations
Chapter 5 Statistical Models in Simulations 5.1 Contents Basic Probability Theory Concepts Discrete Distributions Continuous Distributions Poisson Process Empirical Distributions Useful Statistical Models
More informationPolynomial approximation of mutivariate aggregate claim amounts distribution
Polynomial approximation of mutivariate aggregate claim amounts distribution Applications to reinsurance P.O. Goffard Axa France - Mathematics Institute of Marseille I2M Aix-Marseille University 19 th
More informationSolutions to Homework Discrete Stochastic Processes MIT, Spring 2011
Exercise 1 Solutions to Homework 6 6.262 Discrete Stochastic Processes MIT, Spring 2011 Let {Y n ; n 1} be a sequence of rv s and assume that lim n E[ Y n ] = 0. Show that {Y n ; n 1} converges to 0 in
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More information