Recursive Calculation of Finite Time Ruin Probabilities Under Interest Force
|
|
- Angelica Alicia Daniels
- 5 years ago
- Views:
Transcription
1 Recursive Calculation of Finite Time Ruin Probabilities Under Interest Force Rui M R Cardoso and Howard R Waters Abstract In this paper we consider a classical insurance surplus process affected by a constant interest force. We present numerical algorithms for the calculation of finite time ruin probabilities using a discrete time Markov chain to approximate the risk process. Based on this method, upper and lower bounds are also obtained. 1 Introduction We are interested in the following insurance surplus process: U(t) =ue δt + c s t t 0 e δ(t v) ds(v), t 0, where time is measured in some basic unit, which we shall refer to as a year, and: U(t) is the insurer s surplus at time t, u is the insurer s initial surplus, s t δ is the constant force of interest p.a. earned on the insurer s reserves, c is the rate of premium income p.a. for the insurer, = (exp(δt) 1)/δ, so that c s t is the accumulated premium income in the time interval [0,t], and, Department of Actuarial Mathematics & Statistics, Heriot-Watt University, Edinburgh EH14 4AS, Great Britain Support from CMA and Departamento de Matemática, FCT, Universidade Nova de Lisboa and Fundação para a Ciência e a Tecnologia - FCT/POCTI is gratefully acknowledged 1
2 S(v) is the aggregate claims in [0,v]. We assume that {S(t)} t=0 is a compound Poisson process. The Poisson parameter, and hence the expected number of claims each year, is λ, the (i.i.d.) individual claim amounts have cdf G(x) and pdf g(x) with G(0) = 0, so that all claim amounts are positive. The insurer s rate of premium income can be written c =(1+η)λm 1, where m k is the k th moment about 0 of the individual claim amount distribution and η is the insurer s premium loading factor. We assume, without loss of generality, that: λ =1=m 1 The time to ruin, T, for this process is ined as follows: { inf{t : U(t) < 0} T = if U(t) 0 for all t>0 and the finite time ruin probability, ψ(u, t), is ined as ψ(u, t) =Pr[T t]. (1) Our objective in this paper is to derive recursive numerical algorithms for the (approximate) calculation of, and bounds for, ψ(u, t). In the special case δ =0,{U(t)} t=0 is the classical insurance surplus process and recursive algorithms for the calculation of ψ(u, t) have been developed by De Vylder & Goovaerts (1988) and Dickson & Waters (1991). The key to these algorithms is that the time/amount plane is replaced by a rectangular grid and the continuous (time and amount) surplus process is approximated by a discrete (time and amount) process whose state space is a set of equally spaced points on the amounts axis. The probabilistic properties of this discrete approximating process are easier to derive than the properties of the original process. One of the key simplifying features of this approximating process is that it can be scaled so that in a single unit of time the premium income, and hence the maximum increase in the surplus, is 1. Dickson & Waters (1999) derived recursive algorithms for the calculation of ψ(u, t) for the general model (δ >0) using the discretisation ideas developed by De Vylder & Goovaerts (1988) and Dickson & Waters (1991) for the special case (δ = 0). A complication with this approach was that in any fixed interval of time the maximum increase in the surplus depends not only on the rate of premium income but also on the level of the surplus at the start of the interval. A consequence of this complication was that the resulting algorithms were very calculation-intensive and hence heavy on computer run time. 2
3 In this paper we use a different approach to derive algorithms to calculate ψ(u, t). Our approach approximates or bounds the continuous surplus process by discrete time Markov chains with countable state spaces, an idea which has been used by Dickson & Gray (1984) and Cardoso & Egídio dos Reis (2002) to calculate ruin probabilities for the classical surplus process, i.e. the special case δ =0. Numerical values for the probability of ultimate ruin, i.e. ψ(u, ), for the general model can be calculated using, for example, methods derived by Sundt & Teugels (1995) or by De Vylder (1996). In the special case of exponentially distributed claim amounts, Segerdahl derived an analytic formula for this probability. See Dickson & Waters (1999) for more details of these results. Numerical algorithms for the calculation of values of ψ(u, t) have been presented by Dickson & Waters (1999) and Brekelmans & De Waegenaere (2001). Exact solutions for special cases have been given by Albrecher et al (2001). In Section 2 we derive algorithms for the calculation of upper and lower bounds for ψ(u, t). We also present an algorithm for the (approximate) calculation of this probability. In Section 3 we describe a truncation procedure which can be used to reduce the number of calculations required to produce the approximation or the bounds, while keeping the error introduced by this procedure within specified limits. In Section 4 we give some numerical examples. 2 The Markov Chain Algorithm 2.1 The accumulated aggregate claim distribution Let h be some small unit of time and let S(h) denote the accumulated aggregate claims in the time interval [0,h], with cdf F (x), so that: [ Nh ] F (x) =P[ S(h) x]=p e δ(h Ti) X i x where: i=1 N h X i T i is the number of claims in [0,h], which has a Poisson distribution with mean λh, is the amount of the i th claim, and, is the time of the i th claim. 3
4 We are going to approximate and bound the continuous time surplus process by discrete time processes ined at the time points 0,h,2h,...,sothe smaller the value of h, the better the approximation and the bounds are likely to be. In what follows, it will be convenient for t to be an integer multiple of h, so we will assume that t/h is some positive integer, say, K. For n =0,1,2,..., let p n denote the probability of the event (N h = n). The following result follows from formulae in Karlin & Taylor (1975, pp ) and Ross (1996, pp ): F (x) = p G n n (x) (2) n=0 where G n is the n th convolution of G, and: h G(x) = 1 G ( xe δ(h u)) du h 0 Note that G is the cdf of X i e δ(h τ i), where τ i is uniformly distributed on (0,h). Now let β be some large positive number. We can approximate G by a discrete distribution function, Gd, and bound it by discrete distribution functions G d and G d, so that: G d (x) G(x) G d (x) with all three of these discrete distributions having masses only at the points 0, 1/β, 2/β,... For example, Gd could be constructed using the method of De Vylder & Goovaerts (1988), and the( two bounding distributions could be constructed by concentrating the mass G( n+1 ) G( ) n) at either (n +1)/β β β or n/β, as appropriate. The larger the value of β, the closer the bounding distributions will be, and the closer the approximating distribution is likely to be, to G. Formula (2) shows that F has a compound Poisson distribution. We can now use the recursion formula due to Panjer (1981), in conjunction with formula (2), replacing G by G d, Gd and G d in turn, to calculate cdfs F d, Fd and F d such that: F d (x) F (x) F d (x) and F d (x) F (x) 4
5 This generalised version of Panjer s recursion goes back to Boogaert & De Waegenaere (1990) and was used by Brekelmans & De Waegenaere (2001). From the way in which F d ( F d ) has been constructed, it can be seen that a discrete time surplus process ined at the time points 0,h,2h,..., starting at u, having the same rate of premium income and interest on reserves as the continuous time surplus process and for which the accumulated aggregate claims in a time interval of length h has cdf F d ( F d ) will, with probability 1, always be below (above) the surplus process at these time points. 2.2 Construction of the Markov Chains We are going to construct three discrete time Markov chains with countable state spaces which approximate and bound below and above the continuous time surplus process until ruin occurs. These discrete time processes, denoted {Z n } n=0, {Z n } n=0 and {Z n } n=0, will be ined at times 0,h,2h,..., where h is as in Section 2.1. All three chains start from u, the initial value of the surplus process, so that: P[Z 0 = u] =P[Z 0 =u]=p[z 0 =u] = 1 (3) and for all chains zero is an absorbing state (although they are allowed to start from u = 0), so that for n 1: P[Z n+1 =0 Z n =0]=P[Z n+1 =0 Z n =0]=P[Z n+1 =0 Z n =0]=1 From step 1 onwards, the state space for the two bounding chains is a set of non-negative and increasing numbers {x j },j =0,1,... ined as follows: x 0 = 0 x j+1 = x j e δh + c s h for j =0,1,... Note that if the surplus process were to start at x j, then in a time interval of length h it would reach x j+1 if there were no claims in this interval. For i, n 1 and 1 j i + 1, the transition probabilities for {Z n } n=0 are specified as follows: P[Z n+1 =0 Z n =x i ] = 1 F d (x i+1 ) P[Z n+1 = x j Z n = x i ] = F d ((x i+1 x j 1 ) ) F d ((x i+1 x j ) ) To complete the inition of {Z n } n=0, we need to specify the transition probabilities at the first step. If u = x i for some non-negative integer i, these 5
6 transition probilities are as above with n = 0. Suppose u (x i 1,x i ). Then for 1 j i +1: P[Z 1 =0] = 1 F d ((ue δh + c s h ) ) P[Z 1 = x j ] = F d ((ue δh + c s h x j 1 ) ) F d ((ue δh + c s h x j ) ) Intuitively, this Markov chain behaves like a discrete time surplus process ined at the time points 0,h,2h,..., starting at u, having the same rate of premium income and interest on reserves as the continuous time surplus process and for which the accumulated aggregate claims in a time interval of length h has cdf F d with the extra features that when this process takes a value in an interval (x j 1,x j ], its value is adjusted upwards to x j and that the process stops, taking the value 0, following the first negative value of the surplus. The Markov chain {Z n } n=0 is ined similarly to {Z n } n=0, the differences being that the cdf F d is used in its inition, rather than F d, and at each step values in the interval [x j 1,x j ) are adjusted downwards to x j 1. Formally, this chain is specified as follows. For i, n 1 and 1 j i +1: P[Z n+1 =0 Z n =x i ] = 1 F d (x i+1 x 1 ) P[Z n+1 = x j Z n = x i ] = F d (x i+1 x j ) F d (x i+1 x j+1 ) If u = x i for some non-negative integer i, the first step transition probabilities are as above with n =0. Ifu (x i,x i+1 ), then for 1 j i +1: P[Z 1 = 0] = 1 F d (ue δh + c s h x 1 ) P[Z 1 = x j ] = F d (ue δh + c s h x j ) F d (ue δh + c s h x j+1 ) The Markov chain {Z n } n=0 is ined as follows. For j =1,2,..., let M j =(x j 1 +x j )/2, so that M j is the mid-point of the interval x j 1 x j, and ine M 0 =0. Forn 1the state space for Z n is the set of real numbers 0,M 1,M 2,... Intuitively, this chain will move from M i to M j if, starting from M i, the surplus process time h later has a value in the interval (x j 1,x j ]. More precisely, for i, n 1 and 1 j i: P[Z n+1 = M i+1 Z n = M i ] = F d ((M i e δh + c s h x i ) ) P[Z n+1 = M j Z n = M i ] = F d ((M i e δh + c s h x j 1 ) ) F d ((M i e δh + c s h x j ) ) P[Z n+1 =0 Z n =M i ] = 1 F d ((M i e δh + c s h ) ) Define the non-negative integer γ(u)tobe0ifu= 0 and to be i if u (x i 1,x i ], for some positive integer i. Then the first step transition probabilities are 6
7 ined as follows for 1 j i: P[Z 1 =0] = 1 F d ((ue δh + c s h ) ) P[Z 1 = M γ(u)+1 ] = F d ((ue δh + c s h x γ(u) ) ) P[Z 1 = M j ] = F d ((ue δh + c s h x j 1 ) ) F d ((ue δh + c s h x j ) ) For the Markov chain {Z n } n=0 (respectively, {Z n } n=0 and {Z n } n=0), let P and P (u) (respectively P and P (u) and P and P (u) ) denote the matrix of transition probabilities after the first step and the vector of transition probabilities at the first step. For integers n 1 and j 0, let p (n) (u) j (respectively p(n) (u) j and p (n) (u) j ) denote the probability of the event (Z n = M j ) (respectively (Z n = x j ) and (Z n = x j )). For positive integers m and n, let p (n) ij (respectively p (n) and ij p (n) ij ) denote the probability of the event (Z m+n = M j Z m = M i ) (respectively (Z m+n = x j Z m = x i ) and (Z m+n = x j Z m = x i )). It follows from the inition of x 1 that if at any time τ, U(τ) < 0, then U(τ + h) cannot be greater than or equal to x 1. It then follows from our constructions that: (Z K =0) (U(τ)<0 for some τ,0 <τ t) (Z K =0) (Recall that K = t/h.) Hence: p (K) ψ(u, t) p(k) and: p (K) (u )0 ψ(u, t) The probability p (K) is the first entry in the vector P (u) P K 1. A recursive formula for this probability, in the spirit of a formula due to De Vylder & Goovaerts (1988, formula (9)) can be written as follows: γ(u)+1 p (K) =p + j=1 γ(u)+1 p (u) j p (K 1) j 0 = j=0 p (u) j p (K 1) j 0 (4) Similar recursive formulae can be written down for p (K) and p(k). 2.3 Comments on the method In principle, the vectors P (u), P (u), P (u) and matrices P, P, P have infinite dimension. However, since the Markov chains can move up at most one state 7
8 in a single time interval, only a finite number of entries in each vector can have non-zero values. For example, P (u) can have non-zero values in only its first (γ(u) + 2) entries and P (u) P n can have non-zero entries only in its first (γ(u) n) entries. In essence, our method replaces the continuous time surplus process by discrete time processes constrained to take values within a countable set. Expressed in this way, this is similar to the method of Dickson & Waters (1999). An important difference between the two methods is that the countable set of values used by Dickson & Waters (1999) was equally spaced whereas in our method the intervals (M j M j 1 ) and (x j x j 1 ) are increasing with j, so that the discretisation by amount becomes coarser as the surplus increases. A coarser discretisation results in fewer states for our processes and hence faster computations, whereas a finer discretisation results in greater accuracy. Our method is intuitively appealing in that it uses a finer discretisation for values of the surplus more likely to lead to ruin. The larger the value of β and the smaller the value of h, the more accurate our approximations are likely to be. These two parameters can be set independently, but in our numerical examples we have always set: h = 1 (1 + η)β so that the premium income in a time interval of length h is 1, before accumulating with interest. This relationship means that an increase in β automatically results in a decrease in h. Once the parameter h has been set, the value of x j, j =1,2,..., cannot be greater than: x j 1 e δh + c s h If it were, Z n, starting from x j 1, could never reach x j and so this process could never exceed its starting value. 3 Truncation of the numerical algorithms In their paper, De Vylder & Goovaerts (1988) show how the number of calculations in their recursive formula corresponding to our formula (4) can be reduced by a truncation procedure in such a way that the error introduced can be controlled. We can use these ideas in connection with formula (4). We will do this in two stages. 8
9 3.1 Truncation - Stage 1 Let ɛ (0 <ɛ<1) be fixed. For each positive integer i, ine the integer valued function α(i) as follows: { } j α(i) = min j : p ik ɛ Similarly, for the initial surplus u, we ine α((u)) as follows: { } j α((u)) = min j : p (u) k ɛ We ine ɛ p (u)j and ɛ p ij for i =1,2,... and j =0,1,... as follows: k=0 k=0 ɛp (u) j ɛp (u) j ɛp ij ɛp ij = 0 if j<α((u)) = p (u) j if j α((u)) = 0 if j<α(i) = p ij if j α(i) Intuitively, for each state i, and the initial surplus u, we are setting to zero the probability of reaching a state if the sum of the probabilities of reaching that, or a lower, state is less than ɛ. We now ine ɛ p (n) ij and ɛ p (n) recursively for n =2,3,... as follows: ɛp (n) i 0 ɛp (n) = = i+1 j=α(i) γ(u)+1 j=α((u)) ɛp ij ɛ p (n 1) j 0 ɛp (u) j ɛ p (n 1) j 0 Calculating ɛ p (n) (u) j should be quicker than calculating p(n) (u) j since the summation for the former is, possibly, over a shorter range of values. Our first result tells us that using ɛ p (n) (u) j in place of p(n) (u) j introduces an error which can be quantified and which grows only linearly with n. RESULT 1: For n, i =1,2,... we have: 0 p (n) i 0 ɛp (n) i 0 nɛ (5) 0 p (n) ɛp (n) nɛ (6) 9
10 Proof: Consider first the case n =1. Ifα((u)) = 0, then ɛ p = p. On the other hand, if α((u)) > 0, then p <ɛand ɛ p = 0. In either case, formula (6) holds. The proof of formula (5) for n = 1 is similar. Now suppose the result holds for some positive integer n. We have: ɛp (n+1) = p (n+1) = γ(u)+1 j=α((u)) γ(u)+1 j=0 ɛp (u) j ɛ p (n) j 0 p (u) j p (n) j 0 We have: p ɛ p by construction p (u) j ɛ p (u) j by construction p (n) j 0 ɛ p (n) j 0 by assumption Hence: Also: p (n+1) ɛ p (n+1) = p (n+1) γ(u)+1 j=α((u)) ɛ p (n+1) p (u) j [p (n) j 0 ɛp (n) nɛ + ɛ =(n+1)ɛ j 0 ]+ α((u)) 1 j=0 p (u) j p (n) j 0 Hence, formula (6) holds by induction. The proof of formula (5) is similar. 3.2 Truncation - Stage 2 We can extend the ideas of the previous section by, additionally, setting to zero the probability of going from any state i to 0 if this is less than ɛ. We proceed recursively as follows. For n, i =1,2,...: ɛ p i0 ɛ p i 0 ɛ p (n+1) i 0 ɛ p (n+1) i 0 = ɛ p i 0 if this is ɛ = 0 otherwise i+1 = ɛp ij ɛ p (n) j0 j=α(i) = 0 otherwise if this is ɛ 10
11 Finally, ine: and for n =1,2,...: ɛ p = ɛ p as before ɛ p (n+1) = γ(u)+1 j=α((u)) ɛp (u) j ɛ p (n) j 0 RESULT 2: For n, i =1,2,... we have: 0 p (n) i 0 ɛ p (n) i 0 2nɛ (7) 0 p (n) ɛ p (n) 2nɛ (8) Proof: Consider first the case n = 1. Formula (8) follows from formula (6) since: ɛ p = ɛ p Formula (7) is proved using Result 1 as follows: p i 0 ɛ p i 0 ɛ p i 0 p i 0 ɛ p i 0 = p i 0 ɛ p i 0 + ɛ p i 0 ɛ p i 0 ɛ + ɛ =2ɛ Now suppose the result holds for some positive integer n. We will show that formula (8) holds for n + 1, and hence, by induction, for all n. The proof of formula (7) is similar. p (n+1) = ɛ p (n+1) = γ(u)+1 j=0 γ(u)+1 j=α((u)) p (u) j p (n) j 0 ɛp (u) j ɛ p (n) j 0 We have: p ɛ p by construction p (u) j ɛ p (u) j by construction p (n) j 0 ɛ p (n) j 0 by assumption Hence: p (n+1) ɛ p (n+1) 11
12 Also: p (n+1) ɛ p (n+1) = γ(u)+1 j=α((u)) p (u) j [p (n) j 0 ɛ p (n) 2nɛ + ɛ 2(n +1)ɛ j 0 ]+ α((u)) 1 j=0 p (u) j p (n) j Comments on the truncation In this section we showed how to reduce the number of calculations in the computation of approximations to the ruin probability and, consequently, the computational effort by using a truncation procedure. This procedure can also be applied in the case of the lower and upper bounds obtaining similar results to the ones presented. From the results derived we can also control the error due to this truncation. Furthermore, we conclude that the Markov chain methods presented are stable. In Figures 1 and 2 we compare the approximations of ψ(u, t) without the truncation procedure (solid line) and with the truncation procedure (dashed line). The dotted line gives the values of the latter approximations plus Kɛ. We considered β = 100 and η =0.1. In Figure 1 the individual claim amounts are exponentially distributed with mean 1 and in Figure 2 they follow a Pareto(3,2) distribution. In the first case 30 t 40, u = 5 and δ =0.1. For the second one 35 t 40, u = 10 and δ = For both cases we set ɛ = (so, at most, the error is ). We observe that the two approximations are close. Of course, we can control the closeness of the values, but the point is that the difference between the two approximations is much smaller than the error bound. We found this pattern in other cases studied. 4 Examples We will now present numerical results in cases where the individual claim amount have an exponential distribution with mean 1 or a Pareto(3,2) distribution. In all examples we considered η =0.1 and β = 100. For the case δ = 0 we do not present any values although the approximations computed by the Markov chain method have the same accuracy as the approximations given by the algorithm of Dickson & Waters (1991, section 8). The values obtained for different combinations of u, δ and t are summarized in the following tables. The lower and upper bounds are denoted 12
13 without truncation with truncation truncation plus error Figure 1: Approximations to ψ(5,t) with and without truncation for exponential(1) claim size distribution, with δ =
14 without truncation with truncation truncation plus error Figure 2: Approximations to ψ(10,t) with and without truncation for Pareto(3,2) claim size distribution, with δ =
15 by Lu and Ul, respectively, and the approximations denoted by A. Where it is possible we compare these numbers with the values found in Dickson & Waters (1999), denoted by DW, and in Brekelmans & De Waegenaere (2001). These latter authors also split the time horizon into smaller intervals of equal length and they derived an algorithm to determine a lower and a upper bound, denoted by LB and UB, respectively, for ψ(u, t) by assuming that premium income is received, respectively, at the beginnning and at the end of each interval. Then by averaging these bounds they get an approximation (AVG) to the ruin probability. They also get simulated values (SIM) of the ruin probabilities (see Brekelmans & De Waegenaere (2001, subsection 4.1)) In tables 1 and 2, the values from Brekelmans & De Waegenaere (2001) were calculated with a step size h =0.01 while in tables 4 and 5, h = Note that the values produced by Dickson & Waters (1999) and the ones obtained by our methods were calculated considering β = 100, i.e., a step size h = Table 1: Approximations to, and bounds for, ψ(0,t) - exponential claims δ t =1 t=5 t=10 t=20 t=40 Lu Ul A DW LB UB AVG SIM Lu Ul A DW Lu Ul A DW
16 Table 2: Approximations to, and bounds for, ψ(5,t) - exponential claims δ t =1 t=5 t=10 t=20 t=40 Lu Ul A DW LB UB AVG SIM Lu Ul A DW Lu Ul A DW Table 3: Approximations to, and bounds for, ψ(10,t) - exponential claims δ t =1 t=5 t=10 t=20 t=40 Lu Ul A DW Lu Ul A DW Lu Ul A DW
17 Table 4: Approximations to, and bounds for, ψ(0,t) - Pareto claims δ t =1 t=5 t=10 t=20 t=40 LB UB AVG SIM Lu Ul A DW LB UB AVG SIM Lu Ul A DW LB UB AVG SIM Lu Ul A DW LB UB AVG SIM Lu Ul A DW
18 Table 5: Approximations to, and bounds for, ψ(5,t) - Pareto claims δ t =1 t=5 t=10 t=20 t=40 LB UB AVG SIM Lu Ul A DW LB UB AVG SIM Lu Ul A DW LB UB AVG SIM Lu Ul A DW LB UB AVG SIM Lu Ul A DW Table 6: Approximations to, and bounds for, ψ(10,t) - Pareto claims δ t =1 t=5 t=10 t=20 t=40 Lu Ul A DW Lu Ul A DW Lu Ul A DW Lu Ul A DW
19 Table 7 shows approximations to ψ(u, t) for large values of t considering the individual claim amounts exponentially distributed with mean 1. The values in column (t = ) are exact and provided by Segerdahl s formula. Table 7: Approximations to ψ(u, t) for large values of t - exponential claims u δ t =50 t= 100 t = 200 t = 300 t = 400 t = 500 t = Albrecher et al (2001) investigated when it is suitable to represent the probability of survival, 1 ψ δ (u, t), as gamma series. In particular, they derived exact analytical solutions for exponentially distributed claim sizes with mean 1/θ and for integer values of λ/δ. From their paper, second example in Section 3, we find that, for λ =2δ, with ψ(u, t) =1 b 0 (t) (1 e θu )b 1 (t) (1 e θu θue θu )b 2 (t) (9) b 0 (t) = θ2 c 2 D+e R1t (δ 2 θc + D(δθc + δ 2 ) δ 3 )+e R2t ( δ 2 θc + δ 3 + D(δθc + δ 2 )) D(θ 2 c 2 +2δθc +2δ 2 ) b 1 (t)= δ(2dθc e R1t ( 4δ 2 δθc + Dθc) e R2t (4δ 2 + δθc + Dθc)) D(θ 2 c 2 +2δθc +2δ 2 ) b 2 (t)= δ2 (e R1t (2θc +3δ+D) 2D+e R2t ( 2θc 3δ + D)) D(θ 2 c 2 +2δθc +2δ 2 ) where R 1 = 1 2 (2θc +3δ D), R 1 = 1 2 (2θc +3δ+D) and D = δ(4θc + δ). In Table 8 and Figures 3 to 5 we compare the values from our algorithms with the exact ones provided by (9). In this case we considered δ =0.5 and θ = 1. We can observe that the approximations are close to the exact values of the ruin probability. 19
20 Table 8: Approximations to, and bounds for, ψ(u, t), with δ =0.5 - exponential claims u t =1 t=5 t=10 t=20 Lu Ul A Exact Lu Ul A Exact Lu Ul A Exact L-u A Exact U-l Figure 3: Approximations to, and bounds for, ψ(0,t) with δ =0.5 - exponential(1) claims 20
21 L-u A Exact U-l Figure 4: Approximations to, and bounds for, ψ(5,t) with δ =0.5 - exponential(1) claims 21
22 L-u A Exact U-l Figure 5: Approximations to, and bounds for, ψ(10,t) with δ =0.5 - exponential(1) claims 22
23 5 References Albrecher, H., Teugels, J. L., and Tichy, R. F. (2001). On a gamma series expansion for the time-dependent probability of collective ruin. Insurance: Mathematics and Economics, 29: Boogaert, P. and De Waegenaere, A. (1990). Macro-economic version of a classical formula in risk theory. Insurance: Mathematics and Economics, 9: Brekelmans, R. and De Waegenaere, A. (2001). Approximating the finitetime ruin probability under interest force. Insurance: Mathematics and Economics, 29: Cardoso, R. M. R. and Egídio dos Reis, A. D. (2002). Recursive calculation of time to ruin distributions. Insurance: Mathematics and Economics, 30: De Vylder, F. E. (1996). Advanced risk theory. Editions de l Université de Bruxelles. De Vylder, F. E. and Goovaerts, M. J. (1988). Recursive calculation of finite time survival ruin probabilities. Insurance: Mathematics and Economics, 7:1 8. Dickson, D. C. M. and Gray, J. R. (1984). Approximations to ruin probability in the presence of an upper absorbing barrier. Scandinavian Actuarial Journal, Dickson, D. C. M. and Waters, H. R. (1991). Recursive calculation of finite time survival probabilities. Astin Bulletin, 21: Dickson, D. C. M. and Waters, H. R. (1999). Ruin probabilities with compounding assets. Insurance: Mathematics and Economics, 25: Karlin, S. and Taylor, H. M. (1975). A first course in stochastic processes. Academic Press, San Diego, 2nd edition. Panjer, H. H. (1981). Recursive calculation of a family of compound distributions. Astin Bulletin, 12:
24 Ross, S. M. (1996). Stochatic processes. Wiley, New York, 2nd edition. Sundt, B. and Teugels, J. L. (1995). Ruin estimates under interest force. Insurance: Mathematics and Economics, 16:
A Study of the Impact of a Bonus-Malus System in Finite and Continuous Time Ruin Probabilities in Motor Insurance
A Study of the Impact of a Bonus-Malus System in Finite and Continuous Time Ruin Probabilities in Motor Insurance L.B. Afonso, R.M.R. Cardoso, A.D. Egídio dos Reis, G.R Guerreiro This work was partially
More informationOn the probability of reaching a barrier in an Erlang(2) risk process
Statistics & Operations Research Transactions SORT 29 (2) July-December 25, 235-248 ISSN: 1696-2281 www.idescat.net/sort Statistics & Operations Research c Institut d Estadística de Transactions Catalunya
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 informationA polynomial expansion to approximate ruin probabilities
A polynomial expansion to approximate ruin probabilities P.O. Goffard 1 X. Guerrault 2 S. Loisel 3 D. Pommerêt 4 1 Axa France - Institut de mathématiques de Luminy Université de Aix-Marseille 2 Axa France
More informationThe finite-time Gerber-Shiu penalty function for two classes of risk processes
The finite-time Gerber-Shiu penalty function for two classes of risk processes July 10, 2014 49 th Actuarial Research Conference University of California, Santa Barbara, July 13 July 16, 2014 The finite
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 informationReinsurance and ruin problem: asymptotics in the case of heavy-tailed claims
Reinsurance and ruin problem: asymptotics in the case of heavy-tailed claims Serguei Foss Heriot-Watt University, Edinburgh Karlovasi, Samos 3 June, 2010 (Joint work with Tomasz Rolski and Stan Zachary
More informationarxiv: v1 [math.pr] 19 Aug 2017
Parisian ruin for the dual risk process in discrete-time Zbigniew Palmowski a,, Lewis Ramsden b, and Apostolos D. Papaioannou b, arxiv:1708.06785v1 [math.pr] 19 Aug 2017 a Department of Applied Mathematics
More informationNonlife Actuarial Models. Chapter 5 Ruin Theory
Nonlife Actuarial Models Chapter 5 Ruin Theory Learning Objectives 1. Surplus function, premium rate and loss process 2. Probability of ultimate ruin 3. Probability of ruin before a finite time 4. Adjustment
More informationSome Approximations on the Probability of Ruin and the Inverse Ruin Function
MATIMYÁS MATEMATIKA Journal of the Mathematical Society of the Philippines ISSN 115-6926 Vol. 38 Nos. 1-2 (215) pp. 43-5 Some Approximations on the Probability of Ruin and the Inverse Ruin Function Lu
More informationRuin problems for a discrete time risk model with non-homogeneous conditions. 1 A non-homogeneous discrete time risk model
Ruin problems for a discrete time risk model with non-homogeneous conditions ANNA CASTAÑER a, M. MERCÈ CLARAMUNT a, MAUDE GATHY b, CLAUDE LEFÈVRE b, 1 and MAITE MÁRMOL a a Universitat de Barcelona, Departament
More informationOn Finite-Time Ruin Probabilities in a Risk Model Under Quota Share Reinsurance
Applied Mathematical Sciences, Vol. 11, 217, no. 53, 269-2629 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/ams.217.7824 On Finite-Time Ruin Probabilities in a Risk Model Under Quota Share Reinsurance
More informationRuin probabilities of the Parisian type for small claims
Ruin probabilities of the Parisian type for small claims Angelos Dassios, Shanle Wu October 6, 28 Abstract In this paper, we extend the concept of ruin in risk theory to the Parisian type of ruin. For
More informationEVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION
An. Şt. Univ. Ovidius Constanţa Vol. 92), 2001, 181 192 EVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION Raluca Vernic Dedicated to Professor Mirela Ştefănescu on the occasion of her
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 informationDiscounted probabilities and ruin theory in the compound binomial model
Insurance: Mathematics and Economics 26 (2000) 239 250 Discounted probabilities and ruin theory in the compound binomial model Shixue Cheng a, Hans U. Gerber b,, Elias S.W. Shiu c,1 a School of Information,
More informationON BOUNDS FOR THE DIFFERENCE BETWEEN THE STOP-LOSS TRANSFORMS OF TWO COMPOUND DISTRIBUTIONS
ON BOUNDS FOR THE DIFFERENCE BETWEEN THE STOP-LOSS TRANSFORMS OF TWO COMPOUND DISTRIBUTIONS BJORN SUNDT I AND JAN DHAENE 2 ABSTRACT In the present note we deduce a class of bounds for the difference between
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 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 informationOn Finite-Time Ruin Probabilities for Classical Risk Models
On Finite-Time Ruin Probabilities for Classical Risk Models Claude Lefèvre, Stéphane Loisel To cite this version: Claude Lefèvre, Stéphane Loisel. On Finite-Time Ruin Probabilities for Classical Risk Models.
More informationOptimal stopping of a risk process when claims are covered immediately
Optimal stopping of a risk process when claims are covered immediately Bogdan Muciek Krzysztof Szajowski Abstract The optimal stopping problem for the risk process with interests rates and when claims
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 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 informationType II Bivariate Generalized Power Series Poisson Distribution and its Applications in Risk Analysis
Type II Bivariate Generalized Power Series Poisson Distribution and its Applications in Ris Analysis KK Jose a, and Shalitha Jacob a,b a Department of Statistics, St Thomas College, Pala, Arunapuram, Kerala-686574,
More informationUpper Bounds for the Ruin Probability in a Risk Model With Interest WhosePremiumsDependonthe Backward Recurrence Time Process
Λ4flΛ4» ν ff ff χ Vol.4, No.4 211 8fl ADVANCES IN MATHEMATICS Aug., 211 Upper Bounds for the Ruin Probability in a Risk Model With Interest WhosePremiumsDependonthe Backward Recurrence Time Process HE
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 informationModelling the risk process
Modelling the risk process Krzysztof Burnecki Hugo Steinhaus Center Wroc law University of Technology www.im.pwr.wroc.pl/ hugo Modelling the risk process 1 Risk process If (Ω, F, P) is a probability space
More informationAnalysis of the ruin probability using Laplace transforms and Karamata Tauberian theorems
Analysis of the ruin probability using Laplace transforms and Karamata Tauberian theorems Corina Constantinescu and Enrique Thomann Abstract The classical result of Cramer-Lundberg states that if the rate
More informationTwo binomial methods for evaluating the aggregate claims distribution in De Pril s individual risk model
Two binomial methods for evaluating the aggregate claims distribution in De Pril s individual risk model Sundt, Bjørn Vital Forsikring ASA P.O.Box 250, N-1326 Lysaker, Norway Phone: (+47) 67 83 44 71 Fax:
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 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 informationDEP ARTEMENT TOEGEP ASTE ECONOMISCHE WETENSCHAPPEN
DEP ARTEMENT TOEGEP ASTE ECONOMISCHE WETENSCHAPPEN ONDERZOEKSRAPPORT NR 9513 Some Moment Relations for the Hipp Approximation by JanDHAENE Bj0mSUNDT Nelson DE PRIL Katholieke Universiteit Leuven Naamsestraat
More informationThreshold dividend strategies for a Markov-additive risk model
European Actuarial Journal manuscript No. will be inserted by the editor Threshold dividend strategies for a Markov-additive risk model Lothar Breuer Received: date / Accepted: date Abstract We consider
More information(implicitly assuming time-homogeneity from here on)
Continuous-Time Markov Chains Models Tuesday, November 15, 2011 2:02 PM The fundamental object describing the dynamics of a CTMC (continuous-time Markov chain) is the probability transition (matrix) function:
More informationSimulation methods in ruin models with non-linear dividend barriers
Simulation methods in ruin models with non-linear dividend barriers Hansjörg Albrecher, Reinhold Kainhofer, Robert F. Tichy Department of Mathematics, Graz University of Technology, Steyrergasse 3, A-8
More informationTRANSACTIONS OF SOCIETY OF ACTUARIES 1992 VOL. 44 A PRACTICAL ALGORITHM FOR APPROXIMATING THE PROBABILITY OF RUIN COLIN M.
TRANSACTIONS OF SOCIETY OF ACTUARIES 1992 VOL. 44 A PRACTICAL ALGORITHM FOR APPROXIMATING THE PROBABILITY OF RUIN COLIN M. RAMSAY ABSTRACT For a compound Poisson process, the sequence of record jumps generated
More informationMeasuring the effects of reinsurance by the adjustment coefficient in the Sparre Anderson model
Measuring the effects of reinsurance by the adjustment coefficient in the Sparre Anderson model Centeno, Maria de Lourdes CEMAPRE, ISEG, Technical University of Lisbon and Centre for Actuarial Studies,
More informationMonotonicity and Aging Properties of Random Sums
Monotonicity and Aging Properties of Random Sums Jun Cai and Gordon E. Willmot Department of Statistics and Actuarial Science University of Waterloo Waterloo, Ontario Canada N2L 3G1 E-mail: jcai@uwaterloo.ca,
More informationRuin Theory. A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science at George Mason University
Ruin Theory A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science at George Mason University by Ashley Fehr Bachelor of Science West Virginia University, Spring
More informationA Dynamic Contagion Process with Applications to Finance & Insurance
A Dynamic Contagion Process with Applications to Finance & Insurance Angelos Dassios Department of Statistics London School of Economics Angelos Dassios, Hongbiao Zhao (LSE) A Dynamic Contagion Process
More informationAsymptotics of random sums of heavy-tailed negatively dependent random variables with applications
Asymptotics of random sums of heavy-tailed negatively dependent random variables with applications Remigijus Leipus (with Yang Yang, Yuebao Wang, Jonas Šiaulys) CIRM, Luminy, April 26-30, 2010 1. Preliminaries
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 informationStochastic Areas and Applications in Risk Theory
Stochastic Areas and Applications in Risk Theory July 16th, 214 Zhenyu Cui Department of Mathematics Brooklyn College, City University of New York Zhenyu Cui 49th Actuarial Research Conference 1 Outline
More informationSMSTC (2007/08) Probability.
SMSTC (27/8) Probability www.smstc.ac.uk Contents 12 Markov chains in continuous time 12 1 12.1 Markov property and the Kolmogorov equations.................... 12 2 12.1.1 Finite state space.................................
More informationA conceptual interpretation of the renewal theorem with applications
Risk, Reliability and Societal Safety Aven & Vinnem (eds) 2007 Taylor & Francis Group, London, ISBN 978-0-415-44786-7 A conceptual interpretation of the renewal theorem with applications J.A.M. van der
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 Distribution of the Minimum of Independent Phase Type Random Variables
The Distribution of the Minimum of Independent Phase Type Random Variables Darcy Bermingham, supervised by Yoni Nazarathy The University of Queensland School of Mathematics and Physics December 215 In
More informationPractical approaches to the estimation of the ruin probability in a risk model with additional funds
Modern Stochastics: Theory and Applications (204) 67 80 DOI: 05559/5-VMSTA8 Practical approaches to the estimation of the ruin probability in a risk model with additional funds Yuliya Mishura a Olena Ragulina
More information1 Gambler s Ruin Problem
1 Gambler s Ruin Problem Consider a gambler who starts with an initial fortune of $1 and then on each successive gamble either wins $1 or loses $1 independent of the past with probabilities p and q = 1
More informationA Ruin Model with Compound Poisson Income and Dependence Between Claim Sizes and Claim Intervals
Acta Mathematicae Applicatae Sinica, English Series Vol. 3, No. 2 (25) 445 452 DOI:.7/s255-5-478- http://www.applmath.com.cn & www.springerlink.com Acta Mathema cae Applicatae Sinica, English Series The
More informationThe Compound Poisson Risk Model with a Threshold Dividend Strategy
The Compound Poisson Risk Model with a Threshold Dividend Strategy X. Sheldon Lin Department of Statistics University of Toronto Toronto, ON, M5S 3G3 Canada tel.: (416) 946-5969 fax: (416) 978-5133 e-mail:
More informationStochastic processes. MAS275 Probability Modelling. Introduction and Markov chains. Continuous time. Markov property
Chapter 1: and Markov chains Stochastic processes We study stochastic processes, which are families of random variables describing the evolution of a quantity with time. In some situations, we can treat
More informationHANDBOOK OF APPLICABLE MATHEMATICS
HANDBOOK OF APPLICABLE MATHEMATICS Chief Editor: Walter Ledermann Volume II: Probability Emlyn Lloyd University oflancaster A Wiley-Interscience Publication JOHN WILEY & SONS Chichester - New York - Brisbane
More informationConditional independence of blocked ordered data
Conditional independence of blocked ordered data G. Iliopoulos 1 and N. Balakrishnan 2 Abstract In this paper, we prove that blocks of ordered data formed by some conditioning events are mutually independent.
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 3. Calculaus in Deterministic and Stochastic Environments Steve Yang Stevens Institute of Technology 01/31/2012 Outline 1 Modeling Random Behavior
More informationHSC 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 informationTHREE METHODS TO CALCULATE THE PROBABILITY OF RUIN. University of Lausanne, Switzerland
THREE METHODS TO CALCULATE THE PROBABILITY OF RUIN BY FRANCOIS DUFRESNE AND HANS U. GERBER University of Lausanne, Switzerland ABSTRACT The first method, essentially due to GOOVAERTS and DE VYLDER, uses
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 informationIEOR 3106: Introduction to Operations Research: Stochastic Models. Professor Whitt. SOLUTIONS to Homework Assignment 2
IEOR 316: Introduction to Operations Research: Stochastic Models Professor Whitt SOLUTIONS to Homework Assignment 2 More Probability Review: In the Ross textbook, Introduction to Probability Models, read
More informationRuin, Operational Risk and How Fast Stochastic Processes Mix
Ruin, Operational Risk and How Fast Stochastic Processes Mix Paul Embrechts ETH Zürich Based on joint work with: - Roger Kaufmann (ETH Zürich) - Gennady Samorodnitsky (Cornell University) Basel Committee
More informationTHE FAST FOURIER TRANSFORM ALGORITHM IN RUIN THEORY FOR THE CLASSICAL RISK MODEL
y y THE FST FORIER TRNSFORM LGORITHM IN RIN THEORY FOR THE CLSSICL RISK MODEL Susan M. itts niversity of Cambridge bstract We focus on numerical evaluation of some quantities of interest in ruin theory,
More informationPricing and Risk Analysis of a Long-Term Care Insurance Contract in a non-markov Multi-State Model
Pricing and Risk Analysis of a Long-Term Care Insurance Contract in a non-markov Multi-State Model Quentin Guibert Univ Lyon, Université Claude Bernard Lyon 1, ISFA, Laboratoire SAF EA2429, F-69366, Lyon,
More informationValue Iteration and Action ɛ-approximation of Optimal Policies in Discounted Markov Decision Processes
Value Iteration and Action ɛ-approximation of Optimal Policies in Discounted Markov Decision Processes RAÚL MONTES-DE-OCA Departamento de Matemáticas Universidad Autónoma Metropolitana-Iztapalapa San Rafael
More informationAsymptotics of the Finite-time Ruin Probability for the Sparre Andersen Risk. Model Perturbed by an Inflated Stationary Chi-process
Asymptotics of the Finite-time Ruin Probability for the Sparre Andersen Risk Model Perturbed by an Inflated Stationary Chi-process Enkelejd Hashorva and Lanpeng Ji Abstract: In this paper we consider the
More informationA Comparison: Some Approximations for the Aggregate Claims Distribution
A Comparison: ome Approximations for the Aggregate Claims Distribution K Ranee Thiagarajah ABTRACT everal approximations for the distribution of aggregate claims have been proposed in the literature. In
More informationComonotonicity and Maximal Stop-Loss Premiums
Comonotonicity and Maximal Stop-Loss Premiums Jan Dhaene Shaun Wang Virginia Young Marc J. Goovaerts November 8, 1999 Abstract In this paper, we investigate the relationship between comonotonicity and
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 information4 Branching Processes
4 Branching Processes Organise by generations: Discrete time. If P(no offspring) 0 there is a probability that the process will die out. Let X = number of offspring of an individual p(x) = P(X = x) = offspring
More informationJackknife Euclidean Likelihood-Based Inference for Spearman s Rho
Jackknife Euclidean Likelihood-Based Inference for Spearman s Rho M. de Carvalho and F. J. Marques Abstract We discuss jackknife Euclidean likelihood-based inference methods, with a special focus on the
More informationTopics in Probability Theory and Stochastic Processes Steven R. Dunbar. Waiting Time to Absorption
Steven R. Dunbar Department of Mathematics 203 Avery Hall University of Nebraska-Lincoln Lincoln, NE 6888-030 http://www.math.unl.edu Voice: 402-472-373 Fax: 402-472-8466 Topics in Probability Theory and
More informationIndividual loss reserving with the Multivariate Skew Normal framework
May 22 2013, K. Antonio, KU Leuven and University of Amsterdam 1 / 31 Individual loss reserving with the Multivariate Skew Normal framework Mathieu Pigeon, Katrien Antonio, Michel Denuit ASTIN Colloquium
More informationSurvival probabilities in bivariate risk models, with application to reinsurance
Survival probabilities in bivariate risk models, with application to reinsurance A. Castañer a, M.M. Claramunt a,, C. Lefèvre b a Dept. Matemàtica Econòmica, Financera i Actuarial, Universitat de Barcelona,
More information1 Gambler s Ruin Problem
Coyright c 2017 by Karl Sigman 1 Gambler s Ruin Problem Let N 2 be an integer and let 1 i N 1. Consider a gambler who starts with an initial fortune of $i and then on each successive gamble either wins
More informationDistribution of the first ladder height of a stationary risk process perturbed by α-stable Lévy motion
Insurance: Mathematics and Economics 28 (21) 13 2 Distribution of the first ladder height of a stationary risk process perturbed by α-stable Lévy motion Hanspeter Schmidli Laboratory of Actuarial Mathematics,
More informationApproximating diffusions by piecewise constant parameters
Approximating diffusions by piecewise constant parameters Lothar Breuer Institute of Mathematics Statistics, University of Kent, Canterbury CT2 7NF, UK Abstract We approximate the resolvent of a one-dimensional
More informationTail negative dependence and its applications for aggregate loss modeling
Tail negative dependence and its applications for aggregate loss modeling Lei Hua Division of Statistics Oct 20, 2014, ISU L. Hua (NIU) 1/35 1 Motivation 2 Tail order Elliptical copula Extreme value copula
More informationProblem set 2 The central limit theorem.
Problem set 2 The central limit theorem. Math 22a September 6, 204 Due Sept. 23 The purpose of this problem set is to walk through the proof of the central limit theorem of probability theory. Roughly
More information0.1 Naive formulation of PageRank
PageRank is a ranking system designed to find the best pages on the web. A webpage is considered good if it is endorsed (i.e. linked to) by other good webpages. The more webpages link to it, and the more
More informationMIXED POISSON DISTRIBUTIONS ASSOCIATED WITH HAZARD FUNCTIONS OF EXPONENTIAL MIXTURES
MIXED POISSON DISTRIBUTIONS ASSOCIATED WITH HAZARD FUNCTIONS OF EXPONENTIAL MIXTURES Moses W. Wakoli and Joseph A. M. Ottieno 2 Abstra The hazard funion of an exponential mixture charaerizes an infinitely
More informationDiscrete time Markov chains. Discrete Time Markov Chains, Limiting. Limiting Distribution and Classification. Regular Transition Probability Matrices
Discrete time Markov chains Discrete Time Markov Chains, Limiting Distribution and Classification DTU Informatics 02407 Stochastic Processes 3, September 9 207 Today: Discrete time Markov chains - invariant
More informationESTIMATING THE PROBABILITY OF RUIN FOR VARIABLE PREMIUMS BY SIMULATION. University of Lausanne, Switzerland
ESTIMATING THE PROBABILITY OF RUIN FOR VARIABLE PREMIUMS BY SIMULATION BY FREDERIC MICHAUD University of Lausanne, Switzerland ABSTRACT There is a duality between the surplus process of classical risk
More informationDividend problems in the dual risk model
Dividend problems in the dual risk model Lourdes B. Afonso, Rui M.R. Cardoso & Alfredo D. gídio dos Reis CMA and FCT, New University of Lisbon & CMAPR and ISG, Technical University of Lisbon, Portugal
More informationOn a discrete time risk model with delayed claims and a constant dividend barrier
On a discrete time risk model with delayed claims and a constant dividend barrier Xueyuan Wu, Shuanming Li Centre for Actuarial Studies, Department of Economics The University of Melbourne, Parkville,
More informationModelling Operational Risk Using Bayesian Inference
Pavel V. Shevchenko Modelling Operational Risk Using Bayesian Inference 4y Springer 1 Operational Risk and Basel II 1 1.1 Introduction to Operational Risk 1 1.2 Defining Operational Risk 4 1.3 Basel II
More informationLecture 12. F o s, (1.1) F t := s>t
Lecture 12 1 Brownian motion: the Markov property Let C := C(0, ), R) be the space of continuous functions mapping from 0, ) to R, in which a Brownian motion (B t ) t 0 almost surely takes its value. Let
More informationStatistics 3657 : Moment Generating Functions
Statistics 3657 : Moment Generating Functions A useful tool for studying sums of independent random variables is generating functions. course we consider moment generating functions. In this Definition
More information6.1 Moment Generating and Characteristic Functions
Chapter 6 Limit Theorems The power statistics can mostly be seen when there is a large collection of data points and we are interested in understanding the macro state of the system, e.g., the average,
More informationClaims Reserving under Solvency II
Claims Reserving under Solvency II Mario V. Wüthrich RiskLab, ETH Zurich Swiss Finance Institute Professor joint work with Michael Merz (University of Hamburg) April 21, 2016 European Congress of Actuaries,
More informationGENERALIZED ANNUITIES AND ASSURANCES, AND INTER-RELATIONSHIPS. BY LEIGH ROBERTS, M.Sc., ABSTRACT
GENERALIZED ANNUITIES AND ASSURANCES, AND THEIR INTER-RELATIONSHIPS BY LEIGH ROBERTS, M.Sc., A.I.A ABSTRACT By the definition of generalized assurances and annuities, the relation is shown to be the simplest
More informationstochnotes Page 1
stochnotes110308 Page 1 Kolmogorov forward and backward equations and Poisson process Monday, November 03, 2008 11:58 AM How can we apply the Kolmogorov equations to calculate various statistics of interest?
More informationASTIN Colloquium 1-4 October 2012, Mexico City
ASTIN Colloquium 1-4 October 2012, Mexico City Modelling and calibration for non-life underwriting risk: from empirical data to risk capital evaluation Gian Paolo Clemente Dept. Mathematics, Finance Mathematics
More informationLECTURE 10: REVIEW OF POWER SERIES. 1. Motivation
LECTURE 10: REVIEW OF POWER SERIES By definition, a power series centered at x 0 is a series of the form where a 0, a 1,... and x 0 are constants. For convenience, we shall mostly be concerned with the
More informationAnalysis of a Bivariate Risk Model
Jingyan Chen 1 Jiandong Ren 2 July 23, 2012 1 MSc candidate, Department of Statistical and Actuarial Sciences, The University of Western Ontario, London, Ontario, Canada. 2 Associate Professor, Department
More informationStandard Error of Technical Cost Incorporating Parameter Uncertainty
Standard Error of Technical Cost Incorporating Parameter Uncertainty Christopher Morton Insurance Australia Group This presentation has been prepared for the Actuaries Institute 2012 General Insurance
More informationHeriot-Watt University. M.Sc. in Actuarial Science. Life Insurance Mathematics I. Tutorial 2 Solutions
Heriot-Watt University M.Sc. in Actuarial Science Life Insurance Mathematics I Tutorial 2 Solutions. (a) The recursive relationship between pure enowment policy values is: To prove this, write: (V (t)
More informationLecture 20 : Markov Chains
CSCI 3560 Probability and Computing Instructor: Bogdan Chlebus Lecture 0 : Markov Chains We consider stochastic processes. A process represents a system that evolves through incremental changes called
More informationSTOCHASTIC PROCESSES Basic notions
J. Virtamo 38.3143 Queueing Theory / Stochastic processes 1 STOCHASTIC PROCESSES Basic notions Often the systems we consider evolve in time and we are interested in their dynamic behaviour, usually involving
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationChain ladder with random effects
Greg Fry Consulting Actuaries Sydney Australia Astin Colloquium, The Hague 21-24 May 2013 Overview Fixed effects models Families of chain ladder models Maximum likelihood estimators Random effects models
More informationMath Stochastic Processes & Simulation. Davar Khoshnevisan University of Utah
Math 5040 1 Stochastic Processes & Simulation Davar Khoshnevisan University of Utah Module 1 Generation of Discrete Random Variables Just about every programming language and environment has a randomnumber
More information