Moment-Based Approximation with Mixed Erlang Distributions
|
|
- Mark Day
- 6 years ago
- Views:
Transcription
1 Moment-Based Approximation with Mixed Erlang Distributions by Hélène Cossette, David Landriault, Etienne Marceau, Khouzeima Moutanabbir Abstract Moment-based approximations have been extensively analyzed over the years (see, e.g., Osogami and Harchol-Balter (2006) and references therein). A number of speci c phase-type (and non phasetype) distributions have been considered to tackle the moment-matching problem (see, for instance, Johnson and Taa e (989)). Motivated by the development of more exible moment-based approximation methods, we develop and examine the use of nite mixture of Erlangs with a common rate parameter for the moment-matching problem. This is primarily motivated by Tijms (994) who shows that this class of distributions can approximate any continuous positive distribution to an arbitrary level of accuracy, as well as the tractability of this class of distributions for various problems of interest in quantitative risk management. We consider separately situations where the rate parameter is either known or unknown. For the former case, a direct connection with a discrete moment-matching problem is established. A parallel to the s-convex stochastic order (e.g., Denuit et al. (998)) is also drawn. Numerical examples are considered throughout. Keywords : Risk Theory; Mixed Erlang distributions; Moment-matching; Distribution Fitting; Phase-type approximation. Introduction Mixed Erlang distributions are known to yield analytic solutions to many risk management problems of interest. This is primarily due to the tractable features of this distributional class. Among others, the class of mixed Erlang distributions is closed under various operations such as convolutions and Esscher transformations (e.g., Willmot and Woo (2007) and Willmot and Lin (20)). As such, risk aggregation and ruin problems can more easily be tackled under mixed Erlang assumptions (e.g., Cheung and Woo (204), Cossette et al. (202), and Landriault and Willmot (2009)). Also, Tijms (994) showed that the class of mixed Erlang distributions is dense in the set of all continuous and positive distributions. Therefore, we consider a moment-based approximation method which capitalizes on the aforementioned properties of the mixed Erlang distribution. More precisely, we propose to approximate a distribution with known moments by a moment-matching mixed Erlang distribution. Moment-based approximations have been extensively developed in various research areas including performance evaluation, queueing theory, and risk theory, to name a few. Osogami and Harchol-Balter (2006) identi es the following four criteria to evaluate moment-matching algorithms: () the number of moments matched; (2) the computational e ciency of the algorithm; (3)
2 the generality of the solution; (4) the minimality of the number of parameters (phases). It also seems desirable for the approximation to be in itself a distribution. This is not mentioned in Osogami and Harchol-Balter (2006) for the obvious reason that they consider phase-type distributions as their momentbased approximation class. There exists an extensive literature on the approximation of distributions by a speci c subset of phase-type distributions using moment-based techniques. For instance, Whitt (982) proposed a mixture of two exponential distributions or a generalized Erlang distribution as a moment-based approximation when either the coe cient of variation (CV) is greater than or less than, respectively. Also, both Altiok (985) and Vanden Bosch et al. (2000) proposed an alternative to the moment-based approximation of Whitt (982) when CV> using a Coxian distribution. Alternatively, Johnson and Taa e (989) considered a mixture of Erlangs with a common shape (order) parameter as their moment-based approximation. Most predominantly, there exists a substantial body of literature on the three-moment approximation within the phase-type class of distributions (e.g., Telek and Heindl (2002), Bobbio et al. (2005), and references therein). Matching the rst three moments is often viewed as e ective to provide a reasonable approximation to the underlying system (e.g., Osogami and Harchol-Balter (2006) and references therein). However, as illustrated in this paper and many others, three moments does not always su ce, triggering the development of more exible moment-based approximations. Among others, we mention the work of Johnson and Taa e (989) on mixed Erlang distributions of common order. Also, Dufresne (2007) proposes two approximation techniques based on Jacobi polynomial expansions and the logbeta distribution to t combinations of exponential distributions. This paper is obviously complementary to the aforementioned ones by considering the family of nite mixture of Erlangs with common rate parameter to approximate a distribution on R +, as theoretically justi ed in the continuous case by Tijms (994, Theorem 3.9.). The reader is also referred to Lee and Lin (200) where tting of the same class of distributions is considered using the EM algorithm (which relies on the knowledge of the approximated distribution rather than only its moments). It is worth pointing out that other non-phase type approximation methods have been widely used in actuarial science. A good survey paper on this topic is Chaubey et al. (998). One of these approximation classes are re nements to the Normal approximation such as the Normal Power and the Cornish Fisher approximations (e.g., Ramsay (99), Daykin et al. (994), and Lee and Lin (992)). These approximations are based on the rst few moments. However, the resulting approximation is often not a proper distribution. Other moment-based distributional approximations are the translated gamma distribution (e.g., Seal (977)), translated inverse gaussian distribution (e.g., Chaubey et al. (998)) and the generalized Pareto distribution (e.g., Venter (983)). It should be noted that all these approximation methods are designed to t a speci c number of moments, and thus lack the exibility to match an arbitrary number of moments. The rest of the paper is constructed as follows. In Section 2, a brief review on admissible moments, mixed Erlang distributions and the approximation method of Johnson and Taa e (989) is provided. Section 3 is devoted to our class of nite mixture of Erlangs with common rate parameter. Theoretical and practical considerations related to the approximation method are drawn. Various examples are considered to examine the quality of the resulting approximation. In Section 4, we consider applications of our moment-based approximations of Section 3 when the underlying distribution is of mixed Erlang form with known rate parameter. A parallel is drawn with a discrete moment-matching problem and certain stochastic orderings, notably the s-convex stochastic order (e.g. Denuit et al. (998)). An application of Cossette et al. (2002) will later be examined in more details. 2
3 2 Background 2. Admissible moments Karlin and Studden (966) provide the necessary and su cient conditions for a set of (raw) moments m = ( ; :::; m ) to be from a probability distribution de ned on R +. To state this result, de ne the matrices P k and Q k (k ) as 0 0 k 2 k+ 2 k+ P k = C. A ; Q 2 3 k+2 k = C. A. k k+ 2k k+ k+2 2k+ As stated in Courtois and Denuit (2007), there exists a non-negative random variable (rv) with distribution function (df) F and rst m moments m if and only if the following two conditions are satis ed: det P k > 0, for k = ; :::; b(m )=2c ; det Q k > 0, for k = ; :::; bm=2c ; where bxc holds for the integer part of x. In what follows, we silently assume the moment set m is from a probability distribution on R Mixed Erlang Distribution We now review some known properties of mixed Erlang distributions with common rate parameter. A more elaborate review of this class of distributions can be found in Willmot and Woo (2007), Lee and Lin (200), and Willmot and Lin (20). Let W be a mixed Erlang rv with common rate parameter > 0 and df F W (x) = X k2a l k H(x; k; ); () where A l = f; 2; :::; lg, and f k g l k= is the probability mass function (pmf) of a discrete rv K with support A l for a given l 2 f; 2; :::g [ fg. The Erlang df H is de ned as Xk H(x; k; ) H(x; k; ) = e x i=0 (x) i, x 0. (2) i! where the parameters k and of the Erlang df are known as the shape and rate parameters, respectively. An alternative and useful representation of the mixed Erlang rv W is as a compound rv, i.e. W = P K k= C k where fc k g k are iid exponential rv s with mean =, independent of K. Remark As in e.g., Willmot and Woo (2007), we consider the class of mixed Erlang dfs () rather than the more general class of combinations of Erlangs where some k s are possibly negative. For the latter 3
4 class, additional constraints on f k g l k= exist to ensure that the right-hand side of () is a non-decreasing function in x. This presents additional challenges in the subsequent moment-matching application, challenges which do not arise in the mixed Erlang case. It is well known that the j-th moment of W is given by E W j = j P n k= Qj k i=0 o. (k + i) Of particular importance in actuarial science and quantitative risk management (see, e.g., McNeil et al. (2005) and references therein) are the VaR and TVaR risk measures. For the mixed Erlang rv W, there is in general no closed form expression for VaR (W ) = inf fx 2 R : F W (x) g where 0 <, but its value can be obtained using a routine numerical procedure. As for its TVaR, Lee and Lin (200) showed that TVaR (W ) Z VaR u (W ) du = X k= k k H (VaR (W ) ; k + ; ) : (3) Another quantity of interest is the stop-loss premium de ned as W (b) = E (W b) + with (x)+ = max fx; 0g. For the mixed Erlang df (), we have W (b) = X k= k k H (b; k + ; ) bh (b; k; ) ; b 0. (see also Willmot and Woo (2007, Eq. 3.6) for the higher-order stop-loss moments). Tijms (994) showed that this class of distributions can approximate any continuous positive distribution with an arbitrary level of accuracy. For completeness, the theoretical foundation of this result is given next. Theorem 2 (Tijms (994, Theorem 3.9.)) Let F be the df of a positive rv. For any given h > 0, de ne F h (x) = X (F (kh) F ((k ) h)) H k= Then, lim h!0 F h (x) = F (x) for any continuity point x of F. x; k;. (4) h Note that F h in (4) is a mixed Erlang df of the form () with k = (F (kh) F ((k ) h)) (k = ; 2; :::) and rate parameter = =h. Several approximation methods motivated by Tijms theorem were proposed over the years (see Section for more details). In general, these moment-based approximations propose to work with a speci c subclass of all nite and in nite mixed Erlang distributions. Among them, we recall the method of Johnson and Taa e (989) who will be used later for comparative purposes. 2.3 Method of Johnson and Taa e (989) Johnson and Taa e (989) investigated the use of mixtures of Erlang distributions with common shape parameter for moment-matching purposes. More precisely, mixtures of n (or fewer) Erlangs with common shape parameters are used to match the rst (2n ) moments (whenever the set of moments is within the feasible set). For the three-moment matching problem, Johnson and Taa e (989) generalized the 4
5 approximation of Whitt (982) and Altiok (985) by enlarging the set of feasible moments 3 when CV >. Their method is also valid for some combinations of 3 when CV <. Their three-moment approximation is a mixture of two Erlangs with common shape parameter r (see Theorem 3 of Johnson and Taa e (989)), i.e. F (x) = ph (x; r; ) + ( p) H (x; r; 2 ), (5) where p = r 2 = 2, and f=i g 2 i= are the solutions of As2 + Bs + C = 0 with r + A = r (r + 2) 2 B = r + r 3 r r C = 3 r + r r (r + 2) + r + r r 2 + (n + 2) 2 r r The choice of the shape parameter r is discussed in Johnson and Taa e (989, Proposition 4).! 3 Moment-based approximation with mixed Erlang distribution In this section, we propose to use a di erent subclass of mixed Erlang distributions to examine momentbased approximation techniques. 3. Description of the approach For a given l 2 N +, let ME( m ; A l ) be the set of all nite mixture of Erlangs with df () and rst m moments m. From Section 2.2, this consists in the identi cation of all solutions to the following problem: lx Q j i=0 (k + i) k j = j ; j = ; :::; m, (6) k= under the constraints that > 0 and f k g l k= is a probability measure on A l. Remark 3 For a rv W with df () and rst m moments m, we indi erently write W 2 ME( m ; A l ) or F W 2 ME( m ; A l ). This will also apply to the other distributional classes. Also, let ME res ( m ; A l ) be the (restricted) subset of ME( m ; A l ) with at most m non-zero mixing probabilities f k g l k=. Given that MEres ( m ; A l ) has a nite number of solutions, we propose to use ME res ( m ; A l ) as our approximation class. It is clear that ME res ( m ; A l ) ME res ( m ; A l 0) for l l 0. Note that for a continuous positive distribution with moments m, we know from Theorem 2 that there exists a l large enough such that ME( m ; A l ) is not empty. Even though no formal conclusion can 5
6 be reached for the restricted class ME res ( m ; A l ), all our numerical studies have shown that this set has a large number of distributions (see, for instance, the examples of subsections 3.2. and 3.2.2) for a given m when l is chosen large enough. Distributions in the ME res ( m ; A l ) class are identi ed as follows: for a given set fi k g m k= A l with i < i 2 < ::: < i m l, (6) can be rewritten in matrix form as G m m = M ; where m = ( i ; i2 ; :::; im ) T, M = ( ; 2 2 ; :::; m m ) T, and 0 i i 2 i m i (i + ) i 2 (i 2 + ) i m (i m + ) G m =. B mq mq mq A (i + i) (i 2 + i) (i m + i) i=0 i=0 It follows that m = Gm M under the constraint that T me =, where e is a column vector of s. Note that T me = Gm T M e is a polynomial of degree (at most) m in. Thus, we only consider the real and positive solutions (in ) of Gm T M e =, and complete their mixed Erlang representation with the identi cation of the mixing weights m = Gm M. The procedure is systematically repeated for all l m possible sets of m distinct elements in Al. i=0 l Remark 4 Given that the above procedure is repeated m -times, the computational e ciency of the proposed methodology is mostly driven by this number, and hence the parameters m and l should be chosen accordingly. For a given number of moments m, we observe that: (a) larger values of l result in a more time-consuming numerical procedure; (b) however, l should be chosen large enough for the approximation class ME res ( m ; A l ) to have a reasonable number of members (to legitimally produce a "good" approximation). From our numerical studies, we observe that the selection of l (for a given m) can be problem-speci c, and thus this tradeo in the choice of l should be handled with care. However, as a rule of thumb when m is relatively small (i.e. m 6) which is traditionally in moment-matching exercises, a value of l between 50 and 00 leads to reasonable mixed Erlang approximations. We refer the reader to the numerical illustrations and subsequent remarks in Section 3.2 for a more detailed discussion on this topic. Among all ME res ( m ; A l ) distributions, we propose to use as a criterion of the quality of these approximations the Kolmogorov-Smirnov (KS) distance de ned for two rv s S and W (with respective dfs F S and F W ) as d KS (S; W ) = sup x0 jf S (x) F W (x)j. The KS distance is commonly used in the context of continuous distributions (e.g., Denuit et al. (2005)). Therefore, the chosen mixed Erlang approximation within ME res ( m ; A l ) is the one minimizing the KS distance with the true df F S. We denote by F Wm;l this approximation, i.e. d KS (S; W m;l ) = inf sup jf S (x) F W (x)j, F W 2ME res ( m ;A l ) where W m;l is a rv with df F Wm;l. This requires the calculation of the KS distance for each mixed Erlang distribution in ME res ( m ; A l ) to identify its minimizer W m;l. In general, an explicit expression for this x0 6
7 KS distance does not exist, and hence we propose to numerically nd this value by evaluating the distance between the two dfs over all multiples (up to a given high value) of a small discretization span. Note that other distances such as the stop-loss distance (e.g., Gerber (979)) could have been used to select our approximation distribution. Alternatively, one could have relied on another criterion to identify this approximation distribution (for instance, select the distribution in ME res ( m ; A l ) with the closest subsequent moment to the true distribution). 3.2 Numerical examples We consider a few simple examples to illustrate the quality of the approximation. For comparative purposes, other approximation methods will also be discussed. Some concluding remarks on the mixed Erlang approximation method are later made based on the numerical experiment conducted next Lognormal distribution: Dufresne (2007, Example 5.4) Let S = exp (Z) where Z is a normal rv with mean 0 and variance The rst 5 moments of S are: 5 = (:33; :6487; 3:0802; 7:389; 22:7599). We consider the class of mixed Erlang distributions ME res ( m ; A 70 ) (m = 3; 4; 5) which have a total of 398, and distributions, respectively. The resulting mixed Erlang approximations F Wm;70 (m = 3; 4; 5) are F W3;70 (x) = 0:8209H(x; 6; 6:329) + 0:727H(x; 2; 6:329) + 0:0064H(x; 26; 6:329); F W4;70 (x) = 0:6350H(x; 7; 8:3334)+0:2950H(x; 2; 8:3334)+0:0672H(x; 20; 8:3334)+0:0029H(x; 40; 8:3334); and F W5;70 (x) = 0:6273H(x; 7; 8:3608) + 0:3063H(x; 2; 8:3608) + 0:0609H(x; 20; 8:3608) + 0:0055H(x; 34; 8:3608) + 0:000H(x; 69; 8:3608); for x 0 with respective KS distances of d KS (S; W 3;70 ) = 0:0040, d KS (S; W 4;70 ) = 0:008 and d KS (S; W 5;70 ) = 0:00. Note that the quality of the mixed Erlang approximation (as measured by the KS distance) increases with the number of moments matched. For comparative purposes, the three-moment approximation (5) of Johnson and Taa e (989) is given by F WJT (x) = 0:0087H (x; 4; :2804) + 0:993H (x; 4; 3:5855), In Figure, we compare the density functions of W m;70 (m = 3; 4; 5) ; W JT and S. All three mixed Erlang approximations provide an overall good t to the exact distribution. To further examine the tail t, speci c values of VaR and TVaR for the exact and approximated distributions are provided in Tables and 2, respectively. 7
8 Figure : Density function: Lognormal vs Approximations VaR (W JT ) VaR (W 3;70 ) VaR (W 4;70 ) VaR (W 5;70 ) VaR (S) Table. Values of VaR for W JT, W m;70 (m = 3; 4; 5) and S TVaR (W JT ) TVaR (W 3;70 ) TVaR (W 4;70 ) TVaR (W 5;70 ) TVaR (S) Table 2. Values of TVaR for W JT, W m;70 (m = 3; 4; 5) and S We observe that the VaR and TVaR values of the mixed Erlang approximations compare very well to their lognormal counterparts, especially for the 5-moment approximation. This is particularly true given that the lognormal distribution is known to have a heavier tail than the mixed Erlang distribution. Note that the improvement is indeed not monotone with the number of moments matched, as increasing this number does not necessarily lead to a higher quality approximation in moment-matching techniques. 8
9 3.2.2 Mixture of Two Gamma Distributions: Lee and Lin (200, Section 5, Example ) Let S be a mixture of two Gamma distributions with density f S (s) = 0:2 (3:2s)2:6 e 3:2s s (2:6) + 0:8 (:2s)6:3 e :2s, s 0, s (6:3) and rst 6 moments 6 = (4:3623; 25:7308; 76:9624; 369:8272; 754:249; 0674:454). We consider the class of mixed Erlang distributions ME res ( m ; A 70 ) for m = 3; 4; 5; 6 which are composed of 6000, 83797, and distributions, respectively. The resulting mixed Erlang approximations F Wm;70 (m = 3; 4; 5; 6) are F W3;70 (x) = 0:240H(x; 2; 2:0835) + 0:525H(x; 9; 2:0835) + 0:2645H(x; 5; 2:0835); F W4;70 (x) = 0:2266H(x; 2; 2:0469)+0:4440H(x; 9; 2:0469)+0:2963H(x; 3; 2:0469)+0:0330H(x; 9; 2:0469); and F W5;70 (x) = 0:2023H(x; 3; 3:727) + 0:209H(x; 2; 3:727) + 0:3936H(x; 9; 3:727) + 0:805H(x; 28; 3:727) + 0:045H(x; 42; 3:727); F W6;70 (x) = 0:0768H(x; 2; 3:073) + 0:325H(x; 3; 3:073) + 0:2739H(x; ; 3:073) + 0:3928H(x; 7; 3:073) + 0:88H(x; 25; 3:073) + 0:0052H(x; 37; 3:073); with respective KS distances of d KS (S; W 3;70 ) = 0:048, d KS (S; W 4;70 ) = 0:0050, d KS (S; W 5;70 ) = 0:0035 and d KS (S; W 6;70 ) = 0:0024. Lee and Lin (200) used the EM algorithm to t a mixed Erlang distribution to the same distribution which resulted in the following model: F WEM (x) = 0:2282H(x; 2; :9603) + 0:5430H(x; 9; :9603) + 0:2288H(x; 4; :9603); with KS distance d KS (S; W EM ) = 0:0094. Note that the KS distance for the 3-moment approximation is greater than the KS distance obtained using the EM estimation. We recall that the EM algorithm nds maximum likelihood estimates using the approximated distribution as an input, while our method is based on only partial information on the approximated distribution (e.g., its rst m moments). We observe that the t improves and the KS distance decreases when more moments are included in the approximation. For illustrative purposes, we also provide in Tables 3 and 4 some values of VaR and TVaR for W m;70 (m = 3; 4; 5; 6), W EM and S. VaR (W EM ) VaR (W 3;70 ) VaR (W 4;70 ) VaR (W 5;70 ) VaR (W 6;70 ) VaR (S) Table 3. Values of VaR for W m;70 (m = 3; 4; 5; 6) ; W EM, and S 9
10 TVaR (W EM ) TVaR (W 3;70 ) TVaR (W 4;70 ) TVaR (W 5;70 ) TVaR (W 6;70 ) TVaR (S) Table 4. Values of TVaR for W m;70 (m = 3; 4; 5; 6) ; W EM and S Gompertz distribution The Gompertz distribution with df F S (x) = exp f B (c x ) = ln cg (x > 0) has been extensively applied in various life contingency contexts (e.g., Bowers et al. (997)). Lenart (204) provides an expression for the j-th moment, namely j = j! (ln c) j e B ln c E j B, (7) ln c where Es j (z) = R (ln x) j j! x s e zx dx is the generalized integro-exponential function (see Milgram (985)). Here, we consider an example of Melnikov and Romaniuk (2006) on the USA mortality data of the Human mortality database where the parameters B and c were estimated to B = 6: and c = :0959. Using (7), the rst 5 moments are 5 = (76:3437; 6037:202; :3; ; ). We consider the class of mixed Erlang distributions ME res ( m ; A 90 ). The resulting mixed Erlang approximations F W3;90 and F W4;90 are and F W3;90 (x) = 0:054H(x; 22; :0928) + 0:220H(x; 65; :0928) + 0:7637H(x; 90; :0928); F W4;90 (x) = 0:0347H(x; 35; :0972)+0:0834H(x; 66; :0972)+0:009H(x; 67; :0972)+0:780H(x; 90; :0972), with KS distances d KS (S; W 3;90 ) = 0:088, and d KS (S; W 4;90 ) = 0:0239. In Figure 2, we compare the t of the 3- and 4-moment approximations to the exact distribution by plotting their densities (left) and dfs (right). Overall, we observe that the t is quite reasonable. Note that the KS distance increases from the 3-moment to the 4-moment approximation. We notice that both F W3;90 and F W4;90 use the Erlang-90 df, where 90 is the largest element of A 90. As such, a mixed Erlang approximation with a smaller KS distance can likely be found in both cases by choosing a larger support A l (i.e., l > 90) Some remarks To provide insight on the quality of the moment-based mixed Erlang approximation proposed in this section, we brie y revisit the results of the above three examples. In the rst two (lognormal and mixture of two gammas), the mixed Erlang approximation is easy to implement and provides a very satisfactory t to the true distribution. As none of the resulting mixed Erlang approximations uses the 0
11 Figure 2: Density (left) and df (right): 3- and 4-moment approximations vs Gompertz distribution Erlang-70 df (as l = 70 in the rst two examples), it is unlikely that a better approximation (from the viewpoint of KS distance) can be found by increasing the value of l. Given that the best KS- t is found from a mixture of Erlang distributions with relatively small shape parameter (which corresponds to the parameter k in the Erlang df (2)), the proposed method seems particularly well suited for these two cases. As for any approximation method, limitations can also be found as evidenced by the Gompertz example. Indeed for this example, both the 3-moment and 4-moment mixed Erlang approximations use the Erlang-90 df (recall l = 90 for this example). As mentioned earlier, one can likely reduce the KS distance (if so desired) of the resulting mixed Erlang approximation by increasing the value of l (in light of the comments in Remark 4). This implies that the KS-optimal mixed Erlang approximation would likely involve Erlang distributions with large shape parameters (which have smaller variances for a given rate parameter ). In general, distributions with negative skewness and sharp density peak(s) may require the use of Erlang distributions with large shape parameters to provide a good approximation. Computational time of the proposed mixed Erlang methodology may become a non-negligible issue in these cases (especially as the number of moments matched increases). However, a slight adjustment to the proposed methodology may be considered to address this time-consuming issue. Indeed, one may replace the set A l = f; 2; :::; lg in () by a set of the form fa; 2a; :::; jag for positive integers a and j (note that the two sets coincide when a = and j = l). To illustrate this, we have re-considered the Gompertz example by replacing the set A 90 by the set f5; 0; 5; ::::; 200g. The resulting mixed Erlang approximation, denoted by the rv W m;5:200 when m moments are matched (m = 3; 4; 5), are given by F W3;5:200 (x) = 0:0447H(x; 45; :385) + 0:2573H(x; 85; :385) + 0:6980H(x; 0; :385); F W4;5:200 (x) = 0:068H(x; 40; :6422)+0:097H(x; 85; :6422)+0:3044H(x; 5; :6422)+0:588H(x; 40; :6422), and F W5;5:200 (x) = 0:0038H(x; 20; :9095) + 0:0393H(x; 75; :9095) + 0:262H(x; 0; :9095) + 0:3275H(x; 40; :9095) + 0:5032H(x; 65; :9095);
12 Figure 3: Density (left) and df (right): 5-moment approximation vs Gompertz distribution with KS distances d KS (S; W 3;5:200 ) = 0:07, d KS (S; W 4;5:200 ) = 0:0086, and d KS (S; W 5;5:200 ) = 0:0072. Note that the KS distance in the 4-moment approximation is considerably lower than for d KS (S; W 4;90 ) = 0:0239. In Figure 2, we compare the t of the 5-moment approximation W 5;5:200 to the Gompertz distribution by plotting their densities (left) and dfs (right). We can see that the t is quite acceptable and of a better quality than the two approximations displayed in Figure 2. 4 Moment-based mixed Erlang approximation with known 4. Basic de nitions We consider here a slightly di erent context than the one of Section 3. Instead of approximating a general df F S with known moments m, we rather assume that the df F S is known to be of mixed Erlang form () with given rate parameter > 0, and rst m moments m. However, the mixing weights f k g l k= are assumed unknown or di cult to obtain. Various applications in risk theory and credit risk fall into this context (e.g., Cossette et al. (2002), Lindskog and McNeil (2003), and McNeil et al. (2005); see also the application of Section 4.4). For a given rate parameter > 0, let ME( m ; ) be the set of all mixed Erlang distributions with df () (as l! ) and rst m moments m. Also, de ne ME( m ; A l ; ) to be the subset of ME( m ; ) with df () for a given l 2 N +, and let ME ext ( m ; A l ; ) be a further subset of ME( m ; A l ; ) such that at most (m + ) of the mixing weights f k g l k= are non-zero. Note that a distribution in ME( m; A l ; ) can be expressed as a convex combination of distributions in ME ext ( m ; A l ; ) (see, e.g., De Vylder (996)). For a given function, we consider two approaches to derive bounds and approximations for E [ (S)] when S 2 ME( m ; A l ; ) (in cases when the expectation exists). The rst approach is based on discrete s- convex extremal distributions while the second is based on moment bounds on discrete expected stop-loss transforms. 2
13 Remark 5 Naturally, the set ME( m ; A l ; ) tends to ME( m ; ) as l!. As such, when S 2 ME( m ; ), bounds for risk measures on ME( m ; ) can be approximated by their counterparts in ME( m ; A l ; ) for l reasonably large. For a mixed Erlang rv S with df (), its j-th moment is known to satisfy (6). Using the identity j Y (k + i) = i=0 jx s (j; n) k n, n= where the s (j; n) s are the (signed) Stirling numbers of the rst kind (e.g., Abramowitz and Stegun (972)), (6) becomes jx j j = s (j; n) n, (8) n= for j = ; :::; m where n = P l k= kk n. In matrix form, we have M = s T m; where s = fs (j; n)g m j;n= (with s (j; n) = 0 for n > j) and m = ( ; :::; m ). Isolating m yields m = s M T. (9) n o m From Comtet (974, p.23) we know that s c = ( ) i+j c (i; j) where the c (i; j) s are the i;j= Stirling numbers of the second kind de ned as c (i; j) = (j!) P j k=0 ( )j k k j k i (e.g., Abramowitz and Stegun (972)). Thus, for a given > 0, the class ME( m ; A l ; ) can be found through the identi cation of all discrete distributions with support A l and rst m moments given by the right-hand side of (9). Equivalently, the class ME ext ( m ; A l ; ) can be identi ed by restricting the discrete distributions on A l to have at most (m + ) non-zero mass points. This argument is formalized in the next section. 4.2 Discrete s-convex extremal distributions Let D ( m ; A l ) be all discrete distributions with support A l and rst m moments m = ( ; 2 ; :::; m ). Also, denote by D ext ( m ; A l ) the subset of D ( m ; A l ) with distributions having at most (m + ) nonzero mass points. For a given > 0, each distribution in D ( m ; A l ) (and D ext ( m ; A l )) with m as de ned in (9) corresponds to a mixed Erlang distribution in ME( m ; A l ; ) (and ME ext ( m ; A l ; )). This is a one-to-one correspondence (e.g., De Vylder (996, part 2)). Remark 6 Because of this one-to-one correspondence between D( m ; A l ) and ME( m ; A l ; ), conditions under which ME( m ; A l ; ) is not empty can be found from its discrete counterpart D( m ; A l ). We refer the reader to e.g., De Vylder (996), Marceau (996) or Courtois and Denuit (2009). This allows us to make use of the theory developed in Prékopa (990), Denuit and Lefèvre (997), Denuit, Lefèvre and Mes oui (999), and Courtois et al. (2006) to derive bounds/approximations for E [ (S)] when S 2 ME( m ; A l ; ). 3
14 First, we brie y recall the de nitions of s-convex function and s-convex order introduced by Denuit et al. (998). De nition 7 Denuit, Lefèvre and Shaked (998, 2000). Let C be a subinterval of R or a subset of N, and a function on C. For a positive integer s and x 0 < x < < x s 2 C, we recursively de ne the divided di erences as [x 0; x ; :::; x k ] = [x ; x 2 ; :::; x k ] [x 0 ; x ; :::; x k ] x k x 0 = kx i= (x i ), k = ; 2; :::; s; ky (x i x j ) j=0;j6=i where [x k ] = (x k ) for k = 0; ; :::; s. The function is s-convex if [x 0; x ; :::; x s ] 0 for all x 0 < x < < x s 2 C. We mention that the de nition of s-convex function, which refers to higher-convexity, should not be confused with the one for Schur-convex function, also known as S-convex function. De nition 8 Denuit, Lefèvre and Shaked (2000). For two rv s X and Y de ned on C, X is smaller than Y in the s-convex sense, namely X C s cx Y, if E [ (X)] E [ (Y )] for all s-convex functions (provided the expectation exists). We mention that the -convex order corresponds to the usual stochastic dominance order, and the 2-convex order is the usual convex order (see Müller and Stoyan (2002), Denuit et al. (2005), and Shaked and Shanthikumar (2007) for a review on stochastic orders). Also, as stated in Theorem.6.3 of Müller and Stoyan (2002), the s-convex order can only be used to compare rv s with the same rst (s ) moments (which explains why s is chosen to be m + in what follows). Examples of s-convex functions are (x) = x s+j for j 2 N and (x) = exp (cx) for c 0. For a general treatment of the s-convex order, see, e.g., Denuit et al. (999, 2000) and section.6 of Müller and Stoyan (2002). Let K (m+) min and K (m+) max be the (m + )-extremum rv s on D( m ; A l ), i.e. those which satisfy E K (m+) min E [ (K)] E K(m+) max ; for any (m + )-convex function and K 2 D( m ; A l ). The general distribution forms of K (m+) min and K (m+) max are given in Prékopa (990) (see also Courtois et al. (2006, Section 4)), and are repeated here: m + even o support of K (m+) min nj ; j + ; :::; j m+ ; j m o support of K (m+) max n; j ; j + ; :::; j m ; j m + ; l 2 2 n m + odd o ; j ; j + ; :::; j m ; j m n o j ; j + ; :::; j m ; j m + ; l 2 2 (0) where < j < j + < j 2 < ::: < l. From (0), it is clear that the support of K (m+) most m + elements. min(max) has at 4
15 Let W K = P K j= C j be a mixed Erlang rv, i.e. fc j g j are a sequence of iid exponential rv s with mean =, independent of K. The following result of Denuit, Lefèvre and Utev (999, Property 5.7) relates to the stability of the s-convex order under compounding. Lemma 9 If K A l s cx K0, then W K R+ s cx W K 0: We apply Lemma 9 to de ne the mixed Erlang rv s W K(m+) min and W K(m+) max, which are the (m + )-extremum rv s on ME( m ; A l ; ). It is immediate that, for W 2 ME( m ; A l ; ); W K(m+) min R+ (m+) cx W R+ (m+) cx W K (m+) max ; () For instance, using (), the (m + )-convex functions (x) = x m++j (j 2 N) and (x) = exp (cx) (c 0) yield h i E W m++j K (m+) min E W m++j h i E W m++j K (m+) max, and respectively. h i h i E exp cw K(m+) min E [exp (cw )] E exp cw K(m+) max, 4.3 Moment bounds on discrete expected stop-loss transforms Extrema on the (m + )-convex order yield bounds for E [ (W )] when is (m + )-convex and W 2 ME( m ; A l ; ). However, this approach is not appropriate to derive bounds on TVaR and the stop-loss premium when the number of known moments is greater than 2. It is well known that two rv s with the same mean and variance cannot be compared under the convex order. Consequently, we use an approach inspired from Courtois and Denuit (2009) (see also Hürlimann (2002)) to derive bounds on TVaR and the stop-loss premium. We consider D( m ; A l ) and determine lower and upper bounds for E (K k) + for all k 2 Al. From the lower bound, we de ne the corresponding rv K m down on A l via the df F Km infk2d(m;a down (k) = l ) E (K k) + inf K2D(m;Al ) E (K k ) + ; k = ; 2; :::; l ; k = l., (2) for k 2 A l. Similarly, K m up is de ned as in (2) by replacing inf by sup in the de nition. Given that, for K 2 D( m ; A l ), E (K m down k) + E (K k)+ E (Km up k) + ; for all k 2 A l, it implies that K m down(up) is smaller (larger) than K under the increasing convex order, namely K m down icx K icx K m up (see, e.g., Courtois and Denuit (2009)). Note that K m down and K m up do not belong to D( m ; A l ), but both have rst moment. The increasing convex order is stable under compounding, and thus W Km down icx W K icx W Km up : (3) 5
16 Then, from Denuit et al. (2005, Proposition 3.4.8), it follows that TVaR W Km down TVaR (W K ) TVaR W Km up ; (4) for 2 (0; ). Clearly, the rv s W Km down and W Km up will most likely not belong to ME( m ; A l ; ). Remark 0 Note that, when neither of the two aforementioned approaches are applicable, we propose to derive approximate bounds for E [ (S)] with S 2 ME ext ( m ; A l ; ) by calculating E [ (W )] for all W 2 ME ext ( m ; A l ; ) and choosing E [ (W )] min(max) = inf (sup) E [ (W )]. W 2ME ext ( m ;A l ;) Obviously, E [ (S)] does not necessarily lie between E [ (W )] min and E [ (W )] max. However, the interval estimate [E [ (W )] min ; E [ (W )] max ] may give an idea of the variability of all solutions on ME( m ; A l ; ). 4.4 Example: Portfolio of dependent risks We consider a portfolio of n dependent risks as described in the common mixture model of Cossette and al. (2002). Let S = X + ::: + X n be the aggregate claim amount with X i = B i I i. Conditional on a common mixture rv with pmf a j = P ( = j) for j = ; 2; :::, fi i g n i= are assumed to form a sequence of independent Bernoulli rv s with P (I i = j = j ) = (r i ) j for r i 2 (0; ). As for the B i s, they are assumed to form a sequence of iid exponential rv s of mean, independent of fi i g 20 i= and. In this context, it is clear that S is a two-point mixture of a degenerate rv at 0 and a mixed Erlang rv of the form () with l = n and =, i.e. its Laplace transform is given by E e ts p + pe e ty = X Qn a j i= j= (r i ) j + (r i ) j, t 0; + t where p = P j= a jq n i= (r i) j. We perform the moment-based approximation on the rv Y = (S js > 0) rather than S. For illustrative purposes, we assume n = 20 and a logarithmic pmf for, namely a j = (0:5) j =(j ln 2) for j : Also, the constants r i are set such that the (unconditional) mean of I i is q i = E h(r i ) i with q = ::: = q 0 = 0: and q = ::: = q 20 = 0:02. Under the above assumptions, the rst 5 moments of Y are found to be 5 = (:7999; 6:2270; 3:4785; 208:258; 693:7077). Using the approach on discrete (m + )-convex extremal distributions, the dfs F WK(m+) min and F WK(m+) max for m = 4; 5 are: F WK5 min (x) = 0:5365H (x; ; )+0:227H (x; 2; )+0:2232H (x; 3; )+0:0245H (x; 6; )+0:0030H (x; 7; ), F WK5 max (x) = 0:4860H (x; ; )+0:398H (x; 2; )+0:062H (x; 4; )+0:0537H (x; 5; )+0:0000H (x; 20; ), 6
17 F WK6 min (x) = 0:526H (x; ; ) + 0:379H (x; 2; ) + 0:053H (x; 3; ) + 0:082H (x; 4; ) + 0:0059H (x; 7; ) + 0:0023H (x; 8; ), and F WK6 max (x) = 0:5322H (x; ; ) + 0:2264H (x; 2; ) + 0:20H (x; 3; ) + 0:0004H (x; 5; ) + 0:0299H (x; 6; ) + 0:0000H (x; 20; ) : Let X K(m+) min(max) be a rv with df F XK(m+)- min(max) (x) = p + pf WK(m+)- min(max) (x) for x 0. It follows from Section 4.2 that lower (upper) bounds for the higher-order moments E S j (j = 4; 5; :::) and the exponential premium principle de ned as ' (S) = ln E e S = ( > 0) can be found from their counterparts for X K(m+) min(max). A few numerical values are provided in Tables 5 and 6, respectively. i i j E hx j K5 min E hx j K6 min E S j i i E hx j K6 max E hx j K5 max 4 38: : : : : : :880 29:880 29:880 49: : : : : :676 Table 5. Higher-order moments of S, X K(m+) min and X K(m+) max (m = 4; 5) ' X K5 min ' X K6 min ' (S) ' X K6 max ' X K5 max 0:2 :5545 :5546 :5546 :5548 :5564 0: :3536 :3536 :3536 :3536 :3536 0:0 :237 :237 :237 :237 :237 Table 6. Exponential premiums of S, X K(m+) min and X K(m+) max (m = 4; 5) As expected, the bounds get sharper as the number of moments involved increases. As for the second approach based on moment bounds with discrete expected stop-loss transforms, Table 7 presents the values of TVaR for X Km down(up) (m = 4; 5) with df F XKm-down(up) (x) = p + pf WKm-down(up) (x) for x 0. Exact 4 moments 5 moments TVaR (S) TVaR X K4 down TVaR X K4 up TVaR X K5 down TVaR X K5 up 0:9 5:0696 4:9222 5:2062 4:9800 5:490 0:95 6:224 5:9708 6:4548 6:0594 6:3642 0:99 8:8460 8:2899 9:330 8:4655 9:767 0:995 9:9589 9:2500 0:5629 9:4679 0:3775 0:999 2:5066 :422 3:4854 :7323 3:382 Table 7. Values of TVaR for S, X Km down and X Km up (m = 4; 5) As expected, the inequality (4) is veri ed. Also, we observe that the interval estimate of TVaR (S) shrinks as the number of moments matched increases. Finally, given that neither method is applicable for the VaR risk measure, we make use of the technique 7
18 discussed in Remark 0. We nd the numerical values provided in Table 8. Exact ME ext ( 4 ; A 20 ; ) ME ext ( 5 ; A 20 ; ) VaR (S) VaR (X) min VaR (X) max VaR (X) min VaR (X) max 0:9 3:3965 3:3896 3:403 3:3897 3:400 0:95 4:5704 4:5539 4:6030 4:5584 4:5730 0:99 7:2334 7:282 7:2690 7:225 7:2533 0:995 8:3604 8:343 8:3946 8:3539 8:3925 0:999 0:9388 0:79 0:978 0:9226 0:9536 Table 8. Minimal/Maximal values of VaR within ME ext ( m ; A 20 ; ) (m = 4; 5) vs VaR (S) For this example, we observe that the exact values of VaR (S) are within the minimal and maximal values of their corresponding risk measures among all members of ME ext ( m ; A 20 ; ) for m = 4; 5. Also, the spread between the minimal and maximal values of VaR are reduced when we go from 4 to 5 moments. An identical exercise for the TVaR risk measure resulted in the same conclusions. 5 References Abramowitz, M., Stegun, I. (972). Handbook of mathematical functions. 9th Edition. Washington, DC: US Government Printing O ce. Altiok, T. (985). On the Phase-Type Approximations of General Distributions. IIE Transactions 7(2): 0-6. Bobbio, A., Horváth, A., Telek, M. (2005). Matching Three Moments with Minimal Acyclic Phase Type Distributions. Stochastic Models 2: Bowers, N.L., Gerber, H.U., Hickman, G.C., Jones, D.A., Nesbit, C.J. (997). Actuarial Mathematics. The Society of Actuaries, Schaumburg, Illinois. Chaubey, Y.P., Garrido, J., Trudeau, S. (998). On the computation of aggregate claims distributions: some new approximations. Insurance: Mathematics and Economics 23(3): Cheung, E. C. K., Woo, J. K. (204). On the discounted aggregate claim costs until ruin in dependent Sparre Andersen risk processes. Scandinavian Actuarial Journal. In press. Comtet, L. (974). Advanced Combinatorics. Dordrecht-Holland. Boston. Cossette, H., Mailhot, M., Marceau, E. (202). T-Var based capital allocation for multivariate compound distributions. Insurance: Mathematics and Economics 50(2), Cossette, H., Gaillardetz, P., Marceau, E. (2002). Common mixtures in the individual risk model. Bulletin de l Association suisse des actuaires, Courtois, C., Denuit, M. (2009). Moment Bounds on Discrete Expected Stop-Loss Transforms, with Applications. Methodology and Computing in Applied Probability : Courtois, C., Denuit, M. (2007). Bounds on convex reliability functions with known rst moments. European Journal of Operational Research 77:
19 Courtois, C., Denuit, M., Van Bellegem, S. (2006). Discrete s-convex extremal distributions: Theory and applications. Applied Mathematics Letters 9: Daykin, C.D., Pentikäinen, T., Pesonen, H. (994). Practical Risk Theory for Actuaries. Chapman- Hall: New York. Denuit, M., Dhaene, J., Goovaerts, M.J., Kaas, R. (2005). Actuarial Theory for Dependent Risks Measures, Orders and Models. Wiley. New York. Denuit, M., Lefèvre, M. (997). Some new classes of stochastic order relations among arithmetic random variables, with applications in actuarial sciences. Insurance: Mathematics and Economics 20: Denuit, M., Lefèvre, C, Shaked, M. (998). The s-convex orders among real random variables, with applications. Mathematical Inequalities and Applications : Denuit, M., Lefèvre, M., Mes oui, M. (999). On s-convex stochastic extrema for arithmetic risks. Insurance: Mathematics and Economics 25: Denuit, M., Lefèvre, C., Shaked, M. (2000). Stochastic Convexity of the Poisson Mixture Model. Methodology and Computing in Applied Probability 2(3): Denuit, M., Lefèvre, C., Utev, S. (999). Generalized stochastic convexity and stochastic orderings of mixtures. Probability in the Engineering and Informational Sciences 3: De Vylder, F.E. (996). Advanced Risk Theory. A Self-Contained introduction. Editions de l Université Libre de Bruxelles Swiss Association of Actuaries, Bruxelles. Dufresne, D. (2007). Fitting combinations of exponential to probability distributions. Applied Stochastic Models in Business and Industry 23: Gerber, H. U. (979). An Introduction to Mathematical Risk Theory. S.S. Huebner Foundation. University of Pennsylvania. Philadelphie. Hürlimann, W. (2002). Analytical bounds for two value-at-risk functionals. ASTIN Bulletin 32: Johnson, M.A., Taa e, M.R. (989). Matching moments to phase distributions: Mixtures of Erlang distribution of common order. Stochastic Models 5: Karlin, S., Studden, W.J. (966). Tchebyche systems: with applications in analysis and statistics. John Wiley. New York. Landriault, D., Willmot, G.E. (2009). On the joint distributions of the time to ruin, the surplus prior to ruin and the de cit at ruin in the classical risk model. North American Actuarial Journal 3(2): Lee, S. C., Lin, X. S. (200). Modeling and evaluating insurance losses via mixtures of Erlang distributions. North American Actuarial Journal 4(): Lee, Y. S., Lin, T. K. (992). Higher-order Cornish Fisher expansion. Applied Statistics 4: Lenart, A. (204). The moments of the Gompertz distribution and maximum likelihood estimation of its parameters. Scandinavian Actuarial Journal 204:
20 Lindskog, F., McNeil, A. J. (2003). Common Poisson shock models: applications to insurance and credit risk modelling. ASTIN Bulletin 33(2): Marceau, E. (996). Classical risk theory and Schmitter s problems. PhD Thesis. Université Catholique de Louvain, Louvain-la-Neuve. McNeil, A.J., Frey, R., Embrechts, P. (2005). Quantitative Risk Management. Princeton University Press, Princeton. Melnikov, A., Romaniuk, Y. (2006). Evaluating the performance of Gompertz, Makeham and Lee Carter mortality models for risk management with unit-linked contracts. Insurance: Mathematics and Economics 39: Milgram, M (70): The generalized integro-exponential function. Mathematics of Computation Müller, A., Stoyan, D. (2002). Comparison methods for stochastic models and risks. Wiley, New York. Osogami, T. Harchol-Balter, M. (2006). Closed form solutions for mapping general distributions to quasi-minimal PH distributions. Performance Evaluation 63: Prékopa, A. (990). The discrete moment problem and linear programming. Discrete Applied Mathematics 27: Ramsay, C. (99). A note on the normal power approximation. ASTIN Bulletin 2: Seal, H. L. (977). Approximation to risk theory s F(x,t) by means of the gamma distribution. ASTIN Bulletin 9: Shaked, M., Shanthikumar, J.G., (2007). Stochastic orders and their applications. 2nd edition. Springer Series in Statistics. Springer-Verlag, New York. Telek, M., Heindl, A. (2002). Matching moments for acyclic discrete and continuous phase-type distributions of second order. International Journal of Simulation Systems, Science & Technology 3(3 4): Tijms, H.C. (994). Stochastic Models: An Algorithmic Approach. Chichester: John Wiley. Vanden Bosch, P.M., Dietz, D. C., Pohl, E. A. (2000). Moment matching using a family of phasetype distributions. Communications in Statistics - Stochastic Models 6(3-4): Venter, G. (983). Transformed beta and gamma distributions and aggregate losses. Proc. Cas. Act. Soc., Whitt, W. (982). Approximating a Point Process by a Renewal Process, I:Two Basic Methods. Operation Research 30(): Willmot, G.E., Woo, J.K. (2007). On the class of Erlang mixtures with risk theoretic applications. North American Actuarial Journal (2): Willmot, G.E., Lin, X.S. (20). Risk modelling with the mixed Erlang distribution. Applied Stochastic Models in Business and Industry 27:
Moment-Based Approximation with Mixed Erlang Distributions
Moment-Based Approximation with Mixed Erlang Distributions by Hélène Cossette, David Landriault, Etienne Marceau, and Khouzeima Moutanabbir ABSTRACT Moment-based approximations have been extensively analyzed
More informationConditional Tail Expectations for Multivariate Phase Type Distributions
Conditional Tail Expectations for Multivariate Phase Type Distributions Jun Cai Department of Statistics and Actuarial Science University of Waterloo Waterloo, ON N2L 3G1, Canada Telphone: 1-519-8884567,
More informationStochastic Comparisons of Weighted Sums of Arrangement Increasing Random Variables
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 4-7-2015 Stochastic Comparisons of Weighted
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 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 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 informationExplicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables
Explicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables Walter Schneider July 26, 20 Abstract In this paper an analytic expression is given for the bounds
More informationRemarks on quantiles and distortion risk measures
Remarks on quantiles and distortion risk measures Jan Dhaene Alexander Kukush y Daniël Linders z Qihe Tang x Version: October 7, 202 Abstract Distorted expectations can be expressed as weighted averages
More informationTail Conditional Expectations for Extended Exponential Dispersion Models
American Researc Journal of Matematics Original Article ISSN 378-704 Volume 1 Issue 4 015 Tail Conditional Expectations for Extended Exponential Dispersion Models Ye (Zoe) Ye Qiang Wu and Don Hong 1 Program
More informationCharacterization of Upper Comonotonicity via Tail Convex Order
Characterization of Upper Comonotonicity via Tail Convex Order Hee Seok Nam a,, Qihe Tang a, Fan Yang b a Department of Statistics and Actuarial Science, University of Iowa, 241 Schaeffer Hall, Iowa City,
More informationWidely applicable periodicity results for higher order di erence equations
Widely applicable periodicity results for higher order di erence equations István Gy½ori, László Horváth Department of Mathematics University of Pannonia 800 Veszprém, Egyetem u. 10., Hungary E-mail: gyori@almos.uni-pannon.hu
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 informationRisk Aggregation with Dependence Uncertainty
Introduction Extreme Scenarios Asymptotic Behavior Challenges Risk Aggregation with Dependence Uncertainty Department of Statistics and Actuarial Science University of Waterloo, Canada Seminar at ETH Zurich
More informationOn the discounted penalty function in a discrete time renewal risk model with general interclaim times
On the discounted penalty function in a discrete time renewal risk model with general interclaim times Xueyuan Wu, Shuanming Li Centre for Actuarial Studies, Department of Economics The University of Melbourne,
More informationComparison of Approximations for Compound Poisson Processes
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 0, Dublin Session CPS08) p.490 Comparison of Approximations for Compound Poisson Processes Choirat, Christine University of Navarra, Department
More informationNotes on Mathematical Expectations and Classes of Distributions Introduction to Econometric Theory Econ. 770
Notes on Mathematical Expectations and Classes of Distributions Introduction to Econometric Theory Econ. 77 Jonathan B. Hill Dept. of Economics University of North Carolina - Chapel Hill October 4, 2 MATHEMATICAL
More informationStochastic Processes
Introduction and Techniques Lecture 4 in Financial Mathematics UiO-STK4510 Autumn 2015 Teacher: S. Ortiz-Latorre Stochastic Processes 1 Stochastic Processes De nition 1 Let (E; E) be a measurable space
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 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 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 informationEconomics Bulletin, 2012, Vol. 32 No. 1 pp Introduction. 2. The preliminaries
1. Introduction In this paper we reconsider the problem of axiomatizing scoring rules. Early results on this problem are due to Smith (1973) and Young (1975). They characterized social welfare and social
More informationTime Series Models and Inference. James L. Powell Department of Economics University of California, Berkeley
Time Series Models and Inference James L. Powell Department of Economics University of California, Berkeley Overview In contrast to the classical linear regression model, in which the components of the
More informationStochastic integral. Introduction. Ito integral. References. Appendices Stochastic Calculus I. Geneviève Gauthier.
Ito 8-646-8 Calculus I Geneviève Gauthier HEC Montréal Riemann Ito The Ito The theories of stochastic and stochastic di erential equations have initially been developed by Kiyosi Ito around 194 (one of
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 informationMean-Variance Utility
Mean-Variance Utility Yutaka Nakamura University of Tsukuba Graduate School of Systems and Information Engineering Division of Social Systems and Management -- Tennnoudai, Tsukuba, Ibaraki 305-8573, Japan
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 informationUpper stop-loss bounds for sums of possibly dependent risks with given means and variances
Statistics & Probability Letters 57 (00) 33 4 Upper stop-loss bounds for sums of possibly dependent risks with given means and variances Christian Genest a, Etienne Marceau b, Mhamed Mesoui c a Departement
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 informationECON2285: Mathematical Economics
ECON2285: Mathematical Economics Yulei Luo Economics, HKU September 17, 2018 Luo, Y. (Economics, HKU) ME September 17, 2018 1 / 46 Static Optimization and Extreme Values In this topic, we will study goal
More informationSupermodular ordering of Poisson arrays
Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore
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 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 informationSupplemental Material 1 for On Optimal Inference in the Linear IV Model
Supplemental Material 1 for On Optimal Inference in the Linear IV Model Donald W. K. Andrews Cowles Foundation for Research in Economics Yale University Vadim Marmer Vancouver School of Economics University
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 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 informationVaR bounds in models with partial dependence information on subgroups
VaR bounds in models with partial dependence information on subgroups L. Rüschendorf J. Witting February 23, 2017 Abstract We derive improved estimates for the model risk of risk portfolios when additional
More informationON THE MOMENTS OF ITERATED TAIL
ON THE MOMENTS OF ITERATED TAIL RADU PĂLTĂNEA and GHEORGHIŢĂ ZBĂGANU The classical distribution in ruins theory has the property that the sequence of the first moment of the iterated tails is convergent
More informationTail Mutual Exclusivity and Tail- Var Lower Bounds
Tail Mutual Exclusivity and Tail- Var Lower Bounds Ka Chun Cheung, Michel Denuit, Jan Dhaene AFI_15100 TAIL MUTUAL EXCLUSIVITY AND TAIL-VAR LOWER BOUNDS KA CHUN CHEUNG Department of Statistics and Actuarial
More informationAssessing Microdata Disclosure Risk Using the Poisson-Inverse Gaussian Distribution
Assessing Microdata Disclosure Risk Using the Poisson-Inverse Gaussian Distribution Michael Carlson Department of Statistics, Stockholm University SE-106 91 Stockholm, Sweden. June 28, 2002 Abstract An
More informationCorrections to Theory of Asset Pricing (2008), Pearson, Boston, MA
Theory of Asset Pricing George Pennacchi Corrections to Theory of Asset Pricing (8), Pearson, Boston, MA. Page 7. Revise the Independence Axiom to read: For any two lotteries P and P, P P if and only if
More information1 Review of di erential calculus
Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts
More informationON APPROXIMATIONS IN LIMITED FLUCTUATION CREDIBILITY THEORY VINCENT GOULET. Abstract
ON APPROXIMATIONS IN LIMITED FLUCTUATION CREDIBILITY THEORY VINCENT GOULET Abstract Different approximation methods to find a full credibility level in limited fluctuation credibility are studied, and
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 informationLecture 8: Basic convex analysis
Lecture 8: Basic convex analysis 1 Convex sets Both convex sets and functions have general importance in economic theory, not only in optimization. Given two points x; y 2 R n and 2 [0; 1]; the weighted
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 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 informationRobust Solutions to Multi-Objective Linear Programs with Uncertain Data
Robust Solutions to Multi-Objective Linear Programs with Uncertain Data M.A. Goberna yz V. Jeyakumar x G. Li x J. Vicente-Pérez x Revised Version: October 1, 2014 Abstract In this paper we examine multi-objective
More informationConstruction of asymmetric multivariate copulas
Construction of asymmetric multivariate copulas Eckhard Liebscher University of Applied Sciences Merseburg Department of Computer Sciences and Communication Systems Geusaer Straße 0627 Merseburg Germany
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 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 information1 The Well Ordering Principle, Induction, and Equivalence Relations
1 The Well Ordering Principle, Induction, and Equivalence Relations The set of natural numbers is the set N = f1; 2; 3; : : :g. (Some authors also include the number 0 in the natural numbers, but number
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 informationEstimating VaR in credit risk: Aggregate vs single loss distribution
Estimating VaR in credit risk: Aggregate vs single loss distribution M. Assadsolimani and D. Chetalova arxiv:172.4388v1 [q-fin.cp] 14 Feb 217 Abstract Using Monte Carlo simulation to calculate the Value
More informationA new approach for stochastic ordering of risks
A new approach for stochastic ordering of risks Liang Hong, PhD, FSA Department of Mathematics Robert Morris University Presented at 2014 Actuarial Research Conference UC Santa Barbara July 16, 2014 Liang
More informationMATHEMATICAL PROGRAMMING I
MATHEMATICAL PROGRAMMING I Books There is no single course text, but there are many useful books, some more mathematical, others written at a more applied level. A selection is as follows: Bazaraa, Jarvis
More information4.3 - Linear Combinations and Independence of Vectors
- Linear Combinations and Independence of Vectors De nitions, Theorems, and Examples De nition 1 A vector v in a vector space V is called a linear combination of the vectors u 1, u,,u k in V if v can be
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 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 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 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 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 informationSolving Extensive Form Games
Chapter 8 Solving Extensive Form Games 8.1 The Extensive Form of a Game The extensive form of a game contains the following information: (1) the set of players (2) the order of moves (that is, who moves
More informationMultivariate Normal-Laplace Distribution and Processes
CHAPTER 4 Multivariate Normal-Laplace Distribution and Processes The normal-laplace distribution, which results from the convolution of independent normal and Laplace random variables is introduced by
More informationThe Central Limit Theorem: More of the Story
The Central Limit Theorem: More of the Story Steven Janke November 2015 Steven Janke (Seminar) The Central Limit Theorem:More of the Story November 2015 1 / 33 Central Limit Theorem Theorem (Central Limit
More informationAlvaro Rodrigues-Neto Research School of Economics, Australian National University. ANU Working Papers in Economics and Econometrics # 587
Cycles of length two in monotonic models José Alvaro Rodrigues-Neto Research School of Economics, Australian National University ANU Working Papers in Economics and Econometrics # 587 October 20122 JEL:
More informationOn a compound Markov binomial risk model with time-correlated claims
Mathematica Aeterna, Vol. 5, 2015, no. 3, 431-440 On a compound Markov binomial risk model with time-correlated claims Zhenhua Bao School of Mathematics, Liaoning Normal University, Dalian 116029, China
More informationSharp bounds on the VaR for sums of dependent risks
Paul Embrechts Sharp bounds on the VaR for sums of dependent risks joint work with Giovanni Puccetti (university of Firenze, Italy) and Ludger Rüschendorf (university of Freiburg, Germany) Mathematical
More informationEstimating the Number of Common Factors in Serially Dependent Approximate Factor Models
Estimating the Number of Common Factors in Serially Dependent Approximate Factor Models Ryan Greenaway-McGrevy y Bureau of Economic Analysis Chirok Han Korea University February 7, 202 Donggyu Sul University
More informationThe Degree of the Splitting Field of a Random Polynomial over a Finite Field
The Degree of the Splitting Field of a Random Polynomial over a Finite Field John D. Dixon and Daniel Panario School of Mathematics and Statistics Carleton University, Ottawa, Canada fjdixon,danielg@math.carleton.ca
More informationSum of Two Standard Uniform Random Variables
Sum of Two Standard Uniform Random Variables Ruodu Wang http://sas.uwaterloo.ca/~wang Department of Statistics and Actuarial Science University of Waterloo, Canada Dependence Modeling in Finance, Insurance
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 informationSOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL VALUES GUY KATRIEL PREPRINT NO /4
SOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL VALUES GUY KATRIEL PREPRINT NO. 2 2003/4 1 SOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL
More informationMONOTONICITY OF RATIOS INVOLVING INCOMPLETE GAMMA FUNCTIONS WITH ACTUARIAL APPLICATIONS
MONOTONICITY OF RATIOS INVOLVING INCOMPLETE GAMMA FUNCTIONS WITH ACTUARIAL APPLICATIONS EDWARD FURMAN Department of Mathematics and Statistics York University Toronto, Ontario M3J 1P3, Canada EMail: efurman@mathstat.yorku.ca
More informationOn the Convolution Order with Reliability Applications
Applied Mathematical Sciences, Vol. 3, 2009, no. 16, 767-778 On the Convolution Order with Reliability Applications A. Alzaid and M. Kayid King Saud University, College of Science Dept. of Statistics and
More informationCharacterizations on Heavy-tailed Distributions by Means of Hazard Rate
Acta Mathematicae Applicatae Sinica, English Series Vol. 19, No. 1 (23) 135 142 Characterizations on Heavy-tailed Distributions by Means of Hazard Rate Chun Su 1, Qi-he Tang 2 1 Department of Statistics
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 informationTail Conditional Expectations for. Extended Exponential Dispersion Models
Tail Conditional Expectations for Extended Exponential Dispersion Models A Thesis Presented to the Faculty of the Department of Mathematical Sciences Middle Tennessee State University In Partial Fulfillment
More informationTHE MINIMAL OVERLAP COST SHARING RULE
THE MINIMAL OVERLAP COST SHARING RULE M. J. ALBIZURI*, J. C. SANTOS Abstract. In this paper we introduce a new cost sharing rule-the minimal overlap cost sharing rule-which is associated with the minimal
More informationCalculation of Bayes Premium for Conditional Elliptical Risks
1 Calculation of Bayes Premium for Conditional Elliptical Risks Alfred Kume 1 and Enkelejd Hashorva University of Kent & University of Lausanne February 1, 13 Abstract: In this paper we discuss the calculation
More informationIntroduction: structural econometrics. Jean-Marc Robin
Introduction: structural econometrics Jean-Marc Robin Abstract 1. Descriptive vs structural models 2. Correlation is not causality a. Simultaneity b. Heterogeneity c. Selectivity Descriptive models Consider
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 informationNonlinear Programming (NLP)
Natalia Lazzati Mathematics for Economics (Part I) Note 6: Nonlinear Programming - Unconstrained Optimization Note 6 is based on de la Fuente (2000, Ch. 7), Madden (1986, Ch. 3 and 5) and Simon and Blume
More informationBoundedly complete weak-cauchy basic sequences in Banach spaces with the PCP
Journal of Functional Analysis 253 (2007) 772 781 www.elsevier.com/locate/jfa Note Boundedly complete weak-cauchy basic sequences in Banach spaces with the PCP Haskell Rosenthal Department of Mathematics,
More informationThe main purpose of this chapter is to prove the rst and second fundamental theorem of asset pricing in a so called nite market model.
1 2. Option pricing in a nite market model (February 14, 2012) 1 Introduction The main purpose of this chapter is to prove the rst and second fundamental theorem of asset pricing in a so called nite market
More informationChapter 1. GMM: Basic Concepts
Chapter 1. GMM: Basic Concepts Contents 1 Motivating Examples 1 1.1 Instrumental variable estimator....................... 1 1.2 Estimating parameters in monetary policy rules.............. 2 1.3 Estimating
More informationSome results on Denault s capital allocation rule
Faculty of Economics and Applied Economics Some results on Denault s capital allocation rule Steven Vanduffel and Jan Dhaene DEPARTMENT OF ACCOUNTANCY, FINANCE AND INSURANCE (AFI) AFI 0601 Some Results
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 information1 Basic concepts from probability theory
Basic concepts from probability theory This chapter is devoted to some basic concepts from probability theory.. Random variable Random variables are denoted by capitals, X, Y, etc. The expected value or
More informationAsymptotic Bounds for the Distribution of the Sum of Dependent Random Variables
Asymptotic Bounds for the Distribution of the Sum of Dependent Random Variables Ruodu Wang November 26, 2013 Abstract Suppose X 1,, X n are random variables with the same known marginal distribution F
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationa11 a A = : a 21 a 22
Matrices The study of linear systems is facilitated by introducing matrices. Matrix theory provides a convenient language and notation to express many of the ideas concisely, and complicated formulas are
More informationE cient Simulation and Conditional Functional Limit Theorems for Ruinous Heavy-tailed Random Walks
E cient Simulation and Conditional Functional Limit Theorems for Ruinous Heavy-tailed Random Walks Jose Blanchet and Jingchen Liu y June 14, 21 Abstract The contribution of this paper is to introduce change
More informationProbability Transforms with Elliptical Generators
Probability Transforms with Elliptical Generators Emiliano A. Valdez, PhD, FSA, FIAA School of Actuarial Studies Faculty of Commerce and Economics University of New South Wales Sydney, AUSTRALIA Zurich,
More information2 Interval-valued Probability Measures
Interval-Valued Probability Measures K. David Jamison 1, Weldon A. Lodwick 2 1. Watson Wyatt & Company, 950 17th Street,Suite 1400, Denver, CO 80202, U.S.A 2. Department of Mathematics, Campus Box 170,
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 informationGMM-based inference in the AR(1) panel data model for parameter values where local identi cation fails
GMM-based inference in the AR() panel data model for parameter values where local identi cation fails Edith Madsen entre for Applied Microeconometrics (AM) Department of Economics, University of openhagen,
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 informationThe Kuhn-Tucker Problem
Natalia Lazzati Mathematics for Economics (Part I) Note 8: Nonlinear Programming - The Kuhn-Tucker Problem Note 8 is based on de la Fuente (2000, Ch. 7) and Simon and Blume (1994, Ch. 18 and 19). The Kuhn-Tucker
More informationDistance between multinomial and multivariate normal models
Chapter 9 Distance between multinomial and multivariate normal models SECTION 1 introduces Andrew Carter s recursive procedure for bounding the Le Cam distance between a multinomialmodeland its approximating
More information